[发明专利]将Mike流场结果文件转换成shapefile矢量文件的方法和装置

说明书

本发明实施例提供一种将Mike流场结果文件转换成shapefile矢量文件的方法和装置，属于水环境仿真与地理信息系统技术交叉领域。该方法包括：读取所述Mike流场结果文件，该Mike流场结果文件包含一个或多个保存时间步、该一个或多个保存时间步所对应的计算网格及所述计算网格对应的场变量信息；针对每一保存时间步，遍历该保存时间步所对应的所有计算网格，并针对每一该计算网格，根据该计算网格对应的水面高程与预设淹没水深阈值确定该计算网格是否被包括在淹没区域；以及将针对每一保存时间步所确定的被包括在所述淹没区域的计算网格的坐标及预设场变量信息写入所述shapefile矢量文件。藉此，实现了为用户提供流场数据的在线浏览分析及网络可视化的多元需求。

技术领域

本发明涉及水环境仿真与地理信息系统技术交叉领域，具体地涉及一种将Mike流场结果文件转换成shapefile矢量文件的方法和装置。

背景技术

流动具有明显的空间三维特性，由于流场数据结构复杂，动态多变，为了能更好的表现水体流动的复杂过程，显示局部细节为防洪调度提供决策基础，对流场数据的动态可视化显示提出了更高的要求。目前在没有安装有DHI-Mike软件的终端上查看流场数据，传统的方法主要为从Mike流场结果文件中导出数据，借助于第三方软件如Excel、Tecplot、Matlab等对流场进行可视化分析，这种方法很难把握利用数学模型计算得到的流场结果的合理性和准确性。

发明内容

本发明的目的是提供一种将Mike流场结果文件转换成shapefile矢量文件的方法和装置，其可解决或至少部分解决上述问题。

为了实现上述目的，本发明的一个方面提供一种用于将Mike流场结果文件转换成shapefile矢量文件的方法，该方法包括：读取所述Mike流场结果文件，该Mike流场结果文件包含一个或多个保存时间步、该一个或多个保存时间步所对应的计算网格及所述计算网格对应的场变量信息；针对每一保存时间步，遍历该保存时间步所对应的所有计算网格，并针对每一该计算网格，根据该计算网格对应的水面高程与预设淹没水深阈值确定该计算网格是否被包括在淹没区域；以及将针对每一保存时间步所确定的被包括在所述淹没区域的计算网格的坐标及预设场变量信息写入所述shapefile矢量文件。

可选地，所述针对每一保存时间步遍历该保存时间步所对应的所有计算网格并针对每一该计算网格根据该计算网格对应的水面高程与预设淹没水深阈值确定该计算网格是否被包括在淹没区域以及将针对每一保存时间步所确定的被包括在所述淹没区域的计算网格的坐标及预设场变量信息写入所述shapefile矢量文件包括：针对每一保存时间步，创建以该保存时间步命名的shapefile矢量文件以及GDALDataset类型的对象polygonDS，从polygonDS中创建OGRLayer类型的图层指针对象polygonLayer；针对每一保存时间步所对应的计算网格中的每一计算网格，执行以下操作：(1)判断该计算网格的水面高程是否大于所述预设淹没水深阈值；(2)在该计算网格的水面高程大于所述预设淹没水深阈值的情况下，依据polygonLayer创建OGRFeature类型的指针对象polygonFeature，并将该计算网格对应的所述预设场变量信息写入所述指针对象polygonFeature；创建OGRLinearRing类型的对象ring；将该计算网格的坐标写入所述对象ring中；创建OGRPpolygon类型的对象Polygon，并将所述对象ring加入所述对象Polygon中；将所述对象Polygon设为所述指针对象polygonFeature的几何对象；以及将所述指针对象polygonFeature销毁。

可选地，所述将该计算网格的坐标写入所述对象ring中包括：读取该计算网格中的每一节点的坐标值；判断所读取的坐标值是否需要转换坐标系；在所读取的坐标值需要转换坐标系的情况下，将所读取的坐标值转到预设坐标系下，并将转换后的坐标值写入所述对象ring中；以及在所读取的坐标值不需要转换坐标系的情况下，将所读取的坐标值写入所述对象ring中。

可选地，该方法还包括：接收对所述预设场变量信息的设定。