微波遥感实验报告

1、实习报告撰写的内容与要求实习任务：介绍实习的目的、意义、任务及实习单位的概况等内容。通常以前言或引言形式表述，不单列标题及序号。实习内容：先介绍实习安排概况，包括时间、地点、内容等，然后逐项介绍具体实习流程与实习工作内容，以及专业知识与专业技能在实习过程中的应用。本部分内容应以记叙或白描手法为基调，在完整叙述的基础上，对自己认为有重要意义或需要研究解决的问题进行重点叙述，其它内容则可简述。实习结果：围绕实习任务要求，对实习中发现的问题进行分析、思考，提出解决问题的对策、建议等。分析问题、解决问题要有依据（如有参考文献可在正文后附录）。分析讨论的内容、推理过程及所提出的对策与建议作为实习报告的重

2、要内容之一，是反映或评价实习报告水平的重要依据。实习总结或体会：对实习效果进行综合评价，着重介绍自身的收获与体会，内容较多时可列出小标题，逐一列举。总结或体会的最后部分，应针对实习中发现的自身不足，简要地提出今后学习，努力的方向。将实习日记按照时间顺序以附件形式放在实习报告正文后面。实习报告封皮由学校统一印发，正文一律采用计算机排版、A4纸打印。题目为三号黑体字居中（题目前、后各空一行），正文字体为小四号宋体，要求语句通顺、论述严谨、规范、正确。字数：不少于3000。目录单雷达影像处理TOC o 1-5 h z HYPERLINK l bookmark10 o Current Document

3、 导入数据3 HYPERLINK l bookmark12 o Current Document 影像多视处理4 HYPERLINK l bookmark14 o Current Document 滤波5 HYPERLINK l bookmark16 o Current Document 分析滤波影像7 HYPERLINK l bookmark18 o Current Document 1.5地理编码和辐射定标9 HYPERLINK l bookmark20 o Current Document 配准10InSAR生成DEM11基线估算11 HYPERLINK l bookmark22 o C

4、urrent Document 干涉图生成12 HYPERLINK l bookmark24 o Current Document 去平处理13 HYPERLINK l bookmark26 o Current Document 自适应滤波及相干性计算14 HYPERLINK l bookmark30 o Current Document 相位解缠16 HYPERLINK l bookmark32 o Current Document 选择GCP17 HYPERLINK l bookmark34 o Current Document 轨道精炼和重去平18 HYPERLINK l bookmar

5、k36 o Current Document 相位高程转换19 HYPERLINK l bookmark38 o Current Document 思考题21微波遥感实验报告1.单雷达影像处理1.1导入数据导入武汉市的EnvisatASAR数据，数据文件名：ASA_IMS_lPNDPA20081221_141624_000000162074_00483_35607_5531.Nl操作过程：SARscape-Basic-ImportData-StandardFormatsDatatype选择IMS,然后选择数据的输入和输出路径，路径必须都是英文的。数据的导出，选择JPEG格式，选择输出路径，路径

6、中那个不能有中文。导入的图像显示结果如下：（原图太大，插入之后word没法运行，只能截图）1.2影像多视处理多视的基本原理：在方位向或距离向降低处理器带宽，从而将方位向或距离向的频谱分割成若干部分，称为视窗（look），各视窗分别成像后非相干叠加，以抑制噪声。视窗的个数就称之为视数，它对应于前述的独立样本个数。INSAR原始数据一般是单视数的，即未作过多视处理的数据，多视处理往往在形成干涉图之后。多视处理是在图像的距离向和方位向上的分辨率做了平均，目的是为了抑制SAR图像的斑点噪声。多视的图像提高了辐射分辨率，降低了空间分辨率。选择多视比例系数为5：1（方位向5，距离向1）。操作过程：SARs

7、capeBasicMultilooking5*MultilacikinQJnpuLl11ikLD：*i.E0J&10LOMHFljTSZD&!2c141D:Z01010LOMliHLLLookinc*i.ZOJ010Ltvrl|Ou.tpks|Box-3.rK.v.3lxq-廿匚吕上弋百匸w亡HiltchIjhklD二ILdoIck影像多视化处理的结果如下:1、罰*1Bondl;2Q16101_pwrfilrOrrlnvEnhonrrGQbWindow尊血Ecrallill口1.3滤波操作过程：SARscape-Basic-Filtering-Single(1)Frost滤波Frost滤波的系

