利用kmz文件裁剪SAR图像

我们进行 DInSAR 、 PS-InSAR 和 SBAS-InSAR 等处理时候，经常要对 SLC （单视复数数据）进行裁剪，针对所需工程区进行 InSAR 分析，减少工作量提高效率。利用 Shp 矢量文件（无地理坐标）裁剪斜距坐标和地理坐标的 SAR 图像可参考（ http://blog.sina.com.cn/s/blog_764b1e9d0102xbw8.html ）。

SARscpae5.5 开始支持 kmz 和 kml 文件对 SAR 图像的裁剪，本文介绍直接使用 kmz 文件对斜距坐标的 SAR 图像进行裁剪操作过程， kmz 文件可以非常方便的在 Google Earth 中绘制获取。

注：本文在 SARscpe5.5 中操作完成，第二步也适应于基于地理坐标的 Shp 文件裁剪 SAR 图像的操作。

第一步：制作 kmz 文件

（1） 打开 Google Earth ，定位到需要绘制范围的地点。

（2） 在 工具 栏上，打开“绘制多边形”按钮。

（3） 鼠标左键在影像地图上点击绘制需要裁剪的范围，拖动顶点可以随时修改裁剪范围。

（4） 在对话框中点击确定完成绘制。

（5） 在左边“我的地点”上，刚绘制完成的“裁剪区域”右键选择“将位置另存为”菜单，以字母或数字命名保存 kmz 文件。

注：右键选择“属性”菜单，可对裁剪区域进行编辑。

图：绘制多边形对话框

第二步：利用 kmz 文件裁剪 SAR 图像

（1） 启动 SARscape ，在 ToolBox 中打开 /SARscape/General Tools/Sample Selections/Sample Selection SAR Geometry Data 工具。

（2） 在数据输入面板（ Input File ） 中，选择输入一景或者多景 SAR 图像， SAR 图像可以是 _slc 、 _pwr 、 _slc_list,_gr 。

（3） 可选文件面板（ Optional File ）中。

l   Vector File ：输入第一步中获取的 .kmz 文件。

l   DEM File ：输入带有地理坐标的 DEM 数据提供参考坐标系。

注：如果用地理坐标范围裁剪，在此必须输入带有地理坐标的 DEM 文件，以提供参考坐标系。如果用 SAR 斜距坐标范围裁剪，不用选择 DEM 。

l   Input Reference File ：参考数据，以下情况需要输入：

输入的是 SCANSAR 或者 TOPSAR _slc_list 数据时候，并且基于 SAR 斜距坐标的矢量文件进行裁剪时，或者基于雷达坐标系的坐标进行裁剪时，这里输入待裁剪的该文件的 _pwr 数据。

   本例子使用带有地理坐标的 kmz 文件裁剪，这个不需要输入参考文件。

图：可选文件面板（ Optional File ）

（4） 参数设置面板（ Parameters ）中。

主要参数（ Principal Parameters ）：

l   Make Coregistration ：默认为 False 。若激活该参数，多幅影像裁剪时，会输出裁剪后的配准结果。

l   Coregistration With DEM ：如果上面参数设置为 True ，该参数设置配准时是否考虑 DEM 。

l   Geographical Region ：设置为 Ture ，基于地理坐标系的裁剪，如果输入的是地理坐标或者矢量，该参数设置为 True 。

以下四个参数是基于对角坐标进行裁剪，这里基于 kmz 文件裁剪，选择默认。

l   West/First Column ：最小经度 / 起始列号

l   North/First Row ：最大纬度 / 起始行号

l   East/Last Column ：最大经度 / 最大列号

l   South/Last Row ：最小纬度 / 最大行号

l   Use Min and Max Coordinates ：设置为 True 。如果该参数设置为 True ，会按照矢量文件的最大的坐标范围对数据进行裁剪。如果做 InSAR 处理，建议选择 Ture 。

图：主要参数（ Principal Parameters ）面板

    如果输入裁剪的数据为 TOPSAR/ScanSAR ，切换到 Cut 参数选项，根据裁剪区域的大小，缩小 Min Valid Square Dimension （默认为 500 像素， SAR 图像分块落在裁剪区的像素小于这个值，该分块将会放弃不参与裁剪）阈值。

（5） 数据输出面板（ Output Files ）

输出文件（ Output file list ）：自动将数据输出到 ENVI 默认的数据输出路径下并自动命名，添加 _cut 的标识。

如需要更改输出路径，在某一行数据上右键选择“ Change Output Directories ”菜单。

（6） 点击 Exec 执行 。