1、2018/7/2,Beijing oriental titan technology co., ltd.北京东方泰坦科技股份有限公司,,2018年7月2日9时7分,地理编码及DEM的生成,地理编码,SAR图像的地理编码INSAR图像的地理编码是依据相关基准(datum)将影像数据从雷达成像时的坐标系统(距离-方位)转换到某种较为通用的参考坐标系统。对地图坐标系下的图像编码就是应用几何变换把图像由地图坐标系转换为距离-方位坐标系。,SAR坐标系,UTM坐标系,编码的时候依据地面形状可以分为椭球纠正编码和地形纠正编码，椭球纠正编码不需要额外的DEM，而地形纠正编码需要一个要么是地图投影或者SAR坐
2、标系的DEM。,编码流程图,DEM数据的准备Gamma只支持.hgt的DEM，对于.tif的DEM，可以借助matlab，参考第一次上课的程序。DEM参数文件的产生采用create_dem_parcreate_dem_par N36W116.dem_par要求用户定义定义一些参数DEM投影的选择多种投影类型选择，本次实验采用的是等角投影，故选择EQA，由于默认选择为EQA，所以直接敲回车键。DEM基准、椭球和地图投影的定义要求用户提供DEM数据所在的国家名，然后屏幕出现一系列可用的基准和椭球，本次试验采用WGS72。,其他参数的定义DEM title：自己定义 data format：INTE
3、GER*2DEM height offset(m):0.0 DEM height scale factor:1DEM width:3601 DEM length:3601Posting:-2.777e-04 2.777e-04Northing/easting or latitude/longitude of first DEM sample:37 -116,显示DEM图像disdem_pardisdem_par N36W116.hgt N36W116.dem_par,DEM数据的改正(缺失数据)通过内插可以填补上图中的数据空缺，具体步骤如下：GAMMA软件里面把0作为丢失值，一次必须把DEM数
4、据中的0m的值用一个和它相差不大的值代替，本实验采用1取代。用到的命令：replace_valuesreplace_values N36W116.hgt 0 1 temp_dem 3601 0 4采用同样的命令可以取代数据中的缺失值replace_values temp_dem -32768 0 temp_dem2 3601 0 4采用interp_ad内插填补空缺，根据空缺块的大小，选择合适的内插窗口，块越大，选择的窗口越大。interp_ad temp_dem2 N36W116.dem 3601 16 40 81 2 4显示内插后的图像disdem_par N36W116.dem N36W
5、116.dem_par,图像拼接gamma提供 mosaic 和multi_mosaic对图像进行拼接。 mosaic 6 N28E057.dem N28E057.dem_par N28E058.dem N28E058.dem_par N28E059.dem N28E059.dem_par N29E057.dem N29E057.dem_par N29E058.dem N29E058.dem_par N29E059.dem N29E059.dem_par N36W116.dem N36W116.dem_par 1 3显示拼接后的图像的命令：disdem_par N36W116.dem N36W
6、116.dem_par,地形纠正编码编码步骤：1)DEM/MAP参数文件的生成，采用的create_dem_par，DEM的投影变换，采用命令dem_trans。2)初始编码查询表和模拟的SAR图像的生成，采用命令gc_map。3)模拟SAR强度图像的重采样，采用geocode依据初始查询列表执行前置转换。4)采用create_diff_par ,init_offsetm, offset_pwrm和 offset_fitm精确配准模拟强度图像和真实的SAR强度图像。采用gc_map_fine利用偏差多项式精化查询列表。5)geocode_back执行后置编码把图像从SAR坐标系转化为地图坐标系
7、；geocode执行前置编码把图像从地图坐标系转化为SAR坐标系。,MLI强度图像的生成采用multi_look生成SAR强度图像 multi_look 25394.SLC 25394.SLC.par 25394.mli 25394.mli.par 2 10显示强度图像dispwr 25394.mli 1250,初始编码查询列表的生成一旦用于编码的DEM和SAR图像准备就绪，SAR和地图变换的几何关系就可以初步估计。采用命令gc_map生成初始查询列表。gc_map 25394.mli.par - N36W116.dem_par N36W116.dem N36W116.dem_seg_par
8、N36W116.dem_seg N36W116.eqa.rdc 2 2 N36W116.sim_sar显示查询列表和模拟SAR强度图像：dismph N36W116.eqa.rdc 3800dispwr N36W116.sim_sar 3800,模拟SAR图像的重采样采用geocode把模拟的SAR图像从地图转化为几何关系 geocode N36W116.eqa.rdc N36W116.sim_sar 3806 nsim_sar 1250 900 1 0显示重采样和采样前的SAR模拟图像： dis2pwr nsim_sar N36W116.sim_sar 1250 3806显示采样后和真实的S
9、AR强度图像： dis2pwr nsim_sar 25394.mli 1250 1250,计算偏移量和配准多项式的生成精化编码查询表需要计算真实SAR图像和模拟SAR图像的偏移量。首先要产生包含偏移量和配准多项式系数的DIFF/GEO的参数文件，采用命令：create_diff_parcreate_diff_par 25394.mli.par - 25394.diff_par 1初始估计偏移量：采用命令init_offsetmoffset_pwrm 25394.mli nsim_sar 25394.diff_par offs snr 512 512 offsets 1 24 24 7.0重复使
10、用offset_pwrm / offset_fitm可以改进偏差估计的质量，第一次为了提高效率可以采用少的偏移量和采样数，接着为了求得准确的偏移量和配准多项式系数采用大的偏移值和采样数。,offset_pwrm 25394.mli nsim_sar 25394.diff_par offs snr 512 512 offsets 1 8 8 7.0offset_fitm offs snr 25394.diff_par coffs coffsets 7.0 3offset_pwrm 25394.mli nsim_sar 25394.diff_par offs snr 128 128 offsets
11、 2 24 24 7.0offset_fitm offs snr 25394.diff_par coffs coffsets 7.0 3,精化编码查询列表采用gc_map_finegc_map_fine N36W116.eqa.rdc 3800 25394.diff_par N36W116.neqa.rdc 1,前/后向地理编码：后向地理编码：SAR DEMgeocode_back 25394.mli 1250 N36W116.neqa.rdc 25394.mli.eqa 3800 2654 0 0 0显示编码后的图像dispwr 25394.mli.eqa 3800前向地理编码：DEM SA
12、Rgeocode N36W116.neqa.rdc N36W116.dem_seg 3806 25394.dem 1250 - 1 0生成的DEM转换为EQA投影geocode_back 25394.dem 1250 N36W116.neqa.rdc 25394.dem.eqa 3800 2654 0 0 0显示编码后的DEMdishgt 25394.dem 25394.mli 1250dishgt 25394.dem.eqa 25394.mli.eqa 3806,由解缠相位生成DEM基线的精密估算利用已知高程控制点数据GCPs可以用命令gcp_ras从地图中提取；也可以用extract_gc
13、p从SAR坐标系下的DEM提取。本实例中给定的gcp文件是用第二种方法得到的： 提取控制点处解缠相位：gcp_phase 25394_16242.mod.unw 25394_16242.off N36W116.GCP N36W116.gcp_ph基线的最小二乘估计：base_ls 25394.SLC.par 25394_16242.off N36W116.gcp_ph 25394_16242.base 1 1 1 1 1,将解缠相位转换为高程并生成一个地距文件hgt_map 25394_16242.mod.unw 25394.SLC.par 25394_16242.off 25394_16242.base N36W116.hgt N36W116.gr 1此时的高程指的是斜距向上的的高程，采用res_map将斜距向高程转换为正射高程res_map N36W116.hgt N36W116.gr 25394.mli 25394.slc.par 25394_16242.off N36W116.ghgt N36W116.grdata 7 7 20显示正射高程：dishgt N36W116.ghgt - 2594,
