matlab.zip_eightvce_matlab_matlab photons_upper83t_生物图像

%三波长图像 %正周期傅里叶 figure; I_three=cell(1,1184); % for k=5:1150 for k=6:1189 %img_three=imread(['E:\李鹤\experiment\three-wavelength exper\data\three wavelength\',int2str(k-1),'.JPG']); %m1是一个公共变量会不停的被下一张图片信息覆盖最后是最后一张图片的信息 %img_three=imread(['E:\李鹤\experiment\three-wavelength exper\data2\three wavelength\',int2str(k-1),'.JPG']); %img_three=imread(['E:\李鹤\experiment\three-wavelength exper\data3\three wavelength\',int2str(k-1),'.JPG']); img_three=imread(['E:\李鹤\experiment\three-wavelength exper\data13\three wavelength\',int2str(k-1),'.JPG']); I_three{k}=img_three; %I{1}.....就是每一张图片的信息 end [row,col]=size(I_three{1}); b_fft=zeros(1,1184); for k=6:1189 b_fft(k)=sum(sum(I_three{k})); end bb_fft=b_fft(6:1189); y_b=fft(bb_fft,1184); %对信号进行快速Fourier变换 mag=abs(y_b); %求得Fourier变换后的振幅 plot(mag); axis([0 120 0 1.5*10.^9]); % 设置坐标轴在指定的区间 %title('三波长图像FFT') % %三波长图像 % figure; % I_three=cell(1,1200); % % for k=5:1150 % for k=1:1200 % %img_three=imread(['E:\李鹤\experiment\three-wavelength exper\data\three wavelength\',int2str(k-1),'.JPG']); %m1是一个公共变量会不停的被下一张图片信息覆盖最后是最后一张图片的信息 % %img_three=imread(['E:\李鹤\experiment\three-wavelength exper\data2\three wavelength\',int2str(k-1),'.JPG']); % %img_three=imread(['E:\李鹤\experiment\three-wavelength exper\data3\three wavelength\',int2str(k-1),'.JPG']); % img_three=imread(['E:\李鹤\experiment\three-wavelength exper\data13\three wavelength\',int2str(k-1),'.JPG']); % I_three{k}=img_three; %I{1}.....就是每一张图片的信息 % end % [row,col]=size(I_three{1}); % b_fft=zeros(1,1200); % for k=1:1200 % b_fft(k)=sum(sum(I_three{k})); % end % bb_fft=b_fft(1:1200); % y_b=fft(bb_fft,1200); %对信号进行快速Fourier变换 % mag=abs(y_b); %求得Fourier变换后的振幅 % plot(mag); % axis([0 120 0 1.5*10.^9]); % 设置坐标轴在指定的区间 % %title('三波长图像FFT') % %双波长图像-红外and橙 % I_io=cell(1,1000); % for k=1:1000 % %img_io=imread(['E:\李鹤\experiment\three-wavelength exper\data\orange&infrared\',int2str(k-1),'.JPG']); %m1是一个公共变量会不停的被下一张图片信息覆盖最后是最后一张图片的信息 % img_io=imread(['E:\李鹤\experiment\three-wavelength exper\data2\orange&infrared\',int2str(k-1),'.JPG']); % I_io{k}=img_io; %I{1}.....就是每一张图片的信息 % end % [row,col]=size(I_io{1}); % b_fft=zeros(1,1000); % for k=1:1000 % b_fft(k)=sum(sum(I_io{k})); % end % y_b=fft(b_fft,1000); %对信号进行快速Fourier变换 % mag=abs(y_b); %求得Fourier变换后的振幅 % %双波长图像-红外and紫 % I_ip=cell(1,1143); % for k=54:1196 % %img_ip=imread(['E:\李鹤\experiment\three-wavelength exper\data\infrared&purple\',int2str(k-1),'.JPG']); %m1是一个公共变量会不停的被下一张图片信息覆盖最后是最后一张图片的信息 % img_ip=imread(['E:\李鹤\experiment\three-wavelength exper\data8\infrared&purple\',int2str(k-1),'.JPG']); % I_ip{k}=img_ip; %I{1}.....就是每一张图片的信息 % end % [row,col]=size(I_ip{1}); % b_fft=zeros(1,1143); % for k=54:1196 % b_fft(k)=sum(sum(I_ip{k})); % end % bb_fft=b_fft(54:1196); % y_b=fft(bb_fft,1143); %对信号进行快速Fourier变换 % mag=abs(y_b); %求得Fourier变换后的振幅 % plot(mag); % %双波长图像-橙and紫 % I_op=cell(1,1000); % for k=1:1000 % img_op=imread(['E:\李鹤\experiment\three-wavelength exper\data2\purple&orange\',int2str(k-1),'.JPG']); %m1是一个公共变量会不停的被下一张图片信息覆盖最后是最后一张图片的信息 % I_op{k}=img_op; %I{1}.....就是每一张图片的信息 % end % [row,col]=size(I_op{1}); % b_fft=zeros(1,1000); % for k=1:1000 % b_fft(k)=sum(sum(I_op{k})); % end % y_b=fft(b_fft,1000); %对信号进行快速Fourier变换 % mag=abs(y_b); %求得Fourier变换后的振幅 % %单波长图像——红外 % I_i=cell(1,1200); % for k=1:1200 % %img_i=imread(['E:\李鹤\experiment\three-wavelength exper\data\infrared\',int2str(k-1),'.JPG']); %m1是一个公共变量会不停的被下一张图片信息覆盖最后是最后一张图片的信息 % img_i=imread(['E:\李鹤\experiment\three-wavelength exper\data2\infrared\',int2str(k-1),'.JPG']); % I_i{k}=img_i; %I{1}.....就是每一张图片的信息 % end % [row,col]=size(I_i{1}); % b_fft=zeros(1,1200); % for k=1:1200 % b_fft(k)=sum(sum(I_i{k})); % end % y_b=fft(b_fft,1200); %对信号进行快速Fourier变换 % mag=abs(y_b); %求得Fourier变换后的振幅 % %单波长图像——紫光 % I_p=cell(1,1200); % for k=1:1200 % %img_p=imread(['E:\李鹤\experiment\three-wavelength exper\data\purple\',int2str(k-1),'.JPG']); %m1是一个公共变量会不停的被下一张图片信息覆盖最后是最后一张图片的信息 % %img_p=imread(['E:\李鹤\experiment\three-wavelength exper\data2\purple\',int2str(k-1),'.JPG']); % img_p=imread(['E:\李鹤\experiment\three-wavelength exper\data3\purple\',int2str(k-1),'.JPG']); % I_p{k}=img_p; %I{1}.....就是每一张图片的信息 % end % [row,col]=size(I_p{1}); % b_fft=zeros(1,1200); % for k=1:1200 % b_fft(k)=sum(sum(I_p{k})); % end % y_b=fft(b_fft,1200); %对信号进行快速Fourier变换 % mag=abs(y_b); %求得Fourier变换后的振幅 % %单波长图像——橙光 % I_o=cell(1,1200); % for k=1:1200 % %img_o=imread(['E:\李鹤\experiment\three-wavelength exper\data\orange\',int2str(k-1),'.JPG']); %m1是一个公共变量会不停的被下一张图片信息覆盖最后是最后一张图片的信息 % img_o=imread(['E:\李鹤\experiment\three-wavelength exper\data2\orange\',int2str(k-1),'.JPG']); % % I_o{k}=img_o; %I{1}.....就是每一张图片的信息 % % end % [row,col]=size(I_o{1}); % b_fft=zeros(1,1200); % for k=1:1200 % b_fft(k)=sum(sum(I_o{k})); % end % y_b=fft(b_fft,1200); %对信号进行快速Fourier变换 % mag=abs(y_b); %求得Fourier变换后的振幅 %