SARscape5.2中的数据裁剪功能
分享
-
2016-06-15
数据裁剪是做SAR数据处理常做的一步准备工作，当数据范围远大于研究区的时候，先对数据进行裁剪，可以减少数据量，节省数据处理的时间。
随着ENVI5版本的推出，SARscape的软件界面也发生了变化，采用统一的工程化界面，更加便于操作。本文介绍SARscape最新版本5.2.1中的SAR数据裁剪工具。SARscape5.0之前的版本请参考文档：http://blog.sina.com.cn/s/blog_764b1e9d0101ciz7.html。
SARscape提供了两个裁剪工具
图 SARscape中的裁剪工具
裁剪功能支持坐标范围裁剪和矢量裁剪两种方式。如果使用坐标范围裁剪，直接在裁剪的参数面板输入4个坐标角点即可，如下图所示。
图 输入坐标范围裁剪
如果使用矢量文件裁剪，直接支持shapefile矢量数据，如果没有现成的矢量数据，需要在ENVI中绘制矢量文件，方法为：ENVI主菜单File->New->Vector，新建矢量文件，选择要绘制矢量文件的栅格数据，设置文件名，点击OK，然后鼠标绘制矢量文件，绘制完成后点击鼠标右键，单击"accept"结束绘制，如果要修改矢量文件，点击Create Vector->Edit Vertex修改。
图 矢量编辑功能
SAR坐标系数据的裁剪
Sample Selection SAR Geometry Data，该工具适用于SAR坐标系数据的裁剪，输入的数据是SAR坐标系下的数据，如_slc、_slc_list、_pwr、_gr……。
注：如果输入的是Sentinel数据的_slc_list文件，矢量需要在多视后的_pwr数据上绘制，该强度数据_pwr也需要作为参考文件（reference file）输入。
如果用地理坐标范围裁剪，输入的是地理坐标下的矢量数据或坐标范围，在参数设置面板中"Geographical Region"参数设置为True；如果用斜距坐标范围裁剪，输入的是斜距坐标系下矢量数据或坐标范围，在参数设置面板中"Geographical Region"参数设置为False。
下面以哨兵1A的SLC数据为例，介绍斜距数据的裁剪。
第一步，在ENVI中打开导入后的哨兵1A数据生成的强度数据文件
sentinel1_41_20150325_115732559_IW_SIW1_A_VV_slc_list_pwr，点击File->New->Vector Layer，设置矢量文件名，点击OK。
图 创建矢量文件
图 在强度图上绘制子区域的矢量
在矢量图层点击右键，选择Save As，保存到本地的矢量文件，矢量绘制的不规则也没有关系，可以通过"Use Min and Max Coordinates "参数设置，根据矢量文件的最大坐标进行规则裁剪。
打开/SARscape/General Tools/Sample Selections/Sample Selection SAR Geometry Data工具。
输入SAR坐标系的数据文件，可以是一景，也可以是同一区域的多景数据，默认的数据标识为_slc、_pwr、_slc_list，如果数据不带这个标识，可以在文件筛选处选择*.*。
图 数据输入面板
可选文件面板（Optional File）
1)Vector File：输入上一步绘制的矢量文件。该文件是可选项，如果用坐标范围裁剪，在此不输入。
2)DEM File：输入带有地理坐标的DEM数据提供参考坐标系。该文件是可选项，如果用地理坐标范围裁剪，在此必须输入带有地理坐标的DEM文件，以提供参考坐标系。在此不输入DEM文件。
3)Input Reference File：参考数据，输入多视的强度数据_pwr。在以下情况下，该数据文件为必选项： 这里输入待裁剪的该文件的_pwr数据。
图 可选文件面板
参数设置面板（Parameters）：主要参数（Principal Parameters）
图 主要参数设置
输出文件（Output file list）：自动将数据输出到ENVI默认的数据输出路径下并自动命名，添加_cut的标识。
图 输出面板设置
设置好参数之后，点击Exec执行即可。
图 _slc_list数据裁剪结果的强度图
图 裁剪之后的数据文件
注：如果裁剪的是哨兵1的_slc_list数据，当裁剪范围很小的时候，可能会遇到裁剪结果不全的情况。原因是软件默认的参数中，子区域在每个条带中，覆盖范围大于20%的进行裁剪，否则该条带不参与裁剪。
解决方法为：（1）扩大要裁剪的子区域范围；（2）设置覆盖范围小于20%。在裁剪工具的Parameters参数面板中，选择Cut参数，Perc Valid参数默认为20，把该参数调小。如下图所示。
图 设置Sentinel _slc_list文件默认的子区覆盖范围
地理坐标系数据的裁剪
Sample Selection Geographic Data，该工具适用于地理坐标系数据的裁剪，输入的数据是地理坐标系下的数据，如_geo、_dem、_disp……。启动/SARscape/General Tools/Sample Selections/Sample Selection Geographic Data。
输入带有地理坐标的数据文件，可以是一景，也可以是同一区域的多景数据，默认的数据标识为_geo，如果数据不带这个标识，可以在文件筛选处选择*.*。
输入.shp格式的矢量文件。该文件是可选项，如果用坐标范围裁剪，在此不输入。
通过带有地理坐标的DEM数据/坐标系，确定要裁剪的地理坐标的参考坐标系。该设置是可选项，如果用坐标范围裁剪，在此必须设置参考坐标系，DEM或者坐标系，选择其中一个。
图 地理坐标数据裁剪工具参数面板
输出文件（Output file list）：自动将数据输出到ENVI默认的数据输出路径下并自动命名，添加_cut的标识。
设置好参数之后，点击Exec执行即可。
文章来源：http://blog.sina.com.cn/s/blog_764b1e9d0102xbw8.html
随着ENVI5版本的推出，SARscape的软件界面也发生了变化，采用统一的工程化界面，更加便于操作。本文介绍SARscape最新版本5.2.1中的SAR数据裁剪工具。SARscape5.0之前的版本请参考文档：http://blog.sina.com.cn/s/blog_764b1e9d0101ciz7.html。
SARscape提供了两个裁剪工具
- Sample Selection SAR Geometry Data ——SAR坐标系数据的裁剪工具
- Sample Selection Geographic Data——对地理坐标系数据的裁剪工具
图 SARscape中的裁剪工具
裁剪功能支持坐标范围裁剪和矢量裁剪两种方式。如果使用坐标范围裁剪，直接在裁剪的参数面板输入4个坐标角点即可，如下图所示。
图 输入坐标范围裁剪
如果使用矢量文件裁剪，直接支持shapefile矢量数据，如果没有现成的矢量数据，需要在ENVI中绘制矢量文件，方法为：ENVI主菜单File->New->Vector，新建矢量文件，选择要绘制矢量文件的栅格数据，设置文件名，点击OK，然后鼠标绘制矢量文件，绘制完成后点击鼠标右键，单击"accept"结束绘制，如果要修改矢量文件，点击Create Vector->Edit Vertex修改。
图 矢量编辑功能
SAR坐标系数据的裁剪
Sample Selection SAR Geometry Data，该工具适用于SAR坐标系数据的裁剪，输入的数据是SAR坐标系下的数据，如_slc、_slc_list、_pwr、_gr……。
注：如果输入的是Sentinel数据的_slc_list文件，矢量需要在多视后的_pwr数据上绘制，该强度数据_pwr也需要作为参考文件（reference file）输入。
如果用地理坐标范围裁剪，输入的是地理坐标下的矢量数据或坐标范围，在参数设置面板中"Geographical Region"参数设置为True；如果用斜距坐标范围裁剪，输入的是斜距坐标系下矢量数据或坐标范围，在参数设置面板中"Geographical Region"参数设置为False。
下面以哨兵1A的SLC数据为例，介绍斜距数据的裁剪。
第一步，在ENVI中打开导入后的哨兵1A数据生成的强度数据文件
sentinel1_41_20150325_115732559_IW_SIW1_A_VV_slc_list_pwr，点击File->New->Vector Layer，设置矢量文件名，点击OK。
图 创建矢量文件
图 在强度图上绘制子区域的矢量
在矢量图层点击右键，选择Save As，保存到本地的矢量文件，矢量绘制的不规则也没有关系，可以通过"Use Min and Max Coordinates "参数设置，根据矢量文件的最大坐标进行规则裁剪。
打开/SARscape/General Tools/Sample Selections/Sample Selection SAR Geometry Data工具。
- 数据输入面板（Input File）
输入SAR坐标系的数据文件，可以是一景，也可以是同一区域的多景数据，默认的数据标识为_slc、_pwr、_slc_list，如果数据不带这个标识，可以在文件筛选处选择*.*。
图 数据输入面板
可选文件面板（Optional File）
1)Vector File：输入上一步绘制的矢量文件。该文件是可选项，如果用坐标范围裁剪，在此不输入。
2)DEM File：输入带有地理坐标的DEM数据提供参考坐标系。该文件是可选项，如果用地理坐标范围裁剪，在此必须输入带有地理坐标的DEM文件，以提供参考坐标系。在此不输入DEM文件。
3)Input Reference File：参考数据，输入多视的强度数据_pwr。在以下情况下，该数据文件为必选项： 这里输入待裁剪的该文件的_pwr数据。
- 待裁剪的数据是Sentinel-1的_slc_list文件
- 基于雷达坐标系的矢量文件进行裁剪时
- 基于雷达坐标系的坐标进行裁剪时
图 可选文件面板
参数设置面板（Parameters）：主要参数（Principal Parameters）
- Make Coregistration：默认为False。若激活该参数，多幅影像裁剪时，会输出裁剪后的配准结果。
- Coregistration With DEM：如果上面参数设置为True，该参数设置配准时是否考虑DEM。
- Geographical Region：设置为False，基于SAR坐标系的裁剪，如果输入的是地理坐标或者矢量，该参数设置为True。
- West/First Column：最小经度/起始列号
- North/First Row：最大纬度/起始行号
- East/Last Column：最大经度/最大列号
- South/Last Row：最小纬度/最大行号
- Use Min and Max Coordinates：设置为True。如果该参数设置为True，会按照矢量文件的最大的坐标范围对数据进行裁剪。
图 主要参数设置
- 数据输出面板（Output Files）
输出文件（Output file list）：自动将数据输出到ENVI默认的数据输出路径下并自动命名，添加_cut的标识。
图 输出面板设置
设置好参数之后，点击Exec执行即可。
图 _slc_list数据裁剪结果的强度图
图 裁剪之后的数据文件
注：如果裁剪的是哨兵1的_slc_list数据，当裁剪范围很小的时候，可能会遇到裁剪结果不全的情况。原因是软件默认的参数中，子区域在每个条带中，覆盖范围大于20%的进行裁剪，否则该条带不参与裁剪。
解决方法为：（1）扩大要裁剪的子区域范围；（2）设置覆盖范围小于20%。在裁剪工具的Parameters参数面板中，选择Cut参数，Perc Valid参数默认为20，把该参数调小。如下图所示。
图 设置Sentinel _slc_list文件默认的子区覆盖范围
地理坐标系数据的裁剪
Sample Selection Geographic Data，该工具适用于地理坐标系数据的裁剪，输入的数据是地理坐标系下的数据，如_geo、_dem、_disp……。启动/SARscape/General Tools/Sample Selections/Sample Selection Geographic Data。
- 数据输入面板（Input File）
输入带有地理坐标的数据文件，可以是一景，也可以是同一区域的多景数据，默认的数据标识为_geo，如果数据不带这个标识，可以在文件筛选处选择*.*。
- 可选文件面板（Optional File）
输入.shp格式的矢量文件。该文件是可选项，如果用坐标范围裁剪，在此不输入。
- 参考坐标系设置面板（DEM/Cartographic System）
通过带有地理坐标的DEM数据/坐标系，确定要裁剪的地理坐标的参考坐标系。该设置是可选项，如果用坐标范围裁剪，在此必须设置参考坐标系，DEM或者坐标系，选择其中一个。
- 参数设置面板（Parameters）：主要参数（Principal Parameters）
Make Max Common Area：激活该参数，在输入的所有数据中，裁剪范围内，按照所有数据的最大交集裁剪，常用于升轨-降轨数据的裁剪或者不同传感器数据的裁剪。
Common Dummy Area：地理坐标裁剪，此功能默认为False
Geographical Region：地理坐标裁剪，此功能默认为True
West/First Column：最小经度
North/First Row：最大纬度
East/Last Column：最大经度
South/Last Row：最小纬度
Use Min and Max Coordinates：如果该参数设置为True，会按照矢量文件的最大的坐标范围对数据进行裁剪。图 地理坐标数据裁剪工具参数面板
- 数据输出面板（Output Files）
输出文件（Output file list）：自动将数据输出到ENVI默认的数据输出路径下并自动命名，添加_cut的标识。
设置好参数之后，点击Exec执行即可。
文章来源：http://blog.sina.com.cn/s/blog_764b1e9d0102xbw8.html
