Implement image mask computation and usecase

The mask is useful to filter the regions that are not interesting,
that could make further processing more difficult.

diff --git a/Bubliny.m b/Bubliny.m
index 0f8f3f569fdf97cf3c7510d9d2dc0e3b98d74471..a2f26fa52ed99fa4308c6a3215b1e7a2ad658e10 100755 (executable)
--- a/Bubliny.m
+++ b/Bubliny.m
@@ -110,4 +110,10 @@ for iii=1:nSelectedDirectory  % Cycles for directories with files
        title('Difference between max and min');\r
        imshow(changefig);\r
 \r
+       mask = get_mask(changefig, meanfig);\r
+       figure\r
+       imshow(meanfig .* uint8(mask));\r
+       figure\r
+       plot(mean(double(meanfig), Weights = mask));\r
+\r
 end  % iii, Cycles for directories with files\r

diff --git a/clean_mask.m b/clean_mask.m
new file mode 100644 (file)
index 0000000..29edf46
--- /dev/null
+++ b/clean_mask.m
@@ -0,0 +1,19 @@
+function mask = clean_mask(mask, strength)
+       if nargin < 2
+               strength = 8;
+       end
+       kernel = [1 1 1; 1 1 1; 1 1 1];
+       mask = 1 - mask;
+       for i = 1:strength
+               mask = conv2(mask, kernel, 'same') > 0;
+       end
+       mask = 1 - mask;
+       for i = 1:2*strength
+               mask = conv2(mask, kernel, 'same') > 0;
+       end
+       mask = 1 - mask;
+       for i = 1:strength
+               mask = conv2(mask, kernel, 'same') > 0;
+       end
+       mask = 1 - mask;
+end

diff --git a/get_mask.m b/get_mask.m
new file mode 100644 (file)
index 0000000..28232d9
--- /dev/null
+++ b/get_mask.m
@@ -0,0 +1,14 @@
+function mask = get_mask(changefig, meanfig, threshold, strength)
+       if nargin < 3
+               threshold = 200;
+       end
+
+       mask = (changefig > threshold);
+       if nargin >= 4
+               mask = clean_mask(mask, strength);
+       else
+               mask = clean_mask(mask);
+       end
+       mask = double(mask);
+end
+