一种确定地面沉降和区域静载荷相关性的方法
【技术领域】
[0001]本发明涉及PS-1nSAR技术和遥感建筑指数法技术的综合应用领域，具体而言，涉及一种确定地面沉降和区域静载荷相关性的方法。
【背景技术】
[0002]。目前针对城市化进程中动静载荷与地面沉降的关系问题，国内外学者采用多种传统的技术手段对其展开研究。M.H.Aly通过研究，指出地铁网络和高密度人口也是导致地面不均匀沉降的因素;S.Stramondo认为在原位有效压力条件下，载荷的重量和建成年份持续的时间对地面沉降贡献不同；Kaynia等将列车载荷简化为间距一定的移动载荷，用成层粘弹性半空间上的E-B梁模型对地基振动载荷效应进行了模拟研究。上海的地面沉降研究结果表明，高层建筑对地面沉降量的贡献率可达到30唐益群等借助离心模型试验手段，探索了不同建筑容积率下密集建筑群区对地面沉降的影响规律;何庆成指出地面上的动荷载在一定条件下，将引起土体的压密变形，诱发地面沉降。但是这些研究大多采用的是传统的技术手段或者研究方法，虽然局部特点明显，但由于区域水文地质条件的非均一性，需要区域尺度的整体研究。
[0003]针对地面沉降和区域静载荷相关性问题，严学新等探究了静载荷与地面沉降的关系。使用建筑容积率(即地上总建筑面积与用地面积的比率)来表征地面载荷，沉降信息的获取自历年的监测数据。对于二者关系的分析，则是基于数据的对比分析。
[0004]上述方案存在的问题是，首先获取大范围内、长时间序列的建筑容积率信息繁琐，困难，调查成本高;其次对于沉降信息的获取，也是基于小范围内的监测数据，涉及范围太小，无法反应区域整体沉降特征。再者，在二者关系的确定时，上述方案仅仅做了简单的对照分析，只能定性地确定二者的关系。
【发明内容】
[0005]本发明提供一种确定地面沉降和区域静载荷相关性的方法，用以克服现有技术中存在的至少一个问题。
[0006]为达到上述目的，本发明提供了一种确定地面沉降和区域静载荷相关性的方法，包括以下步骤:
[0007]对TM遥感影像采用设定窗口选取设定数量的采样单元；
[0008]对每个采用单元的图像进行预处理，得到NDBI指数、MNDWI指数和SAVI指数，根据NDBI指数、MNDWI指数和SAVI指数进行IBI指数反演，得到反映建筑载荷时空变化的遥感建筑用地指数IBI;
[0009]采用PS-1nSAR方法对每个采样单元的图像进行处理，根据所获取的该采样单元内的时序沉降数据得到PS点沉降速率；
[0010]将PS点沉降速率作为因变量，遥感建筑用地指数IBI作为自变量，采用统计相关分析方法，获得PS点沉降速率与遥感建筑用地指数IBI的定量关系。[0011 ] 进一步地，所述设定窗口为单个PS点、20 X 20窗口大小和100 X 100窗口大小中的一种或多种。
[0012]进一步地，当所述设定窗口为单个PS点时，采用Spearman秩相关系数法分析单个PS点的PS点沉降速率与遥感建筑用地指数IBI的定量关系。
[0013]进一步地，当分析20X 20窗口大小采用单元内PS沉降速率与遥感建筑用地指数IBI的定量关系时，首先把量纲不统一的遥感建筑用地指数IBI与PS沉降速率进行归一化，并转换为百分率，再采用相关回归分析方法，分析遥感建筑用地指数IBI与PS沉降速率的定量关系。
[0014]进一步地，当采样单元的窗口大小为100X100时，根据其空间分布和PS点的位置情况，共选取了30个100 X 100窗口大小的采样单元，每个采样单元的沉降速率取其所包含的PS点的平均值，对应的遥感建筑用地指数IBI取其栅格单元的平均值，并分析得到遥感建筑用地指数IBI与PS沉降速率的变化关系曲线。
[0015]本发明首先使用PS-1nSAR方法解决了大范围时序沉降信息的问题，覆盖范围大，成本低，准确度高。其次使用遥感植被指数(IBI)来表针地面建筑物的覆盖程度，解决了建筑用地时空变化数据不易获取，调查成本高，操作复杂的问题;再次，在数据分析方面，以PS点沉降速率作为因变量，遥感建筑用地指数IBI作为自变量，采用不同的统计相关分析方法，定量研究两者相关性，解决了现有技术方案中只是简单定性分析的问题。
【附图说明】
[0016]为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0017]图1为本发明一个实施例的确定地面沉降和区域静载荷相关性的方法流程图。
【具体实施方式】
[0018]下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0019]图1为本发明一个实施例的确定地面沉降和区域静载荷相关性的方法流程图；如图所示，该确定地面沉降和区域静载荷相关性的方法包括以下步骤:
[°02°] 对TM(Thematic Mapper，主题地图)遥感影像采用设定窗口选取设定数量的采样单元；
[0021]对每个采用单元的图像进行预处理，得到NDBI指数、MNDWI指数和SAVI指数，根据NDBI指数、MNDWI指数和SAVI指数进行IBI指数反演，得到反映建筑载荷时空变化的遥感建筑用地指数IBI;
[0022]采用PS-1nSAR方法对每个采样单元的图像进行处理，根据所获取的该采样单元内的时序沉降数据得到PS点沉降速率；
[0023]将PS点沉降速率作为因变量，遥感建筑用地指数IBI作为自变量，采用统计相关分析方法，获得PS点沉降速率与遥感建筑用地指数IBI的定量关系。
[0024]其中，采样单元的大小可根据研究区的建筑物分布和PS点分布确定的，其大小是可变的，例如，上述设定窗口可以为单个PS点、20 X 20窗口大小和100 X 100窗口大小中的一种或多种。
[0025]在本发明的一个较佳的实施例中，在三种不同的采样单元内，基于GIS分析平台，提取PS点的平均沉降速率与对应的遥感建筑用地指数IBI，PS点沉降速率作为因变量，