ArcGIS经纬度下算面积是平方度：做距离面积缓冲区前先投影，单位错误结果全错

一、问题的提出：一个让无数GIS新手翻车的"平方度"

如果你在ArcGIS Pro或ArcMap中加载一份WGS84经纬度矢量数据，右键属性表新建一个双精度字段，然后用"计算几何"（Calculate Geometry）求面积，系统会弹出一个单位下拉框——你会发现面积单位列表里赫然出现"平方度"（Square Degrees）这个选项，而且往往就是默认值。点击确定，属性表里填进去的是一串看起来很正常、但完全没有物理意义的数字。更隐蔽的是，这个错误不会报错、不会警告，数据看起来"算出来了"，直到你把结果交给业务方，对方问"这个面积是多少平方公里"时，你才意识到问题的严重性。

这个场景在GIS教学和工程实践中出现的频率极高。笔者在多个技术社区和培训场合观察到，即便是有一两年经验的从业者，也未必能说清楚"平方度"到底错在哪里、为什么错、以及正确的做法是什么。问题的根源不在于软件操作，而在于对地理坐标系（Geographic Coordinate System, GCS）与投影坐标系（Projected Coordinate System, PCS）本质差异的理解缺失。

本文的分析主线非常明确：任何涉及距离、面积、缓冲区的空间分析，都必须先投影，后计算。这条主线将贯穿全文的每一个章节。我们将从大地测量学基础出发，逐层解剖平方度陷阱的数学根源，然后针对缓冲区、距离、面积三大高频场景给出可落地的投影选型策略和操作路径，最后讨论前沿进展与工程质检方法。全文力求做到"理论讲透、操作可复现、结论可检验"。

二、理论根基：地理坐标系与投影坐标系的量纲本质

2.1 地理坐标系：用角度描述地球表面位置

地理坐标系是基于三维椭球面的坐标系统，用经度（Longitude, λ）和纬度（Latitude, φ）两个角度量来定位地球表面任意一点。经度范围通常为-180°至180°，纬度范围为-90°至90°。这里的关键词是"角度"——地理坐标系的坐标单位是度（Degree），不是米，不是千米，更不是任何长度单位。

一个常见的误解是认为"1度等于111公里"。这个说法只在赤道附近、且仅对纬度方向近似成立。实际上，地球是一个两极稍扁、赤道略鼓的椭球体，1度纬度对应的地面距离随纬度变化而变化，1度经度对应的地面距离则随纬度升高而急剧缩短。在赤道，1度经度约111.32公里；在纬度60°处，1度经度只有约55.8公里；到了极点附近，1度经度趋近于0。这种空间变异性，正是地理坐标系不能直接用于距离和面积计算的物理原因。

目前主流的全球地理坐标系包括WGS84（World Geodetic System 1984）和CGCS2000（中国2000国家大地坐标系）。WGS84由美国国防部建立，是全球GPS定位的基准；CGCS2000由中国国家测绘局建立，自2008年起作为中国法定大地坐标系全面启用。两者在椭球参数上极为接近，扁率差异在毫米级，对绝大多数工程应用而言可以认为等价。此外，还有欧洲的ETRS89、日本的JGD2011等区域基准。

2.2 投影坐标系：把曲面"压平"到平面上

投影坐标系是在地理坐标系的基础上，通过数学变换将椭球面上的点映射到平面上，从而获得以长度单位（米、英尺等）表示的平面坐标（X, Y）。这个过程就是"地图投影"。投影坐标系 = 地理坐标系 + 投影方法 + 投影参数。

地图投影的核心矛盾在于：球面不可展。你无法在不产生变形的前提下，把地球表面完整地铺平到一张纸上。因此，任何投影都会在形状、面积、距离、方向四个属性中至少牺牲一个。这就是为什么存在数百种投影方法——每种投影都是针对特定用途、特定区域的权衡方案。

投影类型	代表投影	保持属性	典型用途
等积投影	Albers、Lambert Azimuthal Equal Area	面积	面积统计、土地利用
等距投影	Equidistant Conic、Plate Carrée	特定方向距离	距离量测、缓冲区
等角投影	Mercator、UTM、Web Mercator	角度/形状	导航、地形图
等方位投影	Stereographic、Gnomonic	从中心点出发的方向	航空、地震