8、数分别设置为9*9和21*21,进行了滤波处理(2)Lee滤波:我将Lee滤波的参数设置为9*9,和Frost滤波进行对比:Frost滤波Lee滤波由上面两幅图下个对比可得:经过Lee滤波的图像相比于经过Frost滤波的图像，整体上稍微变亮了，更为平滑，噪声去除较为明显。1.4分析滤波影像问题：用滤波后的影像，分别寻找存在体散射、镜面反射、透视收缩、阴影的地物目标，并对各种现象进行描述。镜面散射：平整的表面更容易在SAR图像上表现为较暗的区域；如图所示即为由水体造成的镜面反射在雷达图像上造成的暗面。体散射：粗糙的表面时常能够得到更高的后向散射。透视收缩：雷达距山底的距离小于距山顶的距离，所以雷

9、达波束先到山的底部，再到山的顶部，成像也是。假设山坡的长度为L,其斜距显示的距离为Lr,很明显，LrBasicGeocodingGeocodingandRadiometrieCalibration。地理编码和辐射定标之后的结果如下:1.6配准覆盖同一地区的多幅雷达影像，如要进行时间序列分析、动态监测、多时相滤波处理等，需要进行图像间的配准处理。配准处理要求是斜距几何，并且各个图像采用相同的接收几何。配准不同于地理编码，地理编码是将每个像素从斜距几何转化为地图投影。操作过程：SARscape-Basic-Coregistration。II1=12.InSAR生成DEM对于要进行干涉的两幅影像，要

10、现将其导入到ENVI中，再进行基线估算和干涉图的生成。2.1基线估算打开基线估算面板：SARscape-Interferometry-BaselineEstimation配准之后的结果如下:3幫1Stroll“(31回Q#1ZcsornE.=j|Bl并圭FiIhEnh-ancaTcwJ.WineJow*1Baidl.A5A_aODG_pwi_rspUSVIilGCF11*ShLftPwiritteirfileFzzQK4noi4!jLm4ntD4.s411ixaEstSrsii1A5A_C005_pvrInputfileListOutputf11ZLItCcoftc：LSitritlan.fL

11、ieGOf4j4CtJDCREE丽FLJhGHeiftstHeichbourBiLinearDEWVZEOIILS2LI3JkAjaQDfijwrT-xj.|ElOFVkLcXICBTiC-kIjtputwMmu*IIT7iisr琳zaK4nc!i4!jnjYlpwrLm4iilIstSr-ii1on.SA5A_COJ5_cvr3rdCublcCaitvoLuticiji0thCotlcCciJioLu.tic!n|gAZLivayvaanput*ch.3f-t?Caapnt*SkiftFu*HTiA-ter-E尸Canpn“1ilflPu-ftiridtreOtJ.y2.2干涉图生成生成干涉

12、图：输入两景SLC数据，输出数据是经过配准和多视的两景数据的干涉图，和主从影像的强度图。对该数据集，距离向为1方位向为4的多视，约可得到15米的地面分辨率。操作过程：SARscape-lnterferometry-lnterferogramgeneration-WithOutDEM打开InterferometryGenerationwithoutDEM面板。完成干涉后的图像结果显示如下:2.3去平处理操作过程：SARscape-lnterferometry-lnterferogramgeneration-1nterferogramFlattering丿DEM/EllipsoidsFlatten

13、ing回|&去平之后的图像处理如下:防寸1处皿II曲U2.4自适应滤波及相干性计算对生成的干涉图（_dint）进行滤波，去掉由平地干涉引起的位相噪声。同时生成干涉的相干图（描述位相质量）和滤波后的主影像强度图。操作过程：SARscape-Interferometry-AdaptiveFilterandCoherenceGeneration-BoxcarWindow，打开AdaptiveFilterandCoherenceGenerationBoxcarWindow面板：AdaptiveFilterandCdKeren匚已GeoerationiGoldsteinTntarferogrsmfiLe

14、D:它N任5口LEELi)Ihock.i1jnamm_liMasterfile：311crijwiveEmnw5she2ngEperimentaselineEstinatiirL，iA5A_jZIDD5_p,i,rSlavefileE:!i11eri?waveEernoteSens3ngEKperimentIaaelineEstina.11an.ASAUUDGpyrOutputRdOtWalTiC:D:VEm50IILS2VlibiCi!1lc,Lilsnam。CokerenceEgWiaduwSize3CahererLcekzWindowSii-Filteringsetting扌Cohere

