1、ICS07.040CCSA75CH中华人民共和国测绘行业标准CH/TXXXXX-XXXX数字水深模型数据生产技术要求Technicalrequirementforproducingdigitalbathymetricmodeldata(点击此处添加与国际标准一致性程度的标识)(报批稿)(本草案完成时间:2024年04月09日)-XX-XX发布XXXXXXXX实施中华人民共和国自然资源部发布目次前言II1范围32规范性引用文件33术语和定义34基本规定34.1数学基础34.2成果构成及要求45生产流程46生产准备56.1 技术设计56.2 数据准备57数字水深模型构建67.1数据预处理67.2特
2、征数据选取67. 3数字水深模型生成67.4 拼接与裁切67.5 相关文件制作78质量控制与成果验收78. 1质量控制79. 2成果验收79成果整理与提交7附录A(资料性)主要内插算法9A.1线性插值三角网法9A.2最小曲率法9A.3改进谢别德法10参考文献12-LX.1刖5本文件按照GB/T1.1-2020标准化工作导则第1部分:标准化文件的结构和起草规则的规定起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由中华人民共和国自然资源部提出。本文件由全国地理信息标准化技术委员会测绘分技术委员会(SACTC230/SC2)归口。本文件起草单位:广东省国土资源
3、测绘院、自然资源部南海调查中心、武汉大学、自然资源部南海发展研究院、国家基础地理信息中心、中国测绘科学研究院、国家测绘产品质量检验测试中心。本文件主要起草人:刘小丁、钟远军、赵建虎、王重平、丁华祥、张惠军、胡海驹、暴景阳、刘文建、白航、朱亚兵、周巨锁、曾纪胜、魏征、周昆、张文峰、梁亚朋、杨帆、王银霞、江林、谢敬谦。数字水深模型数据生产技术要求1范围本文件规定了基础地理信息数字成果数字水深模型数据生产的数学基础、成果构成及要求、生产流程、生产准备、模型构建、质量控制与成果验收、成果整理与提交。本文件适用于基础地理信息数字成果1:500、1:1000、1:2000、1:5000、1:10000、1
4、25000、1:50000、1:100OOO数字水深模型数据的生产作业。2规范性引用文件文件。GB/TGB/TCH/TCH/TCH/TCH/Z183162435610011004XXXX9026下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中,注日期的引用文件,仅该日期对应的版本适用于本文件;不注日期的引用文件,其最新版本(包括所有的修改单)适用于本数字测绘成果质量检查与验收测绘成果质量检查与验收测绘技术总结编写规定测绘技术设计规定数字水深模型数据成果质量检验技术要求基础地理信息数字成果数字水深模型3术语和定义下列术语和定义适用于本文件。3. 1数字水深模型digital
5、bathymetricmodeI利用一组有序水深数字阵列表达区域水下地形的一种方式。来源:测绘学名词(第4版)3. 2空白水域bIankwaterarea未获取或无水深数据的水域。4基本规定4. 1数学基础坐标系统宜采用2000国家大地坐标系(CGCS2000),采用依法批准的相对独立的平面坐标系统时,应与CGCS2000建立联系。地图投影采用高斯-克吕格投影,1:500、1:1000、1:2000数字水深模型宜采用1.5分带,1:5000、1:10000数字水深模型宜采用30分带,1:25000、1:50000、1:100000数字水深模型宜采用6分带。高程基准采用1985国家高程基准。远离
6、大陆的岛、礁,不便于与国家高程基准实施水准连测时,高程基准允许采用当地平均海面。海洋和感潮河段应采用理论最低潮面作为深度基准面,内河航道宜采用航行基准面作为深度基准面,水库宜采用设计水位作为深度基准面,其他内陆水体可采用当地平均水位作为深度基准面。并应明确与1985国家高程基准之间的换算关系。4.2成果构成及要求数字水深模型成果由数字水深模型数据、元数据及其他相关文件构成,其具体成果要求按照CH/Z9026规定执行。数字水深模型成果数据精度按表1执行。表1数字水深模型成果精度指标单位为米比例尺等级中误差0H2020H3030H5050H1001:500一级+0.21+0.25+0.29+0.5
7、20.152+(H/100)2二级+0.34+0.36+0.39+0.58J3z+QH/IOU;21:1000一级+0.21+0.25+0.29+0.520.152+(H/100)2二级+0.34+0.36+0.39+0.58J.32+(H/100)21:2000一级+0.21+0.25+0.29+0.520.152+(H/100)2二级+0.34+0.36+0.39+0.58J.32+(H/100)21:5000一级+0.25+0.28+0.32+0.54-J.22+(H/100)2二级+0.43+0.45+0.47+0.64-J.42+(H/100)21:10000一级+0.25+0.28
8、0.32+0.54JO爰+(H/100)2二级+0.43+0.45+0.47+0.64-J.42+(H/100)21:25000一级+0.34+0.36+0.39+0.58JO32+(H/1UU)2二级+0.52+0.54+0.56+0.71J.52+(H/100)21:50000一级+0.34+0.36+0.39+0.58J.32+(H/100)2二级+0.52+0.54+0.56+0.71-J.52+(H/100)21:100000一级+0.34+0.36+0.39+0.58J.32+(H/100)2二级+0.62+0.63+0.65+0.78-J.62+(H/100)2注:H代表水深。
9、5生产流程根据生产任务要求和数据情况进行生产准备,选择合适的算法,内插生成规则格网数字水深模型,按照规定的要求进行成果质量控制,对符合设计要求的成果进行整理和提交。数字水深模型作业流程如图1所示。图1生产流程图6生产准备6.1 技术设计技术设计主要要求如下:a)技术设计书的编写要求及主要内容应符合CH/T1004的规定;b)根据项目总体要求、资料分析结果等编写技术设计书;c)根据项目具体情况,技术设计书可以与水深数据采集合并编写,也可以单独编写;d)对于技术指标、技术方法与标准要求不一致的,应在技术设计中进行说明。6.2 数据准备6 .2.1数据收集收集数据时应注意数据的时效性,数据比例尺应不
10、小于数字水深模型成果比例尺。数据收集一般包括以下内容。a)测深杆、测深锤、单波束和多波束测深数据。收集水深成果数据及测量过程中的同步水位观测数据、定位数据等。同时,应收集技术设计书、技术总结、质量验收报告、设备检定报告等。b)机载激光测深数据。收集水深成果数据及测量过程中的Pc)S数据、点云数据、影像数据等。同时,应收集技术设计书、技术总结、质量验收报告、设备检定报告等。c)海图数据和水下地形图件。收集海图数据和水下地形图件,包括其编号、图名、基准、投影方式、数据获取时间、范围、比例尺等。纸质图件可通过扫描仪、数字化仪等设备,配合栅格图像矢量化软件,获取其中的水深数据。7 .2.2数据分析对收
11、集的数据进行整理和分析,内容如下:a)水深测量过程及成果资料是否齐全,了解水深数据采集的基本情况,包括设备状况、测量单位、测量时间、测量依据、测量基准、比例尺、测线等;b)对于纸质图件查看图纸成图时间和成图方法,图面是否完好、平整、无折皱、无污渍,色正、线条与注记是否清晰,印刷套合是否良好;c)水深数据采用的坐标系统、高程基准和深度基准是否需要转换;d)水深数据比例尺是否满足生产要求;e)水深数据质量是否合格,水深数据质量不合格不能进行数字水深模型构建;f)水深数据资料不齐全或缺少测量基准等重要测量要素的,不能进行数字水深模型构建。7数字水深模型构建7.1数据预处理7.1.1格式转换将已有数据
12、转换为生产软件支持的数据格式。7.1.2测绘基准统一将已有水深数据基准转换至规定的测绘基准。7.1.3多源水深数据融合多源水深数据融合前,应检查水深数据的接边处,避免存在断层、跳变等异常情况。数据融合应根据其类型、精度以及现势性,在保持基础数据特征的前提下进行。多源水深数据融合应满足以下原则。a)连续性原则:融合的水深数据应当在空间上具有连续性。b)完整性原则:融合的水深数据应该包含尽可能多的水深信息,在进行数据融合时,应尽量填补数据缺失的部分,确保融合结果的完整性。c)现势性原则:优先采用现势性更强的数据。d)高精度原则:融合的水深数据几何表达精度不一致时,应以精度高的数据为准。7. 1.4
13、特殊地物提取独立礁石、沉船、桩柱等特殊地物对应的水深数据应单独存储,不参与数字水深模型构建。7.2 特征数据选取从融合数据中,选取潮沟、航道边坡、堤坝、码头边坡线、水下陡坎和大型沙丘等能突出地形变化的水深数据以及水边线和空白水域边缘的有效水深数据,作为特征数据。7.3 数字水深模型生成7.3.1 3.1内插建模采用经过前期处理的水深数据,选择合适的算法,内插生成规则格网数字水深模型。内插建模应遵循以下要求:a)所有水深数据和特征数据都应参与内插建模;b)数字水深模型格网大小应符合CH/Z9026的规定;c)对预处理后的水深数据内插建模,内插算法可选用最小曲率法、改进谢别德法、线性插值三角网法等
14、常用内插算法见附录A。7.3.2 数字水深模型编辑利用特征数据对数字水深模型进行检查与编辑,使数字水深模型更加符合水下地形变化。7.3.3 空白水域及水深推测区处理空白水域的格网水深值赋予-9999,数字水深模型中达不到规定精度的水域划为水深推测区。空白水域和水深推测区使用矢量面分别进行标识。7.4 拼接与裁切对生成的数字水深模型进行拼接和裁切,相同比例尺数字水深模型拼接与裁切应满足以下要求:a)相邻数字水深模型拼接时,不少于2行(列)重合格网点;相邻数字水深模型拼接时,应对重合格网点水深差值进行检查,当差值不大于2倍数字水深模型成果中误差时,取重合格网点水深平均值作为格网点的最终水深值;当
15、差值大于2倍数字水深模型成果中误差时,需查明原因进行处理直至符合要求,并保留接边精度报告;b)跨投影带拼接时,将数字水深模型投影到同一分带,再进行拼接;c)对数字水深模型进行拼接时,检查有无漏洞、缝隙;d)根据成图比例尺,按照CH/Z9026规定的范围对数字水深模型进行裁切,生成标准分幅的数字水深模型单元;裁切时,可根据设计要求外扩;e)存储数字水深模型数据时,数据的存储格式应符合CH/Z9026的规定。7.5 相关文件制作相关文件制作应符合以下要求:a)从水深数据中提取的特殊地物和特征数据以点状、线状、面状要素文件分别存储;b)对空白水域和水深推测区标识的结果以面状要素文件分别存储;c)元数
16、据的内容、结构和格式应符合CH/Z9026的要求;d)技术总结的编写应符合CH/T1001的要求。8质量控制与成果验收8. 1质量控制检查验收和质量处理应按照GB/T18316和GB/T24356相关章节内容执行。数据成果质量应通过二级检查一级验收,即依次通过测绘单位作业部门的过程检查,测绘单位质量管理部门的最终检查和项目管理单位组织的验收或委托具有资质的质量检验机构进行质量验收。质量控制应有以下方面的内容:a)测绘基准使用的正确性;b)数据分析的合理性,分析记录的齐全性,分析结果准确性;c)水深数据融合的规范性、正确性,融合后水深数据的时效性;d)特殊地物和特征数据选取的准确性和齐全性;e)
17、内插建模算法的合理性;f)模型编辑的正确性;g)模型精度是否符合设计书要求;h)空白水域及水深推测区处理的正确性;i)模型拼接裁剪的正确性和合理性;j)数字水深模型覆盖范围的完整性,格网尺寸的正确性;k)元数据制作的规范性和正确性;D技术总结等文档资料编写的全面性和规范性;m)提供成果资料的完整性。8.2成果验收按照CH/TXXXX的规定对数字水深模型成果进行产品检验,出具验收报告。9成果整理与提交成果资料按CH/Z9026的规定进行整理,经检查无误后提交,主要包括:a)数据文件D原始水深数据;2)控制测量数据、水位资料;3)数字水深模型数据;4)特殊地物文件;5)特征数据文件;6)空白水域;
18、7)水深推测区;8)元数据;9)图幅接合表。b)文档资料D成果清单;2)技术设计书;3)技术总结;4)质量检查报告;5)验收报告。c)其他相关资料。(资料性)主要内插算法A.1线性插值三角网法线性插值三角网法是指在数据点间构建出DeIalmay三角网,拟合点值由所在三角形顶点值线性内插计算得到。该方法适用于中等数量的数据,计算方法见公式(A.1):P-(1-u-v)P1+UP2+VP3(A.1)式中:Pl、P2、P3三角形的三个顶点;P三角形内的拟合点;uPl对P点的权重贡献;VP3对P点的权重贡献。根据线性内插原理求解、V,计算方法见公式(A.2):Pr=(l-w-v)Pl.ywP2.r+v
19、P3.r(A.2)P.v=(1-w-v)P,y+wP2.lvP3,v式中:P.xPl.xP2.、P3xP、Pl、P2、P3点的X轴坐标;P.,、Pl.y、P2.y、P3yp、Pl、P2、P3点的y轴坐标。A.2最小曲率法最小曲率法是指采用全区的数据进行网格化,生成到各样点差值最小的拟合曲面。该方法适用于大量数据,算法原理如下:建立带张力的最小曲率微分方程,具体见公式(A.3):(l-)v2(v2Z)-()v2Z=O(A.3)式中:V2拉普拉斯算子;Ti内部张力系数,取值区间为0,1);Z拟合曲面函数。所满足的自然边界条件见公式(A.4):2ZZ(A.4)(19)。+W)=Oonon5(V2Z)
20、U、n式中:Tb边界张力系数,取值区间为0,1);n边界线法方向。在边界角点满足公式(A.5):2Zxy(A.5)式中:X边界角点X方向;y边界角点y方向。A.3改进谢别德法A.3.1关键函数改进谢别德法是指以反距离加权法为基础,改进了插值过程中的距离权重函数,使其只能在局部范围内起作用,其关键函数见公式(A.6):W%(x)Z(X)/(X)=上匕(A.6)r=1式中:/(%)改进谢别德法插值函数;W1(x)权重函数;q/x)局部近似函数;n全部的采样点个数。A.3.2权重函数各采样点权重函数的大小取决于到插值点的距离。当远于设定的影响半径时,权重函数为0,即只有离插值点近的采样点才会有权重
21、权重函数的计算方法见公式(A.7):p-12=Riidl(x)J,(A.7)Pdl(x)Rw式中:W(x)权重函数;di(%)点天到插值点的距离;Rw影响半径。A.3.3局部近似函数设方为三维空间内已知采样点M的属性值,用二次多项式构建局部近似函数%,使0满足Sa)=/,使用最小二乘法求取S,计算公式见(A.8)、(A.9):minZW2A=Iki(E)-4)%=R:dg04(JW4(怎)叫)(A.8)(A.9)式中:Wki权重系数;dk(xz)点Xk到天之间的距离;R点七附近用于最小二乘拟合的半径。参考文献1 GB/T1.1-2020标准化工作导则第1部分:标准化文件的结构和起草规则2 GB12327-2022海道测量规范3 CH/T1015.2-2007基础地理信息数字产品1:100001:50000生产技术规程第2部分:数字高程模型(DEM)4 CH/T7003内陆水域水下地形测量技术规程5 CH/T9020.2-2013基础地理信息数字产品1:5001:10001:2000生产技术规程第2部分:数字高程模型(DEM)6张立华,贾帅东等.数字水深模型构建理论与方法.北京:测绘出版社,2019年