遥感图像处理图像变换

1、图像变换实验报告一.内容?快速傅里叶变换?主成分分析基本操作?代数运算?彩色变换二.目的在理解原理基础上，熟悉熟练图像的几种常用变换操作步骤三.实验过程傅里叶变换1.打开ENVI5.0,并导入CJ_spot遥感影像。图1:原始影像2.进行傅里叶变换操作川阜ClassificationTransformFifrerjSpectralMapVectorTopographiFilterDefinitionInverseFFT图2:图像的傅里叶变换kitiratinnTnrui+nrnnFillerSpprtralMspVprtarTopngraphi图3:设定滤波器图4:选择阻断滤波器图5:自定义多

2、边形滤波器图6:绘制自定义滤波器后的图像kClassificationTranstorm|filterSpectralMapVectorIopograph图7:逆变换ConvuljiiwiiidtidMorpholDy/iextureAdaptiveFlFTFilteringForiAdidFFTFitpDpfiniTionInverseFFT图8:波段列表窗口注解:HV FFT 仃盯(Bitid l：eK2Ll-piu :(TJ pod ; Il1D fr LB enonZ 10 Di ft 氏n&d Cit Filter耳 13 Mnoryl 1!fT (Krsi ze CEkirui 1

3、tl)CJ spg!is 13-官定.匕w;皮目m骷像Tse率域匡健一南始国量图9:FFT变换后影像图10:变换前后影像对比拓展练习11:滤波器12:滤波器13:滤波器#1FiIterDefinition图14:对于滤波器3修改边界像素值。利用上述方法进行FFT变换，并使RG为原始图像的数据。B分别为flf4,产生四个不同的彩色图像。主成分分析图15:主成分变换菜单图16:王成分刖向变换参数2EigenvalueNlitiber672COv13口2345e7SandNjmber昵1=1StatT一/Flrn=niLt；.D山_ok，最后始枭，夷勉数明y.alt3百Dins:FullScene(

4、3bDLUODciats)特证值Sand2Band3Band4Bandn6-31cO1O11221X1日口二i匚StotaBand1Band2Band3Daud4Band5日ar.d白Daud7Ssm-O957-53002545M2122212Mcaij8?20788133054647439032EQ47327G157633HB12372S3-1225640100abdcv647170149608959.42DOOGLI.02DG13ZZ0?B210333999013C04643Eigenvalue7m417731G5205934S7119429.2730176.15D9D1471231OO；

5、6G2CovarianceBand1in图18:主成分统计结果图19:选择特征值子集QSelectOuipuiPCBarick图20:特征值及累计贡献率TrnsfcrmFiterSpectralMapVectorTopographicRadarWineIiiidqeSlidipniriLjBandRatios|FirKip日ICompQnengForwardPCRotdtioniru/er三己PCRoimtioFIndependentComponents图21:主成分逆变换菜单图22:主成分变换前图25:代数运算对话框图23:主成分变换后代数运算1）整体增强图像的亮度|日444T口Transf

6、ormFiResizeData(Spotlal/Spectral)SubsetDataviaROIsRotarc/FlipDstaLayerStackingConvertDsta(RQ.BTI,EIP)StretchDdt靠StatisticSpdtialStatisticsChangeDelectionMedsureinerit1oolIpjndMaihIspectralMath图24 :代数运算选择界面SegmentationImage图26:代数运算后图像2）利用ENVI功能计算NDVI图27:NDWI生成界面WindowProfilesPolarization SignaturesRe

7、gion OfimerestENVI Color 1 ablestDensity Wlice.JControl RGB Image Planes.Class Color Mapping.图28 :对生心Ma即ingCursorLocatiorVValue.PixelLocatorPointCollection.成图像进行密度分割29:分级数为3图30:密度分割后遥感影像图32:不同格3）利用代数表达式计算NDVI。Bdrid Mdtli ORC-nEHein1图31:输入NDVI表达式日SMemcryS*nEandlathCtflat(b3)/(b4)+float(h8)iQMemoryU口B

8、andKath(float(b4-b3)ffioai(b4+b3)臼SMemory3.Hand任百th(b4-bG)/(M+b2)日爵AAP口Band1nEanl2口Rsnrl3*Tir卜式下的EDVI公式图24:不同表达式ENVI后遥感影像比较比较这三种表达式下的计算结果，会发现（b3-b4）/（b3+b4结果较好彩色变换米图25:将多光谱图像重采样为15图26:选择RGB俞入为543图27:选择全色波段为高分辨率数据将（5,4,3）进行RGB&成显示，窗口为#1,并对scroll窗口进行图像的均衡化。然后，点击上述菜单，选择#1作为输入图30:使用图像数据作为输入.图31:均衡化拉伸后输出

9、图32:彩色正变换选择窗口RGB to HLSRGB to HISV (USGS Munsell)彩色正变换IlnrIFlt-=rpArtvlFflnDC9”aphRairV/1n1nHfIiTidgeSlidiptminqBdndRdtiuiPrincipalComponentsIndependentCompon白ntwMNFRotationlLcdq1qrforms1DecorrelationStretchPhotographicStretch图33:彩色正变换前图34:彩色正变换后色逆变换选择窗口彩色逆变换mjyc Sharpcf lingBandPrincipal CrmponntsIn de pendent ComponentsMNRoQtion| Cui4T u匕f口中in* |Vector Topographic Rada r Window H图35:彩DecorrelationStretchPhotographicStretchSatiJratii-onStretchSyntJwticColorIiridgeNDV1Tass