SARscape-常见问题 - 图文

SARscape Q &A

常见问题

Q. – 要运行SARscape软件，电脑上必须装ENVI并且有ENVI的许可吗？ A. – 是的.

Q. - SARscape支持的操作系统有哪些？

A. - WINDOWS (XP, SP2, Vista, 7) 32 和64位, LINUX 64 位.

Q. – SARscape软件采用的是硬件狗加密的形式，如何确定许可里面有SARscape的哪些模块？ A. – 可以双击SARscape安装包中的\文件, 会显示相关的加密狗的信息 (如：许可期限、可用模块等).

Q. – 要运行SARscape软件，ENVI至少是哪个版本？ A. –ENVI从4.3版本开始，支持所有SARscape的功能

Q. – SARscape软件LINUX版本的，要比WINDOWS版本的好用吗？

A. – 是的，在绝大部分的处理上，在相同的硬件环境下，LINUX平台要比WINDOWS平台运行效率高. 而且，比如说大数据量的处理（好多个G），在WINDOWS平台下可能会报错，在LINIX平台下可以顺利的运行。

Q. – 要运行SARscape，是否必须有IDL环境？ A. – 不需要.

Q. – 在SARscape中生成的数据，可以用ENVI的功能进行处理，然后再在SARscape中处理和使用吗？

A. – 可以，用专门的数据导入功能，该功能可以将开始生成的SARscape结果的信息保存，不过要注意的是，在ENVI中处理之后的数据，栅格数据的参数（如行列号）要保持不变。

Q. – SARscape中可以自定义SAR传感器吗，读取自定义传感器SAR数据和在SARscape和ENVI中处理？

A. – 可以。用户可以根据在“custom”自定义传感器中提供的通用的数据格式来读取用户自己的SAR数据。在做聚焦处理的时候，SLC文件需要准备成SARscape标准的数据格式：浮点型的复数矩阵数据文件，以及两个头文件：.sml文件和.hdr文件。其中，.sml文件是SAR的参数说明（SARscape用来读取SAR数据的参数），.hdr文件是基本的栅格数据说明（ENVI读取数据的参数）。

SARscape头文件（.sml文件）中的字段\是自定义传感器的标识，可以让软件在所有处理时能正确的识别数据类型和进行处理（如辐射定标）.

Q. –如何处理Error信息？ A. –错误由很多原因导致，有可能是输入了错误的数据、不正确的参数、或者软件本身的bug，

当有错误信息的时候，您可以： 1) 查看错误信息的内容；

2) 查看SARscape帮助中提供的该步操作的内容，核对\Values\默认参数中的设置，从而检查是否由于错误的设置导致的错误；

3) 当不能确定错误原因时，将程序自动保存的错误报告发送给当地的SARscape技术支持

Q. – 每次操作的参数，输入/输出的文件列表，在SARscape4.0版本中是存放在输入文目录中，在4.1版本中，找不到这些文件了，是存在别的路径下了吗？

A. – 从4.2版本开始，很多处理的参数都保存到SARscape的默认路径下了，这些临时文件会被不断更新，存放的是最后一次的结果。

Q. – 运行该软件有屏幕分辨率的限制吗？

A. –分辨率至少是1024X768. 再小的话可能一些参数面板和按钮不能正常显示。

Q. – SARscape对硬件环境有哪些要求？

A. – 硬件环境至少是：内存：2G及以上、硬盘空间：50G及以上、处理器：1GHz及以上、当可用内存比较小的时候，软件可以运行，但是处理时间会很长，有一些大数据量的处理会因为内存不足而报错，如相位解缠、ScanSAR的处理等

Q. - SARscape支持64位操作系统吗？

A. –支持：WINDOWS 32和64 位、LINUX 64位。

Q. – SARscape是否能利用多CPU环境？

A. –利用多核资源优化数据处理是一项正在进行和不断改进的工作。有些功能已经能利用多核的优势，提高处理效率。

Q. – 为了保证处理效果，建议选择哪种SAR数据类型/格式？ A. –建议选择带有原始几何信息的（如斜距产品）数据产品，这种数据典型的是单视复数（SLC）产品，这种数据无论是几何信息还是辐射信息，都保留的是最原始的。相反的，应该尽量避免选择地距产品（如ERS PRI, Radarsat-1 SGF, ENVISAT ASAR IMP or APP等），因为很多情况下，这种数据的几何信息和辐射信息已经削减了，而且，一般已经用参考椭球高程做过了局部的辐射校正。此外，SLC数据包含强度和相位信息（而地距产品只包含多视后的强度信息），不仅可以做振幅图像的解译，也可以做干涉测量。

Q. – 在SARscape中处理的数据，会自动的生成8位的TIFF数据，如何可以设置让自动生成32位或者16位的TIFF数据？

