diff options
| author | Pasca <edoardo.pasca@stfc.ac.uk> | 2017-03-29 16:46:28 +0100 |
|---|---|---|
| committer | Pasca <edoardo.pasca@stfc.ac.uk> | 2017-03-29 16:46:28 +0100 |
| commit | b325933591cd1d0d534a90ad5a417c2d03a0c6f3 (patch) | |
| tree | 716a1c5eab60dfe0622b7f98e6e3b8056cdcae51 /supp/add_wedges.m | |
| download | regularization-b325933591cd1d0d534a90ad5a417c2d03a0c6f3.tar.gz regularization-b325933591cd1d0d534a90ad5a417c2d03a0c6f3.tar.bz2 regularization-b325933591cd1d0d534a90ad5a417c2d03a0c6f3.tar.xz regularization-b325933591cd1d0d534a90ad5a417c2d03a0c6f3.zip | |
Initial revision
Diffstat (limited to 'supp/add_wedges.m')
| -rw-r--r-- | supp/add_wedges.m | 30 |
1 files changed, 30 insertions, 0 deletions
diff --git a/supp/add_wedges.m b/supp/add_wedges.m new file mode 100644 index 0000000..8b8f2a7 --- /dev/null +++ b/supp/add_wedges.m @@ -0,0 +1,30 @@ +% create a wedge mask to simulate the missing wedge + +[rows, columns] = size(sino_zing_rings); +grayImage = ones(rows, columns, 'uint8'); +xCoords = [0 360 0]; +yCoords = [35 7 7]; +mask = poly2mask(xCoords, yCoords, rows, columns); +grayImage(mask) = 0; + +xCoords = [727 360 727]; +yCoords = [35 7 7]; +mask = poly2mask(xCoords, yCoords, rows, columns); +grayImage(mask) = 0; + +xCoords = [0 360 0]; +yCoords = [145 173 173]; +mask = poly2mask(xCoords, yCoords, rows, columns); +grayImage(mask) = 0; + +xCoords = [727 360 727]; +yCoords = [145 173 173]; +mask = poly2mask(xCoords, yCoords, rows, columns); +grayImage(mask) = 0; + +grayImage(1:7,:) = 0; +grayImage(173:end,:) = 0; + +%figure; imshow(grayImage, [0 1]); +MW_sino_artifacts = sino_zing_rings.*double(grayImage); +%Dweights = Dweights.*double(grayImage);
\ No newline at end of file |