GAMMA 进行GEO

##下载指定区域小范围内的SRTM

安装elevation 模块

#pip install elevation

具体安装可参照以下网址，直接百度也行

安装了乌版图条件几乎都满足，差不多安装个GDAL就行

elevation · PyPI

例子：

eio --product SRTM1 clip -o ./mosaic.tif --bounds 90.0 30.0 91.0 31.0

--product可选择30m左右的和90m左右的

-o 存放路径，并命名

--bounds 左下角和右上角得经纬度

(GAMMA做GEO时,不需要手动下载SRTM，不需要将hgt转tif,虽然可以再USGS下载tif格式，但是自动还是比较方便）

GAMMA做GEO时，可根据SLC_corners 显示研究区的经纬度范围

1、srtm2dem

2、multi_look

3、gc_map（pixel_area）生成初始查找表/模拟强度图

4、offset_pwrm 计算强度图之间的偏移量

5、offset_fitm 计算拟合系数（检查多项式拟合标准差）

6、gc_map_fine 根据拟合模型精化初始查找表

7、geocode_back

这种地理编码模式是基于地形校正的编码模式，在地形起伏和大区域编码结果较好（个人经验，可仅供参考）。基本思想：根据外部DEM和几何参数模拟的SAR强度图与多视强度图进行配准，精化查找表。

注意：一定要检查final model fit std. dev. (samples) range和number of offsets above correlation threshold，不然对后续结果有影响。通过不断调整第4和第5步骤来得到合适的结果。但也不用过于关注相应的指标，可将SAR数据编码至地理坐标，与光学影像进行对比。

#若是其他的DEM，就不要用第一步得啦，但是基本操作都是这些。

另外特别提醒，对于小区域的平坦地方可能很难达到合适得精度，自己做的时候注意，官方的例子数据也说明了这一点，虽然文档中有说明人工选择控制点，进行GEO，但是自己还没一尝试过。

#如有任何问题，欢迎大家批评指正交流！！！