A. –SARscape中自动的生成8位的TIFF数据，主要用来直接打开数据可视化查看，所以SARscape对TIFF数据做了拉伸，而数据的所有信息（如地理编码和辐射定标之后的结果）保存在数据结果中（如地理编码后的结果_geo），是32位的浮点型数据，和两个头文件，用户可以用ENVI的输出数据功能，另存为TIFF或GeoTIFF格式（16或32位）或者其他数据格式。

Q. –SARscape支持NITF格式吗？

A. – SARscape直接支持的NITF数据格式目前只有SAR-Lupe数据，其他的NITF格式的数据可

以先在ENVI中读取，再通过数据导入成SARscape的数据格式。

Q. –“0-多普勒效应”和“非0-多普勒效应”数据的区别是什么？

A. - “0-多普勒效应”几何的意思是象元的多普勒效应是被压缩的，在聚焦的时候，象元的位置和多普勒抛物线的最高点是一致的；“非0-多普勒效应”（或斜距几何）的意思是，在聚焦的时候，象元的位置是和后项散射系数的最高点一致的。需要注意的是，这两种情况，都保留了原始的相位和光谱信息。

所有的SARscape功能都支持这两种几何形式的数据，不同几何形式的数据也可以组成一对数据来处理，如干涉测量中的主从像对。

Q. – SARscape可以对其他软件中处理的数据结果进行进一步处理吗？比如在ROI_PAC中处理得到的干涉图，能在SARscape中进行后面的处理吗？

A. – 在SARscape中可以读入这样的结果数据，比如一般的二进制格式或者tiff格式的，但是SAR数据原始的参数信息，如轨道信息和其他存在SARscape头文件中的信息，还有相关的信息（比如主-从振幅信息和其他的）都丢失了，所以就不能做后续的处理了。

Q. – 为什么在处理中点击输入数据的时候，没找到数据文件？

A. –在直接找不到输入数据文件的时候，点击文件类型\，只有文件符合相应的扩展名的时候才能自动列出来，在文件类型后面选择*号，可以列出所有的数据文件。

Q. – 在处理中会生成一些不必须的输出文件，能否不将这些文件输出到硬盘里？

A. – 一些处理会生成一些临时文件，如果想在处理结束后自动的删除掉这些临时文件，可以在Default Values>General 面板选中\选项。

Q. – 为什么每一步会自动生成TIFF数据？

A. – 输出的TIFF数据目的是将结果可视化，这些文件比较小，读取起来方便（不论是振幅还是相位数据），如果不想输出TIFF格式的结果，可以在默认参数的面板中把Generate Tiff的功能关闭掉。

Q. – 我正在处理PALSAR FBS和FBD数据，后者的距离向分辨率较低（大概是前者的2倍），想知道SARascape是如何处理不同分辨率的数据的。

A. – 对这种情况，软件在做配准的时候会做处理，配准后的结果分辨率是一致的，基础处理和干涉处理都可以用FBS和FBD的数据，当其中一个作为参考数据的时候，另外一个就会重采样到参考图像上，比如输入FBS或者FBD作为参考数据，重采样意味着对从影像的过采样或者采样不足，在处理面板中的视数设置会根据参考影像的情况来设置。

Q. – 为什么输入TerraSAR-X的SLC数据，输出的数据是两个波段的？

A. – 一些短的复数格式，如TerraSAR-X和Tandem-X，输入的时候会保留原始的数据（16位的实部和16位的虚部），这样存储可以减少输出数据的数据量， ENVI中可以正常打开和显示这样的单波段数据。短的复数数据包含两个波段（分别是实部和虚部），用Tools>Transform Raster Data功能，可以将这样的数据转换成单波段的浮点型的数据（Data Type>Complex 32）。

Q. – 使用双星（Bistatic） 模式的Tandem-X数据有什么优点？

A. - Tandem-X和TerraSAR-X是一对姐妹星，可以获取到精确的TS-X的数据，这一对卫星的

运行模式为：其中一个发射和接收电磁波，另一个只是接收电磁波，即双星模式，这样可以采集到同步信号，不用做时间去相关（对生产DEM非常有优势）。

Q. – 支持的DEM有哪些？如何读取这些DEM数据？

A. – 在数据源方面（如数字化的地图、遥感数据中获取的DEM、航空摄影测量产品），任何DEM都可以。考虑到DEM的可靠性对结果的质量有很大的影响，可以在数据输入的面板选择数据单位“Data Units”，将数据准备成二进制的数据（比如BIL形式的）或者TIFF数据。

Q. –输入TIFF格式的DEM报错的问题

A. – 不是所有的TIFF格式都是SARscape直接支持的。当DEM是以分块的形式存储的TIFF，是不能直接输入的，可以先将分块存储的TIFF数据转成ENVI标准格式（在ENVI中可以做），然后再在SARscape中导入ENVI标准数据，生成SARScape数据的头文件。

