中心圆锥曲线¶
5.0.0 新版功能.
这是圆锥体切线的中心投影 :option:lat_1` `纬度，与 ``conic() 投影自 mapproj R包装。
Classification |
圆锥曲线 |
可用表格 |
正反向球面投影 |
限定区域 |
全球，但最好在标准平行线附近使用 |
Alias |
ccon公司 |
Domain |
二维 |
输入类型 |
大地坐标 |
输出类型 |
投影坐标 |
项目字符串： +proj=ccon +lat_1=52 +lon_0=19¶
使用¶
这种简单的投影很少使用，因为它不是等距的、等面积的，也不是共形的。
使用的一个例子（以及在proj4中实现该投影的主要原因）是波兰的ATPOL地球植物网格，该网格于1970年代在波兰克拉科夫的Jagiellonian大学植物研究所开发 [Zajac1978] . 这个网格最初是手写在纸质地图上，然后再手工复制。在DOS时代的第一个软件开发过程中，投影（连同奇怪的地球半径）被它的创造者选为适合现有地图的折衷方案。许多年后，它仍然是事实上的标准网格在波兰地植物学研究。
ATPOL坐标可通过以下参数获得：
+proj=ccon +lat_1=52 +lon_0=19 +axis=esu +a=6390000 +x_0=330000 +y_0=-350000
有关详细信息，请参阅 [Komsta2016] 和 [Verey2017] .
参数¶
要求的¶
- +lat_1=<value>¶
标准平行投影。
可选的¶
- +lon_0=<value>¶
投影中心的经度。
默认为0.0。
- +x_0=<value>¶
假东距。
默认为0.0。
- +y_0=<value>¶
假北距。
默认为0.0。
数学定义¶
正向投影¶
\[r=\cot\phi\u 0-\tan（\phi-\phi\u 0）\]
\[x=r\sin（\lambda\sin\phi\u 0）\]
\[y=\cot\phi\u 0-r\cos（\lambda\sin\phi\u 0）\]
反投影¶
\[y=\cot\phi\u 0-y\]
\[\phi=\phi\u 0-\tan^{-1}（\sqrt{x^2+y^2}-\cot\phi\u 0）\]
\[\lambda=\frac{\tan^{-1}\sqrt{x^2+y^2}}{\sin\phiu 0}\]
参考值¶
对于ATPOL to WGS84测试，运行以下脚本：
#!/bin/bash
cat << EOF | src/cs2cs -v -f "%E" +proj=ccon +lat_1=52 +lat_0=52 +lon_0=19 +axis=esu +a=6390000 +x_0=330000 +y_0=-350000 +to +proj=longlat
0 0
0 700000
700000 0
700000 700000
330000 350000
EOF
结果应该是
1.384023E+01 5.503040E+01 0.000000E+00
1.451445E+01 4.877385E+01 0.000000E+00
2.478271E+01 5.500352E+01 0.000000E+00
2.402761E+01 4.875048E+01 0.000000E+00
1.900000E+01 5.200000E+01 0.000000E+00
可以运行类似的脚本进行反向测试：
cat << EOF | src/cs2cs -v -f "%E" +proj=longlat +to +proj=ccon +lat_1=52 +lat_0=52 +lon_0=19 +axis=esu +a=6390000 +x_0=330000 +y_0=-350000
24 55
15 49
24 49
19 52
EOF
结果如下：
6.500315E+05 4.106162E+03 0.000000E+00
3.707419E+04 6.768262E+05 0.000000E+00
6.960534E+05 6.722946E+05 0.000000E+00
3.300000E+05 3.500000E+05 0.000000E+00