資源描述:
《遙感數(shù)字圖像處理實(shí)習(xí)報(bào)告含Matlab處理代碼.doc》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫。
1、遼寧工程技術(shù)大學(xué)《數(shù)字圖像處理》上機(jī)實(shí)習(xí)報(bào)告教學(xué)單位遼寧工程技術(shù)大學(xué)專業(yè)攝影測量與遙感實(shí)習(xí)名稱遙感數(shù)字圖像處理班級測繪研11-3班學(xué)生姓名路聚峰學(xué)號(hào)指導(dǎo)教師孫華生實(shí)習(xí)1讀取BIP、BIL、BSQ文件一、實(shí)驗(yàn)?zāi)康挠肕atlab讀取BIP、BIL、BSQ文件,并將結(jié)果顯示出來。遙感圖像包括多個(gè)波段,有多種存儲(chǔ)格式,但基本的通用格式有3種,即BSQ、BIL和BIP格式。通過這三種格式,遙感圖像處理系統(tǒng)可以對不同傳感器獲取的圖像數(shù)據(jù)進(jìn)行轉(zhuǎn)換。BSQ是像素按波段順序依次排列的數(shù)據(jù)格式。BIL格式中,像素先以行為單位塊,在每個(gè)塊內(nèi),按照波段順序排列像素。BIP格式中,以像素為核心,
2、像素的各個(gè)波段數(shù)據(jù)保存在一起,打破了像素空間位置的連續(xù)性。用Matlab讀取各個(gè)格式的遙感數(shù)據(jù),是圖像處理的前提條件,只有將圖像讀入Matlab工作空間,才能進(jìn)行后續(xù)的圖像處理工作。二、算法描述1.調(diào)用fopen函數(shù)用指定的方式打開文件。2.在for循環(huán)中調(diào)用fread函數(shù),用指定的格式讀取各個(gè)像素。3.用reshape函數(shù),重置圖像的行數(shù)列數(shù)。4.用imadjust函數(shù)調(diào)整像素的范圍,使其有一定對比度。5.用imshow顯示讀取的圖像。三、Matlab源代碼1.讀取BSQ的源代碼:clearallclclines=400;samples=640;N=6;img=fop
3、en('D:sample_BSQ','rb');fori=1:Nbi=fread(img,lines*samples,'uint8');band_cov=reshape(bi,samples,lines);band_cov2=band_cov';band_uint8=uint8(band_cov2);tif=imadjust(band_uint8);mkdir('D:MATLAB','tifbands1')name=['D:MATLABtifbands1tif',int2str(i),'.tif'];imwrite(tif,name,'tif');tilt=[
4、'波段',int2str(i)];subplot(3,2,i),imshow(tif);title(tilt);endfclose(img);2.讀取BIP源代碼clearallclclines=400;samples=640;N=6;fori=1:Nimg=fopen('D:MATLABsample_BIP','rb');b0=fread(img,i-1,'uint8');b=fread(img,lines*samples,'uint8',(N-1));band_cov=reshape(b,samples,lines);band_cov2=band_cov';%×
5、a??band_uint8=uint8(band_cov2);tif=imadjust(band_uint8);mkdir('E:MATLAB','tifbands')name=['E:MATLABtifbandstif',int2str(i),'.tif'];imwrite(tif,name,'tif');%imwrite(A,filename,fmt)tilt=['波段',int2str(i)];subplot(3,2,i),imshow(tif);title(tilt);fclose(img);end3.讀取BIL的源代碼clearallclclines=
6、400;samples=640;N=6;fori=1:Nbi=zeros(lines,samples);forj=1:samplesimg=fopen('D:MATLABsample_BIL','rb');bb=fread(img,(i-1)*640,'uint8');b0=fread(img,1*(j-1),'uint8');bandi_linej=fread(img,lines,'uint8',1*(N*samples-1));fclose(img);bi(:,j)=bandi_linej;endband_uint8=uint8(bi);tif=imadjust
7、(band_uint8);mkdir('D:MATLAB','tifbands')name=['D:MATLABtifbandstif',int2str(i),'.tif'];imwrite(tif,name,'tif');tilt=['2¨??',int2str(i)];subplot(3,2,i),imshow(tif);title(tilt);end四、運(yùn)行結(jié)果圖1:讀取文件的六個(gè)波段圖實(shí)習(xí)2均值/中值濾波、邊緣信息提取一、實(shí)驗(yàn)?zāi)康呐c原理各種圖像濾波算子可以實(shí)現(xiàn)圖像的增強(qiáng),去噪,邊緣提取等。圖像增強(qiáng)的目的在于:1.采用一系列