ENVI对SAR数据的预处理过程(详细版).docx

一、数据的导入：(1) 在Toolbox中，选择SARscape->Basic->ImportData->StandardFormats->ALOSPALSAR⑵在打开的面板中，数据类型(DataType)：JAXA-FBDLevel1.1注：这些信息可以从数据文件名中推导而来⑶单击Leader/Paramfile，选择d1300816-005-ALPSRP246750820-H1.1_A\LED-ALPSRP246750820-H1.1__A文件⑷点击Datalist，选择d1300816-005-ALPSRP246750820-H1.1_A\IMG-HH-ALPSRP246750820-H1.1__A文件⑷单击Outputfile，选择输出路径注：软件会在输入文件名的基础上增加几个标识字母，如这里增加“_SLC”⑸单击Start执行，最后输出结果是ENVI的slc文件，sml格式的元数据文件，hdr格式的头文件等⑹可在ENVI中打开导入生成的以slc为后缀的SAR图像文件二、多视单视复数(SLC)SAR图像产品包含很多的斑点噪声，为了得到最高空间分辨率的SAR图像，SAR信号处理器使用完整的合成孔径和所有的信号数据。

多视处理是在图像的距离向和方位向上的分辨率做了平均，目的是为了抑制SAR图像的斑点噪声多视的图像提高了辐射分辨率，降低了空间分辨率1) 在Toolbox中，选择SARscape->Basic->Multilooking2) 单击Inputfile按钮，选择一景SLC数据(前面导入生成的ALOSPALSAR数据)注意：文件选择框的文件类型默认是*_slc，就是文件名以_slc结尾的文件，如不是，可选择*.*⑶设置：方位向视数(AzimuthLooks)：5,距离向视数(RangeLooks)：1注：详细的计算方法如下所述另外，单击Look按钮可以估算视数4) BorderResize选项，选择此项，会对检测结果边缘中的无效值，进而重新计算输出图像的大小这里不选择5) 输出路径会依据软件默认参数设置自动添加或自行修改，单击Start按钮执行回word格式-可编辑-感谢下载支持(6)计算完之后在Display中显示结果，可以看到图像的斑点噪声得到的抑制，但是降低了空间分辨率◎Multiookng注：方位向视数(AzimuthLooks)；距离向视数(RangeLooks)视数跟距离向分辨率、方位向分辨率以及中心入射角有关。

精确的计算方法如下地面分辨率=pixelspacingslantrange/sin(incidenceangle)• pixelspacingazimuth=3.169m(方向分辨率)• pixelspacingslantrange=9.368m(距离分辨率)• incidenceanglescenecentre=34.3°•地距的距离分辨率：9.368/sin(34.3°)=16.62m,距离向视数为1• 地距的方位向分辨率经过多视后保持与地距的距离分辨率一致• 16.62/3.169=5,因此方位向视数为5三、滤波从连贯SAR传感器中获取的图像都有斑点噪声,SARscape提供两大类滤波,用于单波段雷达图像的滤波和多时相雷达图像滤波1) Toolbox中,选择SARscape->Basic->Filtering->SingleImage2) 打开FilteringSingleImage面板,单击Inputfile按钮,选择前面得到的多视处理结果⑶有8中滤波供选择，选Frost,同时会打开FrostLee/RefinedLee对话框4) 在FrostLee/RefinedLee对话框中,分别设置AzimuthWindowSize和RangeWindowSize：9(5) 单击Start执行。

I回P冨StoreBa.1匚hCancelcliffuaianMean-Median.閘口日已EPSFrostInjutfilelistE:5C-HH-ALE噫RP24&T£[iE总D-HL.L_k_pvrOutputfilelistE:MMG-HH-ALFSRFZ4B75QBEa-HLL_l_pvr_f:l440FilteringSingle-Image四、地理编码和辐射定标SAR系统是测量发射和返回脉冲的功率比，这个比值(就是后向散射)被投影为斜距几何由于不同SAR传感器或者不同接收模式，为了更好的对比SAR图像几何和辐射特征，需要将SAR数据从斜距或地距投影转换为地理坐标投影(制图参考系)1)Toolbox中，选择SARscape->Basic->Geocoding->GeocodingandRadiometricCalibration2) 在GeocodingandRadiometricCalibration面板中，单击Inputfile按钮，选择前面Frost滤波做的结果⑶CartographicSystem,设置输出投影参数⑸像元大小(GRIDSIZE)：x：25,y：25。

