地理编码（Radarsat2 + SRTM 30M）

1. 新建文件夹，放 DEM 进去
   这里用的是 SRTM 30m，下载地址：
   https://step.esa.int/auxdata/dem/SRTMGL1/
   

2. 生成 DEM 参数文件

create_dem_par N18E108.dem_par
EQA  # 与默认一致，直接回车
WGS84
1  # 4326
N18E108  # 名字随意
INTEGER*2  # 整型
0.0  # 与默认一致，直接回车
1.0  # 与默认一致，直接回车
3601  # 像素列数
3601  # 像素行数
-2.7777778e-04  2.7777778e-04  # 自动计算，直接回车
19 108  # DEM 左上角的经纬度

disdem_par N18E109.hgt N18E109.dem_par  # 查看 DEM

3. 镶嵌

mk_tab . hgt dem_par hgt_tab  # 生成 hgt_tab 文件
multi_mosaic hgt_tab sanya.dem sanya.dem_par 0 3  # 镶嵌
replace_values sanya.dem 0 1 sanya.dem 7201 0 4 1  # 小于等于 0 的数改为 1，防止重采样有振铃效应。7201 是 sanya.dem 的 width。4 指 DEM 是 short integer 格式
dem_trans sanya.dem_par sanya.dem sanya2.dem_par sanya2.dem 2 2 - - - 2 -  # 重采样。DEM 和多视图分辨率相差太大则需要重采样。
disdem_par sanya2.dem sanya2.dem_par  # 查看结果

4. SAR 影像多视（Radarsat2，分辨率3米）

SLC_copy 20241130_VV.slc 20241130_VV.slc.par 20241130_VV.rslc 20241130_VV.rslc.par - - 2830 16535 25891  # 裁剪 SAR 影像，最后四个数字：距离向起始，距离向大小，方位向起始，方位向大小
multi_look 20241130_VV.rslc 20241130_VV.rslc.par 20241130_VV.mli 20241130_VV.mli.par 1 2  # 多视（使其和 DEM 分辨率一致），距离向：方位向 = 1:2
raspwr 20241130_VV.mli 16535  # 生成 bmp
eog 20241130_VV.mli.bmp  # 查看

5. 初始地理编码

gc_map 20241130_VV.mli.par - sanya2.dem_par sanya2.dem DEM_seg_par DEM_seg lookup_table 1 1 sim_sar  # 生成查找表、WGS84坐标系的模拟SAR影像
geocode lookup_table sim_sar 7555 nsim_sar 16535 8939  # 模拟SAR影像地理编码至 SAR 影像坐标系。7555 是 DEM_seg_par 的 width，16535 8939 是多视图的 range、azimuth
dis2pwr nsim_sar 20241130_VV.mli 16535 16535  # 查看模拟和真实 SAR 影像。16535 是两者的 range

6. 精化查找表

rm 20241130.diff_par
create_diff_par 20241130_VV.mli.par - 20241130_VV.diff_par 1 0  # 生成参数文件
offset_pwrm nsim_sar 20241130_VV.mli 20241130_VV.diff_par offs ccp 128 128 offsets 1 - - 0.2  # 计算偏移量。128 128 是窗口大小，2 是过采样因子（2的倍数），24 24 是偏移估计数，0.2 是互相关阈值
offset_fitm offs ccp 20241130_VV.diff_par coffs coffsets 0.25 6  # 计算偏移多项式。6 是阶数。0.25 是参与计算多项式的同名点相干性阈值。距离向、方位向标准差优于 1 像元视为成功。
# final model fit std. dev. (samples) range: 1.1493   azimuth: 1.1901
gc_map_fine lookup_table 16535 20241130_VV.diff_par lookup_table.fine 1  # 优化查找表。16535 是多视图的 range

7. 坐标系转换

geocode lookup_table.fine DEM_seg 7555 sanya20241130.dem 16535 8939  # DEM 由 WGS84 坐标系转为 SAR 坐标系。7555 是 DEM_seg_par 的 width，16535 8939 是多视图的 range、azimuth
# dis2dt_pwr sanya20241130.dem sanya20241130.dem 20241130.mli 1653 1653 - - - - - - 2  # 查看。2 指自动识别最大最小值缩放
rasdt_pwr sanya20241130.dem 20241130_VV.mli 16535 - - - - - - 2  # 生成 bmp。16535 是多视图的 range。2 指自动识别最大最小值缩放
eog sanya20241130.dem.bmp

结束

*8. 反向地理编码

geocode_back 20241130_VV.mli 16535 lookup_table.fine 20241130_VV.mli.eqa 7555  # 多视图反向地理编码至 WGS84 坐标系。16535 是多视图的 range，7555 是 DEM_seg_par 的 width
# dispwr 20241130_VV.mli.eqa 7555
raspwr 20241130_VV.mli.eqa 7555
eog 20241130_VV.mli.eqa.bmp
data2geotiff DEM_seg_par 20241130_VV.mli.eqa 2 20241130_VV.mli.eqa.tif  # 生成 tif 文件