SARscape5.7中环境减灾二号05星SAR数据处理方法

环境减灾二号05卫星(别名:5米S-SAR 01星)于2022年10月13日发射成功,与2023年8月9日发射成功的环境减灾二号06星进行组网运行。

缺少环境减灾二号05卫星产品说明书,本文缺少对该卫星数据的详细说明,文件命令为常见的命名方式,如下所示:

HJ2E_KSC_STRIP_006412_E119.2_N35.3_20231207_SLC_HV_L10000141614

能得到该数据为条带成像模式(STRIP)的HV极化的L1级SLC产品。

最新的SARsacpe5.7版本已经开始支持环境减灾二号05雷达卫星的数据处理。本文是在SARscape5.7下完成,包括数据导入、多视处理、滤波处理、地理编码。

注:不出意外,环境减灾二号06也可以采用相同的处理方法。

软件试用方式:

1、安装SARscape5.7:https://envi.geoscene.cn/sar_license

2、下载补丁并安装:

https://pan.baidu.com/s/1RbB-V4PZRpSIXSg8f-Zcrg?pwd=ENVI

1数据导入

本文档使用的是环境减灾二号05 L1级产品,它是单视复数据(SLC),条带模式(STRIP)成像模式,标称分辨率为5米。

在处理之前,我们选择一套默认参数,打开Toolbox/SARscape/Preferences/Preferences specific,在打开的界面中,选择Load Preferences->VHR(better than 6m)。再点击General parameters选项,设置Cartographic Grid Size(m):5。

(1)在Toolbox中,选择SARscape/Import Data/SAR Spaceborne/Single Sensor/HJ-2。

(2)在打开的面板中,

  • 数据输入面板(Input Files)
  • 输入文件(Input File List):输入过滤的.xml文件。支持多景输入批量输入。
  • 参数设置面板(Parameters):主要参数(Principal Parameters)
  • 对数据重命名(Rename the File Using Parameters):True。软件会自动在输入文件名的基础上增加几个标识字母,如增加“_HH_slc”。

注:如果输入多景同一天成像数据,输出的文件名会重复,这里设置False。

  • 数据输出面板(Output Files)

输出文件(Output file list):自动读取ENVI默认的数据输出目录以及输入面板中的数据文件名。

注:1.如果要修改输出的路径,在右边单击文件夹图标选择输出文件夹目录。

    2、如果要修改输出的路径,在输出文件名右键选择Edit菜单。

(3)单击Exec按钮开始执行。

图:数据导入

图:导入后图像(水域区域有部分像素为0值,自动透明显示)

生成的结果除了图像文件外,还包括Shapefile和KML格式的图像轮廓线。

2流程化工具处理

SARscape提供单时相SAR流程化处理工具和分步处理工具。本文采用的是单时相SAR流程化处理工具。

注:分步工具可以从/SARscape/Basic/Intensity Processing中选择对应的多视、滤波、地理编码工具,分布式处理工具一次可以处理N景影像。

启动流程化的SAR处理工具:双击/SARscape/Basic/Intensity Single Image Workflow,打开Intensity Single Image流程化工具面板。流程化的处理工具界面左侧是工作流的各个步骤,右侧是相应步骤的输入输出及参数设置,下方是步骤控制按钮Back和Next按钮,在处理过程中可点击Back和Next按钮切换至中间某一步查看参数或调整参数重新处理。

第一步:数据输入(Input)

  (1) 数据输入(Input File)

  • Area of interest in Geographic Coordinates:输入裁剪区域。经纬度坐标的shp文件、KMZ、KML文件格式。
  • Input File:输入上一步导入得到的文件。

(2)参考坐标系设置(DEM/Cartographic System),输入参考高程。

Reference下拉菜单中提供了两种方式:

  • Input DEM:输入已有的参考DEM文件。
  • DEM Download:下载相应区域的DEM数据。

本例选择DEM Download方式,单击Type选择ALOS World 3D 30m。

(3)参数设置(Parameters):Grid Size:5m。

(4)单击Next按钮。

第二步:输入原始SAR数据(Import Generic SAR Data)

支持一些传感器直接输入原始文件。本例已经导入过数据,这一步不用设置,直接单击Next按钮。

第三步:多视(Multilooking)

(1)主要参数设置(Principal Parameters):

  • 多视方法(Multilooking Method):Time domain。
  • 软件根据元数据文件和设置的制图分辨率(Grid Size for Suggested Looks:5m)自动算出距离向视数(Range Looks):1,方位向视数(Azimuth Looks):2。

(2)单击Next按钮执行下一步。

注:设置Global->Generate Quick Look:False,不输出快视图,能提高一点处理速度。

图:多视处理步骤

第四步:滤波(Single Image Filtering)

(1)主要参数设置(Principal Parameters):

  • 滤波方法(Filter Method):提供了7种滤波方法,本例选择Lee方法。
  • 等值视数(Equivalent Number of Looks):-1

注:1、需要设置每个滤波方法的参数,可选择Filters parameters选项进行参数设置。

2、根据需求可设置Global->Generate Quick Look:False,不输出快视图,能提高一点处理速度。

(2)单击Next按钮执行下一步。

第五步:地理编码与定标(Geocoding and Radiometric Calibration)

(1)主要参数(Principal Parameters):

  • 像元大小(X Grid Size):5。
  • 像元大小(Y Grid Size):5。
  • 辐射定标(Radiometric Calibration):False。

注:环境减灾二号05星目前元数据中没有提供定标所需的参数,后面几个参数不需要设置。

  • 散射面积(Scattering Area):Local Incidence Angle。
  • 辐射归一化(Radiometric Normalization):False。

注:辐射归一化后得到的为Gamma计量单位。

  • 局部入射角校正(Local Incidence Angle):False。
  • 输出类型(Output type):Linear。

注:输出类型linear和Db的关系为:∂^0 (Db)=10*log∂^0

(2)单击Next按钮执行下一步。

图:地理编码与定标参数设置

最后一步:输出(Output)

(1)选择输出结果路径。

(2)Delete Temporary Files:勾选删除中间结果。

(3)单击Finish按钮。

3浏览结果

本实例数据覆盖区域海边区域,如下为HV极化图像,可以看到图像纹理清晰,能分辨浅海海底部分信息,海面上的养殖场依稀可见。

图:局部效果图

叠加到天地图卫星底图上,缩放到SAR合适尺度上,肉眼几乎区分不出几何偏差,几何定位非常精确。

图:叠加天地图卫星图像

posted @ 2024-07-08 15:09  ENVI-IDL技术殿堂  阅读(474)  评论(0编辑  收藏  举报