Q. – 可以用ASTER的DEM吗？这种DEM数据可以用于SARscape中的数据处理吗？

A. –SARscape可以支持任何的DEM数据。我们做过测试，对ASTER DEM和SRTM-3的DEM，发现在绝大多数情况下，SRTM都要优于ASTER的DEM（在ASTER DEM上发现了很多人工制作的地方）。而两个数据的区别在于：

- SRTM 的产品可以通过工具，利用SAR数据的覆盖范围自动的提取和镶嵌；

- ASTER 的产品, 和其他的DEM产品不能利用\Elevation Model Extraction\功能得到，必须用通用数据输入的功能可以生成。

Q. –下载路径有默认的位置吗？比如SRTM-3或GTOPO30数据的下载。

A. – 内置下载链接的每种版本的DEM，在下载的时候会自动生成相应的路径，为\SA\\SARscape x.x.xxx\\work\，可以在设置里面修改默认的路径。

Q. –有无默认的存储精确轨道参数的路径？软件可以自动读取和使用

A. – 用户提供的精确的轨道参数，当软件在读取数据的时候可以自动提取和使用的，可以在Default Values>Directories下面设置路径。该功能只是针对ENVISAT ASAR数据的。

Q. –可以在图像的某个子区域使用更新了的轨道参数吗？

A. – 是的，更新轨道参数的功能可以对整景影像也可以对一部分图像做，如果要对一个子区域做的话，必须事先用Sample Selection功能提取出子区域。

Q. – 在处理SAR数据的时候可以用更精确的轨道参数吗？有什么主要的优点？ A. – ERS, ENVISAT ASAR 和 RADARSAT-2的数据产品支持更新的轨道参数。对于ERS，SARscape有功能更新原始的轨道参数；对于ENVISAT ASAR，可以在1级数据的导入或0级数据聚焦处理的时候输入更精确的轨道参数；对于RADARSAT-2数据，当输入数据的时候可以输入精确的轨道信息。

利用精确的轨道信息的优点主要体现在干涉测量的时候，特别是对去除轨道相关的干涉条纹，典型的是在输出的轨道位置上有错误的情况（即使是约几厘米的非常小的错误）。需要注意的是，使用精确的轨道参数可以提高位置向量的精度，但不会对那些影像定位精度参数有什么修改，如获取时间、斜距、脉冲频率（PRF）、载波频率等，当这些参数有误的时候，可以通过已知的地面控制点（在地理编码和辐射定标、干涉图去平、轨道精炼的时候）来修正，或是用轨道校正工具（Orbit Correction）自动计算和修正

Q. – 在多视的时候自动计算出来的系数得到的图像分辨率要比原始象元单视的分辨率粗，多视的系数是怎么计算的？

A. –“look”的按钮可以计算出最合适的方位向和距离向的系数，是用原始数据的距离向/方位向和入射角，以及在默认参数Default Values>General里面设置的制图分辨率“Cartographic Grid Size”一起计算出的。制图分辨率“Cartographic Grid Size”设置的越大，计算出的距离向和方位向的视数就越大，多视后的图像分辨率就越粗。

Q. – 聚束式模式获取的数据，多视的计算标准会不会变化？因为这种模式是要得到高分辨率的小范围的数据。

A. – 聚束模式获取的数据是在方位向进行了聚束，对多视的计算没有影响，因为入射角（和相应的距离向分辨率）无论是对聚束式还是扫描式的获取模式都是一样计算的。需要注意的是，在距离向上入射角的变化（对任何一种获取方式而言），会引起地距分辨率的变化，对该数据（不管是在几何上还是辐射上）最佳的计算方法是用Geocoding>Optimal Resolution resampling工具，对单视SL数据做。

Q. – 软件在处理条带式（Stripmap）和聚束式（Spotlight）的数据时有无区别？在处理过程中有需要区分考虑的情况吗？

A. – 从数据处理的角度来看，条带式和聚束式这两种获取数据的方式最主要的区别在于方位向上多普勒图心的可变性（对前者来说影响很小，但对后者有很大的影响），数据读取的时候会自动的设别数据的获取模式，然后将所有处理与多普勒图心变化的参数信息添加到SARscape头文件中。有了这些注释信息，数据就会用相同的方式正常的处理，如使用多普勒图心的信息来调节方位向的内插滤波器（在配准的时候），或者是做局部带宽滤波的时候 总之，所有的处理过程对用户来说都是易懂的，没有因为不同的模式而做特殊的处理。

Q. – 可以对同一个传感器获取的有重叠区域的两景数据做图像配准吗？ A. – 当然可以

Q. – 在做配准的时候只是用轨道参数来做的吗？

