作者: 中科图新技术团队 标签: #GIS #测绘 #DEM #三维可视化 阅读时间: 11 分钟

用地形数据的时候,经常听到一句「这份地形数据是 10 米精度」。这句话里其实压着三件不同的事:点多密、文件里点与点隔多远、整体高程偏差有多大。三者混着用,就会得出「重采样之后精度提高了」「10 米格网能看清 10 米宽的路」这类站不住的结论。本文只做一件事:把这三个数值口径分开讲清,再给一份用地形/DEM数据的核对清单。

标称分辨率:发布方声明的「点多密」

标称分辨率是数据发布时写在规格里的一个数,指的是格网间距的声明值,回答的是「每隔多远放一个高程点」,而不是「每个点有多准」。它属于产品规格,属于交付前给出的说明,与你这块测区的实际偏差之间没有直接换算关系。

哥白尼 DEM 数据集说明页的写法很典型:Digital Surface Model (DSM) provided in 3 different resolutions (90m, 30m, 10m),同一页的规格里,GLO-30 的纬度方向格网间距写作 1.0″,GLO-90 写作 3.0″。能看到「分辨率」与「格网间距」在这里是同一条目的两种写法——30 米与 1.0 角秒说的是同一个间距,只是单位不同(角秒折合的长度随纬度略有差异)。这条还提醒另一件事:它是 DSM,含建筑、基础设施与植被顶面,与经编辑去掉建筑和植被的 DTM 不是一回事。也就是说,「多密」和「代表什么表面」是两个维度,看到前者不能顺带推断后者。

所以看到「30 米分辨率」,能确定的只有采样间隔。至于高程能不能对到 ±2 米,规格里通常另有说法,要去中误差那一栏找;找不到,就只能当成未知,不能当成「差不多」。

格网间距:文件里量出来的那个距离

标称是声明,格网间距是事实。同一份数据经过裁剪、投影转换、重采样、拼接之后,文件里相邻两点的实际距离可能与标称一致,也可能已经被改过。用工具读一下文件头里的仿射参数,或者直接量两个相邻像元中心之间的距离,几秒钟就能确认。

这里容易出现一类不易察觉的错位:文件名还叫「30m」,内容已经被重采样到 5 米;或者原数据是地理坐标(单位为度),被当成投影坐标(单位为米)读取,算出来的格网间距差了几个数量级。前一种情况数据量涨了 36 倍,信息量并没有跟着涨;后一种情况连坐标都不对,任何面积与土方量算都不可用。

实务上的判断顺序建议是:先读坐标系与仿射参数,算出实际格网间距;再与标称值比对,差异超出合理舍入范围就追问处理过程;最后把结论记下来,而不是记「应该是 30 米」。这一步不花时间,但能挡住后面不少返工。

高程中误差:一个统计量,不是单点上限

高程中误差是对照高精度检查点算出来的整体偏差水平,通常是各检查点高程差的平方和取均值再开方,反映的是「这一批点平均偏多少」。它描述一个数据集、一个图幅或一个测区的总体表现,不是任意单点的误差上限。

测绘行业标准《基础地理信息 三维模型产品规范》(征求意见稿)里有一句把三者并列的表述:地形精度可通过成图比例尺、网格间距以及高程中误差三者的综合角度来表达。同一份文件还写,中误差的两倍值为最大误差,并且 DEM 内插点的高程中误差按上表规定放宽 0.2 倍。这几句放在一起,恰好说明三件事各自承担不同角色:一个是采样密度,一个是统计口径,一个是抽样条件下的边界换算。标准里把三件事放在同一句话里表达,本身就说明它们不能互相替代。

需要留意这个统计量的两个边界:检查点的数量、分布与测量精度,决定了这个数可不可信;统计范围与你要用的范围是否重合,决定了这个数能不能外推。一个全测区 1.5 米的中误差,不等于测区边缘那 200 平方米也是 1.5 米——如果那片正好落在系统性偏移区,单点偏差可能明显大于整体平均值。

