java 弧度 转经纬度_GPS 经纬度转换为 经过旋转后的平面坐标详解
/**
* 地图工具
*
* @author dxm
*
*/
public class MapUtil {
// 参考点
private double[] d34 = null;
private double[] d33 = null;
private double[] gs34 = null;
private double[] gs33 = null;
// 参考点位置
private double x34, y34;
// 计算参数
private double a, b, c, bl;
/**
* 构造器
*
* @param ds34 参考点经纬度
* @param ds33 参考点经纬度
* @param length 两个参考点在本坐标系长度
*/
public MapUtil(String[] ds34, String[] ds33, double x34, double y34, double length) {
this.x34 = x34;
this.y34 = y34;
d34 = convertGPS(ds34[0], ds34[1]);
d33 = convertGPS(ds33[0], ds33[1]);
gs34 = BLToGauss(d34[0], d34[1]);
gs33 = BLToGauss(d33[0], d33[1]);
// 点33到点34的高斯横向距离
a = gs34[0] - gs33[0];
// 点33到点34的高斯纵向距离
b = gs34[1] - gs33[1];
// 点33到点34的高斯直线距离
c = Math.sqrt(a * a + b * b);
// 高斯距离与本坐标系的比例
bl = length / c;
}
/**
* 根据经纬度算XY
*
* @param lng
* @param lat
* @return
*/
public int[] getXY(double lng, double lat) {
// 与参考点比较经度计算坐标
if (lng > d34[0]) {
return getXY1(lng, lat);
} else {
return getXY2(lng, lat);
}
}
/**
* 已有定点为基准算XY坐标
*/
private int[] getXY1(double lng, double lat) {
// 将经纬度转换为高斯投影坐标
double[] gs = BLToGauss(lng, lat);
// 与点34的横向高斯距离
double d = Math.abs(gs[0] - gs34[0]);
double e = d * c / b;
double f = a * e / c;
// 与点34的纵向高斯距离
double gf = Math.abs(gs[1] - gs34[1]);
double g = gf - f;
double h = a * g / c;
double i = b * g / c;
double x = x34 - i * bl;
double y = y34 + (h + e) * bl;
return new int[] { (int) Math.rint(x), (int) Math.rint(y) };
}
/**
* 已有定点为基准算XY坐标
*/
private int[] getXY2(double lng, double lat) {
// 将经纬度转换为高斯投影坐标
double[] gs = BLToGauss(lng, lat);
// 与点34的横向高斯距离
double d = Math.abs(gs[0] - gs34[0]);
double e = d * c / a;
double f = e * b / c;
// 与点34的纵向高斯距离
double gf = Math.abs(gs[1] - gs34[1]);
double g = gf - f;
double h = a * g / c;
double i = b * g / c;
double x = x34 - (e + i) * bl;
double y = y34 + h * bl;
return new int[] { (int) Math.rint(x), (int) Math.rint(y) };
}
/**
* 经纬度度分秒转带小数点的度
*
* @param jw
* @return
*/
private double[] convertGPS(String lng, String lat) {
lng = lng.replaceAll("N", "");
lng = lng.replaceAll("S", "");
lng = lng.replaceAll("E", "");
lng = lng.replaceAll("W", "");
lat = lat.replaceAll("N", "");
lat = lat.replaceAll("S", "");
lat = lat.replaceAll("E", "");
lat = lat.replaceAll("W", "");
double lngD = 0;
String[] fs = lng.split("°");
if (fs.length > 0) {
lngD = Double.parseDouble(fs[0]);
}
fs = fs[1].split("′");
if (fs.length > 0) {
double fen = Double.parseDouble(fs[0]);
fen = fen / 60;
lngD += fen;
}
fs = fs[1].split("″");
if (fs.length > 0) {
double mi = Double.parseDouble(fs[0]);
mi = mi / (60 * 60);
lngD += mi;
}
double latD = 0;
fs = lat.split("°");
if (fs.length > 0) {
latD = Double.parseDouble(fs[0]);
}
fs = fs[1].split("′");
if (fs.length > 0) {
double fen = Double.parseDouble(fs[0]);
fen = fen / 60;
latD += fen;
}
fs = fs[1].split("″");
if (fs.length > 0) {
double mi = Double.parseDouble(fs[0]);
mi = mi / (60 * 60);
latD += mi;
}
return new double[] { lngD, latD };
}
/**
* 由高斯投影坐标反算成经纬度
*
* @param X
* @param Y
* @return
*/
@SuppressWarnings("unused")
private double[] GaussToBL(double X, double Y) {
int ProjNo;
int ZoneWide; // //带宽
double[] output = new double[2];
double longitude1, latitude1, longitude0, X0, Y0, xval, yval;// latitude0,
double e1, e2, f, a, ee, NN, T, C, M, D, R, u, fai, iPI;
iPI = 0.0174532925199433; // //3.1415926535898/180.0;
// a = 6378245.0; f = 1.0/298.3; //54年北京坐标系参数
a = 6378140.0;
f = 1 / 298.257; // 80年西安坐标系参数
ZoneWide = 6; // //6度带宽
ProjNo = (int) (X / 1000000L); // 查找带号
longitude0 = (ProjNo - 1) * ZoneWide + ZoneWide / 2;
longitude0 = longitude0 * iPI; // 中央经线
X0 = ProjNo * 1000000L + 500000L;
Y0 = 0;
xval = X - X0;
yval = Y - Y0; // 带内大地坐标
e2 = 2 * f - f * f;
e1 = (1.0 - Math.sqrt(1 - e2)) / (1.0 + Math.sqrt(1 - e2));
ee = e2 / (1 - e2);
M = yval;
u = M / (a * (1 - e2 / 4 - 3 * e2 * e2 / 64 - 5 * e2 * e2 * e2 / 256));
fai = u + (3 * e1 / 2 - 27 * e1 * e1 * e1 / 32) * Math.sin(2 * u) + (21 * e1 * e1 / 16 - 55 * e1 * e1 * e1 * e1 / 32) * Math.sin(4 * u) + (151 * e1 * e1 * e1 / 96) * Math.sin(6 * u) + (1097 * e1 * e1 * e1 * e1 / 512) * Math.sin(8 * u);
C = ee * Math.cos(fai) * Math.cos(fai);
T = Math.tan(fai) * Math.tan(fai);
NN = a / Math.sqrt(1.0 - e2 * Math.sin(fai) * Math.sin(fai));
R = a * (1 - e2) / Math.sqrt((1 - e2 * Math.sin(fai) * Math.sin(fai)) * (1 - e2 * Math.sin(fai) * Math.sin(fai)) * (1 - e2 * Math.sin(fai) * Math.sin(fai)));
D = xval / NN;
// 计算经度(Longitude) 纬度(Latitude)
longitude1 = longitude0 + (D - (1 + 2 * T + C) * D * D * D / 6 + (5 - 2 * C + 28 * T - 3 * C * C + 8 * ee + 24 * T * T) * D * D * D * D * D / 120) / Math.cos(fai);
latitude1 = fai - (NN * Math.tan(fai) / R) * (D * D / 2 - (5 + 3 * T + 10 * C - 4 * C * C - 9 * ee) * D * D * D * D / 24 + (61 + 90 * T + 298 * C + 45 * T * T - 256 * ee - 3 * C * C) * D * D * D * D * D * D / 720);
// 转换为度 DD
output[0] = longitude1 / iPI;
output[1] = latitude1 / iPI;
return output;
}
/**
* 由经纬度反算成高斯投影坐标
*
* @param longitude
* @param latitude
*/
private double[] BLToGauss(double longitude, double latitude) {
int ProjNo = 0;
// 带宽
int ZoneWide = 6;
double longitude1, latitude1, longitude0, X0, Y0, xval, yval;
double a, f, e2, ee, NN, T, C, A, M, iPI;
// 3.1415926535898/180.0;
iPI = 0.0174532925199433;
// 54年北京坐标系参数
a = 6378245.0;
f = 1.0 / 298.3;
// 80年西安坐标系参数
// a=6378140.0;
// f=1/298.257;
ProjNo = (int) (longitude / ZoneWide);
longitude0 = ProjNo * ZoneWide + ZoneWide / 2;
longitude0 = longitude0 * iPI;
// 经度转换为弧度
longitude1 = longitude * iPI;
// 纬度转换为弧度
latitude1 = latitude * iPI;
e2 = 2 * f - f * f;
ee = e2 * (1.0 - e2);
NN = a / Math.sqrt(1.0 - e2 * Math.sin(latitude1) * Math.sin(latitude1));
T = Math.tan(latitude1) * Math.tan(latitude1);
C = ee * Math.cos(latitude1) * Math.cos(latitude1);
A = (longitude1 - longitude0) * Math.cos(latitude1);
M = a * ((1 - e2 / 4 - 3 * e2 * e2 / 64 - 5 * e2 * e2 * e2 / 256) * latitude1 - (3 * e2 / 8 + 3 * e2 * e2 / 32 + 45 * e2 * e2 * e2 / 1024) * Math.sin(2 * latitude1) + (15 * e2 * e2 / 256 + 45 * e2 * e2 * e2 / 1024) * Math.sin(4 * latitude1) - (35 * e2 * e2 * e2 / 3072) * Math.sin(6 * latitude1));
xval = NN * (A + (1 - T + C) * A * A * A / 6 + (5 - 18 * T + T * T + 72 * C - 58 * ee) * A * A * A * A * A / 120);
yval = M + NN * Math.tan(latitude1) * (A * A / 2 + (5 - T + 9 * C + 4 * C * C) * A * A * A * A / 24 + (61 - 58 * T + T * T + 600 * C - 330 * ee) * A * A * A * A * A * A / 720);
X0 = 1000000L * (ProjNo + 1) + 500000L;
Y0 = 0;
xval = xval + X0;
yval = yval + Y0;
return new double[] { xval, yval };
}
}
java 弧度 转经纬度_GPS 经纬度转换为 经过旋转后的平面坐标详解相关推荐
- java 墨卡托转经纬度_经纬度,墨卡托等坐标转换
1.经纬度转墨卡托 1.1Java: public static double[] lngLat2Mercator(double lng, double lat) { double[] xy = ne ...
- Java中常见RuntimeException与其他异常表及Exception逻辑关系详解
Java中常见RuntimeException与其他异常表及Exception逻辑关系详解 前言 常见`RuntimeException` 其他错误类型 `Error`类 `Exception`类 E ...
- java中北大学ppt总结+课后习题第四章(小宇特详解)
java中北大学ppt总结+课后习题第四章(小宇特详解) 继承 子类与父类 继承是根据现有类创建新的类的机制,由继承而得到的新类称为子类(subclass)或派生类(derived class),被继 ...
- java调用webservice_笃学私教:Java开发网站架构演变过程-从单体应用到微服务架构详解...
原标题:笃学私教:Java开发网站架构演变过程-从单体应用到微服务架构详解 Java开发网站架构演变过程,到目前为止,大致分为5个阶段,分别为单体架构.集群架构.分布式架构.SOA架构和微服务架构.下 ...
- Java调用SMSLib用单口短信猫发送短信详解
技术园地 当前位置:短信猫网站主页 > 技术园地 > [转载]Java调用SMSLib用单口短信猫发送短信详解 发布时间:2017/02/09 点击量:620 SMSLib是Apache的 ...
- Java 集成高德地图(源码复制即用)实现地址转经纬度,经纬度转地址,计算两点间距离
步骤: (更改方法返回值类型为基本类型可省去这步) copy 封装好的返回对象 导入工具类所需枚举类(避免出现魔法值) 导入工具类(将高德应用Key改为自己申请的key) 导入 controller ...
- java 墨卡托转经纬度,坐标转换:墨卡托 经纬度
坐标转换:墨卡托 & 经纬度 var lonlat={x:114.514864,y:38.04232}; var mercator={x:12747736.346966475,y:458540 ...
- 米转经纬度;经纬度算距离
//米转地理距离度 private static double meterToDegree(double len){double ra = Math.pow((Math.pow(len, -1)) * ...
- java中文字符怎么保证出现正确_JAVA中文字符编码问题详解
JAVA中文字符编码问题详解 JAVA的中文字符乱码问题一直很让人头疼.特别是在WEB应用中.网上的分析文章和解决方案都很多,但总是针对某些特定情况的.很多次遇到乱码问题后,经过极为辛苦的调试和搜索资 ...
最新文章
- VS2013+OpenCV3.1.0配置教程
- 11.14PMP试题每日一题
- 【边缘检测】BDCN:Bi-Directional Cascade Network for Perceptual Edge Detection
- ubuntu下vscode使用cmake编译运行c++配置文件
- scrollbarStyle属性
- mysql 转int_MySQL索引凭什么能让查询效率提高这么多?
- linux屏幕截图之滚动截屏
- C语言实现万年历记事本,简单的日历记事本jQuery插件e-calendar(带样式美化)
- ## Android Studio 开发(四)--蓝牙通信
- 将1自动补位为01_苹果图标自动补位怎么关闭
- echart地图飞线图
- @Value,@ConfigurationProperties,@EnableConfigurationProperties,@PropertySource,@PropertySources
- 安霸入伙,“开放视觉”ADAS 俱乐部再迎新成员
- LintCode 597.具有最大平均数的子树
- 大一寒假培训(七)——队列与优先队列
- 利用Linux的/proc/stat获取指定进程的cpu占有率及内存
- outlook的archive pst丢失后
- AIoT方案|机智云NFC “碰一碰”智慧家庭一碰配网解决方案
- 英语口语 8 级是这么炼成的!你也可以!
- KONICA MINOLTA China | 柯尼卡美能达-SMB扫描问题