A. – 如果只想用轨道参数来校正，必须关闭 \和\默认参数，初始化参数必须选上 \项。需要注意的是，如果选择了\Accuracy\选项，程序会用9窗口的振幅交叉检验方法自动执行轨道的可靠性检查。如果\选项关闭, 就会基于轨道的偏移量，在距离向和方位向的常数项来进行轨道精炼。

Q. – 用基础模块中的配准功能和干涉模块里的工具生成的配准的强度图（_pwr 和_rsp）有什么区别？

A. - 基础模块中的配准功能和干涉模块里相应的工具是非常类似的，但是也不完全一样。 干涉叠加模块中考虑了振幅的相关性和相位（对slc数据来说），以得到非常精确的偏移参数，配准精度必须保证达到1/10象元，这样可以生成高质量的干涉图。在基础模块中的Basic Module>Feature Extraction>Coherence的配准功能，也是同样的方法。反过来说，基础模块中对于强度数据的配准，半个象元的精度已经足够了，所以基础模块中的配准功能只考虑了振幅的相关性，没有考虑相位。

此外，在干涉模块，在slc数据配准的时候使用了光谱移动滤波（spectral shift filter），

为了优化相干性，这就把一些对于干涉处理中不重要的信号“屏蔽”掉了，这会影响生成的强度数据和相应的辐射定标结果。所以，不能对干涉模块中生成的强度数据进行“地理编码和辐射定标”（\）。

Q. – 椭球体高程（Ellipsoidal Height）是什么意思？

A. - 椭球体高程的意思就是相对于参考椭球体的垂直距离，是椭球体表面上的点到标准椭球体上的量测距离。

Q. - 正高（Orthometric Height）是什么意思？

A. – 正交-或等位-高程，代表到大地水准面的垂直距离，是沿着力作用线到大地水准面表面的距离。

Q. –什么是大地水准面（Geoid）？

A. – 大地水准面是地球表面，是以地球质心为圆心的一个椭球体，代表地球的形状，假设地球表面是均匀的，也无地形起伏。

Q. – 用带有椭球体高程的DEM或者者用大地水准面，得到的结果不一样，采用哪一种更好呢？

A. – 带有椭球体高程的DEM无疑是更好的，可以避免在特定的地理位置用大地水准面带来的误差，采用椭球体高程保证了和现有的大多数卫星系统的一致性，这些系统的量测和参数都是基于地球中心的。下图表示了大地水准面和椭球体高度的区别。

Q. – 在SARscape里输入的DEM都使用的什么大地水准面？ A. – 使用的是以下两个：

-GTOPO30 和SRTM-3 DEM使用的是ACE 的EGM96 大地水准面. - RAMP DEM使用的是OSU91A大地水准面

Q. – 基准面偏移（Datum Shift）的意义是什么？

A. – 基准面偏移参数是将椭球体的圆心（ellipsoid's origin）转换到地球的中心（Earth's centre）。

Q. –好几处SARscape的面板（如配准、地理编码、投影转换等），能使用不同的重采样/内插方法，在什么处理情况下或者输入数据应该选择什么方法，有没有依据？ A. –一般来说，需要注意：

- 最近邻法是只考虑到了临近的象元，最适用于对分类数据采样.

-双线性内插法，三次卷积法、四次卷积法、考虑到了周围临近的4、8、16个象元，一般情

况下，内插的次数越高得到的结果越精确，不过有的时候，4次卷积会得到一些负值的结果. - Sinc采样方法是一种更完美的内插方法，考虑了周围256个象元，不过在后向散射很强的地方还是会得到一些负值象元，所以对强度数据要谨慎使用该采样方法.

- 最佳分辨率（The Optimal Resolution）采样方法是地理编码的功能特有的，适用于SLC数据，可以最大程度的保留辐射和几何特性，特别是在丘陵地带或山区。

Q. – 为什么制图栅格大小是25米，而输入的数据（RADSAT-2）分辨率高于5米呢？

A. – 输出的制图分辨率大小是在默认的Default Values>General>Cartographic Grid Size选项下设置的. 这个值（如25米）就是所有处理的时候输出的默认的栅格大小，（如 Basic>Geocoding and Radiometric Calibration, Interferometry>Phase to Height Conversion and Geocoding, Interferometry>Phase to Displacement Conversion and Geocoding and others），用户可以在处理的时候在参数面板上设置输出网格大小用于输出的结果。

Q. – 怎么输入地理坐标：度、分、秒，或是小数形式的度？ A. – 地理坐标是以小数形式的度单位输入的，要用小数点（如29.3），不能是逗号（如29,30）

Q. – 在SARscape处理过程中哪里会用到地面控制点，是如何使用的？

A. – 在一些功能会用到，如地理编码、干涉去平、PS和SBAS处理中的\GCP\文件，还有其他的一些功能用GCP来修正在方位向的定位误差（如错误的方位向开始时间），和/或距离向上（如错误的斜距范围）。这种校正不是用仿射多项式变换，而是用严格的多普勒算法，所以一个精确的GCP点就可以用来做校正，（如果在一个处理中有更多的控制点，最终的偏移量是从所有控制点的平均值计算出来的）

