[发明专利]矿区地表移动变形3D WebGIS预计分析方法及系统

权利要求书

1.一种矿区地表移动变形3D WebGIS预计分析方法，其特征在于，包括以下步骤：

步骤1，矿区地表三维数据获取及格式转换：

步骤11，数据获取：

获取矿区地表多视倾斜立体影像，结合外业实测的像控点坐标，基于SfM方法和密集匹配方法经由解析空中三角测量、密集匹配、不规则三角网构建及纹理映射生成osgb格式的矿区地表三维实景模型数据；

步骤12，数据转换：

将osgb格式的矿区地表三维实景模型数据转换为Cesium支持的3DTiles数据格式的矿区地表三维实景模型数据；

步骤2，单体化矢量数据提取及属性数据挂接

步骤21，根据步骤1得到的矿区地表三维实景模型数据采用模型切割方法、ID单体化方法或动态单体化方法分离并选中模型单体，得到采煤工作面及关键地物要素的二维矢量数据；

步骤22，属性数据挂接，将属性数据挂接在工作面及关键地物要素的二维矢量数据上，属性数据包括工作面沉陷预计参数结构表、电力设施数据结构表、道路数据结构表、建筑物数据结构表；

步骤3，对3DTiles数据格式的矿区地表三维实景模型数据进行服务化管理；

步骤4，利用MySQL关系型数据库存储采煤工作面及关键地物要素的二维矢量数据及其挂接的属性数据；

步骤5，3D WebGIS分析系统单元构建，3D WebGIS分析系统单元包括前端和后端，前端包括图层控制及场景显示模块、开采沉陷3D WebGIS空间分析模块、损害等级云图绘制模块及统计报表输出模块；后端包括服务请求响应模块、开采沉陷预计模块及GIS空间分析功能函数封装调用模块及数据库管理模块；

步骤51，后端实现

步骤511，开采沉陷预计模块，采用概率积分法预计模型并结合工作面沉陷预计参数结构表对地下采煤引起的地表沉陷情况进行预计，内插相应类别的地表变形等值线；

步骤512，GIS空间分析功能函数封装调用模块，通过将地理空间目标划分为点、线、面不同的类型，获得这些不同类型目标的形态结构；将空间目标的空间数据和属性数据结合起来，进行特定任务的空间计算与分析；地理空间实体抽象为点、线、面三种；在GIS中，上述实体间的空间关系有点-面、线-面及面-面三种实体间的空间拓扑关系；将封闭的地表变形等值线视为面状类型数据，继而定量研究分析其与代表矿区不同地物类型的点、线、面三种实体间的空间拓扑关系；

步骤513，数据库管理模块，建立采煤工作面及关键地物要素地理空间数据库和属性数据库；

步骤52，前端实现

步骤521，前端框架建立，采用Cesium开源三维引擎框架，实现在浏览器中展示三维虚拟地球的功能；

步骤522，图层控制及场景显示模块：

步骤5221，配置图层目录树，使各关键地物要素及采煤工作面矢量图层和矿区场景模型图层按顺序分层显示，并通过配置显示参数，控制图层的显示效果；

步骤5222，基于各关键地物要素及采煤工作面矢量数据，利用Cesium的贴地面矢量绘制功能实现其三维绘制与渲染；

步骤523，开采沉陷3D WebGIS空间分析模块：

步骤5231，调用Cesium前端空间分析模块并通过交互式地选择一个或多个所要预计的工作面；

步骤5232，前端请求并调用后端工作面沉陷预计参数结构表数据库及开采沉陷预计模块，获得开采沉陷地表下沉及地表变形等值线；

步骤5233，将开采沉陷影响范围传至后端，联合后端采煤工作面及关键地物要素地理空间数据库并同时调用后端的GIS空间分析功能函数封装调用模块执行空间分析操作，将分析结果返回至前端进行解析、显示；

步骤524，损害等级云图绘制模块，依据开采沉陷预计模块分析结果和不同损害等级的边界值，按轻微损害、轻度损害、中度损害和严重损害绘制四个等级的损害区域云图并用颜色区分；

步骤525，统计报表输出模块，依据开采沉陷预计模块分析结果统计采煤工作面的走向、倾向和/或用户指定方向上的最大最小值、沉陷区域的面积和体积、指定深度的沉陷区域面积和体积，并以表格的形式导出。

2.根据权利要求1所述矿区地表移动变形3D WebGIS预计分析方法，其特征在于：步骤11中采用倾斜摄影方法获取矿区地表多视倾斜立体影像。