地理编码及DEM生成.pptVIP

* Beijing oriental titan technology co., ltd. 北京东方泰坦科技股份有限公司 * 地理编码及DEM的生成 地理编码 SAR图像的地理编码 INSAR图像的地理编码是依据相关基准(datum)将影像数据从雷达成像时的坐标系统(距离-方位)转换到某种较为通用的参考坐标系统。对地图坐标系下的图像编码就是应用几何变换把图像由地图坐标系转换为距离-方位坐标系。 SAR坐标系 UTM坐标系 编码的时候依据地面形状可以分为椭球纠正编码和地形纠正编码，椭球纠正编码不需要额外的DEM，而地形纠正编码需要一个要么是地图投影或者SAR坐标系的DEM。 编码流程图 DEM数据的准备 Gamma只支持.hgt的DEM，对于.tif的DEM，可以借助matlab，参考第一次上 课的程序。 DEM参数文件的产生 采用create_dem_par create_dem_par N36W116.dem_par 要求用户定义定义一些参数 DEM投影的选择 多种投影类型选择，本次实验采用的是等角投影，故选择EQA， 由于默认选择为EQA，所以直接敲回车键。 DEM基准、椭球和地图投影的定义 要求用户提供DEM数据所在的国家名，然后屏幕出现一系列可用的基准和椭球， 本次试验采用WGS72。 其他参数的定义 DEM title：自己定义 data format：INTEGER*2 DEM height offset(m):0.0 DEM height scale factor:1 DEM width:3601 DEM length:3601 Posting:-2.777e-04 2.777e-04 Northing/easting or latitude/longitude of first DEM sample:37 -116 显示DEM图像disdem_par disdem_par N36W116.hgt N36W116.dem_par DEM数据的改正(缺失数据) 通过内插可以填补上图中的数据空缺，具体步骤如下： GAMMA软件里面把0作为丢失值，一次必须把DEM数据中的0m的值用一个和它相差不大的值代替，本实验采用1取代。用到的命令：replace_values replace_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 N36W116.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 N36W116.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坐标系转化为地图坐标系； geocode执行前置编码把图像从地图坐标系转化为SAR坐标系。