当然还有一些误差，通过上面说的方法没有纠正过来的，如原始数据的脉冲重复频率（PRF）误差，就会用另外的方法来校正这种错误，用至少两个控制点来校正方位向的距离。 不过大的轨道误差引起的宏观的场景旋转不能被校正。

还有一点就是，当提示\geolocation error\时（如在地理编码的时候没有用到任何GCP），是在视觉上很难觉察到（1个像素或者是亚像素），所以很难找到这么精确的GCP来校正这么小的误差。

还有另外一种GCP，叫做轨道GCP（\GCP\），用在干涉功能中的轨道精炼和重去平的步骤（InSAR、DInSAR、PS、SBAS、ScanSAR），在这种情况下，会用到较多的GCP，来校正即使是很小的轨道误差。

Q. –如何交互式的编辑SARscape生成的结果？

A. –相位解缠后的相位编辑功能，可以用于其他SARscape生成的数据的编辑。

Q. – SARscape是如何处理用ENVI的功能处理和生成的数据的？

A. – 任何数据要在SARscape中处理，需要事先进行数据的导入，有以下两种可能：

- ENVI的数据是来自于之前SARscape处理过的结果，其栅格文件的信息没有改变，这种情况下，可以直接使用之前的SARscape头文件（.sml），用Import->ENVI->SARscape Original 功能导入数据.

- ENVI的数据是不是由SARscape处理过的，或者在处理的过程中数据的栅格文件的信息改变了，这种情况下，用Import->ENVI->ENVI Original 功能导入数据.

Q. – 用UPS和球面投影，有什么区别？

A. - UPS 是球心投影，这是专门为极地地区采用的投影，投影平面与球体的表面相切，投影的点在球体的中心；球面投影是方位正形投影，适用于不同的纬度，不同的地理位置的参数可以自己定义。

Q. – 用SARscape相关工具得到的SAR数据的定标结果进行三波段彩色合成，彩色合成的结果上，高亮显示的红色区域，即使不把整体后向散射高的分配给红色，还是会高亮显示红色，该现象如何解释？