上表梳理了四大类投影的属性保持特征。需要强调的是，没有一种投影能同时保持面积和距离都准确。等积投影保面积但距离变形，等距投影保特定方向距离但面积变形。这意味着，做面积分析要选等积投影，做距离/缓冲区分析要选等距投影或局部变形极小的投影（如UTM）。选错投影类型，即使单位是"米"，结果依然可能错得离谱。

2.3 中国常用坐标系速查

在中国做GIS项目，以下几套坐标系是必须熟悉的：

• WGS84 (EPSG:4326)：GPS原始数据坐标系，全球通用，经纬度单位。
• CGCS2000 (EPSG:4490)：中国法定大地坐标系，经纬度单位，与WGS84在工程精度上可互操作。
• CGCS2000 3度带高斯-克吕格 (EPSG:4513-4534等)：中国大比例尺地形图标准，投影坐标，单位为米。
• Web Mercator (EPSG:3857)：Web地图事实标准，投影坐标，单位为"米"，但高纬度地区面积严重变形。
• UTM Zone 48N-53N (EPSG:32648-32653)：中国跨多个UTM带，投影坐标，单位为米，适合中小范围分析。

关于中国高斯-克吕格投影的分带规则，3度带从东经1.5°开始，每3°一带，带号n与中央经线L的关系为L = 3n；6度带从东经0°开始，每6°一带，中央经线L = 6n - 3。中国领土范围大致覆盖3度带的25带至45带，6度带的13带至23带。选带错误会导致坐标值异常，但ArcGIS通常不会主动提示。

三、平方度陷阱的数学解剖：为什么角度相乘没有物理意义

3.1 从量纲分析说起

量纲分析（Dimensional Analysis）是物理学的基本工具。面积在物理上的量纲是长度的平方（L²），单位是平方米（m²）、平方公里（km²）等。而"平方度"的量纲是角度的平方（θ²），它描述的是球面上的一块区域对应的立体角投影，与地面实际面积之间不存在固定的换算系数。

举个具体例子：一个1°×1°的经纬度网格，在赤道附近对应的地面面积约为111.32 km × 110.57 km ≈ 12308 km²；而在纬度60°附近，同样的1°×1°网格对应的面积约为111.32 km × 55.8 km ≈ 6212 km²，只有赤道处的一半左右。两者在"平方度"单位下都是1平方度，但实际面积相差一倍。这就是平方度没有物理意义的直接证据。

3.2 ArcGIS的"计算几何"到底做了什么

当你在ArcGIS中对经纬度数据执行"计算几何"求面积时，软件内部执行的是平面几何的鞋带公式（Shoelace Formula），直接对经纬度坐标对进行多边形面积计算。这个公式本身没有错，错的是输入——它把角度当成了平面直角坐标。

鞋带公式的原理是：对于多边形顶点(x₁,y₁), (x₂,y₂), ..., (xₙ,yₙ)，面积A = ½|Σ(xᵢ·yᵢ₊₁ - xᵢ₊₁·yᵢ)|。当x和y都是角度时，算出来的A单位就是"度²"，即平方度。ArcGIS忠实地报告了这个单位，但用户往往忽略了它。

值得注意的是，ArcGIS Pro在较新版本中会对地理坐标系下的面积计算给出更明显的提示，部分版本甚至会在单位下拉框中将"平方度"标记为不推荐。但ArcMap及早期版本基本没有这类提示。这种"沉默的错误"是GIS工程中最危险的一类问题。

3.3 平方度与真实面积的换算：为什么没有通用公式

有人会问：能不能给一个"1平方度等于多少平方公里"的换算系数？答案是：不能。原因在于，平方度到真实面积的换算依赖于该区域所处的纬度，而且对于非规则形状，还依赖于区域的形状和位置。

对于球面上的一个微小区域，其面积微元可以表示为dA = R²·cos(φ)·dφ·dλ，其中R为地球半径，φ为纬度，λ为经度。这个公式清楚地表明，面积微元与cos(φ)成正比——纬度越高，同样的经纬度跨度对应的实际面积越小。因此，任何试图用单一系数换算平方度的做法都是错误的。

正确的做法只有一条：先将数据投影到合适的投影坐标系，再计算面积。投影坐标系下的坐标已经是长度单位，鞋带公式算出来的就是平方米或平方公里，量纲正确，物理意义明确。

四、投影选型方法论：不同分析场景的坐标系决策树

