地理编码（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 temp_dem 3601 0 4 1  # 0 改为 1
# replace_values temp_dem -32768 0 temp_dem2 3601 0 4  # -32768 改为 0
# interp_ad temp_dem2 sanya.dem 3601 16 40 81 2 4  # 内插
# rm temp_dem temp_dem2  # 删除临时文件
disdem_par sanya.dem sanya.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.mli 20241130.mli.par 10 20  # 多视（使其和 DEM 分辨率一致），距离向：方位向 = 10:20
raspwr 20241130.mli 1653  # 生成 bmp
eog 20241130.mli.bmp  # 查看

5. 初始地理编码

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

6. 精化查找表

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

7. 坐标系转换

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

结束

*8. 反向地理编码

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