A. – RGB彩色合成是对图像解译很敏感的一种方法。取决于如何将原始的标准化数据从原始的浮点型格式分配到24位（8位是一个颜色通道）的RGB tiff数据，在生成彩色合成的功能\Color Composite\用每个通道不同的比例(default setting)或者是相同的比例(\来进行3个图层的彩色合成，最后一个选项是能最好的保留原始信息的方法，常用于处理后向散射变化的情况，不仅是对表面粗糙度的局部变化（如洪水区域、植被生长、森林砍伐等），还有介电常数的变化在图像的大部分区域。

举个例子，如果我们想监测由于洪水导致的后向散射的变化，假设RGB图像是由指定的洪水发生图像为红和绿通道，这样可以看到在图上的很多区域有明显的黄色（表面粗糙度由于洪水的影响呈现不平状，由介电常数导致的后向散射增加，是由土壤、植被和空气中湿度的增加导致的）。反之，可以得到明显的蓝色，在有洪水覆盖的地区，在红色和绿色的后向散射信号几乎是0。为了解决该目视解译的问题，生成RGB彩色合成图像时，可以手动的做直方图拉伸（可以分别对三个通道做），直到没有发生变化的区域呈现不明显的颜色。可靠的解译和变化监测是用真实的后向散射系数（从浮点型的标准产品计算出来的时间信号）来进行的，而不是用从RGB合成图像上的一般视觉上的评估

Q. – 数据处理的时候可以用掩膜吗？ A. – 通常数据处理都是可以用掩膜的。

Q. – 很多软件在处理的时候，其他处理功能就不可用了，能有什么办法在软件处理的时候，其他功能还可以用？

A. –SARscape在执行批处理程序而不是在面板上直接start执行的时候，可以运行ENVI或SARscape的任何功能，需要注意的是，如果一台机器上同时执行好几个批处理的程序，要把每个批处理设置成不同的工作路径。

Q. – 在处理一个连续的批处理程序的时候，选择文件后缀（扩展名）是比较困难的？

A. – 当用户对SARscape和其文件命名熟悉之后，这些会容易些。对输出的文件扩展后缀名熟悉后对于文件的选择是很有帮助的。

Q. – 我们想从相干性图上得到森林的变化，轮廓清晰的范围比植被和树木范围有更高的相干性，我们知道相干性受地形和后向散射强度的影响，如何减少由地形引起的相干性？ A. – 首先是光谱位移滤波（spectral shift filter），一般在默认的时候都选中了，就是一个处理该问题的有效参数，该功能是在干涉生成的时候执行，这一步用已知的DEM，将局部的坡进行滤波，当没有DEM输入的时候，就被认为是平地没有坡度的影响。此外，减少坡度影响可以用去除残余相（默认参数里的Remove Residual Phase Frequency），这个参数会在滤波和相干性生成（\）的时候发挥作用。最后，适当应用相干性来生成（或建模）森林参数，一个有用的参数是入射角图像（_lia），在地理编码的

时候选中有关选项会生成该结果。经过地理编码的入射角图像可以转换为斜距图像，利用地距转斜距的功能（Map to Slant Range Image Transformation）。

Q. – 利用SAR做海事上的应用，有相关的介绍资料吗？

A. – 将会推出SAR海洋用户手册，能从海洋的SAR图像中获取各种信息，以及海洋SAR数据的分析方法，为那些想将SAR数据用于他们研究的非专业人员但是科研工作者提供的。

干涉叠加模块

Q. – 使用该模块就能利用数据的全分辨率（Full Resolution）进行处理，这样的说法对吗？ A. – 所有的模块都可以利用数据的全空间分辨率，之所以PS方法得到的形变精度比InSAR高，是因为PS的位移监测（差分相位）目标是稳定的散射体。

Q. – 在处理流程中，有没有对属于同一时间序列的数据，进行过特殊的处理？让其靠近其多普勒质心。 Does the processing chain applies any procedure specifically intended to make closer the Doppler Centroids of the acquisitions belonging to the same input temporal series? A. - The interferometric processing takes properly care of the data Doppler Centroids either during the design of the interpolation filters or when the common azimuth bandwidth filter is performed (as it is for example done in range direction for the baseline de-correlation/spectral shift filtering). These filtering steps can be activated by setting the relevant flag in the Default Values. Specifically, the common Doppler bandwidth (i.e. Doppler Filter) and the Spectral Shift Filter are never performed (Default Values flags always off) within the PS module due to the nature of the algorithm and the type of targets that are considered, which are actually point targets; for these objects there is not spectral shift, baseline decorrelation or decorrelation of the Doppler bands, thus none of the filters in range and azimuth direction is necessary, moreover the activation of these filters cause the loss of much information for real point targets and it eliminates the advantage of exploiting large baselines for obtaining a very precise estimate of the PS height.

On the other hand the SBAS module, which focuses on distributed targets where common Doppler bandwidth and spectral shift filtering are meaningful, normally activates these filters.

Q. – 有没有对比过PS和SBAS的处理结果？

A. – 对比过，对日本沉降监测的结果，还采集了实测数据。两种算法的对比结果非常吻合而且和实测数据也是一致的（可以参考PS和SBAS的详细文档对比）。 两种方法的区别有以下几点：

◆ PS (永久散射体) – 通常对于局部目标的形变监测更好

◆ SBAS (短基线) – 形变信息是均匀的，而且可以监测非线性的形变 不过在这里是用了一个PS应用的简单例子，因为： 1、 用到的数据数量比较多，用了34景数据 2、 监测的沉降是在长时间序列上呈线性形变的 3、 在检测区有明显的高相干的散射体

当要监测的情况不满足上面的某个条件时，SBAS方法要更加可靠。

Q.-PS方法监测形变的精度可以达到多高？

A. – 要考虑到图像的个数和形变的速率和PS要求的匹配程度，能达到的精度是： ? 能得到小于1厘米/每年的位移速率

? 在每个输入的时相相对于参考时相，PS点的相对位移的时期都约是厘米级的。 参考时相，对于假设没有形变“0位移”的参考时相，不一定是在输入面板的“Reference file”里（手动或自动），而是有可能是固定的一个最早的图像，在输入列表（input file list）里面，在相关的面板（default value）里面检查设置。

Q. – 当PS点的密度小于每平方公里100个点的时候，结果的精确度就会有问题，对这种情况下有什么建议？

A. – 当PS点密度低的时候，最重要的问题就是大气干扰的估算和去除，这也是为什么PS方法适合于城市地区，但是在郊外和农村地区，会常常得到错误的结果，这种情况下，我们建议用SBAS方法。

Q. - 我有19个时相的用来做干涉的图，是在略城市化的农村区域，建议用什么方法？ A. – 实际上，当数据量小于最小数据量要求（20景），以及所在的区域可能获取不到足够的PS点的时候，最好的方法是用SBAS方法。不过也可以尝试用PS方法，减少默认的PS点密度，这也是在略城市化的区域一种常见的设置。

Q. – 在线的文件里面提到过，在干涉模块用到的一些干涉处理算法，在干涉叠加模块的常规处理中也用到了，在干涉叠加处理的时候，有没有用到什么特殊的干涉模型方法？

A. – 干涉计算模型，在干涉叠加常规处理中是必须的，执行的时候会自动调用，而PS/SBAS的输出不需要在这个模块调用任何功能。

Q. –在其他默认的面板里面设置的基线阈值的意义是什么？

A. – 默认的阈值（500米）相当于临界基线（用基线估算算出来的临界基线）的5倍，实际上，在干涉处理时和基线相关的约束对分析分布目标是适用的（干涉模块和干涉叠加模块），相反的，如果要分析的是局部/点目标（用PS的方法），那么临界基线就不再有任何限制了。

Q. – 在PS处理时，对相干性阈值的设置或者调整有没有什么特定的标准？ A. – 用默认的设置进行第一次迭代处理之后，看找到的PS点的情况，如果在找到的PS点里面很多可能噪声像素，可以增加相干性阈值，如果找到的PS点太少，可以减少阈值。如果不选择“Generate Slant Range Products”，选择“Generate Geocoded Products”选项的话，第二次的迭代是在生成最终地理编码产品的时候执行的，这种情况下，第一次和最长的一次处理，相干性阈值改变的话，对第一次和时间最长的一次处理，不会有影响，而这两步就跳过了。

Q. –参考数据的选择有什么依据？PS处理中精度估计的结果包含了什么信息？

A. –参考数据是自动选择出来的，是输入数据中平均基线最小的，这就说明这个数据是： 拥有相对于其他数据的完美的空间和时间位置，这有利于图像配准以及所有的处理会变得容易，还有提供较高的相干性（短的基线对大量的去相关是不太敏感的）。选择参考图像的时候，要避免选择大气影响大的图像。

精度估计 \功能可以得到对测量精度的估计的初步的结果，对所有时相的

数据和所有PS点都是有效的。在PS处理中的精度因子的计算方法和GPS系统的精度因子一样，考虑到时序图像的几何，基线估算和相应的2π高程模糊，还会用到默认的PS点密度。

最后，当PS处理结束后，会计算出每个PS点的实际精度，基于象元相干性和实际的PS密度

Q. -如何知道哪些数据是因为配准失败而被丢弃的？

A. –当配准失败的时候，在输出路径下会生成相应的文件： \，文件中会将由于配准问题而丢弃的数据列出来。

Q. –在PS分析完成之后，有没有哪些软件可以检查的关键要素，确保结果的可靠性？ A. – 其中一项就是检查相关的配准过程：\文件包含了那些没有配准的图像。另一项要检查的是大气影响的去除，根据经验来说，对比\（大气影像去除之前）和\（大气影响去除之后），多时间相干性增加，说明大气校正是成功的。

> 对比前后的cc文件，可以更好的确定如何设置滤波器参数（低通，以米为单位；高通，以天为单位），以避免相干性的损失，来更好的估算形变量。

Q. - SBAS 的首字母分别地表什么？ A. – 代表 “small-baseline subset”，该技术最早起源于：Berardino, Fornaro, Lanari and Sansosti (参看参看文献).

Q. – 在SBAS连接图生成的时候，最大最小的空间和时间基线的确定有什么标准？

A. – 设置最小和最大基线一般是根据数据处理的目的来确定的，看是要用SBAS做形变监测还是提取DEM。如果是要做形变监测，典型的一个设置就是能得到所有的连接图，反之，如果要生成DEM，同样是不连贯的像对可以得到高精度的DEM。

基于上面的叙述，如果在生成DEM的时候选择了“允许不连贯的区域”（Allow Disconnected Blocks），最大最小基线就可以设置为可去除最小像对基线的值，比如临界基线的20%的像对就去除了，这些是对高程测量无用的数据。

另一方面，不管是提取DEM还是测量形变，最好避免选择太大基线的像对，这样的数据相干性很低，这种情况下，空间基线的上限可以设置为临界基线的50%，时间基线的上限也可以参考这样的设置，因为相干性和基线距离是成反比的，影响时间失相干的因素和地表覆盖有关的，如植被、土壤含水量、积雪覆盖等。

在处理中会考虑到这些阈值，在处理过程中是公开的，在干涉生成这一步之后，分析结果以确定在连接表编辑的时候去除哪些像对。

Q. – 在SBAS反算的时候，如何设置大气的低通（Atmosphere Low Pass）和高通（High Pass）滤波参数？

A. – 窗口大小越小，滤波器就越强，滤波器中引进的平滑可以评估，对比\中的时间信号（没有进行大气去除的结果）和最终的形变结果。这有助于理解是否大气滤波平滑去

除了重要的形变信息。这可以发生在那些大气模式（atmospheric patterns）很少的情况下（如山区范围）因此滤波大小必须要相应的调整，这可以通过在SBAS反算之前观察干涉图评估出来。然而，如果形变模式（displacement patterns）窗口就很小的话（等于或小于大气的），这种情况下不能过分的减少窗口，以防止将形变作为大气信息去除掉了。

Q. –软件如何处理数据集的低相干性的散射区域？

A. - The example below shows a typical case of scattered coherence areas. The area of interest, where subsidence phenomena have been reported, is in the lower right corner.

This area is well represented and the coherence is good, but the surrounding low coherence zones are can affect the SBAS inversion process and eventually the measurement accuracy. However, if there are enough interferograms (theoretically at least five per acquisition) the program is able to \the missing information by means of the 3-dimensions unwrapping approach.

Q. – 干涉图系列的解缠方法选择的是什么？

A. - Actually the unwrapping execution depends on the Interferometric Stacking approach which is adopted:

◆ in the Persistent Scatterers, based on the original publication of Ferretti et al., the unwrapping is not performed for the estimation of the displacement rate and height corrections, since a pixel-wise spectral analysis approach in the time-baseline plane is exploited; this approach has the advantage of avoiding the need of unwrapping by working on the complex data only.

◆ in the SBAS Inversion, it is possible to choose between two main methods: Region Growing or Minimum Cost Flow, this last one either with a square or with an irregular triangulated - Delaunay - grid (Default Values>Interferometry). In case the Delaunay method is adopted, a 2- or 3-dimensions unwrapping (the third dimension being represented by the time) can be selected. At this regards it must be noted that the 3-dimensions approach provides superior results when there are disconnected areas (typically due to low coherence), as it exploits the high coherence

interferograms (third dimension) to estimate how to create new connections in \coherence) interferograms; the disadvantage of this method is that it is much more costly in terms of processing time. The unwrapping, in the SBAS processing chain, is carried out two times: once before a first estimate of average displacement rate and height correction and once more afterwards to refine the first results.