4.1 决策树总览

投影选型不是拍脑袋，而是有章可循的决策过程。笔者总结了一套决策树，按分析类型、研究区范围、精度要求三个维度进行判断：

分析类型？
├─ 面积统计
│   ├─ 研究区 < 一个省 → 高斯-克吕格3度带 / UTM
│   ├─ 研究区 = 全国/大洲 → Albers等积投影（中国常用Albers中央经线105°E，双标准纬线25°N和47°N）
│   └─ 研究区 = 全球 → Equal Earth / Mollweide / Robinson
├─ 距离量测 / 缓冲区
│   ├─ 研究区 < 一个城市 → 高斯-克吕格3度带 / UTM（局部变形<0.1%）
│   ├─ 研究区 = 一个省 → UTM / 高斯-克吕格6度带
│   └─ 研究区 = 全国 → 分段投影 / 测地线计算（Geodesic）
└─ 导航 / 可视化
    └─ Web Mercator（仅可视化，不做量测）

这棵决策树的核心逻辑是：研究区范围越小，投影变形越小，可选投影越多；研究区范围越大，越需要选择专门针对该范围优化的投影。对于中国全国尺度的面积统计，Albers等积投影是行业标准选择，其双标准纬线通常设为25°N和47°N，中央经线105°E，这套参数能在中国范围内将面积变形控制在1%以内。

4.2 面积分析：等积投影优先

面积分析的第一原则是选择等积投影。等积投影保证投影前后面积比例不变，虽然形状会有变形，但对于面积统计而言，形状变形无关紧要。

中国常用的等积投影包括：

• Albers Equal Area Conic：最适合中国全国尺度的面积统计，EPSG:102025（Esri自定义）或按需自定义参数。
• Lambert Azimuthal Equal Area：适合大洲尺度，欧洲环境署（EEA）推荐用于欧洲面积统计。
• Equal Earth：2018年提出的新型全球等积投影，视觉效果好，适合全球专题图。

需要注意的是，高斯-克吕格和UTM都是等角投影，不是等积投影。在小范围（如一个县、一个市）内，等角投影的面积变形可以忽略不计，因此用于面积统计问题不大。但在省级以上尺度，等角投影的面积变形会显著增大，此时应改用等积投影。

4.3 距离与缓冲区：等距或局部等角投影

距离量测和缓冲区分析对投影的要求与面积分析不同。缓冲区分析的本质是"以某点为中心，向外扩展固定距离"，因此最关心的是投影在局部范围内的距离保真度。

UTM和高斯-克吕格投影在中央经线附近的距离变形极小（UTM中央经线处比例因子0.9996，即1000米实际距离在投影上为999.6米，误差0.04%），在带边缘（±3°）处变形约为0.1%。对于大多数工程应用，这个精度足够。但如果研究区跨越多个UTM带，就需要考虑分带处理或改用其他投影。

对于跨带或大范围的距离量测，更严谨的做法是使用测地线（Geodesic）计算。ArcGIS Pro提供了"测地线缓冲区"（Geodesic Buffer）工具，可以直接在地理坐标系下生成基于椭球体测地线的缓冲区，无需投影。这是目前精度最高的方案，但计算开销也更大。

4.4 一张表看懂投影选型

分析场景	推荐投影	单位	精度说明
县市级面积统计	CGCS2000 3度带	平方米	变形<0.01%
省级面积统计	Albers (105°E, 25°N, 47°N)	平方米	面积变形<1%
城市缓冲区	UTM / 高斯3度带	米	距离变形<0.1%
跨省缓冲区	测地线缓冲区	米	椭球体精确
全球面积对比	Equal Earth / Mollweide	平方米	面积精确

五、缓冲区分析实战：从经纬度到米的全流程操作

5.1 错误示范：直接对经纬度数据做缓冲区

假设你有一份WGS84的点数据，想在每个点周围生成500米缓冲区。如果你直接打开"缓冲区"工具，距离输入500，单位选择"米"，ArcGIS会怎么处理？

答案是：ArcGIS会尝试将500米换算成"度"。换算的方式是基于数据坐标系的原点纬度或某个参考纬度进行粗略估算。这个换算是不精确的，而且随纬度变化。在高纬度地区，误差可能达到数倍。更糟糕的是，如果单位选择不当，软件可能直接把500当成"度"来处理，生成一个覆盖半个地球的巨型缓冲区。

