利用米勒投影法实现经纬度和平面坐标的相互转换

常用的地球经纬度与平面坐标的转换方法有米勒投影、墨卡托投影、横轴墨卡托投影（也叫UTM投影，百度地图api就用的是这个）、高斯-克吕格投影、Lambert等角正割圆锥投影等

因为我这里需求是地球经纬度坐标转平面笛卡尔坐标，所以选和墨卡托投影方式类似的米勒投影。这种投影简单是说，就是假设有一个和赤道垂直的圆柱套在地球上，然后在地心点亮一盏灯，灯光将地球各个点投影在圆柱上，在把圆柱展开，就得到地球的平面投影了，示意图如下：

使用这种方式得到的投影地图在两极会拉长，如图所示：

米勒投影和墨卡托投影类似，只是在几点区域面积变形不如后者。使用米勒投影法实现地球经纬度坐标与笛卡尔平面坐标互相转换的C++源码如下：

#include <math.h>
#include <iostream>
using namespace std;
#define M_PI       3.14159265358979323846double* 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];
}

利用米勒投影法实现经纬度和平面坐标的相互转换相关推荐

投影法字符分割matlab,利用投影法进行字符分割
1.概述字符分割有很多方法,但并不是每一种方法是万能的,那么就需要根据自己的需要来分析.例如:我现在项目的需求是将一串编号给切分开来.查了网上的资料和文献,大致适合项目的有两种方法:投影分割法和连通 ...
【OpenCV】利用投影法进行字符分割
1.概述字符分割有很多方法,但并不是每一种方法是万能的,那么就需要根据自己的需要来分析.例如:我现在项目的需求是将一串编号给切分开来.查了网上的资料和文献,大致适合项目的有两种方法:投影分割法和连通 ...
# 畸变矫正_边缘投影法对文本图像矫正——python
图像矫正就是通过一些已知的参考点,即无失真图象的某些象素点和畸变图象相应象素的坐标间对应关系,拟合出映射关系中的未知系数,并作为恢复其它象素的基础. 1.矫正原理在进行图片矫正时,有些图片具有小角度 ...
墨卡托投影法将经纬度转换为平面坐标
地球经纬度转换为平面坐标,想必是数学建模里面经常会遇到的问题. 由于地球是一个近似椭圆,因此用经纬度确定的坐标,无法直接使用平面几何的计算公式计算距离等数据.使用墨卡托投影法将经纬度坐标投影为平面坐标 ...
python 图像倾斜校正_边缘投影法对文本图像矫正——python
图像矫正就是通过一些已知的参考点,即无失真图象的某些象素点和畸变图象相应象素的坐标间对应关系,拟合出映射关系中的未知系数,并作为恢复其它象素的基础. 1.矫正原理在进行图片矫正时,有些图片具有小角度 ...
基于灰度投影法的图片偏移测量
目录前言一.灰度投影法简介二.代码实现总结前言前几天导师给了个需求,要求给出经过运动放大后的视频和原视频对应帧图片的整体偏移量,上网搜索,似乎光流法非常流行,但是还有一种灰度投影法看大家评 ...
曲面积分的投影法_在家学|第一类曲面积分与第二类曲面积分的计算
利用投影法计算第一类曲面积分设函数为定义在曲面上的连续函数.曲面的方程为.具有对和的连续偏导数,即此曲面是光滑的,且其在平面上的投影为可求面积的.则如果曲面由方程给出,在坐标面上的投影区域为,函数 ...
曲面积分的投影法_曲线曲面积分与重积分知识点汇总
本文源自扬哥去年生日发的推送, 主要梳理曲线曲面积分与重积分的各种计算方法以及对应的一些联系, 题目多数来自每日一题与裴礼文, 还有部分为华师大课本例题. 计算题中, 对称性要放在战略的高度. 其次, ...
曲面积分的投影法_大学数学：第一、二型曲面积分：投影、高斯公式补面挖点怎么用...
两类线面积分一直都是考查重点,主要考查各类线面积分的计算,这部分重点是第二类积分的计算.平面曲线积分与路径无关的判定和相关计算与证明,要熟练掌握与格林公式.高斯公式相关的各类题型及思路.方法.技巧等. ...

利用米勒投影法实现经纬度和平面坐标的相互转换

利用米勒投影法实现经纬度和平面坐标的相互转换相关推荐

最新文章

热门文章