SBAS矢量结果说明

SBAS处理得到的形变速率结果可以选择输出矢量,在SARscape5.4以上的版本中,可在SBAS Geocoding参数界面选择Make Geocoded Shape参数为True,即可输出SBAS矢量结果。

图 SBAS选择输出矢量结果

在SARscape5.2.1及更低版本中。用栅格转矢量的工具可得到SBAS的矢量结果,工具为:/SARscape/Interferometric Stacking/SBAS/Raster to Shape Conversion。

工具启动工具后,在Input File中输入工程文件auxiliary.sml:


图 输入工程文件面板

切换到Option Files中,这里可以输入矢量文件,用来确定栅格转矢量的范围,如果不输入,则会对整个区域进行转换。


图 可以输入shp文件定义要转换的范围,如果不输入,会整景范围做转换

切换到Parameters面板,面板中参数如下:

West、East、North、South这四个可以确定范围,如果没有输入矢量范围文件,也可以在这里输入坐标范围来确定栅格矢量转换的范围。如果不输入,默认的-9999,代表没有子区范围,所有范围都会进行栅格矢量转换。

Product Coherence Threshold最小相干性阈值,大于该阈值的进行转换。

Shape Max Nr of Points一个shp文件中最大点数量,如果一个shp文件放不下,会生成若干个shp文件,带有 _01; _02等标识。建议这里不要超过100000,否则显示起来会很困难。

KML Max Nr of Points一个kml文件中最大点数量,如果一个kml放不下,会生成若干个kml文件,带有_01,_02等标识,建议这里不要超过50000个点,否则数据量太大,显示会有问题,为了显示顺畅,建议这里点个数设置少一些,比如10000个点。

Generate Shape生成shp文件,该选项设置为ture

Generate Kml生成kml文件,该选项设置为ture

Scaling RangeKML文件的图例范围

Generate Shape Time Series生成矢量文件的位移时序

Generate Kml Time Series生成kml文件的位移时序


图 参数设置面板

生成的结果位置在生成的结果在_SBAS_processing\vector\coherence_thr_路径下,有shp和/或.kml结果,如下图是样例区域生成的4个shp文件


图 生成若干个shp文件,文件名带有_0,_1等标识

打开shp文件的属性表,可以查看属性字段:

 

 

属性字段代表的意义为:

字段

意义

Velocity

该点的形变速率,单位mm/y

Coherence

该点的相干性,值的范围是:相干性阈值~1

H_Precisio

高程测量平均精度估算值,单位为米

V_Precisio

速率测量平均精度的估算值,单位为mm/y

L1Norm

所有时相的位移量绝对值的总和,单位是mm。

ChiSqr

拟合和反演质量的量度。这个值越高,拟合和反演质量就越差

Lat/Lon/X/Y/Z

纬度/经度/高程

LOS_Az

Azimuth Line of Sight:视线方位角。从北顺时针方向为正;从北逆时针方向为负。

LOS_In

Incidence angle of the Line of Sight:视线入射角,测量的角度是在椭球面上的视线和垂直方向上的夹角。

Height_cor

高程校正值,单位为米。

Term0

反演多项式的常数项,单位为mm。如果选择 “无位移模型”时候,不生成该值。

Term1

反演多项式的一次方的系数单位为mm(以毫米为单位)。果选择 “无位移模型”时候,不生成该值。

Tot_Displ

所有时相累积的位移量,单位为mm

D_YYYYMMDD

每个时相相对于第一时相的位移量,第一个时相为0,单位为mm

Shp矢量文件包含了模型计算出的所有信息,每个点的相干性、历史形变量等。可以像PS一样,先将矢量文件在ENVI中显示出来,再用/SARscape/General Tools/Time Series Analyzer/Vector工具进行可视化分析。

形变点的颜色渲染:在TS Analisy面板上点击Apply


多个shp的情况下,可以在矢量图层上点击右键,选择Set As Actve Layer,



多个shp文件根据速率值进行彩色渲染结果:


图 速率彩色渲染图

形变点的历史形变过程绘制:选中一个点,点击Plot Time Series,可以绘制出该点的历史形变图。



双击.kml文件,打开google显示形变点。

注:如果google earth版本太高,可能图例不能正常显示。


图 栅格转矢量生成的kml在Google Earth显示

注:处理区域较大时,kml文件较大,不便于显示,SARscape5.4版本已经取消了该功能,要将形变结果叠加在google earths上显示请参考SARscape得到的形变结果显示与分析

posted @ 2022-07-28 10:19  ENVI-IDL技术殿堂  阅读(4829)  评论(1编辑  收藏  举报