笔者在实际项目中见过最离谱的案例：某团队对全国POI数据做1公里缓冲区，结果生成的多边形面积总和超过了中国国土面积。排查后发现，数据是WGS84经纬度，缓冲区距离单位被误设为"度"，1公里被当成了1度。

5.2 正确流程：投影 → 缓冲 → 可选反投影

正确的缓冲区分析流程分为三步：

第一步：投影。使用"投影"（Project）工具，将经纬度数据转换到合适的投影坐标系。对于城市级分析，推荐CGCS2000 3度带或UTM。投影后，数据的坐标单位变为米。

第二步：缓冲。在投影坐标系下执行"缓冲区"（Buffer）工具，距离输入500，单位选择"米"。此时软件直接按平面几何计算，结果准确。

第三步：可选反投影。如果最终成果需要以经纬度形式交付，可以用"投影"工具将缓冲区结果转回WGS84。但要注意，反投影后的缓冲区边界在经纬度下看起来会有轻微变形，这是正常的。

5.3 缓冲区距离单位的陷阱

即使投影正确，缓冲区工具的距离单位设置仍有陷阱。ArcGIS的缓冲区工具允许用户选择距离单位，包括米、千米、英里、度等。如果数据是投影坐标系（单位为米），而你选择了"千米"，软件会自动换算；但如果数据是地理坐标系，选择"米"时软件会进行粗略的角度换算，精度不可靠。

笔者的建议是：永远在投影坐标系下做缓冲区，且距离单位与数据坐标单位保持一致。如果数据是米，距离就用米；如果数据是英尺，距离就用英尺。不要依赖软件的自动换算。

5.4 测地线缓冲区：跨带场景的终极方案

对于跨越多个UTM带或大范围的研究区，投影缓冲区会遇到"带边缘变形"问题。此时可以使用ArcGIS Pro的"测地线缓冲区"工具。该工具直接在地理坐标系下工作，基于椭球体测地线计算缓冲区边界，无需投影。

测地线缓冲区的原理是：以中心点为起点，沿各个方位角计算测地线，在测地线上取距离等于缓冲半径的点，连接这些点形成缓冲区边界。这种方法在数学上是精确的，但计算量比平面缓冲区大得多。对于大数据量，建议先做投影缓冲区，只在关键区域使用测地线缓冲区做精度验证。

六、距离与面积量测：工具选择、单位校验与误差控制

6.1 ArcGIS中的距离量测工具对比

ArcGIS提供了多种距离量测工具，适用场景各不相同：

工具	所在工具箱	适用坐标系	精度
计算几何	属性表右键	投影坐标系	平面几何，高
Near	分析工具 → 邻域分析	投影坐标系	平面几何，高
Point Distance	分析工具 → 邻域分析	投影坐标系	平面几何，高
测地线距离	数据管理 → 字段计算	地理坐标系	椭球体，最高
测量工具	交互式工具条	任意	依赖坐标系，中

从表中可以看出，投影坐标系下的平面几何工具精度高且计算快，是大多数场景的首选。测地线距离精度最高，但计算慢，适合做精度验证或小范围高精度分析。

6.2 面积量测的三种方法及其误差

面积量测有三种常见方法，误差特征不同：

方法一：投影坐标系下计算几何。这是最常用的方法。在合适的投影下，误差主要来自投影变形，通常小于1%。

方法二：椭球体面积计算。ArcGIS Pro的"计算几何"工具在数据为地理坐标系时，可以选择"椭球体"面积计算方式。这种方法基于椭球体参数直接计算，精度最高，但计算速度较慢。

方法三：平方度（错误方法）。如前所述，这是错误的方法，结果无物理意义。

为了直观展示三种方法的差异，笔者用一份模拟数据进行了对比。数据为一个约0.1°×0.1°的矩形，中心位于北纬40°。三种方法的结果如下：

方法	结果	单位	相对误差
平方度法	0.01	平方度	无法比较
投影法（UTM 50N）	约 72,800,000	平方米	基准
椭球体法	约 72,850,000	平方米	+0.07%

6.3 单位校验清单

每次完成距离或面积计算后，建议执行以下校验：

1. 检查数据坐标系的单位是"度"还是"米"。
2. 检查计算结果字段的单位标注是否正确。
3. 用已知尺寸的要素（如一个标准矩形）验证计算结果。
4. 将结果与公开数据（如行政区面积）交叉比对。
5. 对于大范围数据，检查是否跨越多个投影带。