Q. – 设置分解级别的时候，如何设置最优的级别，有无特定的依据？ A. –这里没有特定的规则，因为该参数的最优设置取决于区域的相干性和解缠的方法。 There are not specific rules since the optimal setting of this parameter depends strongly on the scene coherence and also on the unwrapping method adopted. 大多数情况下，尤其是选择了默认的方法，分解等级为1是最好的设置。有时候，将等级设置为1（1A – 干涉生成工作流这一步）增加到2（在第一次反算这一步），等级的增加会得到更好的结果。

Q. - After the SBAS inversion, some of the Unwrapped Interferograms are still affected by Directional Slopes. Is this due to a problem during the Refinement and Re-Flattening step? What is the best way to address such issues?

A. - The first thing to do it is to visualize and verify the unwrapped data (_upha_list_meta), after the Interferometric Workflow step, in order to understand where the GCPs must be located for the next Refinement and Re-Flattening. The worst pairs (i.e. very low or scattered coherence, which typically causes a bad unwrapping) shall be removed by means of the appropriate tool; if possible before the Refinement and Re-Flattening, otherwise after the First Inversion step.

After the execution of the Refinement and Re-Flattening, in order to ensure that the GCPs have been properly selected (in terms of position, distribution and quantity), another visual analysis is needed to verify that major residual phase ramps are not in; if they are still present, a higher number of GCPs (20 or more) is probably required.

