#include <math.h>
#include <iostream>
using namespace std;
#define M_PI 3.14159265358979323846
double* MillierConvertion(double lat, double lon)
{
double L = 6381372 * M_PI * 2;//地球周长
double W = L;// 平面展开后，x轴等于周长
double H = L / 2;// y轴约等于周长一半
double mill = 2.3;// 米勒投影中的一个常数，范围大约在正负2.3之间
double x = lon * M_PI / 180;// 将经度从度数转换为弧度
double y = lat * M_PI / 180;// 将纬度从度数转换为弧度
y = 1.25 * log(tan(0.25 * M_PI + 0.4 * y));// 米勒投影的转换
// 弧度转为实际距离
x = (W / 2) + (W / (2 * M_PI)) * x;
y = (H / 2) - (H / (2 * mill)) * y;
double* result = new double[2];
result[0] = (int)x;
result[1] = (int)y;
return result;
}
double* MillierConvertion1(double x, double y)
{
double L = 6381372 * M_PI * 2;//地球周长
double W = L;// 平面展开后，x轴等于周长
double H = L / 2;// y轴约等于周长一半
double mill = 2.3;// 米勒投影中的一个常数，范围大约在正负2.3之间
double lat;
lat = ((H / 2 - y) * 2 * mill) / (1.25 * H);
lat = ((atan(exp(lat)) - 0.25 * M_PI) * 180) / (0.4 * M_PI);
double lon;
lon = (x - W / 2) * 360 / W;
double* result = new double[2];
result[0] = lon;
result[1] = lat;
return result;
}
void main()
{
double a, b;
double *test;
double *tet;
a = 0;
b = 0;
tet = MillierConvertion(a, b);
cout << tet[0] << endl << tet[1];
cout << endl;
test = MillierConvertion1(tet[0], tet[1]);
cout << test[0] << endl << test[1];
}
利用米勒投影法实现经纬度和平面坐标的相互转换
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 一. 投影：从不规则梨形球体到平面 我们日常所用的地图是平面地图。但地球实际上是一个赤道略宽两极略扁的不规则的梨形...
- 转自https://www.cnblogs.com/arxive/p/6017260.html 一、基本概念 首先...