⑹重采样方法(RESAMPLE)：有5中方法供选择，从左到右精度提高，但是速度越慢这里选择OptimalResolution⑺选择辐射定标(RadiometricCalibration)(8) 其他可选项：辐射归一化(Radometric)，局部入射角校正(LocalIncidenceAngle)、叠掩/阴影处理(Layover/Shadow)(9) 单击Start执行五、镶嵌工具在Toolbox中启动/Mosaicking/SeamlessMosaic1)点击左上的加号(如图)添加需要镶嵌的影像数据word格式-可编辑-感谢下载支持JSeamle5-sMosaicSeamlessMosaicNogaicScenesIntaASingleRa5terFile£election嗨h嵌PreviewtInputFlL*(e)7LS14W30OK2C»a£>73jIITLJIFileInftirmaliDiLFiibishL51450曲_p曲2QOT&9EL_|TLjWulti£R1314J023200B(IT3L11Sp，CilorCwrrSwimAll.dtiorLS1嗣030JJ磁期前3LJ1TLJI口ti百計.E45O3Oi_O3O2(KS0T22JTL_TLermEEL5145OES_0Z3Z00i391L_llTL_IhsrmsiCmttl9也EkSif卜比勺野FFIr滋"亏田冲4rr*議:醫兰理4理社疋證00址口1■丄TLE1叩旺dLo[jL511303Qi_03a2()OT(M13_MTL_Ihrn(2) 添加进来之后可以看到数据的位置和重叠关系和影像轮廓线。

3) 勾选ShowPreview，可以预览镶嵌效果在DataIgnoreValue一览输入透明值，这里输入0如下图为0值透明的效果匀色处理：提供的匀色方法是直方图匹配1)在CorlorCorrection选项中，勾选HistogramMatching：• "OverlapAreaOnly：重叠区直方图匹配• -EntireScene：整景影像直方图匹配(2)在Main选项中，ColorMatchingAction中右键设置参考(Reference)和校正(Adjust)，根据预览效果确定参考图像接边线与羽化接边线包括自动和手动绘制两种方法，也可以结合起来使用1)选择下拉菜单Seamlines->AutoGenerateSeamlines，自动绘制接边线，如下图所示，自动裁剪掉TM边缘“锯齿”2) 自动生成的接边线比较规整，可以明显看到由于颜色不同而显露的接边线下拉菜单Seamlines->Starteditingseamlines，编辑接边线，可以在接边处绘制多边形，之后自动将绘制的多边形作为新的接边线输出结果切换到Export选项，这是输出文件名、路径、格式、波段、背景值、重采样方法等信息。

eviehSeamlessMosaicItasai.cScentsI并忆ASingle盼馱就□rd*rMa.in|Gi&lerConrS•飆11ib»fFt：.thOpen,将SPOT4数据bldr_sp.img和LandsatTM数据TM-30m.img分别打开2) 在Toolbox中，打开/ImageSharpening/Gram-SchmidtPanSharpening在文件选择框中分别选择bldr_sp.img作为低分辨率影像(LowSpatial)和bldr_sp.img作为高分辨率影像(HighSpatial),单击OK。

打开PanSharpeningParameters面板3) 在PanSharpeningParameters面板中，选择传感器类型(Sensor)：Unknown，重采样方法(Resampling)：CubicConvolution，输出格式为：ENVI4) 选择输出路径及文件名，单击OK执行融合处理注：进度条显示在右下角相同传感器图像融合对于高分辨率影像，同样可以Gram-SchmidtPansharping融合方法达到更好的效果，下面以QuickBird影像为例介绍这种融合方法，数据存放在"…〉8.图像融合\数据\Quikbird"中ENVI5.1中,pansharping融合方法就是Gram-SchmidtPanSharpening选项1) File>Open，打开影像文件qb_boulder_msi.img和qb_boulder_pan.img2) 在Toolbox中，打开/ImageSharpening/Gram-SchmidtPanSharpening在文件选择框中分别选择qb_boulder_msi.img作为低分辨率影像(LowSpatial)和qb_boulder_pan.img作为高分辨率影像(HighSpatial)，单击OK。

打开PanSharpeningParameters面板3) 在PanSharpeningParameters面板中，选择传感器类型(Sensor)：QuickBird，重采样方法(Resampling)：CubicConvolution，输出格式为：ENVI注：传感器类型(Sensor)中还包括：GeoEyeT、Goktruk-2、IKONOS、landsat8_oli、landsat8_tirs、NPPVIIRS、PleiadesTA/B、QuickBird、UI:GSS：Sensorrasat、Spot-6、。