offset_pwrm 25394.mli nsim_sar 25394.diff_par offs snr 128 128 offsets 2 24 24 7.0
offset_fitm offs snr 25394.diff_par coffs coffsets 7.0 3 精化编码查询列表 采用gc_map_fine

初始编码查询列表的生成 一旦用于编码的DEM和SAR图像准备就绪，SAR和地图变换的几 何关系就可以
初步估计。采用命令gc_map生成初始查询列表。
gc_map 25394.mli.par - N36W116.dem_par N36W116.dem N36W116.dem_seg_par N36W116.dem_seg N36W116.eqa.rdc 2 2 N36W116.sim_sar 显示查询列表和模拟SAR强度图像： dismph N36W116.eqa.rdc 3800 dispwr N36W116.sim_sar 3800
SAR坐标系
UTM坐标系

编码的时候依据地面形状可以分为椭球纠正编码和地形纠正编码，椭 球纠正编码不需要额外的DEM，而地形纠正编码需要一个要么是地图 投影或者SAR坐标系的DEM。
编码流程图

DEM数据的准备
Gamma只支持.hgt的DEM，对于.tif的DEM，可以借助matlab，参考第一次上

offset_pwrm 25394.mli nsim_sar 25394.diff_par offs snr 512 512 offsets 1 8 8 7.0 offset_fitm offs snr 25394.diff_par coffs coffsets 7.0 3
offset_pwrm 25394.mli nsim_sar 25394.diff_par offs snr 512 512 offsets 1 24 24 7.0
重复使用offset_pwrm / offset_fitm可以改进偏差估计的质量，第一次为了提 高效率可以采用少的偏移量和采样数，接着为了求得准确的偏移量和配准 多项式系数采用大的偏移值和采样数。

将解缠相位转换为高程并生成一个地距文件 hgt_map 25394_16242.mod.unw 25394.SLC.par 25394_16242.off 25394_16242.base N36W116.hgt N36W116.gr 1 此时的高程指的是斜距向上的的高程，采用res_map将斜距向高程转换为 正射高程

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
地理编码及DEM的生成
Beijing oriental titan technology co., ltd.
2016/2/7 2016年2月7日6时55分
北京东方泰坦科技股份有限公司

地理编码
SAR图像的地理编码
INSAR图像的地理编码是依据相关基准(datum)将影像数据从雷达成像时的 坐标系统(距离-方位)转换到某种较为通用的参考坐标系统。对地图坐标系 下的图像编码就是应用几何变换把图像由地图坐标系转换为距离-方位坐标 系。
gc_map_fine N36W116.eqa.rdc 3800 25394.diff_par N36W116.向地理编码：SAR → DEM geocode_back 25394.mli 1250 N36W116.neqa.rdc 25394.mli.eqa 3800 2654 000 显示编码后的图像 dispwr 25394.mli.eqa 3800 前向地理编码：DEM → SAR geocode 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 显示编码后的DEM dishgt 25394.dem 25394.mli 1250
要求用户提供DEM数据所在的国家名，然后屏幕出现一系列可用的基准和椭球，
本次试验采用WGS72。

其他参数的定义 DEM title：自己定义 DEM height offset(m):0.0 DEM width: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 data format：INTEGER*2 DEM height scale factor:1 DEM length:3601
像和真实的SAR强度图像。采用gc_map_fine利用偏差多项式精化查询列表。
5)geocode_back执行后置编码把图像从SAR坐标系转化为地图坐标系； 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参数文件的产生 采用create_dem_par create_dem_par N36W116.dem_par 要求用户定义定义一些参数 DEM投影的选择 多种投影类型选择，本次实验采用的是等角投影，故选择EQA， 由于默认选择为EQA，所以直接敲回车键。 DEM基准、椭球和地图投影的定义

地形纠正编码
编码步骤： 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图像和模拟SAR图像的偏移量。首先要 产生包含偏移量和配准多项式系数的DIFF/GEO的参数文件，采用命令： create_diff_par create_diff_par 25394.mli.par - 25394.diff_par 1 初始估计偏移量：采用命令init_offsetm

图像拼接 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
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


dishgt 25394.dem.eqa 25394.mli.eqa 3806
由解缠相位生成DEM 基线的精密估算——利用已知高程控制点数据 GCPs可以用命令gcp_ras从地图中提取；也可以用extract_gcp从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 1111

模拟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 显示采样后和真实的SAR强度图像： dis2pwr nsim_sar 25394.mli 1250 1250