七、典型误用案例剖析：Web Mercator、UTM分区与跨带问题

7.1 Web Mercator：可视化利器，量测毒药

Web Mercator（EPSG:3857）是Google Maps、OpenStreetMap、高德、百度等在线地图的事实标准。它的优点是全球覆盖、计算简单、瓦片金字塔结构规整。但它有一个致命缺陷：面积变形随纬度急剧增大。

Web Mercator采用球体近似（而非椭球体），其面积变形因子为1/cos²(φ)。在赤道，变形因子为1；在纬度60°，变形因子为4；在纬度80°，变形因子约为33。这意味着，在Web Mercator下，格陵兰岛看起来比非洲还大，而实际上非洲面积是格陵兰的14倍。

因此，Web Mercator绝对不能用于面积统计和距离量测。它只适合可视化展示。如果你从在线地图下载了Web Mercator数据，做面积分析前必须重新投影到等积投影。

7.2 UTM分区与跨带问题

UTM将全球划分为60个6度带，每个带独立投影。中国领土覆盖UTM 43N至53N带（约东经72°至138°）。跨带数据在做缓冲区或距离分析时，如果统一投影到某一个带，远离中央经线的部分变形会增大。

以UTM 50N（中央经线117°E）为例，东经126°处的距离变形约为0.15%，东经132°处约为0.35%。对于精度要求不高的项目，这个变形可以接受；但对于高精度工程，就需要分带处理或使用测地线方法。

分带处理的流程是：按UTM带将数据分组，每组投影到对应带，分别做分析，最后将结果合并。这个过程可以用ArcGIS的模型构建器或Python脚本自动化。

7.3 高斯-克吕格与UTM的差异

中国的高斯-克吕格投影与UTM都是横轴墨卡托投影的变种，但有两个关键差异：

• 比例因子不同：UTM在中央经线处比例因子为0.9996，高斯-克吕格为1.0。这意味着UTM在中央经线处距离略短，但带边缘变形更均匀。
• 分带方式不同：UTM为6度带，高斯-克吕格有3度带和6度带两种。3度带精度更高，适合大比例尺。

在实际项目中，如果数据已经是CGCS2000 3度带，直接用即可，不必转成UTM。两者在小范围内精度相当，转换反而可能引入不必要的误差。

八、前沿进展：动态投影、椭球体测地线与自动化质检

8.1 动态投影与按需投影

ArcGIS Pro引入了"动态投影"（On-the-fly Projection）机制，允许不同坐标系的数据在同一个地图中叠加显示。这极大方便了数据浏览，但也带来新的风险：地图显示正确不代表分析正确。动态投影只是视觉上的坐标转换，分析工具仍然按数据自身的坐标系计算。

近年来，一些GIS平台开始探索"按需投影"（Projection on Demand）机制，即在分析工具执行时自动检测坐标系并选择合适的投影。Esri在ArcGIS Pro 3.x中增强了这方面的提示，但尚未实现全自动。笔者认为，完全自动化的投影选型在短期内难以实现，因为投影选型涉及精度要求、研究区范围、分析类型等多重因素，需要人工判断。

8.2 椭球体测地线计算的普及

随着计算能力的提升，椭球体测地线计算（Geodesic Computation）正在从专业领域走向普及。GEOS（Geometry Engine Open Source）库自3.9版本起提供了完整的测地线计算支持，PostGIS 3.0+也集成了相关功能。这意味着，未来的GIS分析可以越来越多地直接在地理坐标系下进行，无需投影，同时保证精度。

但需要注意的是，测地线计算的计算开销远大于平面计算。对于大数据量的缓冲区分析，测地线方法可能慢几个数量级。因此，短期内"投影+平面计算"仍是主流方案，测地线计算主要用于精度验证和小范围高精度分析。

8.3 自动化质检与坐标系感知

针对平方度这类"沉默错误"，业界正在发展自动化质检工具。这些工具的核心思路是"坐标系感知"：在分析流程中嵌入坐标系检查节点，自动判断当前坐标系是否适合当前分析类型，并给出警告或建议。

