本文主要是介绍java代码实现CGCS2000大地坐标系XY值转化为对应经纬度,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
以下的方法可以实现CGCS2000大地坐标系XY值转化为对应经纬度信息,
L0参数为中央子午线的经线值,应用方法前需要确定L0的值,否则得出的值会有很大的偏差。
L0我用的是117
/**** @param X 点位坐标X* @param Y 点位坐标Y* @param L0 L0参数为中央子午线的经线值* @return*/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);}
关于中央子午线延申
全球分为二十四个时区,以能够被15整除的经度作为该区域的中央子午线,每一时区占经度15度。在该时区中央子午线以东的地区,时间要加,以西的地区,时间要减,一度4分钟。
中国共分五个时区:
(1)中原时区:以东经120度为中央子午线。
(2)陇蜀时区:以东经105度为中央子午线。
(3)新藏时区:以东经90度为中央子午线。
(4)昆仑时区:以东经75(82.5)度为中央子午线。
(5)长白时区:以东经135(127.5)度为中央子午线。
一个时区的“标准时”,只是一个大地区的统一时间,大家共同遵守的“人工”时间而已,并不是该时区内每个地点的“本地时间(lmt)”——真正的经度时。要用出生地的经度与出生大地区的标准时来加减,全球任何地点都用这个原则。
例如:中原时区包括内蒙古、辽宁、河北、山西、山东、河南、安徽、江苏、湖北、湖南、江西、浙江、福建、广东、海南、香港、澳门、台湾。这个大地区当时钟敲定正午12点时,只有位于东经120度线上的地点才是12点,其它的地方是少于或多于12点。如香港位于东经114度10分,比东经120度偏西5度50分,其真正经度时是11时36分40秒。
这篇关于java代码实现CGCS2000大地坐标系XY值转化为对应经纬度的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!