c++高斯投影正反算_论文推荐丨闫瑾，杨绚，李妮，龚光红：地图投影变形球面大圆弧的度量指标...

《测绘学报》

构建与学术的桥梁 拉近与权威的距离

复制链接，关注《测绘学报》抖音！

【测绘学报的个人主页】长按复制此条消息，长按复制打开抖音查看TA的更多作品##7NsBSynuc88##[抖音口令]

地图投影变形球面大圆弧的度量指标

摘要：地图投影是地图学的重要研究内容。任何地图投影都不可避免地存在变形问题。针对地图投影的变形，本文提出球面大圆弧和互补比率均值相结合的地图投影面积变形与形状变形指标。通过算例验证和相关性分析，大圆指标一方面简化了小圆指标(即互补比率均值)的计算过程，并能与小圆指标的结果保持一致；另一方面，大圆指标与微分指标之间也具有较高的一致性(形状变形指标的皮尔森积矩相关系数大于0.988)。由于大圆指标不依赖于微分计算，且计算简捷，因此大圆指标更具通用性。本文进一步采用回归分析对大圆指标进行分析，结果表明，大圆指标与微分指标具有较好的线性关系(线性回归的平均误差小于1.10‰)。为了降低采样点数量和解决采样点不统一问题，本文还提出了基于随机采样的指标计算方法，并对随机方法进行了验证和分析。依据大圆指标与微分指标的一致性和线性关系，可以认为使用大圆指标能够有效地评估地图投影的变形情况。

关键词：地图投影 几何变形 互补比率均值 大圆指标 回归分析 质量评估

引文格式：闫瑾, 杨绚, 李妮, 等. 地图投影变形球面大圆弧的度量指标. 测绘学报，2020，49(6)：711-723. DOI: 10.11947/j.AGCS.2020.20190101.

阅读全文：http://xb.sinomaps.com/article/2020/1001-1595/2020-6-711.htm

全文概述

利用地图投影的微分性质能够研究地图投影的变形^[23]。基于微分的方法依据地图投影正解公式的偏导数，计算得到长度、面积或形状变形系数^[24]。微分指标是一种传统的方法，但该指标计算过程相对复杂(需要进行偏导数的计算)，且要求被评估的地图投影具有连续可微的正解公式。对于直接球面剖分(direct spherical subdivision)方法，其构造是通过递归的方式进行的，不存在任意等级格网的通用构造公式，即不存在连续的地图投影公式。文献[25]给出的结果也表明直接球面剖分的互补比率均值指标并不连续，因此微分指标在使用中也存在应用受限的问题。面积变形系数也可以根据球面微圆的边界参数方程，对球面微圆所在区域进行数值积分求得^[25]，但计算过程也较复杂，且积分运算一般不存在统一的形式。

本文拟提出一种不依赖于微分或积分运算的大圆面积指标，用于计算地图投影的面积变形系数，同时在文献[25]中提出的小圆指标(即互补比率均值指标)基础上简化得到大圆形状指标，用于计算地图投影的形状变形系数。基于大圆面积指标和大圆形状指标，结合随机采样方法，本文将针对不同类型的地图投影，开展评估、比较和分析，验证指标的有效性，并分析地图投影的质量。分析结果表明，大圆指标能够有效降低小圆指标和微分指标的计算量，并能够有效地评估地图投影质量。

1 不依赖微分计算的地图变形指标

1.1 互补比率均值

互补比率均值(averaged ratio between complementary profiles)^[25]ρ_c^s(文献[25]中以ρ_CP表示)是一种评价球面格网形状变形的方法，该指标由两项构成。为了区分小圆指标^[25](带有上标s)与本文中的大圆意义下指标(带有上标g)，本文中以ρ_c^s表示文献[25]中的ρ_CP，其计算公式为

(1)

式中，P_i、P和P_j这3点的球面小圆弧长；P_iPP_j表示在小圆弧P点的单位切向量；点P为采样区域中心点对应到球面上的点(采样区域中心点的定义见第2节)；点P₁和P₃为采样区域某一条对角线的两端点对应到球面上的点；点P₂和P₄为采样区域另一条对角线的两端点对应到球面上的点。如图 1、图 2所示，图 1(a)中的球面四边形与图 2(c)中的正方形区域相对应，图 1中红色星号表示点P，点P与图 2(c)中心的黑色星号相对应，图 1中洋红色圆形表示大圆弧a)中小圆弧通过点P；图 1(c)为图 1(b)的局部放大图。图 1(b)中的P₁P₃、P₂P₄、P₁PP₃和P₂PP₄并非长度为1的单位向量，只是为了示意各个向量的方向，计算互补比率均值时需要使用单位向量。

图 1 球面大圆与球面小圆(即互补比率均值)指标 Fig. 1 Diagram of great circle arcs and small circle arcs (i.e., average ratio between complementary profiles^[25]) based metrics

图选项

图 2 彭纳投影的底索指标和采样点 Fig. 2 Tissot's indicatrices of Bonne projection and sample points of Bonne projection

图选项

小圆指标^[25]需要计算两个球面小圆弧长P₁PP₃与P₂PP₄)，其中P₁PP₃计算过程由式(2)给出^[25]

(2)

式中，p、p₁和p₃分别为点P、P₁和P₃的笛卡儿坐标；·表示向量的点积；×表示向量的叉积；|*|表示向量的模长。计算P₂PP₄和p₁和p₃替换为p₂和p₄(为了计算上的便利，式(1)和式(2)中采用半径为1的单位球，下文中的式(3)—式(6)、式(9)—式(12)同样使用单位球进行计算)。

文献[25]利用互补比率均值，着重分析和比较了球面离散格网(而非地图投影)的质量。本文基于互补比率均值和球面大圆弧，适当地简化互补比率均值指标，同时提出一种不依赖于微分计算的面积变形指标(见1.2节)，并检验形状变形指标和面积变形指标对评价地图质量的有效性。

1.2 基于大圆弧的形状变形指标与面积变形指标

等角特征和等积特征是评估地图投影变形的两个独立的方面。本节将分别给出用于评估地图投影形状变形和面积变形的指标，关于两个指标的算例分析将在第2节中给出。

首先，本文针对文献[25]提出的小圆形状变形指标(式(1))进行简化，由式(3)给出球面大圆意义下的互补比率均值定义

(3)

式中，