例如，可以编写一个ArcGIS Python工具箱（Python Toolbox），在面积计算前自动检查数据坐标系：如果是地理坐标系，则提示用户先投影；如果是投影坐标系，则检查投影类型是否适合面积统计。这类工具的实现并不复杂，但能有效防止低级错误。

# 坐标系检查伪代码示例
import arcpy

def check_cs_for_area(feature_class):
    desc = arcpy.Describe(feature_class)
    sr = desc.spatialReference
    if sr.type == "Geographic":
        return "警告：地理坐标系，面积计算将得到平方度，请先投影"
    elif sr.type == "Projected":
        if "Equal_Area" in sr.projectionName or "Albers" in sr.projectionName:
            return "通过：等积投影，适合面积计算"
        else:
            return "注意：非等积投影，大范围面积计算可能有偏差"
    else:
        return "未知坐标系类型"

这段伪代码展示了坐标系感知质检的基本逻辑。在实际项目中，可以将其扩展为完整的质检流程，覆盖投影、缓冲区、距离、面积等所有关键环节。

九、工程检查清单与最佳实践总结

9.1 分析前检查清单

1. 确认数据坐标系类型（地理/投影）。
2. 确认数据坐标单位（度/米/英尺）。
3. 确认研究区范围（县/省/全国/全球）。
4. 确认分析类型（面积/距离/缓冲区/可视化）。
5. 根据决策树选择合适投影。
6. 检查是否需要分带处理。

9.2 分析后校验清单

1. 检查结果单位是否正确（平方米/平方公里/米）。
2. 用已知要素验证结果合理性。
3. 与公开数据交叉比对。
4. 检查是否有异常值（如面积为负、距离过大）。
5. 记录使用的投影参数，便于复现。

9.3 最佳实践总结

回到文章开头的问题：ArcGIS经纬度下算面积为什么是平方度？因为地理坐标系的坐标单位是角度，角度相乘得到平方度。这个答案简单，但背后的投影学知识体系却很深厚。希望本文能帮助读者不仅知其然，更知其所以然，在实际项目中避开这个陷阱，做出经得起检验的分析成果。

十、参考文献与拓展资源

10.1 主要参考文献

1. Snyder, J. P. (1987). Map Projections: A Working Manual. US Geological Survey Professional Paper 1395. （经典投影学手册，本文评述：至今仍是投影参数查询的首选参考）
2. Iliffe, J., & Lott, R. (2008). Datums and Map Projections for Remote Sensing, GIS and Surveying (2nd ed.). CRC Press.
3. Esri. (2024). ArcGIS Pro Documentation: Coordinate Systems and Projections. Esri Press. https://pro.arcgis.com/en/pro-app/latest/help/mapping/properties/coordinate-systems-and-projections.htm
4. Battersby, S. E., Finn, M. P., Usery, E. L., & Yamamoto, K. H. (2014). Implications of Web Mercator and its use in online mapping. Cartographica, 49(2), 85-101.
5. Šavrič, B., Patterson, T., & Jenny, B. (2018). The Equal Earth map projection. International Journal of Geographical Information Science, 33(3), 454-465.
6. Karney, C. F. F. (2013). Algorithms for geodesics. Journal of Geodesy, 87(1), 43-55.
7. CGCS2000国家大地坐标系推广使用技术指南. (2008). 中国国家测绘局.
8. GEOS Development Team. (2023). GEOS: Geometry Engine Open Source Documentation. https://libgeos.org/
9. PostGIS Development Team. (2024). PostGIS Manual: Geography Type and Geodesic Calculations. https://postgis.net/docs/

10.2 拓展学习资源

• Esri官方投影教程：https://pro.arcgis.com/en/pro-app/latest/help/mapping/properties/projection-basics.htm
• QGIS投影文档：https://docs.qgis.org/latest/en/docs/user_manual/working_with_projections/working_with_projections.html
• EPSG.io坐标系查询：https://epsg.io/
• ArcGIS Pro测地线缓冲区教程：https://pro.arcgis.com/en/pro-app/latest/tool-reference/analysis/geodesic-buffer.htm
• 空间参考与投影视频课程（Esri Academy）：https://www.esri.com/training/

10.3 数据集说明

本文涉及的模拟数据均基于WGS84椭球参数（长半轴6378137m，扁率1/298.257223563）和UTM投影参数计算，用于说明方法原理，非实测数据。实际项目中的误差大小取决于研究区位置、范围和投影参数，建议以实测数据为准。