Once the re-flattened data have been checked, the First Inversion step can performed, which generates a new set of the unwrapped data (these are stored in the \After this step, the remaining \inversion.

Note that the program is implemented in a way that, even when some small residual ramps remain (of course the less the better...), the inversion process is not notably affected.

Q. – 有没有生成轨道GCP文件的简单方法？

A. – GCP点当需要输入\时，该参数是最重要的参数，其选择标准非常重要。在干涉测量模块的QA里面已经提到过了，而在干涉叠加模块这里，需要增加说明的是： it is that the same set of points must be used to \have been created after the coregistration onto the \not impossible...) to find the best point configuration for all pairs at once, the suggested

approach it is to try a point distribution which is good for the majority of the pairs (for this purpose it becomes useful to visualize the stack of images by means of the meta file generated automatically by the program); in several cases the use of the Default Values>Flattening>Refinement and Re-flattening>Refinement Method>Residual Phase, which allows the use of the less precise but most robust phase correction approach, shall be adopted. It must be also taken into account that some points, which are inside the imaged area for some pairs, can fall outside in others and thus it can be required to enter more points; in other cases, for instance when the original orbits are all very much accurate and there is not any visual evidence of residual phase ramps, it can be sufficient to simply remove a phase offset (constant value) by choosing a Default Values>Flattening>Refinement and Re-flattening>Refinement Method>Residual Phase Poly Degree equal to 1.

In any case, also relying on very precise orbits, the use of the \the reflattening process must be always executed in order to correct the phase offset (i.e. constant phase removal).

Q. – 对SBAS结果最好的分析方法是什么？

A. – 得到的结果是以meta文件分组存放的 ，为了方便做多时相的时间序列分析。Meta文件可以用任意栅格分析的工具，或者转换成矢量文件后，用矢量数据的分析工具进行分析。

本文来源：https://www.bwwdw.com/article/9e26.html