等高线

在哪找:元数据、说明文件与成果报告

三项指标很少整齐地写在同一个地方,通常是一处一项,需要拼起来看。下表是按文件类型分的清点方式。

查的地方 通常能查到 常常缺
元数据(数据头或随附 xml) 坐标系、覆盖范围、实际格网间距、发布方声明的分辨率、产品类型是 DSM 还是 DTM 检查点数量与中误差来源
说明文件(产品说明、readme、数据手册) 采集方式、生产年份、高程基准、精度指标的定义口径 本项目测区的实际偏差
成果报告(生产技术报告、验收报告) 检查点数量与分布、中误差、最大误差、检查日期与方法 文件外那次重采样与裁剪

找不到就问来源方要,把问题问具体:中误差对应多少个检查点、检查点怎么布、覆盖不覆盖我要用的范围、这份交付件在交付后有没有再处理过。猜一个数写进方案,比空着风险大得多。

下面是拿到一份地形数据后的核对路径,可以先照这个顺序走一遍。

地形数据三个数值口径的核对流程与台账判定(架构图)

核对完这三项,落到软件里还有一处要跟着确认:数据的坐标系(含投影高),设错了整体会被抬起来或压下去。

数据导入时设置高程偏移量的对话框(界面示意)

分辨率与比例尺:一个换算,两个方向

三个数讲清楚之后,还有一个问题几乎每次评审都会被问到:这份地形数据,能出多大的图?

测绘成果规范里对这个是有对应的。本文引用的那份征求意见稿,在 5.4 地形精度要求之后紧接一张表,把成图比例尺、网格间距与各级地形类别的高程中误差并排列在一处:

成图比例尺 1:500 1:1000 1:2000 1:5000 1:10000 1:25000 1:50000 1:100000
网格间距(米) 0.5 1 2 2.5 5 10 25 50

两行相除,比值大致稳定:格网间距 ≈ 成图比例尺分母 ÷ 2000。反算同样成立:成图比例尺分母 ≈ 格网间距 × 2000。表里 1:500 至 1:2000 这三级是按分母 ÷1000 取值的,1:25000 一级按 ÷2500,其余四级都接近 ÷2000;用 ÷2000 粗估前三级会把格网估稀一半,落到具体某一级时,回表核一遍。

这个近似可以两头用。

从格网看能出多大的图。 手上是 30 米格网的大范围地形数据,30 × 2000 得 60000,落在表中 1:50000 与 1:100000 之间。这样一句话,比「这份数据是 10 米精度」更能对上评审要问的事。反过来算:10 米格网折算约 1:20000,5 米格网约 1:10000。

从比例尺看该多密的格网。 要出 1:2000 的图,表里对应 2 米格网。拿 30 米数据去出 1:2000 的图:30 米 ÷ 2000 得 15 毫米,图上 1.5 厘米才落一个高程点,坡面走向、沟渠位置都画不出来,图面上只剩一团起伏。

不查表也有一个自查办法:把格网间距除以比例尺分母,看它在图上有多长。 表里各级算下来,格网在图上大致是 0.4 至 1 毫米——1:500 约 1 毫米,1:10000 与 1:50000 约 0.5 毫米。细过 0.1 毫米,已经超出图面的表达下限,再密也画不出来;粗到几毫米,图面上就是一个挨一个的方格,地形要靠想象去补。反方向同样成立:比例尺分母乘 0.1 毫米,是这张图能画出的地面尺寸下限,1:500 是 5 厘米,1:10000 是 1 米,1:50000 是 5 米。

这里要说清一件事:上面换算对上的是「密度」,对不上「准不准」。 30 米格网折算 1:60000,说的是这个间距在 1:60000 的图上刚好画得出来,不等于这份高程在 1:60000 的精度要求上合格;反过来,一份 5 米格网的数据,若原始观测只有稀疏检查点,它的中误差也可能比 30 米的还大。密度管「画不画得下」,中误差管「画上去的数可不可信」,两个问题分开问。

