clc; % clears the Command Window display\r
\r
%% vstupni data T:\00-Orvalho\Travnickova\Bubble column\r
-% \r
+%\r
\r
%%\r
% User dialog for the data directory used for files\r
\r
list={DDirectory.name};\r
[Selected,OK] = listdlg('ListString',list,'ListSize',[200 500]);\r
-SelectedDirectory = list(Selected)'; \r
-if isempty(SelectedDirectory) % handles the error, when no spreadsheet is selected \r
- error('No experiment selected.')\r
+SelectedDirectory = list(Selected)';\r
+if isempty(SelectedDirectory) % handles the error, when no spreadsheet is selected\r
+ error('No experiment selected.')\r
end\r
nSelectedDirectory = numel(SelectedDirectory); % Counts directories\r
\r
%%\r
% read names of all files\r
\r
- FileType = '*.png'; % extension of image file\r
- FileName=([DirectoryName '/' cell2mat(SelectedDirectory(iii)) '/' FileType]);\r
- FFiles = dir(FileName); % Loads files to struct FFiles(1).name, FFiles(1).folder, FFiles(1).date - i s casem\r
- nFFiles = numel(FFiles); % Counts files\r
- \r
+ FileType = '*.png'; % extension of image file\r
+ FileName=([DirectoryName '/' cell2mat(SelectedDirectory(iii)) '/' FileType]);\r
+ FFiles = dir(FileName); % Loads files to struct FFiles(1).name, FFiles(1).folder, FFiles(1).date - i s casem\r
+ nFFiles = numel(FFiles); % Counts files\r
+\r
%% variable allocation\r
\r
- LevelNum = zeros(1,nFFiles);\r
- Heihgt = zeros(1,nFFiles);\r
- Metr_fit_H = zeros(1,nFFiles);\r
- Botom = zeros(1,nFFiles);\r
- Metr_fit_B = zeros(1,nFFiles);\r
+ LevelNum = zeros(1,nFFiles);\r
+ Heihgt = zeros(1,nFFiles);\r
+ Metr_fit_H = zeros(1,nFFiles);\r
+ Botom = zeros(1,nFFiles);\r
+ Metr_fit_B = zeros(1,nFFiles);\r
\r
%%\r
% load of all files in cycles\r
\r
- step=1;\r
- startFile= 40; %945; %200 (20. )\r
- endFile= 110; %nFFiles; %1000; %325 (20.)\r
- nFFiles=endFile; % Attention!!!\r
- for jjj=[1 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
+ step=1;\r
+ startFile= 40; %945; %200 (20. )\r
+ endFile= 110; %nFFiles; %1000; %325 (20.)\r
+ nFFiles=endFile; % Attention!!!\r
+ for jjj=[1 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
- IntensitaAverFin{jjj}= mean(ContrastFigureFin);\r
- \r
- % figure\r
- % imshow(ContrastFigure1)\r
- % figure\r
- % imshow(ContrastFigure2)\r
- % figure\r
- % imshow(ContrastFigure3)\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
+ IntensitaAverFin{jjj}= mean(ContrastFigureFin);\r
+\r
+ % figure\r
+ % imshow(ContrastFigure1)\r
+ % figure\r
+ % imshow(ContrastFigure2)\r
+ % figure\r
+ % imshow(ContrastFigure3)\r
\r
\r
%%\r
% Load the dimensions from the 1st image (p&q in pxl)\r
- if jjj == 1\r
+ if jjj == 1\r
\r
%%\r
-% Meter \r
- D_num=str2double(SelectedDirectory{iii}(21)); % position of camera\r
-\r
- reply = input('Is this a NEW EXPERIMENT that does not yet have calculated parameters for the cm scale? (y/n): ','s'); % y, n Selects where to find parameters for the cm scale. Either in the array or guess them. \r
- if strcmp(reply,'n')\r
- AxisIntensArray=[22,22,22,22,25,27,7,7,7,7,7,23,3,20,11,2,4,19,6,25,3]; % this needs to be guessed\r
- AxisIntensUpArray=[60,60,65,65,34,5,29,29,29,29,29,22,26,2,27,2,28,2,29,2,29]; % this needs to be guessed\r
- AxisEmptyArray=[190,190,190,190,200,190,190,190,190,190,190,190,190,190,190,190,190,190,190,190,190]; % this needs to be guessed 190\r
- LevelMeterArray=[0.95,0.95,0.95,1,0.11,0.8,0.2,0.2,0.2,0.1,0.29,0.9,0.39,0.9,0.39,0.6,0.29,0.9,0.25,0.7,0.27];\r
- AxisIntens=AxisIntensArray(Selected-13); % The number of all files minus 9 in which there is no data\r
- AxisIntensUp=AxisIntensUpArray(Selected-13); \r
- AxisEmpty=AxisEmptyArray(Selected-13); \r
- LevelMeter=LevelMeterArray(Selected-13); \r
- elseif strcmp(reply,'y')\r
- AxisIntens=25; % 1.- 4. =22; 5. - 29 6.-27 7.8.9.11.11.-7 12.=23 13.=3 % Right this needs to be guessed\r
- AxisIntensUp=2; % 3. = 65; 1. a 2. = 60; 5. - 34 6. a 7.8.9.10.11. = 29 12.=22 13.=26 % Left this needs to be guessed\r
- AxisEmpty=190; % 5. = 200; 6. 7.8.9.10.11.12.13 = 190 % this needs to be guessed 190\r
- LevelMeter=0.7; % 5. = 0.11, 4. = 1; 1., 2. a 3. = 0.95, 6. =0.9; 7.8.9.=0.2 10.=0.1 11.=0.29 12.=0.9 13.=0.39 % number on the meter so that the left dash is above the surface and the right is below the surface\r
- end\r
- \r
- Meter=ContrastFigure(floor(p*3/4-p*1/16+p*1/64):floor(p*3/4+p*1/64), 1:q);\r
- ContrastMeter = imadjust(Meter);\r
- [r,s]=size(ContrastMeter);\r
- ContrastMeterMid=ContrastMeter(floor(end/2),:);\r
- \r
- hf1=figure ('Name','Contrasted meter','NumberTitle','off','PaperSize',[8.5 6],'PaperUnits','centimeters'); \r
- imshow(ContrastMeter)\r
- hold on\r
- plot(ContrastMeterMid,'-m','linewidth',2)\r
- hold on\r
-\r
- for kkk=2:1:q \r
- if kkk<(1-LevelMeter)*q && ContrastMeterMid(kkk)>AxisIntensUp % looks for intensity minima for dashes on the axis above the level, Left\r
- ContrastMeterMid(kkk)=0;\r
- elseif kkk>=(1-LevelMeter)*q && ContrastMeterMid(kkk)>AxisIntens % looks for intensity minima for dashes on the axis under the level, Right\r
- ContrastMeterMid(kkk)=0;\r
- else\r
- ContrastMeterMid(kkk)=1;\r
- end\r
- if ContrastMeterMid(kkk) == ContrastMeterMid(kkk-1) % looking for the right side of the peak\r
- ContrastMeterMid(kkk-1)=0;\r
- end\r
- end\r
- for mmm=q:-1:AxisEmpty \r
- if ContrastMeterMid(mmm)==1 % erases the peaks (written numbers) to the left of the dashes on the axis\r
- ContrastMeterMid((mmm-AxisEmpty):mmm-1)=0;\r
- mmm=mmm-AxisEmpty;\r
- end\r
- end\r
- ooo=1;\r
- for nnn=1:1:q \r
- if ContrastMeterMid(nnn)==1 % Writes x places of dashes on the axis\r
- MeterScale(ooo)=nnn;\r
- ooo=ooo+1;\r
- end\r
- end\r
-\r
- if D_num==1\r
- plot(MeterScale, floor(r/2) * ones(10, 1),'oc','linewidth',2)\r
- elseif D_num==2\r
- plot(MeterScale(1:end-1), floor(r/2) * ones(8, 1),'oc','linewidth',2)\r
- end\r
+% Meter\r
+ D_num=str2double(SelectedDirectory{iii}(21)); % position of camera\r
+\r
+ reply = input('Is this a NEW EXPERIMENT that does not yet have calculated parameters for the cm scale? (y/n): ','s'); % y, n Selects where to find parameters for the cm scale. Either in the array or guess them.\r
+ if strcmp(reply,'n')\r
+ AxisIntensArray=[22,22,22,22,25,27,7,7,7,7,7,23,3,20,11,2,4,19,6,25,3]; % this needs to be guessed\r
+ AxisIntensUpArray=[60,60,65,65,34,5,29,29,29,29,29,22,26,2,27,2,28,2,29,2,29]; % this needs to be guessed\r
+ AxisEmptyArray=[190,190,190,190,200,190,190,190,190,190,190,190,190,190,190,190,190,190,190,190,190]; % this needs to be guessed 190\r
+ LevelMeterArray=[0.95,0.95,0.95,1,0.11,0.8,0.2,0.2,0.2,0.1,0.29,0.9,0.39,0.9,0.39,0.6,0.29,0.9,0.25,0.7,0.27];\r
+ AxisIntens=AxisIntensArray(Selected-13); % The number of all files minus 9 in which there is no data\r
+ AxisIntensUp=AxisIntensUpArray(Selected-13);\r
+ AxisEmpty=AxisEmptyArray(Selected-13);\r
+ LevelMeter=LevelMeterArray(Selected-13);\r
+ elseif strcmp(reply,'y')\r
+ AxisIntens=25; % 1.- 4. =22; 5. - 29 6.-27 7.8.9.11.11.-7 12.=23 13.=3 % Right this needs to be guessed\r
+ AxisIntensUp=2; % 3. = 65; 1. a 2. = 60; 5. - 34 6. a 7.8.9.10.11. = 29 12.=22 13.=26 % Left this needs to be guessed\r
+ AxisEmpty=190; % 5. = 200; 6. 7.8.9.10.11.12.13 = 190 % this needs to be guessed 190\r
+ LevelMeter=0.7; % 5. = 0.11, 4. = 1; 1., 2. a 3. = 0.95, 6. =0.9; 7.8.9.=0.2 10.=0.1 11.=0.29 12.=0.9 13.=0.39 % number on the meter so that the left dash is above the surface and the right is below the surface\r
+ end\r
+\r
+ Meter=ContrastFigure(floor(p*3/4-p*1/16+p*1/64):floor(p*3/4+p*1/64), 1:q);\r
+ ContrastMeter = imadjust(Meter);\r
+ [r,s]=size(ContrastMeter);\r
+ ContrastMeterMid=ContrastMeter(floor(end/2),:);\r
+\r
+ hf1=figure ('Name','Contrasted meter','NumberTitle','off','PaperSize',[8.5 6],'PaperUnits','centimeters');\r
+ imshow(ContrastMeter)\r
+ hold on\r
+ plot(ContrastMeterMid,'-m','linewidth',2)\r
+ hold on\r
+\r
+ for kkk=2:1:q\r
+ if kkk<(1-LevelMeter)*q && ContrastMeterMid(kkk)>AxisIntensUp % looks for intensity minima for dashes on the axis above the level, Left\r
+ ContrastMeterMid(kkk)=0;\r
+ elseif kkk>=(1-LevelMeter)*q && ContrastMeterMid(kkk)>AxisIntens % looks for intensity minima for dashes on the axis under the level, Right\r
+ ContrastMeterMid(kkk)=0;\r
+ else\r
+ ContrastMeterMid(kkk)=1;\r
+ end\r
+ if ContrastMeterMid(kkk) == ContrastMeterMid(kkk-1) % looking for the right side of the peak\r
+ ContrastMeterMid(kkk-1)=0;\r
+ end\r
+ end\r
+ for mmm=q:-1:AxisEmpty\r
+ if ContrastMeterMid(mmm)==1 % erases the peaks (written numbers) to the left of the dashes on the axis\r
+ ContrastMeterMid((mmm-AxisEmpty):mmm-1)=0;\r
+ mmm=mmm-AxisEmpty;\r
+ end\r
+ end\r
+ ooo=1;\r
+ for nnn=1:1:q\r
+ if ContrastMeterMid(nnn)==1 % Writes x places of dashes on the axis\r
+ MeterScale(ooo)=nnn;\r
+ ooo=ooo+1;\r
+ end\r
+ end\r
+\r
+ if D_num==1\r
+ plot(MeterScale, floor(r/2) * ones(10, 1),'oc','linewidth',2)\r
+ elseif D_num==2\r
+ plot(MeterScale(1:end-1), floor(r/2) * ones(8, 1),'oc','linewidth',2)\r
+ end\r
%--------------------------------------------------------------------------\r
- print ('-dpng', '-r1200', ['Contrasted_meter.png'])\r
+ print ('-dpng', '-r1200', ['Contrasted_meter.png'])\r
%--------------------------------------------------------------------------\r
\r
\r
- hf2=figure ('Name','Intensity minima for dashes on the axis','NumberTitle','off','PaperSize',[8.5 6],'PaperUnits','centimeters'); \r
- plot(ContrastMeterMid,'-m','linewidth',2)\r
- \r
- xlim([1,q]);\r
- ylim([0,1]);\r
+ hf2=figure ('Name','Intensity minima for dashes on the axis','NumberTitle','off','PaperSize',[8.5 6],'PaperUnits','centimeters');\r
+ plot(ContrastMeterMid,'-m','linewidth',2)\r
+\r
+ xlim([1,q]);\r
+ ylim([0,1]);\r
%--------------------------------------------------------------------------\r
- print ('-dpng', '-r1200', ['Intensity_places_dashes.png'])\r
+ print ('-dpng', '-r1200', ['Intensity_places_dashes.png'])\r
%--------------------------------------------------------------------------\r
clear Meter\r
clear ContrastMeter\r
clear ContrastMeterMid\r
\r
% Calculation of the actual position in mm using the scale of the meter\r
- \r
- hf3=figure ('Name','Dashes scale','NumberTitle','off','PaperSize',[8.5 6],'PaperUnits','centimeters'); \r
- if D_num==1\r
- plot(MeterScale,[1000;900;800;700;600;500;400;300;200;100],'ob','linewidth',2)\r
- hold on\r
- polypar_Metr = polyfit(MeterScale,[1000;900;800;700;600;500;400;300;200;100],2);\r
- elseif D_num==2\r
- plot(MeterScale(1:end-1),[1800;1700;1600;1500;1400;1300;1200;1100],'ob','linewidth',2)\r
- hold on\r
- polypar_Metr = polyfit(MeterScale(1:end-1),[1800;1700;1600;1500;1400;1300;1200;1100],2);\r
- end\r
- no_layer_fit=1: 1 : q;\r
- Metr_fit=polyval(polypar_Metr,no_layer_fit);\r
- plot(no_layer_fit,Metr_fit,'-r','linewidth',2);\r
- hold on\r
- if D_num==1\r
- xlim([0,q]);\r
- ylim([0,1200]);\r
- elseif D_num==2\r
- xlim([0,q]);\r
- ylim([0,2000]);\r
- end\r
+\r
+ hf3=figure ('Name','Dashes scale','NumberTitle','off','PaperSize',[8.5 6],'PaperUnits','centimeters');\r
+ if D_num==1\r
+ plot(MeterScale,[1000;900;800;700;600;500;400;300;200;100],'ob','linewidth',2)\r
+ hold on\r
+ polypar_Metr = polyfit(MeterScale,[1000;900;800;700;600;500;400;300;200;100],2);\r
+ elseif D_num==2\r
+ plot(MeterScale(1:end-1),[1800;1700;1600;1500;1400;1300;1200;1100],'ob','linewidth',2)\r
+ hold on\r
+ polypar_Metr = polyfit(MeterScale(1:end-1),[1800;1700;1600;1500;1400;1300;1200;1100],2);\r
+ end\r
+ no_layer_fit=1: 1 : q;\r
+ Metr_fit=polyval(polypar_Metr,no_layer_fit);\r
+ plot(no_layer_fit,Metr_fit,'-r','linewidth',2);\r
+ hold on\r
+ if D_num==1\r
+ xlim([0,q]);\r
+ ylim([0,1200]);\r
+ elseif D_num==2\r
+ xlim([0,q]);\r
+ ylim([0,2000]);\r
+ end\r
%--------------------------------------------------------------------------\r
- print ('-dpng', '-r1200', 'Dashes_scale.png')\r
+ print ('-dpng', '-r1200', 'Dashes_scale.png')\r
%--------------------------------------------------------------------------\r
\r
- xlswrite(ResFile,{'Contrasted meter'},['Met ' cell2mat(SelectedDirectory(iii))],'A2'); %[cell2mat(SelectedDirectory(iii)) ': Meter scale'],'A1'); ['Meter scale ' num2str(iii)]\r
- %%xlswritefig(hf1, ResFile, ['Met ' cell2mat(SelectedDirectory(iii))] , 'A4')\r
- xlswrite(ResFile,{'Intensity: places of dashes on the axis'},['Met ' cell2mat(SelectedDirectory(iii))],'A6');\r
- %%xlswritefig(hf2, ResFile, ['Met ' cell2mat(SelectedDirectory(iii))], 'B8')\r
- xlswrite(ResFile,{'Dashes scale'},['Met ' cell2mat(SelectedDirectory(iii))],'K6'); \r
- %%xlswritefig(hf3, ResFile, ['Met ' cell2mat(SelectedDirectory(iii))], 'L8')\r
- xlswrite(ResFile,[{'Meter Scale'} {'Dashes scale'}],['Met ' cell2mat(SelectedDirectory(iii))],'A28'); \r
- xlswrite(ResFile,[{'AxisIntens'} {'AxisIntensUp'} {'AxisEmpty'} {'LevelMeter'}],['Met ' cell2mat(SelectedDirectory(iii))],'D28'); \r
- if D_num==1\r
- xlswrite(ResFile,[MeterScale' [1000;900;800;700;600;500;400;300;200;100]],['Met ' cell2mat(SelectedDirectory(iii))],'A29'); \r
- elseif D_num==2\r
- xlswrite(ResFile,[MeterScale(1:end-1)' [1800;1700;1600;1500;1400;1300;1200;1100]],['Met ' cell2mat(SelectedDirectory(iii))],'A29'); \r
- end\r
- xlswrite(ResFile,[{AxisIntens} {AxisIntensUp} {AxisEmpty} {LevelMeter}],['Met ' cell2mat(SelectedDirectory(iii))],'D29'); \r
-\r
- end\r
+ xlswrite(ResFile,{'Contrasted meter'},['Met ' cell2mat(SelectedDirectory(iii))],'A2'); %[cell2mat(SelectedDirectory(iii)) ': Meter scale'],'A1'); ['Meter scale ' num2str(iii)]\r
+ %%xlswritefig(hf1, ResFile, ['Met ' cell2mat(SelectedDirectory(iii))] , 'A4')\r
+ xlswrite(ResFile,{'Intensity: places of dashes on the axis'},['Met ' cell2mat(SelectedDirectory(iii))],'A6');\r
+ %%xlswritefig(hf2, ResFile, ['Met ' cell2mat(SelectedDirectory(iii))], 'B8')\r
+ xlswrite(ResFile,{'Dashes scale'},['Met ' cell2mat(SelectedDirectory(iii))],'K6');\r
+ %%xlswritefig(hf3, ResFile, ['Met ' cell2mat(SelectedDirectory(iii))], 'L8')\r
+ xlswrite(ResFile,[{'Meter Scale'} {'Dashes scale'}],['Met ' cell2mat(SelectedDirectory(iii))],'A28');\r
+ xlswrite(ResFile,[{'AxisIntens'} {'AxisIntensUp'} {'AxisEmpty'} {'LevelMeter'}],['Met ' cell2mat(SelectedDirectory(iii))],'D28');\r
+ if D_num==1\r
+ xlswrite(ResFile,[MeterScale' [1000;900;800;700;600;500;400;300;200;100]],['Met ' cell2mat(SelectedDirectory(iii))],'A29');\r
+ elseif D_num==2\r
+ xlswrite(ResFile,[MeterScale(1:end-1)' [1800;1700;1600;1500;1400;1300;1200;1100]],['Met ' cell2mat(SelectedDirectory(iii))],'A29');\r
+ end\r
+ xlswrite(ResFile,[{AxisIntens} {AxisIntensUp} {AxisEmpty} {LevelMeter}],['Met ' cell2mat(SelectedDirectory(iii))],'D29');\r
+\r
+ end\r
%%\r
% Figure with average intensity, without the bands with hole and meter\r
- figure\r
- imshow(ContrastFigure)\r
- hold on\r
- plot(IntensitaAverFin{jjj},'-y','linewidth',2)\r
- hold on\r
- if jjj>=startFile+step %subtracts the intensities\r
- plot(IntensitaAverFin{jjj-step},'Color',[0.9290 0.6940 0.1250],'linewidth',2)\r
- hold on\r
- IntensitaAverFin_2{jjj}=IntensitaAverFin{jjj}-IntensitaAverFin{jjj-step};\r
- waytrend=mean(IntensitaAverFin_2{jjj},2); % waytrend will tell whether to calculate the minimum or maximum green intensity IntensitaAverFin_2{jjj}\r
- if waytrend<=0;\r
- plot(-IntensitaAverFin_2{jjj},'-g','linewidth',2)\r
- else\r
- plot(IntensitaAverFin_2{jjj},'-g','linewidth',2)\r
- end\r
- hold on\r
- end\r
-\r
- infmt = 'hh:mm:ss.SS';\r
- Time(jjj)=duration([(FFiles(jjj).name(8:9)) ':' (FFiles(jjj).name(10:11)) ':' (FFiles(jjj).name(12:16))],'Format',infmt)-duration([(FFiles(1).name(8:9)) ':' (FFiles(1).name(10:11)) ':' (FFiles(1).name(12:16))],'Format',infmt);\r
- infmd = ['yy' 'MM' 'dd'];\r
- Date(jjj) = datetime(FFiles(jjj).name(1:6),'InputFormat',infmd);\r
- Level(jjj,1:3)=SelectedDirectory{iii}(7:9);\r
- LevelNum(1,jjj)=str2double(Level(jjj,:))*1000;\r
-\r
- % % plot(lk,abs(pk),'oc','linewidth',2)\r
- if D_num==1\r
- LevelPixels=MeterScale(1+(1000-LevelNum(1,jjj))/100); % x level in pixels in the original image\r
- [pk,lk] = findpeaks(-IntensitaAverFin{jjj}(1:(LevelPixels+30))); %10\r
- [Heihgt(1,jjj),pkH]=max(pk);\r
- qH=lk(pkH);\r
- Metr_fit_H(1,jjj)=polyval(polypar_Metr,qH);\r
- plot(qH,-Heihgt(jjj),'or','linewidth',2)\r
- if jjj>=startFile+step \r
- LevelPixels=MeterScale(1+(1000-LevelNum(1,jjj))/100); % x level in pixels in the original image\r
- if waytrend<=0;\r
- [pk,lk] = findpeaks(-IntensitaAverFin_2{jjj}(1:(LevelPixels+30))); %10\r
- else\r
- [pk,lk] = findpeaks(IntensitaAverFin_2{jjj}(1:(LevelPixels+30))); %10\r
- end\r
- [Heihgt(1,jjj),pkH]=max(pk);\r
- if max(pk)>=10; % error cleaning %30\r
- qH=lk(pkH);\r
- Metr_fit_H(1,jjj)=polyval(polypar_Metr,qH);\r
- end\r
- plot(qH,Heihgt(jjj),'om','linewidth',2)\r
- end\r
- \r
- elseif D_num==2\r
- if jjj<startFile+step \r
- LevelPixels=MeterScale(1); % x level in pixels in the original image\r
- [pk,lk] = findpeaks(-IntensitaAverFin{jjj}(LevelPixels:end)); %10\r
- [Heihgt(1,jjj),pkH]=max(pk);\r
- qH=LevelPixels+lk(pkH);\r
- Metr_fit_H(1,jjj)=polyval(polypar_Metr,qH);\r
- plot(qH,-Heihgt(jjj),'or','linewidth',2)\r
- elseif jjj>=startFile+step \r
- LevelPixels=MeterScale(1); % x level in pixels in the original image\r
- if waytrend<=0;\r
- [pk,lk] = findpeaks(-IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
- else\r
- [pk,lk] = findpeaks(IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
- end\r
- [Heihgt(1,jjj),pkH]=max(pk);\r
- qH=LevelPixels+lk(pkH);\r
- Metr_fit_H(1,jjj)=polyval(polypar_Metr,qH);\r
- if max(pk)<30; % error cleaning %30\r
- qH=NaN;\r
- Metr_fit_H(1,jjj)=NaN;\r
- end\r
- plot(qH,Heihgt(jjj),'om','linewidth',2)\r
- end\r
- end\r
-\r
- if D_num==1\r
- if jjj<startFile+step \r
- [pk2,lk2] = findpeaks(-IntensitaAverFin{jjj},'MinPeakProminence',5); % Small peaks at the bubble/liquid interface\r
- \r
- [pk3,lk3] = findpeaks(-IntensitaAverFin{jjj},'MinPeakProminence',5,'MinPeakDistance',30); % Throw out small random peaks in the liquid, peaks far apart cca 1cm\r
- [one,nqq]=size(lk3);\r
- [one,npp]=size(lk2);\r
- \r
- if npp==0\r
- pk2=NaN;\r
- lk2=NaN;\r
- elseif nqq==0\r
- pk3=NaN;\r
- lk3=NaN;\r
- end\r
- \r
- for qqq=1:1:nqq\r
- for ppp=1:1:npp\r
- if lk3(qqq)==lk2(ppp)\r
- pk2(ppp)=NaN;\r
- lk2(ppp)=NaN;\r
- end\r
- end\r
- end\r
- \r
- % % plot(lk2,abs(pk2),'om','linewidth',2)\r
- hold on\r
- \r
- [Botom(1,jjj),pk2B]=max(abs(pk2)); % For 1-4 files\r
- qB=lk2(pk2B); % For 1-4 files\r
- \r
- Metr_fit_B(1,jjj)=polyval(polypar_Metr,qB);\r
- if Metr_fit_B(1,jjj)>=(Metr_fit_H(1,jjj)-2) % Throw out the bubble/liquid interface if it is above the surface, 2cm 2D interface\r
- Metr_fit_B(1,jjj)=NaN;\r
- qB=NaN;\r
- elseif Botom(1,jjj)<=(max(abs(pk))+min(abs(pk)))*0.62 % 0.65 is an estimate of how much intensity the bubbles need\r
- Metr_fit_B(1,jjj)=NaN; % For 1-4 files\r
- qB=NaN;\r
- end\r
- plot(qB,Botom(jjj),'ob','linewidth',2)\r
- end\r
- \r
- if jjj>=startFile+step \r
- LevelPixels=MeterScale(1+(1000-LevelNum(1,jjj))/100); % x level in pixels in the original image\r
- if waytrend<=0;\r
- [pk2,lk2] = findpeaks(-IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
- else\r
- [pk2,lk2] = findpeaks(IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
- end \r
- [Botom(1,jjj),pk2B]=max(pk2);\r
- qB=LevelPixels+lk2(pk2B); \r
- Metr_fit_B(1,jjj)=polyval(polypar_Metr,qB);\r
- if max(pk2)<=45; % error cleaning\r
- qB=NaN;\r
- Metr_fit_B(1,jjj)=NaN;\r
- end\r
- plot(qB,Botom(jjj),'oc','linewidth',2)\r
- end\r
- \r
- elseif D_num==2\r
- if jjj<startFile+step \r
- [pk2,lk2] = findpeaks(-IntensitaAverFin{jjj},'MinPeakProminence',5); % Small peaks at the bubble/liquid interface\r
- \r
- [pk3,lk3] = findpeaks(-IntensitaAverFin{jjj},'MinPeakProminence',5,'MinPeakDistance',30); % Throw out small random peaks in the liquid, peaks far apart cca 1cm\r
- [one,nqq]=size(lk3);\r
- [one,npp]=size(lk2);\r
- \r
- if npp==0\r
- pk2=NaN;\r
- lk2=NaN;\r
- elseif nqq==0\r
- pk3=NaN;\r
- lk3=NaN;\r
- end\r
- \r
- for qqq=1:1:nqq\r
- for ppp=1:1:npp\r
- if lk3(qqq)==lk2(ppp)\r
- pk2(ppp)=NaN;\r
- lk2(ppp)=NaN;\r
- end\r
- end\r
- end\r
- \r
- % % plot(lk2,abs(pk2),'om','linewidth',2)\r
- hold on\r
- \r
- [Botom(1,jjj),pk2B]=max(abs(pk2)); % For 1-4 files\r
- qB=lk2(pk2B); % For 1-4 files\r
- \r
- % [Botom(1,jjj),pk2B]=max(abs(pk)); % For 5. files\r
- % qB=lk(pk2B); % For 5. files\r
- \r
- Metr_fit_B(1,jjj)=polyval(polypar_Metr,qB);\r
- if Metr_fit_B(1,jjj)>=(Metr_fit_H(1,jjj)-2) % Throw out the bubble/liquid interface if it is above the surface, 2cm 2D interface\r
- Metr_fit_B(1,jjj)=NaN;\r
- qB=NaN;\r
- elseif Botom(1,jjj)<=(max(abs(pk))+min(abs(pk)))*0.62 % 0.65 is an estimate of how much intensity the bubbles need\r
- Metr_fit_B(1,jjj)=NaN; % For 1-4 files\r
- qB=NaN;\r
- end\r
- plot(qB,Botom(jjj),'ob','linewidth',2)\r
- end\r
- \r
- if jjj>=startFile+step \r
- LevelPixels=MeterScale(1); % x level in pixels in the original image\r
- if waytrend<=0;\r
- [pk2,lk2] = findpeaks(-IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
- else\r
- [pk2,lk2] = findpeaks(IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
- end\r
- [Botom(1,jjj),pk2B]=max(pk2);\r
- qB=LevelPixels+lk2(pk2B); \r
- Metr_fit_B(1,jjj)=polyval(polypar_Metr,qB);\r
- if max(pk2)<=20; % error cleaning\r
- qB=NaN;\r
- Metr_fit_B(1,jjj)=NaN;\r
- end\r
- plot(qB,Botom(jjj),'oc','linewidth',2)\r
- end\r
- end\r
-\r
- xlswrite(ResFile,[{'ImageName'} {'Date'} {'Time [s]'} {'Intensita Heihgt [-]'} {'Intensita Botom [-]'} {'Heihgt [cm]'} {'Botom [cm]'} {'Level [cm]'}],cell2mat(SelectedDirectory(iii)),'A1:H1'); \r
- xlswrite(ResFile,[{ImageName} {string(Date(1,jjj))} {seconds(Time(1,jjj))} {Heihgt(1,jjj)} {Botom(1,jjj)} {Metr_fit_H(1,jjj)} {Metr_fit_B(1,jjj)} {str2double(SelectedDirectory{iii}(7:9))}], cell2mat(SelectedDirectory(iii)),['A' num2str(jjj+1) ':H' num2str(jjj+1)]); \r
-\r
-\r
- clear Figure\r
- clear ContrastFigure\r
- clear ContrastFigure1\r
- clear ContrastFigure2\r
- clear ContrastFigure3\r
- clear ContrastFigureFin\r
- % clear IntensitaAverFin\r
-\r
- end % jjj, Cycles for figures files\r
- \r
- hf4=figure ('Name','Velocity','NumberTitle','off','PaperSize',[8.5 6],'PaperUnits','centimeters'); \r
- hp1=plot(seconds(Time(startFile:step:end)), Metr_fit_H(startFile:step:nFFiles),'-or','linewidth',2); % [1 startFile:step:nFFiles]\r
- hold on\r
- hp2=plot(seconds(Time(startFile:step:end)), Metr_fit_B(startFile:step:nFFiles),'-ob','linewidth',2);\r
- hold on\r
- hp3=plot(seconds(Time(startFile:step:end)),LevelNum(1,startFile:step:nFFiles)','-og','linewidth',2);\r
- \r
- xlabel('time [s]','FontWeight','bold');\r
- ylabel('position [cm]','FontWeight','bold');\r
- grid('on');\r
- [legend_h,object_h,plot_h,text_strings] = legend([hp1,hp2,hp3],'Heihgt', 'Botom', 'Level','Location','southeast' );\r
- set(legend_h, 'Box', 'off', 'FontSize',7)\r
+ figure\r
+ imshow(ContrastFigure)\r
+ hold on\r
+ plot(IntensitaAverFin{jjj},'-y','linewidth',2)\r
+ hold on\r
+ if jjj>=startFile+step %subtracts the intensities\r
+ plot(IntensitaAverFin{jjj-step},'Color',[0.9290 0.6940 0.1250],'linewidth',2)\r
+ hold on\r
+ IntensitaAverFin_2{jjj}=IntensitaAverFin{jjj}-IntensitaAverFin{jjj-step};\r
+ waytrend=mean(IntensitaAverFin_2{jjj},2); % waytrend will tell whether to calculate the minimum or maximum green intensity IntensitaAverFin_2{jjj}\r
+ if waytrend<=0;\r
+ plot(-IntensitaAverFin_2{jjj},'-g','linewidth',2)\r
+ else\r
+ plot(IntensitaAverFin_2{jjj},'-g','linewidth',2)\r
+ end\r
+ hold on\r
+ end\r
+\r
+ infmt = 'hh:mm:ss.SS';\r
+ Time(jjj)=duration([(FFiles(jjj).name(8:9)) ':' (FFiles(jjj).name(10:11)) ':' (FFiles(jjj).name(12:16))],'Format',infmt)-duration([(FFiles(1).name(8:9)) ':' (FFiles(1).name(10:11)) ':' (FFiles(1).name(12:16))],'Format',infmt);\r
+ infmd = ['yy' 'MM' 'dd'];\r
+ Date(jjj) = datetime(FFiles(jjj).name(1:6),'InputFormat',infmd);\r
+ Level(jjj,1:3)=SelectedDirectory{iii}(7:9);\r
+ LevelNum(1,jjj)=str2double(Level(jjj,:))*1000;\r
+\r
+ % % plot(lk,abs(pk),'oc','linewidth',2)\r
+ if D_num==1\r
+ LevelPixels=MeterScale(1+(1000-LevelNum(1,jjj))/100); % x level in pixels in the original image\r
+ [pk,lk] = findpeaks(-IntensitaAverFin{jjj}(1:(LevelPixels+30))); %10\r
+ [Heihgt(1,jjj),pkH]=max(pk);\r
+ qH=lk(pkH);\r
+ Metr_fit_H(1,jjj)=polyval(polypar_Metr,qH);\r
+ plot(qH,-Heihgt(jjj),'or','linewidth',2)\r
+ if jjj>=startFile+step\r
+ LevelPixels=MeterScale(1+(1000-LevelNum(1,jjj))/100); % x level in pixels in the original image\r
+ if waytrend<=0;\r
+ [pk,lk] = findpeaks(-IntensitaAverFin_2{jjj}(1:(LevelPixels+30))); %10\r
+ else\r
+ [pk,lk] = findpeaks(IntensitaAverFin_2{jjj}(1:(LevelPixels+30))); %10\r
+ end\r
+ [Heihgt(1,jjj),pkH]=max(pk);\r
+ if max(pk)>=10; % error cleaning %30\r
+ qH=lk(pkH);\r
+ Metr_fit_H(1,jjj)=polyval(polypar_Metr,qH);\r
+ end\r
+ plot(qH,Heihgt(jjj),'om','linewidth',2)\r
+ end\r
+\r
+ elseif D_num==2\r
+ if jjj<startFile+step\r
+ LevelPixels=MeterScale(1); % x level in pixels in the original image\r
+ [pk,lk] = findpeaks(-IntensitaAverFin{jjj}(LevelPixels:end)); %10\r
+ [Heihgt(1,jjj),pkH]=max(pk);\r
+ qH=LevelPixels+lk(pkH);\r
+ Metr_fit_H(1,jjj)=polyval(polypar_Metr,qH);\r
+ plot(qH,-Heihgt(jjj),'or','linewidth',2)\r
+ elseif jjj>=startFile+step\r
+ LevelPixels=MeterScale(1); % x level in pixels in the original image\r
+ if waytrend<=0;\r
+ [pk,lk] = findpeaks(-IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
+ else\r
+ [pk,lk] = findpeaks(IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
+ end\r
+ [Heihgt(1,jjj),pkH]=max(pk);\r
+ qH=LevelPixels+lk(pkH);\r
+ Metr_fit_H(1,jjj)=polyval(polypar_Metr,qH);\r
+ if max(pk)<30; % error cleaning %30\r
+ qH=NaN;\r
+ Metr_fit_H(1,jjj)=NaN;\r
+ end\r
+ plot(qH,Heihgt(jjj),'om','linewidth',2)\r
+ end\r
+ end\r
+\r
+ if D_num==1\r
+ if jjj<startFile+step\r
+ [pk2,lk2] = findpeaks(-IntensitaAverFin{jjj},'MinPeakProminence',5); % Small peaks at the bubble/liquid interface\r
+\r
+ [pk3,lk3] = findpeaks(-IntensitaAverFin{jjj},'MinPeakProminence',5,'MinPeakDistance',30); % Throw out small random peaks in the liquid, peaks far apart cca 1cm\r
+ [one,nqq]=size(lk3);\r
+ [one,npp]=size(lk2);\r
+\r
+ if npp==0\r
+ pk2=NaN;\r
+ lk2=NaN;\r
+ elseif nqq==0\r
+ pk3=NaN;\r
+ lk3=NaN;\r
+ end\r
+\r
+ for qqq=1:1:nqq\r
+ for ppp=1:1:npp\r
+ if lk3(qqq)==lk2(ppp)\r
+ pk2(ppp)=NaN;\r
+ lk2(ppp)=NaN;\r
+ end\r
+ end\r
+ end\r
+\r
+ % % plot(lk2,abs(pk2),'om','linewidth',2)\r
+ hold on\r
+\r
+ [Botom(1,jjj),pk2B]=max(abs(pk2)); % For 1-4 files\r
+ qB=lk2(pk2B); % For 1-4 files\r
+\r
+ Metr_fit_B(1,jjj)=polyval(polypar_Metr,qB);\r
+ if Metr_fit_B(1,jjj)>=(Metr_fit_H(1,jjj)-2) % Throw out the bubble/liquid interface if it is above the surface, 2cm 2D interface\r
+ Metr_fit_B(1,jjj)=NaN;\r
+ qB=NaN;\r
+ elseif Botom(1,jjj)<=(max(abs(pk))+min(abs(pk)))*0.62 % 0.65 is an estimate of how much intensity the bubbles need\r
+ Metr_fit_B(1,jjj)=NaN; % For 1-4 files\r
+ qB=NaN;\r
+ end\r
+ plot(qB,Botom(jjj),'ob','linewidth',2)\r
+ end\r
+\r
+ if jjj>=startFile+step\r
+ LevelPixels=MeterScale(1+(1000-LevelNum(1,jjj))/100); % x level in pixels in the original image\r
+ if waytrend<=0;\r
+ [pk2,lk2] = findpeaks(-IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
+ else\r
+ [pk2,lk2] = findpeaks(IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
+ end\r
+ [Botom(1,jjj),pk2B]=max(pk2);\r
+ qB=LevelPixels+lk2(pk2B);\r
+ Metr_fit_B(1,jjj)=polyval(polypar_Metr,qB);\r
+ if max(pk2)<=45; % error cleaning\r
+ qB=NaN;\r
+ Metr_fit_B(1,jjj)=NaN;\r
+ end\r
+ plot(qB,Botom(jjj),'oc','linewidth',2)\r
+ end\r
+\r
+ elseif D_num==2\r
+ if jjj<startFile+step\r
+ [pk2,lk2] = findpeaks(-IntensitaAverFin{jjj},'MinPeakProminence',5); % Small peaks at the bubble/liquid interface\r
+\r
+ [pk3,lk3] = findpeaks(-IntensitaAverFin{jjj},'MinPeakProminence',5,'MinPeakDistance',30); % Throw out small random peaks in the liquid, peaks far apart cca 1cm\r
+ [one,nqq]=size(lk3);\r
+ [one,npp]=size(lk2);\r
+\r
+ if npp==0\r
+ pk2=NaN;\r
+ lk2=NaN;\r
+ elseif nqq==0\r
+ pk3=NaN;\r
+ lk3=NaN;\r
+ end\r
+\r
+ for qqq=1:1:nqq\r
+ for ppp=1:1:npp\r
+ if lk3(qqq)==lk2(ppp)\r
+ pk2(ppp)=NaN;\r
+ lk2(ppp)=NaN;\r
+ end\r
+ end\r
+ end\r
+\r
+ % % plot(lk2,abs(pk2),'om','linewidth',2)\r
+ hold on\r
+\r
+ [Botom(1,jjj),pk2B]=max(abs(pk2)); % For 1-4 files\r
+ qB=lk2(pk2B); % For 1-4 files\r
+\r
+ % [Botom(1,jjj),pk2B]=max(abs(pk)); % For 5. files\r
+ % qB=lk(pk2B); % For 5. files\r
+\r
+ Metr_fit_B(1,jjj)=polyval(polypar_Metr,qB);\r
+ if Metr_fit_B(1,jjj)>=(Metr_fit_H(1,jjj)-2) % Throw out the bubble/liquid interface if it is above the surface, 2cm 2D interface\r
+ Metr_fit_B(1,jjj)=NaN;\r
+ qB=NaN;\r
+ elseif Botom(1,jjj)<=(max(abs(pk))+min(abs(pk)))*0.62 % 0.65 is an estimate of how much intensity the bubbles need\r
+ Metr_fit_B(1,jjj)=NaN; % For 1-4 files\r
+ qB=NaN;\r
+ end\r
+ plot(qB,Botom(jjj),'ob','linewidth',2)\r
+ end\r
+\r
+ if jjj>=startFile+step\r
+ LevelPixels=MeterScale(1); % x level in pixels in the original image\r
+ if waytrend<=0;\r
+ [pk2,lk2] = findpeaks(-IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
+ else\r
+ [pk2,lk2] = findpeaks(IntensitaAverFin_2{jjj}(LevelPixels:end)); %10\r
+ end\r
+ [Botom(1,jjj),pk2B]=max(pk2);\r
+ qB=LevelPixels+lk2(pk2B);\r
+ Metr_fit_B(1,jjj)=polyval(polypar_Metr,qB);\r
+ if max(pk2)<=20; % error cleaning\r
+ qB=NaN;\r
+ Metr_fit_B(1,jjj)=NaN;\r
+ end\r
+ plot(qB,Botom(jjj),'oc','linewidth',2)\r
+ end\r
+ end\r
+\r
+ xlswrite(ResFile,[{'ImageName'} {'Date'} {'Time [s]'} {'Intensita Heihgt [-]'} {'Intensita Botom [-]'} {'Heihgt [cm]'} {'Botom [cm]'} {'Level [cm]'}],cell2mat(SelectedDirectory(iii)),'A1:H1');\r
+ xlswrite(ResFile,[{ImageName} {string(Date(1,jjj))} {seconds(Time(1,jjj))} {Heihgt(1,jjj)} {Botom(1,jjj)} {Metr_fit_H(1,jjj)} {Metr_fit_B(1,jjj)} {str2double(SelectedDirectory{iii}(7:9))}], cell2mat(SelectedDirectory(iii)),['A' num2str(jjj+1) ':H' num2str(jjj+1)]);\r
+\r
+\r
+ clear Figure\r
+ clear ContrastFigure\r
+ clear ContrastFigure1\r
+ clear ContrastFigure2\r
+ clear ContrastFigure3\r
+ clear ContrastFigureFin\r
+ % clear IntensitaAverFin\r
+\r
+ end % jjj, Cycles for figures files\r
+\r
+ hf4=figure ('Name','Velocity','NumberTitle','off','PaperSize',[8.5 6],'PaperUnits','centimeters');\r
+ hp1=plot(seconds(Time(startFile:step:end)), Metr_fit_H(startFile:step:nFFiles),'-or','linewidth',2); % [1 startFile:step:nFFiles]\r
+ hold on\r
+ hp2=plot(seconds(Time(startFile:step:end)), Metr_fit_B(startFile:step:nFFiles),'-ob','linewidth',2);\r
+ hold on\r
+ hp3=plot(seconds(Time(startFile:step:end)),LevelNum(1,startFile:step:nFFiles)','-og','linewidth',2);\r
+\r
+ xlabel('time [s]','FontWeight','bold');\r
+ ylabel('position [cm]','FontWeight','bold');\r
+ grid('on');\r
+ [legend_h,object_h,plot_h,text_strings] = legend([hp1,hp2,hp3],'Heihgt', 'Botom', 'Level','Location','southeast' );\r
+ set(legend_h, 'Box', 'off', 'FontSize',7)\r
\r
%--------------------------------------------------------------------------\r
- print ('-dpng', '-r1200', 'Velocity.png')\r
+ print ('-dpng', '-r1200', 'Velocity.png')\r
%--------------------------------------------------------------------------\r
\r
\r
- %%xlswritefig(hf4, ResFile, cell2mat(SelectedDirectory(iii)), 'J2')\r
+ %%xlswritefig(hf4, ResFile, cell2mat(SelectedDirectory(iii)), 'J2')\r
\r
end % iii, Cycles for directories with files\r
-\r
projects / Bubble_column.git / commitdiff