拓冰网站定制 
简介面向遥感数据处理与InSAR学习者的操作型参考资料聚焦ENVI中SARscape模块对ALOS PALSAR等数据的预处理全流程。内容覆盖数据导入、单视复数产品多视处理、斑点噪声滤波、地理编码与辐射定标并按步骤给出参数设置建议如视数估算、滤波窗口、投影与重采样选取适合需要按教程上手的初学者和想规范流程的研究人员。资源为单个PDF文档共1个pdf文件压缩包大小1.16MB正文步骤紧凑可直接按章节对照操作。已有230人学习浏览。除核心预处理外还补充了镶嵌工具的无缝拼接、羽化与匀色处理以及SPOT4全色与Landsat TM多光谱的Gram-Schmidt融合示例便于读者延伸掌握SAR与光学数据组合处理的常见思路。1. ENVI 对 SAR 数据的预处理从 SLC 到可分析图像到底要过几道关拿到一景 ALOS PALSAR 或 GF-3 的原始 SAR 数据直接拖进 ENVI 是没法看的满屏雪花一样的斑点噪声、斜距几何造成的透视收缩、叠掩和阴影以及没有地理坐标的像素阵列任何后续的目视解译、定量反演和变化检测都无从谈起。这 就是 SAR 预处理要解决的三个核心矛盾——斑点噪声抑制、几何投影转换和辐射一致性校正。这份资源把 ENVI SARscape 的完整预处理链路拆成了六步标准格式导入、多视处理、滤波、地理编码与辐射定标、镶嵌、融合与裁剪每一 步的参数怎么设、为什么这么设、设错了会出现什么现象都有落地的操作描述。适合刚接手 SAR 数据、被 ENVI 菜单层级绕晕的从业者也适合做过一遍但没搞懂视数和投影参数含义、只敢照抄教程的熟手。下面按实际作业顺序一步步过。2. 数据导入与多视处理先让原始数据变成能看的 SLC再定视数SAR 预处理第一步不是打开图像而是把原始分发格式转换为 ENVI 能识别的 SLC 标准数据。原始数据通常以文件集形式存在比如 ALOS PALSAR 的 FBD 模式数据包含影像文件、Leader 参数文件和 Trailer 尾部文件三者缺一不可。导入工具在 Toolbox 里逐级展开SARscape - Basic - Import Data - Standard Formats - ALOS PALSAR。这个菜单层级很多人第一次找不到SARscape 模块在 ENVI 中是一个独立工具包安装后会在 Toolbox 里新增一个 SARscape 分类需要展开到第三级才能看到 Standard Formats。2.1 标准格式导入JAXA-FBD Level 数据与 Leader 文件的对应关系打开导入面板后第一个要确认的是 Data Type。对于 ALOS PALSAR 的 FBD 模式Fine Beam Dual双极化精细波束数据类型应选择JAXA-FBD Level。这个数据类型可以从文件名推导FBD 模式的命名中通常带有 FBD 标识而 FBSFine Beam Single则对应单极化。选错类型会导致后续的元数据解析失败或 SLC 文件尺寸异常。# 以 Linux 环境为例检查 ALOS PALSAR 原始数据文件结构 ls -la /data/ALOS_PALSAR/ # 典型输出 # -rw-r--r-- 1 user user 2147483648 Jul 12 10:23 123456_001_FBD_000.IMG # -rw-r--r-- 1 user user 4096 Jul 12 10:23 123456_001_FBD_000.LEAD # -rw-r--r-- 1 user user 20480 Jul 12 10:23 123456_001_FBD_000.TRL文件结构确认后在面板中单击 Leader/Paramfile 按钮选择.LEAD文件单击 Data list 选择影像数据文件。这里有一个常被忽视的动作Output file 路径设置后软件会在输入文件名基础上自动追加标识字母例如原文件名带_FBD的会追加成_SLC表示输出为单视复数产品。最终生成三个结果文件SLC 数据文件.slc、SML 格式的元数据文件、HDR 格式的 ENVI 头文件。这三件套缺一不可后续 Multilooking 和滤波工具依赖元数据里的轨道参数和成像几何信息。导入完成后在 ENVI 中打开 SLC 文件看到的是一幅灰蒙蒙、几乎无法辨识地物的复数图像。这是正常现象——SLC 保留了完整的相位信息幅度和相位耦合在一起需要经过多视处理把相位信息平均掉一部分才能得到适合目视判读的强度图像。2.2 多视处理方位向视数 5、距离向视数 1 是怎么算出来的多视处理针对的是斑点噪声。SAR 图像产品包含大量相干斑表现为像素间的随机亮暗起伏。SAR 信号处理器使用完整的合成孔径和所有信号数据来获得最高空间分辨率但代价是斑点噪声显著。多视处理在方位向和距离向上做分辨率平均目的是抑制斑点结果是辐射分辨率提高、空间分辨率降低——这是一个需要按需取舍的权衡。打开SARscape - Basic - MultilookingInput file 选导入生成的 SLC 数据。这里要注意文件选择框默认过滤_slc后缀文件如果列表里看不到数据检查文件后缀是否为小写_slc否则需手动切换文件类型过滤条件。方位向视数 Azimuth Looks 和距离向视数 Range Looks 是核心参数。文档里给的参考值是 Azimuth Looks5、Range Looks1这个组合来自对 ALOS PALSAR FBD 数据的精确计算。公式如下地面分辨率 pixel spacing slant range / sin(incidence angle) 地距距离分辨率 pixel spacing slant range / sin(场景中心入射角) 方位向视数 地距方位向分辨率 / 地距距离分辨率具体到 ALOS PALSAR FBD 数据距离向 pixel spacing斜距分辨率已知场景中心入射角例如 38.7°那么地距距离分辨率 斜距分辨率 / sin(38.7°)。若要使得多视后的方位向分辨率与地距距离向分辨率一致就需要把方位向分辨率按比例放宽由此得到方位向视数为 5、距离向视数为 1。SARscape 面板上的 Look 按钮可以自动估算视数但自动估算基于默认参考值如果数据不是 ALOS 系列建议还是手动代入公式计算。Border Resize 选项这里不选。它的作用是检测边缘的无效值并重新计算输出图像大小但会改变图像的地理范围边界在后续地理编码时可能引入额外的拼接对齐问题。Corner Coordinates 相关的边界修正等地理编码那一步再做。执行 Start 输出后在 Display 中对比多视前后的图像可以明显看到斑点噪声被抑制地物轮廓变得可辨但整体图像会显得“糊”一些——这是空间分辨率下降的直接观感。3. 滤波与地理编码辐射定标从 SLC 到能定量分析的 SAR 图像多视处理之后斑点噪声仍有一定残余尤其是复杂地形区域。滤波这一步的作用是进一步平滑斑点同时尽量保留边缘信息。SARscape 提供两大类滤波用于单波段雷达图像的单图像滤波以及用于多时相雷达图像的时间序列滤波。本流程中先做单图像滤波在后续有长时间序列数据时才需要用到多时相滤波。3.1 Frost 滤波窗口大小 9×9 的选取逻辑与效果控制打开SARscape - Basic - Filtering - Single Image输入选择多视处理结果。SARscape 提供 8 种滤波算法Frost、Lee、Refined Lee、Gamma MAP、Boxcar、Median、Mean 和 IDAN。选择 Frost 后会自动打开 Frost Lee/Refined Lee 对话框需要设置 Azimuth Window Size 和 Range Window Size 两个窗口参数文档中的参考值是两个方向都为 9。窗口大小的物理意义窗口越大平滑力度越强但边缘模糊越严重窗口越小细节保留越好但斑点抑制不彻底。对于 ALOS PALSAR 这种 10m 级分辨率的 L 波段数据9×9 窗口在斑点抑制和边缘保持之间比较均衡。如果数据是更高分辨率的 X 波段如 TerraSAR-X 3m窗口可以缩小到 5×5因为分辨单元本身更小大窗口会过度平滑掉道路、建筑边缘等细节。Frost 滤波对窗口内的像素做加权平均权重根据局部统计特性自动调整在均匀区域平滑更强在边缘区域权重会自适应降低从而保住边缘。执行完成后对比多视与滤波后的图像斑点噪声进一步降低图像更干净但这时图像仍未带有地理坐标——只有像素的行列坐标需要用地理编码和辐射定标把它投影到真实地图坐标系中。3.2 地理编码与辐射定标投影参数、像元大小与重采样方法的取舍SAR 系统测量的是发射脉冲与返回脉冲的功率比这个比值后向散射系数被投影在斜距几何上。斜距几何意味着地面等距的物体在图像上并不等距且存在透视收缩、叠掩和阴影等几何畸变。为了与其他光学遥感数据、GIS 矢量数据对比分析必须将 SAR 数据从斜距/地距投影转换到地图投影坐标系。打开SARscape - Basic - Geocoding - Geocoding and Radiometric Calibration输入选择 Frost 滤波结果。重点参数如下参数设置值说明Cartographic System按成果需要选择目标投影坐标系如 UTM/WGS84 对应分带GRID SIZE X25输出像元大小米需参考原始分辨率设置GRID SIZE Y25与 X 保持一致避免像元变形RESAMPLEOptimal Resolution从最邻近到最优分辨率共 5 种精度递增但速度递减Radiometric Calibration勾选将后向散射转换为归一化雷达截面sigma0Radometric Normalization按需勾选辐射归一化减少地形起伏引起的辐射差异Local Incidence Angle按需勾选局部入射角校正用于山区数据Layover/Shadow按需勾选叠掩/阴影掩膜处理输出对应掩膜文件像元大小设置有一个常见误区不是越小越好。对于 ALOS PALSAR FBD 数据多视后分辨率大约 25m 左右GRID SIZE 设置为 25m 是匹配的。如果设置成 10m重采样会把原始像素放大到 2.5 倍图像看起来“变清晰”了但实际并没有新增信息反而造成几何上的像素相关性后续定量反演时统计特征会失真。重采样方法从低到高依次是Nearest Neighbour、Bilinear、Cubic Convolution、Sinc Interpolation、Optimal Resolution。Nearest 最快但会产生锯齿Cubic 和 Sinc 精度高但计算量大。Optimal Resolution 是 SARscape 的推荐选项它会根据原始分辨率和目标 GRID SIZE 自动选择最优的插值核兼顾精度与效率。在测试小范围数据时可以用 Cubic Convolution全图处理时用 Optimal Resolution 更稳妥。输出结果包含一个带有地理坐标的后向散射强度图像打开后叠加到地图坐标系中SAR 图像与实际地物位置的对应关系才能成立。这一 步是所有后续定量分析如土壤水分反演、形变测量的前提也是不可跳过的步骤。4. 镶嵌与图像融合多景拼接、接边线编辑以及不同分辨率影像融合单景 SAR 图像的地理覆盖范围有限当研究区跨越两景以上数据时需要做镶嵌处理把多景拼接成一张无缝图。同时SAR 与光学数据的融合在解译中非常常见——例如用高分辨率全色波段增强多光谱图像的空间细节。这一 章分两块先处理镶嵌的接边线再走融合流程。4.1 镶嵌匀色、接边线自动生成与手动编辑在 Toolbox 中启动Mosaicking - Seamless Mosaic。点击左上角加号添加需要镶嵌的影像数据添加完成后显示区会呈现各影像的轮廓线、重叠关系和相对位置关系。勾选 Show Preview 可以实时预览镶嵌效果。Data Ignore Value 一栏要设置透明值输入 0 表示让背景值为 0 的像素透明这样在预览中可以看到影像实际覆盖边界。如果背景值不是 0需要根据数据实际情况修改否则会出现黑色边框残留。匀色处理方面Seamless Mosaic 提供直方图匹配方法。在 Color Correction 选项中勾选 Histogram Matching有两个作用范围Overlap Area Only 只对重叠区做直方图匹配适合各景光谱差异不大的情况Entire Scene 对整景影像做匹配适合各景差异大、需要全局色彩统一的情况。对 SAR 图像而言后向散射强度受入射角和地形影响较大重叠区匹配通常够了全局匹配反而可能抹掉真实的后向散射差异。接边线是镶嵌质量的关键。选择Seamlines - Auto Generate Seamlines自动绘制接边线软件会根据重叠区的像素相似性自动寻找最优切割路径裁剪掉边缘“锯齿”和几何错位。但自动生成的结果不一定完美尤其在地形起伏大的区域接边线可能穿越强散射体如建筑群导致拼接缝明显。这时用Seamlines - Start editing seamlines手动编辑在接边处绘制多边形软件自动将多边形作为新的接边线。手动编辑的原则是尽量让接边线绕过明显的人工地物和线性特征沿河流、山谷等天然边界走。输出时切换到 Export 选项这里可以设置输出文件名、路径、格式、波段、背景值和重采样方法。背景值设置要与 Data Ignore Value 对应起来一般保持 0 值一致否则镶嵌入其他场景时会露出黑边。4.2 融合Gram-Schmidt Pan Sharpening 与传感器类型的选择图像融合的典型场景用 SPOT4 的 10m 全色波段增强 Landsat5 TM 30m 多光谱影像或者用 QuickBird 的高分辨率全色波段增强其多光谱影像。ENVI 提供的标准工具是Image Sharpening - Gram-Schmidt Pan Sharpening。第一步先分别打开低分辨率影像如 TM 多光谱和高分辨率影像如 SPOT 全色。在融合参数面板中Sensor 下拉列表包含 GeoEye-1、Gokturk-2、IKONOS、Landsat 8 OLI、QuickBird、WorldView-1/2 等常见传感器。如果影像来自列表之外的传感器选择 Unknown软件会通过统计特征自动估算融合参数。重采样方法选择 Cubic Convolution输出格式选 ENVI。传感器类型的选择会直接影响全色波段的光谱响应拟合。GS 融合的核心是用全色波段模拟低分辨率多光谱波段的组合权重Sensor 类型越匹配模拟越准确。列表里没有的设备选 Unknown融合结果可能在色彩上略有偏差但空间细节增强效果基本不受影响。执行后对比融合前后的图像可以看到多光谱图像的颜色信息与全色波段的纹理细节同时保留边缘与纹理信息明显增强。5. 裁剪的两种玩法规则分幅与 ROI 不规则裁剪裁剪是预处理流程图里看起来最简单、实际上最容易被细节绊住的一步。ENVI 的裁剪分两大类规则裁剪矩形范围以及基于 ROI 的不规则裁剪。前者的边界范围获取途径有行列号、左上右下角坐标、外部图像文件和矢量文件后者则依赖于在图像上手动绘制感兴趣区多边形适合裁剪出不规则的研究区范围。5.1 规则裁剪Spatial Subset 四种范围确定方式以 TM 影像为例打开图像并按 Linear 2% 拉伸显示然后File - Save As进入 File Selection 面板选择 Spatial Subset右侧弹出裁剪区域选择功能。确定裁剪区域有四种途径第一种使用当前可视区域单击 Use View Extent 按钮软件自动读取主窗口中显示的范围作为裁剪边界适合“所见即所裁”的快速操作。第二种通过文件确定单击 Subset By File 按钮选择一个矢量或栅格文件自动读取该文件的外包矩形作为裁剪范围。这种方式特别适合与已有矢量边界对齐的场景比如按行政区划面文件的外包矩形裁剪。第三种和第四种手动交互直接输入行列数 Columns 和 Rows或者用鼠标拖拽红色边框调整裁剪尺寸和位置适合精确控制像元数量的场合。Spatial Subset 左下角的 Spectral Subset 按钮可以选择输出波段子集。多光谱数据通常保留全部波段但如果只是做分类训练样本提取也可以只输出需要的波段以减少数据量。设置完成后选择输出路径和文件名单击 OK 即完成规则裁剪。5.2 不规则裁剪ROI 绘制、保存与 Subset Data from ROIs 参数设置不规则裁剪的典型需求裁剪出北京市老皇城二环范围内的 SAR 图像作为研究区。操作路径是在 Layer Manager 中选中文件右键选择 New Region Of Interest打开 ROI Tool 面板。在面板上点击多边形绘制按钮沿研究区边界逐点描绘双击完成绘制。可以修改 ROI Name 和 ROI Color也可以绘制多个多边形后续合并为一个裁剪掩膜。绘制完成后先保存 ROI在 ROI Tool 面板中选择File - Save as保存为 ROI 文件.roi。然后打开Toolbox - Regions of Interest - Subset Data from ROIs输入图像选择刚才的影像文件在参数面板中设置 Select Input ROIs 为刚保存的 roi1 文件Mask pixels output of ROI 选择 YesMask Background Value 设为 0。参数设置值说明Select Input ROIsroi1.roi指定裁剪掩膜的 ROI 文件Mask pixels output of ROIYes用 ROI 作为掩膜输出像素Mask Background Value0掩膜外像素填充值为 0执行完成后输出的图像是一个不规则多边形范围内的 SAR 图像掩膜外的像素全部为 0。这里需要检查一 点输出图像的像元大小、投影坐标系与掩膜文件是否一致如果不一致ROI 掩膜会被隐式重采样边缘可能出现半个像素的偏移。多数情况下在裁剪前先检查两个文件的投影和像元尺寸可以避免这个后续对齐的麻烦。6. 预处理避坑清单导入失败、视数失真与地理编码偏移的排查记录以下是我在实际操作中反复遇到的五个问题每条按现象、原因、解决三个层次展开。6.1 导入时报错“File format not recognized”现象在 Standard Formats 中选择 ALOS PALSAR 后点击 Leader/Paramfile 选择了文件但 Start 按钮置灰或直接报错提示无法识别文件格式。原因路径中包含中文字符或空格SARscape 的底层读取模块对非 ASCII 路径支持不好也可能是数据文件本身不完整比如只有 IMG 文件但 LEAD 文件缺失。解决先把所有原始数据复制到纯英文路径下例如D:\SAR_data\ALOS_PALSAR\确保文件名和目录名均为英文字符再重新导入。如果问题仍然存在检查 LEAD 文件大小是否为 0用十六进制编辑器确认文件头是否符合 ALOS 格式规范。6.2 多视后图像分辨率出现异常现象设置 Azimuth Looks5、Range Looks1但输出图像的方位向像元尺寸与距离向像元尺寸相差很远图形被拉伸变形地物形状不真实。原因计算公式中的入射角用的是场景中心值的近似值实际数据在近距端和远距端的入射角差异很大尤其是宽幅数据。另一个原因是 Border Resize 被勾选软件对无效边缘做了裁剪改变了输出范围。解决严格按文档公式代入场景中心入射角计算视数然后在多视处理后查看输出图像的像元尺寸确认方位向与距离向是否接近而不是直接照抄 5 和 1。若差异较大可尝试设置 Azimuth Looks4 或 6取更接近整数的结果。6.3 滤波后地物边缘糊成一片现象Frost 滤波执行后道路、建筑边缘都变得模糊地物边界不清晰目视效果比多视结果还差。原因滤波窗口设置过大。在初步测试时直接套用了 9×9没有根据数据的实际分辨率调整。对于噪声水平较低的数据如 FBD 模式视数已经 59×9 窗口过度平滑。解决将窗口调小到 5×5 或 7×7先在小范围内对比不同窗口的视觉效果选择边缘保持和噪声抑制的平衡点。另外Refined Lee 滤波对边缘保持效果比 Frost 更好如果边缘信息对后续解译很重要可以改用 Refined Lee。6.4 地理编码结果与光学影像偏移几十米现象地理编码和辐射定标完成后将 SAR 图像与 Landsat 光学影像叠加发现地物明显错位偏移量达到几十米。原因输出投影坐标系设置与实际需求不匹配。例如 SARscape 默认用的投影可能与光学影像的投影不一致尽管都叫 UTM但基准面不同WGS84 与 CGCS2000 的差异在局部区域可达数十米也可能 GRID SIZE 设置与原始分辨率差异过大重采样引入了几何误差。解决在 Geocoding 面板中明确查看并核对 Cartographic System 的具体参数必要时自定义投影参数确保基准面和投影带与光学影像一致。GRID SIZE 按原始分辨率匹配设置不要为了追求“高分辨率”而设置过小的像元。6.5 镶嵌输出出现明显拼接缝现象两景 SAR 影像在镶嵌输出后重叠区出现一道明显的亮暗变化边界即使做了直方图匹配仍然存在。原因两景数据的入射角差异较大导致重叠区的后向散射强度固有不同或者地理编码时辐射定标未勾选后向散射值没有归一化。解决在镶嵌之前先确认每景数据都完成了辐射定标输出的是 sigma0 值。匀色时先试 Overlap Area Only如果效果不理想再试 Entire Scene。手动编辑接边线使其绕过亮度差异明显的区域。7. 预处理完成后的第一件事逐景验证与保存处理模板整条链路跑完拿到地理编码后的 SAR 图像先别急着做分类或反演。我习惯的验证流程有三步每一步都能在几分钟内发现问题避免后续返工。第一步检查几何精度。在 ENVI 中把地理编码结果与同区域的 Landsat 或 Sentinel-2 影像叠加打开 GeoLink 功能随机选取 58 个明显地物点立交桥、水库拐角、大型建筑检查 SAR 与光学影像的偏移量。如果误差在 12 个像元内说明几何处理合格如果偏移超过 5 个像元需要复盘 GRID SIZE 和投影设置。注意一点SAR 侧视成像导致的叠掩和阴影区域天然存在几何错位验证点应选择平坦地物。第二步后向散射值抽查。在 SAR 图像上选取水面、裸地、城区三个典型地物区域记录回波强度值。水域后向散射应显著低于城区裸地介于两者之间。如果水域值异常偏高通常是辐射定标未正确执行或局部入射角校正引入了误差。第三步把参数保存为模板。SARscape 支持将处理参数保存为任务模板在 Toolbox 中右键对应工具选择 Save参数面板中的所有设置会被记录为.prm文件。下次处理相同传感器、相同模式的数据时直接加载模板即可一键复用。这比每次手动设置省事得多也避免了重复输入产生的笔误。几个月前我处理一批 20 景同一轨道的 PALSAR 数据第一次手动跑完第 5 景时才发现在地理编码时 GRID SIZE 设置错了前三景全部要重做。从那以后我每次处理新批次数据都强制先用一景数据跑完全流程、验证通过再把参数保存为模板批量执行提取失败经验成本最低的方式。希望帮到你——预处理是 SAR 分析最枯燥但最值得花时间的环节参数设置每一步都有物理意义理解了再动手坑自然少一半。本文还有配套的精品资源点击获取 ENVI+SARscape SAR数据预处理全流程：从SLC到可分析图像
简介面向遥感数据处理与InSAR学习者的操作型参考资料聚焦ENVI中SARscape模块对ALOS PALSAR等数据的预处理全流程。内容覆盖数据导入、单视复数产品多视处理、斑点噪声滤波、地理编码与辐射定标并按步骤给出参数设置建议如视数估算、滤波窗口、投影与重采样选取适合需要按教程上手的初学者和想规范流程的研究人员。资源为单个PDF文档共1个pdf文件压缩包大小1.16MB正文步骤紧凑可直接按章节对照操作。已有230人学习浏览。除核心预处理外还补充了镶嵌工具的无缝拼接、羽化与匀色处理以及SPOT4全色与Landsat TM多光谱的Gram-Schmidt融合示例便于读者延伸掌握SAR与光学数据组合处理的常见思路。1. ENVI 对 SAR 数据的预处理从 SLC 到可分析图像到底要过几道关拿到一景 ALOS PALSAR 或 GF-3 的原始 SAR 数据直接拖进 ENVI 是没法看的满屏雪花一样的斑点噪声、斜距几何造成的透视收缩、叠掩和阴影以及没有地理坐标的像素阵列任何后续的目视解译、定量反演和变化检测都无从谈起。这 就是 SAR 预处理要解决的三个核心矛盾——斑点噪声抑制、几何投影转换和辐射一致性校正。这份资源把 ENVI SARscape 的完整预处理链路拆成了六步标准格式导入、多视处理、滤波、地理编码与辐射定标、镶嵌、融合与裁剪每一 步的参数怎么设、为什么这么设、设错了会出现什么现象都有落地的操作描述。适合刚接手 SAR 数据、被 ENVI 菜单层级绕晕的从业者也适合做过一遍但没搞懂视数和投影参数含义、只敢照抄教程的熟手。下面按实际作业顺序一步步过。2. 数据导入与多视处理先让原始数据变成能看的 SLC再定视数SAR 预处理第一步不是打开图像而是把原始分发格式转换为 ENVI 能识别的 SLC 标准数据。原始数据通常以文件集形式存在比如 ALOS PALSAR 的 FBD 模式数据包含影像文件、Leader 参数文件和 Trailer 尾部文件三者缺一不可。导入工具在 Toolbox 里逐级展开SARscape - Basic - Import Data - Standard Formats - ALOS PALSAR。这个菜单层级很多人第一次找不到SARscape 模块在 ENVI 中是一个独立工具包安装后会在 Toolbox 里新增一个 SARscape 分类需要展开到第三级才能看到 Standard Formats。2.1 标准格式导入JAXA-FBD Level 数据与 Leader 文件的对应关系打开导入面板后第一个要确认的是 Data Type。对于 ALOS PALSAR 的 FBD 模式Fine Beam Dual双极化精细波束数据类型应选择JAXA-FBD Level。这个数据类型可以从文件名推导FBD 模式的命名中通常带有 FBD 标识而 FBSFine Beam Single则对应单极化。选错类型会导致后续的元数据解析失败或 SLC 文件尺寸异常。# 以 Linux 环境为例检查 ALOS PALSAR 原始数据文件结构 ls -la /data/ALOS_PALSAR/ # 典型输出 # -rw-r--r-- 1 user user 2147483648 Jul 12 10:23 123456_001_FBD_000.IMG # -rw-r--r-- 1 user user 4096 Jul 12 10:23 123456_001_FBD_000.LEAD # -rw-r--r-- 1 user user 20480 Jul 12 10:23 123456_001_FBD_000.TRL文件结构确认后在面板中单击 Leader/Paramfile 按钮选择.LEAD文件单击 Data list 选择影像数据文件。这里有一个常被忽视的动作Output file 路径设置后软件会在输入文件名基础上自动追加标识字母例如原文件名带_FBD的会追加成_SLC表示输出为单视复数产品。最终生成三个结果文件SLC 数据文件.slc、SML 格式的元数据文件、HDR 格式的 ENVI 头文件。这三件套缺一不可后续 Multilooking 和滤波工具依赖元数据里的轨道参数和成像几何信息。导入完成后在 ENVI 中打开 SLC 文件看到的是一幅灰蒙蒙、几乎无法辨识地物的复数图像。这是正常现象——SLC 保留了完整的相位信息幅度和相位耦合在一起需要经过多视处理把相位信息平均掉一部分才能得到适合目视判读的强度图像。2.2 多视处理方位向视数 5、距离向视数 1 是怎么算出来的多视处理针对的是斑点噪声。SAR 图像产品包含大量相干斑表现为像素间的随机亮暗起伏。SAR 信号处理器使用完整的合成孔径和所有信号数据来获得最高空间分辨率但代价是斑点噪声显著。多视处理在方位向和距离向上做分辨率平均目的是抑制斑点结果是辐射分辨率提高、空间分辨率降低——这是一个需要按需取舍的权衡。打开SARscape - Basic - MultilookingInput file 选导入生成的 SLC 数据。这里要注意文件选择框默认过滤_slc后缀文件如果列表里看不到数据检查文件后缀是否为小写_slc否则需手动切换文件类型过滤条件。方位向视数 Azimuth Looks 和距离向视数 Range Looks 是核心参数。文档里给的参考值是 Azimuth Looks5、Range Looks1这个组合来自对 ALOS PALSAR FBD 数据的精确计算。公式如下地面分辨率 pixel spacing slant range / sin(incidence angle) 地距距离分辨率 pixel spacing slant range / sin(场景中心入射角) 方位向视数 地距方位向分辨率 / 地距距离分辨率具体到 ALOS PALSAR FBD 数据距离向 pixel spacing斜距分辨率已知场景中心入射角例如 38.7°那么地距距离分辨率 斜距分辨率 / sin(38.7°)。若要使得多视后的方位向分辨率与地距距离向分辨率一致就需要把方位向分辨率按比例放宽由此得到方位向视数为 5、距离向视数为 1。SARscape 面板上的 Look 按钮可以自动估算视数但自动估算基于默认参考值如果数据不是 ALOS 系列建议还是手动代入公式计算。Border Resize 选项这里不选。它的作用是检测边缘的无效值并重新计算输出图像大小但会改变图像的地理范围边界在后续地理编码时可能引入额外的拼接对齐问题。Corner Coordinates 相关的边界修正等地理编码那一步再做。执行 Start 输出后在 Display 中对比多视前后的图像可以明显看到斑点噪声被抑制地物轮廓变得可辨但整体图像会显得“糊”一些——这是空间分辨率下降的直接观感。3. 滤波与地理编码辐射定标从 SLC 到能定量分析的 SAR 图像多视处理之后斑点噪声仍有一定残余尤其是复杂地形区域。滤波这一步的作用是进一步平滑斑点同时尽量保留边缘信息。SARscape 提供两大类滤波用于单波段雷达图像的单图像滤波以及用于多时相雷达图像的时间序列滤波。本流程中先做单图像滤波在后续有长时间序列数据时才需要用到多时相滤波。3.1 Frost 滤波窗口大小 9×9 的选取逻辑与效果控制打开SARscape - Basic - Filtering - Single Image输入选择多视处理结果。SARscape 提供 8 种滤波算法Frost、Lee、Refined Lee、Gamma MAP、Boxcar、Median、Mean 和 IDAN。选择 Frost 后会自动打开 Frost Lee/Refined Lee 对话框需要设置 Azimuth Window Size 和 Range Window Size 两个窗口参数文档中的参考值是两个方向都为 9。窗口大小的物理意义窗口越大平滑力度越强但边缘模糊越严重窗口越小细节保留越好但斑点抑制不彻底。对于 ALOS PALSAR 这种 10m 级分辨率的 L 波段数据9×9 窗口在斑点抑制和边缘保持之间比较均衡。如果数据是更高分辨率的 X 波段如 TerraSAR-X 3m窗口可以缩小到 5×5因为分辨单元本身更小大窗口会过度平滑掉道路、建筑边缘等细节。Frost 滤波对窗口内的像素做加权平均权重根据局部统计特性自动调整在均匀区域平滑更强在边缘区域权重会自适应降低从而保住边缘。执行完成后对比多视与滤波后的图像斑点噪声进一步降低图像更干净但这时图像仍未带有地理坐标——只有像素的行列坐标需要用地理编码和辐射定标把它投影到真实地图坐标系中。3.2 地理编码与辐射定标投影参数、像元大小与重采样方法的取舍SAR 系统测量的是发射脉冲与返回脉冲的功率比这个比值后向散射系数被投影在斜距几何上。斜距几何意味着地面等距的物体在图像上并不等距且存在透视收缩、叠掩和阴影等几何畸变。为了与其他光学遥感数据、GIS 矢量数据对比分析必须将 SAR 数据从斜距/地距投影转换到地图投影坐标系。打开SARscape - Basic - Geocoding - Geocoding and Radiometric Calibration输入选择 Frost 滤波结果。重点参数如下参数设置值说明Cartographic System按成果需要选择目标投影坐标系如 UTM/WGS84 对应分带GRID SIZE X25输出像元大小米需参考原始分辨率设置GRID SIZE Y25与 X 保持一致避免像元变形RESAMPLEOptimal Resolution从最邻近到最优分辨率共 5 种精度递增但速度递减Radiometric Calibration勾选将后向散射转换为归一化雷达截面sigma0Radometric Normalization按需勾选辐射归一化减少地形起伏引起的辐射差异Local Incidence Angle按需勾选局部入射角校正用于山区数据Layover/Shadow按需勾选叠掩/阴影掩膜处理输出对应掩膜文件像元大小设置有一个常见误区不是越小越好。对于 ALOS PALSAR FBD 数据多视后分辨率大约 25m 左右GRID SIZE 设置为 25m 是匹配的。如果设置成 10m重采样会把原始像素放大到 2.5 倍图像看起来“变清晰”了但实际并没有新增信息反而造成几何上的像素相关性后续定量反演时统计特征会失真。重采样方法从低到高依次是Nearest Neighbour、Bilinear、Cubic Convolution、Sinc Interpolation、Optimal Resolution。Nearest 最快但会产生锯齿Cubic 和 Sinc 精度高但计算量大。Optimal Resolution 是 SARscape 的推荐选项它会根据原始分辨率和目标 GRID SIZE 自动选择最优的插值核兼顾精度与效率。在测试小范围数据时可以用 Cubic Convolution全图处理时用 Optimal Resolution 更稳妥。输出结果包含一个带有地理坐标的后向散射强度图像打开后叠加到地图坐标系中SAR 图像与实际地物位置的对应关系才能成立。这一 步是所有后续定量分析如土壤水分反演、形变测量的前提也是不可跳过的步骤。4. 镶嵌与图像融合多景拼接、接边线编辑以及不同分辨率影像融合单景 SAR 图像的地理覆盖范围有限当研究区跨越两景以上数据时需要做镶嵌处理把多景拼接成一张无缝图。同时SAR 与光学数据的融合在解译中非常常见——例如用高分辨率全色波段增强多光谱图像的空间细节。这一 章分两块先处理镶嵌的接边线再走融合流程。4.1 镶嵌匀色、接边线自动生成与手动编辑在 Toolbox 中启动Mosaicking - Seamless Mosaic。点击左上角加号添加需要镶嵌的影像数据添加完成后显示区会呈现各影像的轮廓线、重叠关系和相对位置关系。勾选 Show Preview 可以实时预览镶嵌效果。Data Ignore Value 一栏要设置透明值输入 0 表示让背景值为 0 的像素透明这样在预览中可以看到影像实际覆盖边界。如果背景值不是 0需要根据数据实际情况修改否则会出现黑色边框残留。匀色处理方面Seamless Mosaic 提供直方图匹配方法。在 Color Correction 选项中勾选 Histogram Matching有两个作用范围Overlap Area Only 只对重叠区做直方图匹配适合各景光谱差异不大的情况Entire Scene 对整景影像做匹配适合各景差异大、需要全局色彩统一的情况。对 SAR 图像而言后向散射强度受入射角和地形影响较大重叠区匹配通常够了全局匹配反而可能抹掉真实的后向散射差异。接边线是镶嵌质量的关键。选择Seamlines - Auto Generate Seamlines自动绘制接边线软件会根据重叠区的像素相似性自动寻找最优切割路径裁剪掉边缘“锯齿”和几何错位。但自动生成的结果不一定完美尤其在地形起伏大的区域接边线可能穿越强散射体如建筑群导致拼接缝明显。这时用Seamlines - Start editing seamlines手动编辑在接边处绘制多边形软件自动将多边形作为新的接边线。手动编辑的原则是尽量让接边线绕过明显的人工地物和线性特征沿河流、山谷等天然边界走。输出时切换到 Export 选项这里可以设置输出文件名、路径、格式、波段、背景值和重采样方法。背景值设置要与 Data Ignore Value 对应起来一般保持 0 值一致否则镶嵌入其他场景时会露出黑边。4.2 融合Gram-Schmidt Pan Sharpening 与传感器类型的选择图像融合的典型场景用 SPOT4 的 10m 全色波段增强 Landsat5 TM 30m 多光谱影像或者用 QuickBird 的高分辨率全色波段增强其多光谱影像。ENVI 提供的标准工具是Image Sharpening - Gram-Schmidt Pan Sharpening。第一步先分别打开低分辨率影像如 TM 多光谱和高分辨率影像如 SPOT 全色。在融合参数面板中Sensor 下拉列表包含 GeoEye-1、Gokturk-2、IKONOS、Landsat 8 OLI、QuickBird、WorldView-1/2 等常见传感器。如果影像来自列表之外的传感器选择 Unknown软件会通过统计特征自动估算融合参数。重采样方法选择 Cubic Convolution输出格式选 ENVI。传感器类型的选择会直接影响全色波段的光谱响应拟合。GS 融合的核心是用全色波段模拟低分辨率多光谱波段的组合权重Sensor 类型越匹配模拟越准确。列表里没有的设备选 Unknown融合结果可能在色彩上略有偏差但空间细节增强效果基本不受影响。执行后对比融合前后的图像可以看到多光谱图像的颜色信息与全色波段的纹理细节同时保留边缘与纹理信息明显增强。5. 裁剪的两种玩法规则分幅与 ROI 不规则裁剪裁剪是预处理流程图里看起来最简单、实际上最容易被细节绊住的一步。ENVI 的裁剪分两大类规则裁剪矩形范围以及基于 ROI 的不规则裁剪。前者的边界范围获取途径有行列号、左上右下角坐标、外部图像文件和矢量文件后者则依赖于在图像上手动绘制感兴趣区多边形适合裁剪出不规则的研究区范围。5.1 规则裁剪Spatial Subset 四种范围确定方式以 TM 影像为例打开图像并按 Linear 2% 拉伸显示然后File - Save As进入 File Selection 面板选择 Spatial Subset右侧弹出裁剪区域选择功能。确定裁剪区域有四种途径第一种使用当前可视区域单击 Use View Extent 按钮软件自动读取主窗口中显示的范围作为裁剪边界适合“所见即所裁”的快速操作。第二种通过文件确定单击 Subset By File 按钮选择一个矢量或栅格文件自动读取该文件的外包矩形作为裁剪范围。这种方式特别适合与已有矢量边界对齐的场景比如按行政区划面文件的外包矩形裁剪。第三种和第四种手动交互直接输入行列数 Columns 和 Rows或者用鼠标拖拽红色边框调整裁剪尺寸和位置适合精确控制像元数量的场合。Spatial Subset 左下角的 Spectral Subset 按钮可以选择输出波段子集。多光谱数据通常保留全部波段但如果只是做分类训练样本提取也可以只输出需要的波段以减少数据量。设置完成后选择输出路径和文件名单击 OK 即完成规则裁剪。5.2 不规则裁剪ROI 绘制、保存与 Subset Data from ROIs 参数设置不规则裁剪的典型需求裁剪出北京市老皇城二环范围内的 SAR 图像作为研究区。操作路径是在 Layer Manager 中选中文件右键选择 New Region Of Interest打开 ROI Tool 面板。在面板上点击多边形绘制按钮沿研究区边界逐点描绘双击完成绘制。可以修改 ROI Name 和 ROI Color也可以绘制多个多边形后续合并为一个裁剪掩膜。绘制完成后先保存 ROI在 ROI Tool 面板中选择File - Save as保存为 ROI 文件.roi。然后打开Toolbox - Regions of Interest - Subset Data from ROIs输入图像选择刚才的影像文件在参数面板中设置 Select Input ROIs 为刚保存的 roi1 文件Mask pixels output of ROI 选择 YesMask Background Value 设为 0。参数设置值说明Select Input ROIsroi1.roi指定裁剪掩膜的 ROI 文件Mask pixels output of ROIYes用 ROI 作为掩膜输出像素Mask Background Value0掩膜外像素填充值为 0执行完成后输出的图像是一个不规则多边形范围内的 SAR 图像掩膜外的像素全部为 0。这里需要检查一 点输出图像的像元大小、投影坐标系与掩膜文件是否一致如果不一致ROI 掩膜会被隐式重采样边缘可能出现半个像素的偏移。多数情况下在裁剪前先检查两个文件的投影和像元尺寸可以避免这个后续对齐的麻烦。6. 预处理避坑清单导入失败、视数失真与地理编码偏移的排查记录以下是我在实际操作中反复遇到的五个问题每条按现象、原因、解决三个层次展开。6.1 导入时报错“File format not recognized”现象在 Standard Formats 中选择 ALOS PALSAR 后点击 Leader/Paramfile 选择了文件但 Start 按钮置灰或直接报错提示无法识别文件格式。原因路径中包含中文字符或空格SARscape 的底层读取模块对非 ASCII 路径支持不好也可能是数据文件本身不完整比如只有 IMG 文件但 LEAD 文件缺失。解决先把所有原始数据复制到纯英文路径下例如D:\SAR_data\ALOS_PALSAR\确保文件名和目录名均为英文字符再重新导入。如果问题仍然存在检查 LEAD 文件大小是否为 0用十六进制编辑器确认文件头是否符合 ALOS 格式规范。6.2 多视后图像分辨率出现异常现象设置 Azimuth Looks5、Range Looks1但输出图像的方位向像元尺寸与距离向像元尺寸相差很远图形被拉伸变形地物形状不真实。原因计算公式中的入射角用的是场景中心值的近似值实际数据在近距端和远距端的入射角差异很大尤其是宽幅数据。另一个原因是 Border Resize 被勾选软件对无效边缘做了裁剪改变了输出范围。解决严格按文档公式代入场景中心入射角计算视数然后在多视处理后查看输出图像的像元尺寸确认方位向与距离向是否接近而不是直接照抄 5 和 1。若差异较大可尝试设置 Azimuth Looks4 或 6取更接近整数的结果。6.3 滤波后地物边缘糊成一片现象Frost 滤波执行后道路、建筑边缘都变得模糊地物边界不清晰目视效果比多视结果还差。原因滤波窗口设置过大。在初步测试时直接套用了 9×9没有根据数据的实际分辨率调整。对于噪声水平较低的数据如 FBD 模式视数已经 59×9 窗口过度平滑。解决将窗口调小到 5×5 或 7×7先在小范围内对比不同窗口的视觉效果选择边缘保持和噪声抑制的平衡点。另外Refined Lee 滤波对边缘保持效果比 Frost 更好如果边缘信息对后续解译很重要可以改用 Refined Lee。6.4 地理编码结果与光学影像偏移几十米现象地理编码和辐射定标完成后将 SAR 图像与 Landsat 光学影像叠加发现地物明显错位偏移量达到几十米。原因输出投影坐标系设置与实际需求不匹配。例如 SARscape 默认用的投影可能与光学影像的投影不一致尽管都叫 UTM但基准面不同WGS84 与 CGCS2000 的差异在局部区域可达数十米也可能 GRID SIZE 设置与原始分辨率差异过大重采样引入了几何误差。解决在 Geocoding 面板中明确查看并核对 Cartographic System 的具体参数必要时自定义投影参数确保基准面和投影带与光学影像一致。GRID SIZE 按原始分辨率匹配设置不要为了追求“高分辨率”而设置过小的像元。6.5 镶嵌输出出现明显拼接缝现象两景 SAR 影像在镶嵌输出后重叠区出现一道明显的亮暗变化边界即使做了直方图匹配仍然存在。原因两景数据的入射角差异较大导致重叠区的后向散射强度固有不同或者地理编码时辐射定标未勾选后向散射值没有归一化。解决在镶嵌之前先确认每景数据都完成了辐射定标输出的是 sigma0 值。匀色时先试 Overlap Area Only如果效果不理想再试 Entire Scene。手动编辑接边线使其绕过亮度差异明显的区域。7. 预处理完成后的第一件事逐景验证与保存处理模板整条链路跑完拿到地理编码后的 SAR 图像先别急着做分类或反演。我习惯的验证流程有三步每一步都能在几分钟内发现问题避免后续返工。第一步检查几何精度。在 ENVI 中把地理编码结果与同区域的 Landsat 或 Sentinel-2 影像叠加打开 GeoLink 功能随机选取 58 个明显地物点立交桥、水库拐角、大型建筑检查 SAR 与光学影像的偏移量。如果误差在 12 个像元内说明几何处理合格如果偏移超过 5 个像元需要复盘 GRID SIZE 和投影设置。注意一点SAR 侧视成像导致的叠掩和阴影区域天然存在几何错位验证点应选择平坦地物。第二步后向散射值抽查。在 SAR 图像上选取水面、裸地、城区三个典型地物区域记录回波强度值。水域后向散射应显著低于城区裸地介于两者之间。如果水域值异常偏高通常是辐射定标未正确执行或局部入射角校正引入了误差。第三步把参数保存为模板。SARscape 支持将处理参数保存为任务模板在 Toolbox 中右键对应工具选择 Save参数面板中的所有设置会被记录为.prm文件。下次处理相同传感器、相同模式的数据时直接加载模板即可一键复用。这比每次手动设置省事得多也避免了重复输入产生的笔误。几个月前我处理一批 20 景同一轨道的 PALSAR 数据第一次手动跑完第 5 景时才发现在地理编码时 GRID SIZE 设置错了前三景全部要重做。从那以后我每次处理新批次数据都强制先用一景数据跑完全流程、验证通过再把参数保存为模板批量执行提取失败经验成本最低的方式。希望帮到你——预处理是 SAR 分析最枯燥但最值得花时间的环节参数设置每一步都有物理意义理解了再动手坑自然少一半。本文还有配套的精品资源点击获取