java代码实现CGCS2000大地坐标系XY值转化为对应经纬度

以下的方法可以实现CGCS2000大地坐标系XY值转化为对应经纬度信息，

注意输入的XY值与给定的值相反，则参数X为已知的Y，参数Y为已知的X。

得出的结果为 [纬度，经度]，不要应用错误。

L0参数为中央子午线的经线值，应用方法前需要确定L0的值，否则得出的值会有很大的偏差。


public static double [] xytolatlon(double X, double Y ,double L0) {double lat ,lon;Y-=500000;double []  result  = new double[2];double iPI = 0.0174532925199433;//pi/180double a = 6378137.0; //长半轴 mdouble b = 6356752.31414; //短半轴 mdouble f = 1/298.257222101;//扁率 a-b/adouble e = 0.0818191910428; //第一偏心率 Math.sqrt(5)double ee = Math.sqrt(a*a-b*b)/b; //第二偏心率double bf = 0; //底点纬度double a0 = 1+(3*e*e/4) + (45*e*e*e*e/64) + (175*e*e*e*e*e*e/256) + (11025*e*e*e*e*e*e*e*e/16384) + (43659*e*e*e*e*e*e*e*e*e*e/65536);double b0 = X/(a*(1-e*e)*a0);double c1 = 3*e*e/8 +3*e*e*e*e/16 + 213*e*e*e*e*e*e/2048 + 255*e*e*e*e*e*e*e*e/4096;double c2 = 21*e*e*e*e/256 + 21*e*e*e*e*e*e/256 + 533*e*e*e*e*e*e*e*e/8192;double c3 = 151*e*e*e*e*e*e*e*e/6144 + 151*e*e*e*e*e*e*e*e/4096;double c4 = 1097*e*e*e*e*e*e*e*e/131072;bf = b0 + c1*Math.sin(2*b0) + c2*Math.sin(4*b0) +c3*Math.sin(6*b0) + c4*Math.sin(8*b0); // bf =b0+c1*sin2b0 + c2*sin4b0 + c3*sin6b0 +c4*sin8b0 +...double tf = Math.tan(bf);double n2 = ee*ee*Math.cos(bf)*Math.cos(bf); //第二偏心率平方成bf余弦平方double c = a*a/b;double v=Math.sqrt(1+ ee*ee*Math.cos(bf)*Math.cos(bf));double mf = c/(v*v*v); //子午圈半径double nf = c/v;//卯酉圈半径//纬度计算lat=bf-(tf/(2*mf)*Y)*(Y/nf) * (1-1/12*(5+3*tf*tf+n2-9*n2*tf*tf)*(Y*Y/(nf*nf))+1/360*(61+90*tf*tf+45*tf*tf*tf*tf)*(Y*Y*Y*Y/(nf*nf*nf*nf)));//经度偏差lon=1/(nf*Math.cos(bf))*Y -(1/(6*nf*nf*nf*Math.cos(bf)))*(1+2*tf*tf +n2)*Y*Y*Y + (1/(120*nf*nf*nf*nf*nf*Math.cos(bf)))*(5+28*tf*tf+24*tf*tf*tf*tf)*Y*Y*Y*Y*Y;result[0] =retain6(lat/iPI);result[1] =retain6(L0+lon/iPI);//System.out.println(result[1]+","+result[0]);return result;}

private static double retain6(double num) {String result = String.format("%.6f", num);return Double.valueOf(result);}

java代码实现CGCS2000大地坐标系XY值转化为对应经纬度相关推荐

java content-length_为什么Content-Length HTTP头字段使用的值不是Java代码中给出的值？...
我有一段 Java代码将字节数组传输到HTTP服务器: HttpURLConnection connection = (HttpURLConnection) url.openConnection(); ...
Hadoop实例之Java代码实现利用MapReduce求π值
需求:假如有一个边长为1的正方形.以正方形的一个端点为圆心,以1为半径,画一个圆弧,于是在正方形内就有了一个直角扇形.在正方形里随机生成若干的点,则有些点是在扇形内,有些点是在扇形外.正方形的面积是1 ...
java代码实现十进制到二进制的转化
前言十进制与二进制之间的转化是小学就学过的内容,但是我相信还是有一部分同学忘了如何"反向求余",代码用python和matlab等编程语言特别好实现,对于刚接触java语言的同学 ...
为什么数据库可以查询到值，java代码里面返回为null值
为什么我的查询用#不行要改成$ 一直查询不到我需要的值改成$可以查询 #{} 这种是防止注入输入进来的数据都会带有' '引号 $这种是不防止注入你输入进来的数据怎么进来的就怎么放进去不会自带引 ...
js中变量和jsp中java代码中变量互相访问解决方案
1.js变量获取jsp页面中java代码的变量值. 方法:var JS变量名 = <%=JAVA变量名 %> 2.java代码获取js变量的值. 说明:在JSP中:Java部分是在 ...
js变量和java变量相等,js中变量和jsp中java代码中变量互相访问解决方案
1.js变量获取jsp页面中java代码的变量值. 方法:var JS变量名 = 我们常常会将js文件和jsp文件分开写,在js文件中,上面的方法似乎不管用了. 也可以通过变通的方法来解决: a.js ...
【运筹优化】AFSA人工鱼群算法求解无约束多元函数最值(Java代码实现)
文章目录前言优化目标求解结果迭代结果可视化算法流程 Java代码可视化代码前言本文以求解二元函数最小值为例,如果需要求解多元函数,只需要修改以下变量即可: varNum:变量维度数 u ...
【智能算法】FA萤火虫算法求解无约束多元函数最值(Java代码实现)
文章目录前言优化目标求解结果搜索过程可视化 Java代码可视化代码前言本文以求解二元函数最小值为例,如果需要求解多元函数,只需要修改以下变量即可: varNum:变量维度数 ub和lb: ...
【运筹优化】SMO蜘蛛猴优化算法求解无约束多元函数最值(Java代码实现)
文章目录前言优化目标优化结果迭代过程可视化 Java代码可视化代码优化流程(图太大了,所以放最后...) 前言本文以求解二元函数最小值为例,如果需要求解多元函数,只需要修改以下变量即可: ...
【运筹优化】CSO蟑螂算法求解无约束多元函数最值(Java代码实现)
文章目录一.前言二.优化目标三.求解结果四.迭代结果可视化五.算法流程图六.Java代码实现 6.1 算法部分代码 6.2 可视化部分代码一.前言本文以求解二元函数最小值为例,如果需要 ...

java代码实现CGCS2000大地坐标系XY值转化为对应经纬度

java代码实现CGCS2000大地坐标系XY值转化为对应经纬度相关推荐

最新文章

热门文章