15、ncegeneralioilJAdaptivefilterVCoherencefromFintSt：=LftStoreBatchICarLcelHelp自适应滤波之后的图像为:*1BandFileEnhanceTc-dIiWindDv加载相干性图，打开CursorLocation/Value查看像元值，相干性系数分布在0-1,值越大说明该区域的相干性越高，值越小，表明该区域在两个时相上发生了变化。0#1Bandtrlnle_name_cc|匕叫回|坯|RleOvertayEnhance-TodHWindow亮色区域相干系数白色区域的相干系数黑色区域的相干系数在InSAR数据处理中，相关系数是衡

16、量图像对的重要指标，它可以用来选择数据对，也可用来评价图像是否已经配准。由上面两幅图对比可得，在较黑的区域相干性趋近于0,在较亮的区域相干性趋近于1。这是因为在黑色区域，雷达回波信号弱，信噪比非常低，故相干系数较小，而在亮色区域却恰好相反。2.5相位解缠干涉相位只能以2n为模，所以只要相位变化超过了2n，就会重新开始和循环。相位解缠是对去平和滤波后的位相进行相位解缠，解决2n模糊的问题。操作过程：SARscape-Interferometry-PhaseUnwrapping，打开PhaseUnwrapping面板相位解缠后的图像如下:2.6选择GCP加载并显示ASA_0006_pwr图像，利用

17、ENVI的矢量工具选择控制点，菜单VectorCreateNewVectorLayerUseRasterImageFile，选择ASA_0006_pwr文件作为选择控制点的基准文件，单击OK,输入矢量图层的名字以及输出路径。5el4CtInputTile:lilw_nQ_Lnl11il14二nsuph中amdm!jw_rapE*L1iDLrormtl：TifLIaC：Sf1Lcr-3iB!iiV434nE；LiiGMj4mDims-EJTEk67KkJJSG5icb:PLc-b.13dlFPjjl门L?9,30Dbyt.FileTnAddNewVectors;在VectorParameters

18、面板，选择ModePoint为NewVectorsType;然后进行选点操作；VectorParameterE:AddNe.|=|El|右图中蓝色框中的红色点即为已经选择采集好的控制点下面是将矢量控制点导入SARscape，操作过程：SARscapeToolsGenerateGroundControlPointFile百tart|StDiTQBztch|Cajim.HeId2.7轨道精炼和重去平操作过程：SARscape-Interferometry-RefinementandReflattening,打开轨道精炼和重去平面板RjefinementandR&-flatfcening|SMast

19、erfileE:kMi：i-:-wo.veHriot亡EerLSin.gl,.EMp亡iti.neiLt1-.!neEsiini-ititlASA_lJJU5uwrSlavefileE:Vilicrowavelen口t已5eninglExperin已丁lasel】neEetinatichlA5AJZ叮CIbjwfIrftsrstgi-irifils：t-.T=：rTrirj：?.iii.：-.rr.Jrfii=祢=_*卅SyrLtketicPhaseI:EHVI5O1IIL92liT?hDokfilc_nanic_3intUnwrippedFhasa:EmiSO1.JDLaZ.h1.hookf

20、ile_n：urieuphaCohet-cncefileI:：DLE21ilAcoAil；_c:51：aritRazgeDEMfileI.lEffVnCCIiLtZililAduklirilum：!：diOrbitalGCFfileE:iMicrowavelen已3ensingVEKperin已丁LtVvector66OTit.zmlBENE辽eDu切utRoottfaiTieE:VilicrowaveflenDt已5engingAEKperine丁ltchoragqapi口glfi1e_najie_reflatRinementSettmg*Re1attuningCartographicSys

21、tenGEO-GLOBALGEOWGSB4重去平后的输出结果_reflat_dint和_reflat_upha加载到波段列表中:dR-EfinementandRe-flatfeEningOrbitalRi血亡mattLttlKg：OrBitaishiftinKdirectaqnn)OrbitsshiftinYdirectioilih)OrbitcLLshiftleiZd.Lract3:-tl(n)A-j.ric.riachievableMISIn=1133132952=-D.4353273077=-1.5332335236-0.4E0ILIB115lependencyoftheshiEtisi

22、Idirsctifrcnth*DependeuczoftheshlftinTdiiectiomfrcathe1ApaiiRDEN&nd5JLR.DENjnjazi.nutJiposiisoitai1Tiatlpositioiliil口thpocitionIE胡TH眈銀ti/pisd)lii/pLKel)(njpi=0.00=6327787=a.0D34,9i：ii：i431=0.flEZ365B37=B1T7L5BSE7Z重去平之后的图像为:2.8相位高程转换这一步是将经过绝对校准和解缠的实际位相，结合合成位相，转换为DEM并进行地理编码。生成进行了地理编码的DEM文件、相干图像，还有在数据镶嵌中会用到的两个图像