还有一层边界:这张表是建库与出图环节的密度参照,不是分析用途的上限。做汇水、坡面、挖填这类量算时,格网够不够要按地形起伏与沟谷宽度单独判断——一条几米宽的冲沟,在 50 米格网上留不下痕迹。

三处常见混用

把格网当成分辨能力

「10 米格网」说的是每 10 米有一个高程采样点,不等于能看清 10 米宽的物体。采样间隔只约束「能不能落到点」,地物的可识别尺寸还受成像方式、点云密度、滤波与建模策略影响。错在这里,判断就会跟着偏:方案里写成「可识别 10 米以上地物」,于是拿它去定 5 米宽沟渠的断面位置、去提取 3 米宽的道路边界,成果会是一条被平滑掉的曲线。稳妥的写法是把用途限定在坡面、汇水、大范围地形起伏这类尺度明显大于格网的量算上。

把重采样当成提精度

把 30 米数据重采样到 5 米,改的是采样密度,不是信息量。新格网上的高程大多由邻域插值得到,原始观测还是那些点。错在这里,判断就会偏:以为文件变小格网就等于误差缩小,于是在 5 米版本上做土方量算并承诺更小的偏差。实际做法建议是:量算仍以原始格网为准,重采样件只用于对齐、配色或渲染;同一块区域用两个版本各算一遍,差异超出预期就回去查插值方法。

把统计量当成单点保证

中误差是对整批检查点算的,把它当成「任意一点都不会超过这个数」,是把统计口径用成了承诺。错在这里,判断就会偏:在两百平方米的基坑范围里写「高程偏差不超过 ±0.5 米」,而该范围的实际系统性偏移可能到 1 米以上,验核时对不上。可行的做法是:小范围用途单独布检查点实测,用那批点自己算偏差;在没有本区实测点之前,只写数据来源与统计口径,不写单点指标。

一处做法:三项抄进同一张台账

与其在评审前临时回忆,不如在数据进来的时候就把三项抄到同一张表上,一行一份数据,三列各管一个口径,再加一列写来源。表不用复杂,能查、能追、能交接就够。

数据名 标称分辨率 实际格网 中误差来源
某区域 DSM 成果(2024 年生产) 30 米 30 米,读仿射参数确认 验收报告,42 个检查点
甲方提供的裁剪件 文件名标注 5 米 5 米,由 30 米重采样 待补

三项不齐的先标「待补」,不要让空栏在方案里悄悄变成默认值。台账跟着数据走,谁接手谁能看懂,比口头交接可靠。

总结

  • 标称分辨率说的是采样间隔,属于发布方声明。它回答「点多密」,回答不了「多准」,也不说明代表的是地表、地表加地物还是裸地。
  • 格网间距要在文件里量出来。重采样只改采样密度,原始观测点没有增加,量算仍以原始格网为准。
  • 格网间距与成图比例尺可以互折:分母 ÷2000 量级得格网,格网 ×2000 量级得比例尺;折出来的是密度对应,不是精度结论,落到具体某一级回表核一遍。
  • 高程中误差是统计量,两倍值只是抽样条件下的边界换算,不能外推到任意小范围,也不能替代本区实测。
  • 三项各有出处,抄进同一张台账并注明检查点来源,比记在脑子里可靠;缺的那项标待补,别猜。

延伸阅读方面,想看标准原文怎么表述这三者的关系,可以按文末链接读那份征求意见稿;想看数据集规格是怎么写「分辨率」的,读哥白尼那份数据集说明页即可。

如对文中的三维地形数据加载与量测感兴趣,可了解图新说在这方面的支持,详见产品文档。


中科图新 | 让每一寸空间都有数字灵魂

参考来源: - Copernicus DEM 数据集说明(Copernicus Data Space Ecosystem) - 《基础地理信息 三维模型产品规范》(征求意见稿)