GAMMA处理数据（四）

之前不管三七二十一，以为所有的影像尺寸会是一样的，就统一对所有多视数据设置了同样大小尺寸进行了bmp文件生成，然后发现有三个bmp图片怪怪的，经请教，发现是MLI具有不同的尺寸，至于原因，我也不是很清楚。

1、解决因多视强度数据尺寸大小不一样导致bmp图片异常

批量从vv.mli.par文件中取出range_samples值并查看，因为没有配准，此时尺寸是不一样的

run_all.pl dates 'get_value $1.vv.mli.par range_smaples >> tmp1'

有三个多视强度数据的尺寸稍微小一点

重新生成了对应尺寸的bmp多视强度图片，可以正常显示

2、地理编码

（1）首先对参考影像进行多视，这里选择中间的影像作为参考影像:20210311.vv.mli作为参考影像，SLC_corners得到脚点坐标

（2）导入DEM，dem_import，之前没注意到还要选择geoid，geoid_par这两个参数，看了官方给的manual，里面有注释：the geoid is used to transform the DEM height from an altitude above sea level(i..e. above the geoid) to an altitude above the ellipsid. 意思就是说：大地水准面用来将正高转换为大地高。DEM采用正高模型进行建模，通过卫星或航空传感器收集影像时，通常在地球的椭球体模型之上测量位置信息（尤其是高度信息），因此高度值会有所不同，所以需要转换。

在arcgis的说明文档中指出：

GAMMA软件提供了两种大地水准面模型：

至于这两种模型的差别：那我选择精度更高的EGM2008模型

DEM par文件：

（3）显示DEM的晕渲图（a shaded relief），大概是可以看出地形的意思？？？至于看效果，还不知道怎么看

（4）生成初始查找表，将DEM数据从地理坐标系转换到SAR坐标系，这里的查找表后面还要进行更新

使用命令：gc_map2，输入一系列参数

输入参数依次为：多视后的par文件、DEM par文件、DEM数据

lat_ovr、lon_ovr是什么过采样因子？？？不明白这个因子的作用是什么，决定生成像素尺寸的大小？

输出文件：dem segment文件以及它的par文件，这些文件是什么我还不清楚，DEM分段文件？？？

初始查找表，类型是fcomplec类型，在原始DEM几何下，我不知道这个类型是什么，百度了一下

ls_map：叠掩和阴影图（原始的DEM几何下）

ls_map_rdc：叠掩和阴影图（SAR几何下）

inc：局部入射角图，然后其余的输出文件感觉用不太上，就没有仔细研究了

输出结果、还有一些输出文件：

 （5）模拟SAR强度影像

使用命令：piexl_area，基于DEM的模拟后向散射图，pix_gamma0为斜距几何结构

输入：多视后的par文件、DEMpar文件、DEM数据、地理编码初始查找表、ls_map以及inc_map、其余参数：nstep、area_fact不明白什么意思，暂时默认

输出文件：sigma0、gamma0，百度了一下，这里应该是在做辐射定标，应该是辐射定标后的输出文件

然后可以显示一下：利用dis2pwr显示模拟强度图和MLI图像

（6）根据模拟影像和真实MLI影像，可以计算出查找表的校正。计算gamma0与多视影像之间的偏移量。

        1.使用命令：create_diff_par

        2.使用命令：offset_pwrm

        3.使用命令：offser_fitm

offset_fitm计算出精化多项式是一个常数，有了高质量的轨道数据，不需要更高阶的多项式，精化质量可以通过这个命令显示在屏幕上

（7）精化查找表

gc_map_fine，这里有个参数不知道怎么设置：width，让我输入查找表的宽度，既然是将DEM变换到SAR坐标系下，是不是说明应该输入DEM的宽度。是在之前生成的20210311_seg.dem_par文件中查找width

生成的hgt文件，显示如下：

（8）计算更新之后的sigma0、gamma0

使用piexl_area命令，跟上面的步骤一样

（9）计算SAR坐标系下的DEM

使用geocode：

在生成hgt如果行数输入错误，就会在配准步骤弹出以下报错信息： 

 改正的话，按照它报出的行数重新生成hgt文件就行，然后进行重新配准，就不会报错了。

（10）地理编码-后向编码

用于将SAR坐标系转到地理坐标系下，输出的是地理坐标系下的影像png

3、配准

（1）选择位于中间的影像作为主影像

比如有25个日期，我可能会选择第12个日期为主影像，让所有影像向它配准

（2）准备批量配准组合

第一个组合为主影像前后的日期

第二个组合为主影像前后的日期与它们自己前后日期的组合

（3）配准最初相邻两景影像（配准组合1）、配准其他数据（配准组合2）

使用命令：ScanSAR_coreg.py，配准的生成文件是以.rslc结尾的文件，之前因为改了名字，导致干涉的时候用错了文件，耽误了好长时间

关于配准的总结：

首先是对参考主影像前后两个时间，与参考主影像进行配准

后面是将主影像与其他的所有影像进行配准，但是这里存在一个问题，比如：主影像与其他影像相隔时间太长的话，配准可能会失败，所以要有已经与主影像配准好的影像作为“辅助”影像，所以在第二个命令中需要一个RSLC_tab3的文件

（5）相位去斜：不明白这个操作的意义是什么？但是在我这个实验中，相位去斜是不必须的

官方对于这个命令的描述如下：

DESCRIPTION
S1_deramp_TOPS_reference is used to deramp the burst SLC data of a S1 TOPS reference scene. For each burst the azimuth phase ramp caused by the Doppler Centroid variation is calculated and subtracted. Both the deramped data (as burst SLC, i.e. with the individual files for each sub-swath) and the phase ramps subtracted (in the same format, but using float values) are generated.

Deramping of the S1 TOPS data may be used to apply oversampling in offset tracking, apply common band filtering in azimuth direction or to apply split-beam interferometry, or to use the spectral diversity criteria in the IPTA for the point target candidate identification (program sp_stat). For differential interferoemtry, on the other hand, deramping of the S1 IWS data is not necessary.

使用命令：

run_all.pl dates 'S1_deramp_TOPS_slave $1_rslc_tab $1 20210311_rslc_tab 8 2 1'

（6）裁减研究区

虽然这次实验可以不进行研究区的裁减，但是还是要了解裁减的命令是什么

run_all.pl 'SLC_copy $1.rslc.deramp $1.rslc.deramp.par $1.rslc.deramp.crop $1.rslc.deramp.crop.par - - x x x x'

这里参数应该是研究区位于哪一range samples? 看到示例参数输入了四个数字，但是我怎么知道我的研究区在哪一range samples？

到这里，数据预处理以及配准步骤就做完了，下一步就是基线估计、偏移量估计、干涉图生成、地形相位去除、滤波等处理