raws-reading-and-fusing.zip_光谱图像 matlab_多光谱融合_多波段融合_读取raw_高光谱

fid=fopen('x1.raw','r'); pixel=fread(fid,inf,'uchar'); tmp=pixel(1:SUM);%SUM为总像素数 m=reshape(tmp,sample,line);%sample为列数 line为行数 n1=(uint8(m))'; fid=fopen('x2.raw','r'); pixel=fread(fid,inf,'uchar'); tmp=pixel(1:SUM); m=reshape(tmp,sample,line); n2=(uint8(m))'; fid=fopen('x3.raw','r'); pixel=fread(fid,inf,'uchar'); tmp=pixel(1:SUM); m=reshape(tmp,sample,line); n3=(uint8(m))'; subplot(1,3,1),imshow(n1); subplot(1,3,2),imshow(n2); subplot(1,3,3),imshow(n3); figure,subplot(1,3,1) N=(n1+n2+n3)/3; imshow(N),title('三波段合成'); fid=fopen('Y.raw','r'); pixel=fread(fid,inf,'uchar'); tmp=pixel(1:SUM); m=reshape(tmp,sample,line); M=(uint8(m))'; % figure,subplot(1,3,1); imshow((N+M)/2),title('直接融合'); % F1=fft2(N); F2=fft2(M); N2=abs(ifft2(F1+F2)/2); %两幅图像在频域内相加后的傅立叶逆变换 subplot(1,3,2),image(N2),title('傅里叶变换后融合'); % [C1,L1]=wavedec2(N,2,'sym4'); [C2,L2]=wavedec2(M,2,'sym4');%分别对两幅图像进行小波分解 C=C1+C2;%融合分解系数 N3=waverec2(C,L1,'sym4');%对融合后的信号进行图像重构 subplot(1,3,3),imshow(N3/2),title('小波变换后融合'); %HIS变换 f=imresize(N,0.5); hsi=rgb2hsi(f); i=hsi(:,:,3);s=hsi(:,:,2);h=hsi(:,:,1);%它首先将RGB颜色空间的3个波段的多光谱图像转化到HIS空间的3个量，即色调(H),亮度(I)、饱和度(S)。亮度(I)代表空间信息，色调(H)代表光谱信息。 hspec=imhist(i); i2=histeq(Y,hspec);%将高空间分辨率图像进行对比度拉伸，使它和亮度分量(I)有相同的直方图 s1=imresize(s,2); h1=imresize(h,2); hsi1=cat(3,h1,s1,mat2gray(i2));%用拉伸后的高空间分辨率图像代替亮度分量(I)，把它同色调(H),饱和度(S)进行HIS逆变换得到融合图像 Z=hsi2rgb(hsi1); imshow(Z),title('HIS变换融合');