Refactor a small function out
authorLukáš Jiřiště <jiriste@icpf.cas.cz>
Mon, 29 Jul 2024 14:05:42 +0000 (16:05 +0200)
committerLukáš Jiřiště <jiriste@icpf.cas.cz>
Mon, 29 Jul 2024 14:05:42 +0000 (16:05 +0200)
The get_representative_intensities is small but I think it may get
expanded significantly. If not this is still more descriptive I think.

Bubliny.m
get_representative_intensities.m [new file with mode: 0644]

index 01a304de4d457ec7562569a9219d1ec130808985..2249ad20d498d24f380ee1a4a62778a967f63270 100755 (executable)
--- a/Bubliny.m
+++ b/Bubliny.m
@@ -76,32 +76,8 @@ for iii=1:1:nSelectedDirectory  % Cycles for directories with files
        for jjj=startFile:step:nFFiles                                      % Cycles for figures files\r
                ImageName = FFiles(jjj).name;\r
                Figure = imread([FFiles(jjj).folder '/' ImageName]);\r
-               [p,q]=size(Figure);\r
                ContrastFigure = imadjust(Figure);                                  % Adjust the contrast of the image using imadjust.\r
-\r
-%%\r
-% Figure bands for average compution intensity, without the bands with hole and meter\r
-\r
-               ContrastFigure1=ContrastFigure(floor(p*1/16):floor(p/8+p*1/16-3), 1:q);\r
-               [p1,q]=size(ContrastFigure1);\r
-               ContrastFigure2=ContrastFigure(floor(p/2+p/32):floor(p*3/4-p*1/16-5), 1:q);\r
-               [p2,q]=size(ContrastFigure2);\r
-               ContrastFigure3=ContrastFigure(floor(p*3/4+p*1/16):floor(p-p*1/16), 1:q);\r
-               [p3,q]=size(ContrastFigure3);\r
-\r
-               ContrastFigureFin(1:p1,1:q)=ContrastFigure1;\r
-               ContrastFigureFin(p1+1:p1+p2,1:q)=ContrastFigure2;\r
-               ContrastFigureFin(p1+p2+1:p1+p2+p3,1:q)=ContrastFigure3;\r
-\r
-               IntensitiesAverFin{jjj}= mean(ContrastFigureFin);\r
-\r
-               % figure\r
-               % imshow(ContrastFigure1)\r
-               % figure\r
-               % imshow(ContrastFigure2)\r
-               % figure\r
-               % imshow(ContrastFigure3)\r
-\r
+               IntensitiesAverFin{jjj} = get_representative_intensities(ContrastFigure);\r
 \r
 % Figure with average intensity, without the bands with hole and meter\r
                figure\r
@@ -254,9 +230,6 @@ for iii=1:1:nSelectedDirectory  % Cycles for directories with files
 \r
                clear Figure\r
                clear ContrastFigure\r
-               clear ContrastFigure1\r
-               clear ContrastFigure2\r
-               clear ContrastFigure3\r
                clear ContrastFigureFin\r
                % clear IntensitiesAverFin\r
 \r
diff --git a/get_representative_intensities.m b/get_representative_intensities.m
new file mode 100644 (file)
index 0000000..1ca91b9
--- /dev/null
@@ -0,0 +1,20 @@
+function intensities = get_representative_intensities(Figure)
+
+               [p,q]=size(Figure);
+               ContrastFigure = imadjust(Figure);                                  % Adjust the contrast of the image using imadjust.
+
+% Figure bands for average compution intensity, without the bands with hole and meter
+
+               ContrastFigure1=ContrastFigure(floor(p*1/16):floor(p/8+p*1/16-3), 1:q);
+               [p1,q]=size(ContrastFigure1);
+               ContrastFigure2=ContrastFigure(floor(p/2+p/32):floor(p*3/4-p*1/16-5), 1:q);
+               [p2,q]=size(ContrastFigure2);
+               ContrastFigure3=ContrastFigure(floor(p*3/4+p*1/16):floor(p-p*1/16), 1:q);
+               [p3,q]=size(ContrastFigure3);
+
+               ContrastFigureFin(1:p1,1:q)=ContrastFigure1;
+               ContrastFigureFin(p1+1:p1+p2,1:q)=ContrastFigure2;
+               ContrastFigureFin(p1+p2+1:p1+p2+p3,1:q)=ContrastFigure3;
+
+               intensities = mean(ContrastFigureFin);
+end