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
\r
clear Figure\r
clear ContrastFigure\r
- clear ContrastFigure1\r
- clear ContrastFigure2\r
- clear ContrastFigure3\r
clear ContrastFigureFin\r
% clear IntensitiesAverFin\r
\r
--- /dev/null
+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