當前位置：首頁 > 编程资源 > 编程问答 >内容正文

编程问答

城建坐标与经纬度转换工具

發布時間：2023/12/9 编程问答 54 豆豆

生活随笔收集整理的這篇文章主要介紹了城建坐标与经纬度转换工具小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

部分底層核心源碼

需要引入的Jar包為: CoodinateCover.jar可在csdn自行搜索

CoodinateCover.java

package PkgCover;import java.io.PrintStream;public class CoodinateCover {private GaussProjection _GP;private double _X;private double _Y;private double _B;private double _L;private String _ErrMsg;private enumLocalCoordinateSys _lcs;public static void main(String[] args){CoodinateCover CC = new CoodinateCover(3, enumEllipseSys.WGS84, enumLocalCoordinateSys.BeiJingCJ);if (!CC.FunCover(39.976732200000001D, 116.41657351000001D).booleanValue()) {System.out.println("出現錯誤！" + CC.get_ErrMsg());}System.out.println("經緯度坐標（40.2321880683303, 116.523973513395）轉換為城建坐標（" + CC.get_X() + ",Y坐標：" + CC.get_Y() + "）");if (!CC.FunAntiCover(CC.get_X(), CC.get_Y()).booleanValue()) {System.out.println("出現錯誤！" + CC.get_ErrMsg());}System.out.println("反轉后為：緯度為：" + CC.get_B() + "，經度為：" + CC.get_L());}public CoodinateCover(int zoneWide, enumEllipseSys ellipse, enumLocalCoordinateSys lcs){this._lcs = lcs;this._GP = new GaussProjection(ellipse, zoneWide);}public Boolean FunCover(double B, double L){if (!this._GP.Projection(B, L, 0.0D).booleanValue()){this._ErrMsg = ("高斯投影失敗：" + this._GP.get_ErrMsg());return Boolean.valueOf(false);}FPData oFPData;FPData oFPData;switch (this._lcs){case BeiJingCJ: LocalCoordinateSys LCS = new BeiJingCJ();oFPData = LCS.getFPData(true, this._GP.get_Maridian());break;default: oFPData = null;}if (oFPData == null){this._ErrMsg = "沒有實例化FPData對象";return Boolean.valueOf(false);}this._X = (oFPData.get_X() + oFPData.get_K() * this._GP.get_X() * Math.cos(oFPData.get_T()) - oFPData.get_K() * this._GP.get_Y() * Math.sin(oFPData.get_T()));this._Y = (oFPData.get_Y() + oFPData.get_K() * this._GP.get_X() * Math.sin(oFPData.get_T()) + oFPData.get_K() * this._GP.get_Y() * Math.cos(oFPData.get_T()));return Boolean.valueOf(true);}public Boolean FunAntiCover(double x, double y){Boolean bRtn = Boolean.valueOf(true);switch (this._lcs){case BeiJingCJ: bRtn = Boolean.valueOf(BeiJintCJToWGS84(x, y));break;default: bRtn = Boolean.valueOf(false);}return bRtn;}private boolean BeiJintCJToWGS84(double x, double y){boolean bRtn = true;LocalCoordinateSys LCS = new BeiJingCJ();FPData oFPData = LCS.getFPData(false, 117);this._X = (oFPData.get_X() + oFPData.get_K() * x * Math.cos(oFPData.get_T()) - oFPData.get_K() * y * Math.sin(oFPData.get_T()));this._Y = (oFPData.get_Y() + oFPData.get_K() * x * Math.sin(oFPData.get_T()) + oFPData.get_K() * y * Math.cos(oFPData.get_T()));if (!this._GP.AntiProjection(this._X, this._Y, 117).booleanValue()){this._ErrMsg = ("高斯投影失敗：" + this._GP.get_ErrMsg());bRtn = false;return bRtn;}if ((this._GP.get_L() < 1.5D) && (this._GP.get_L() > -1.5D)){this._L = (117.0D + this._GP.get_L());bRtn = true;}else if (this._GP.get_L() > 1.5D){this._B = this._GP.get_Lat();this._L = (117.0D + this._GP.get_L());this._ErrMsg = "超出投影范圍，可能不準確";bRtn = false;}else{FPData oFPData2 = LCS.getFPData(false, 114);this._X = (oFPData2.get_X() + oFPData2.get_K() * x * Math.cos(oFPData2.get_T()) - oFPData2.get_K() * y * Math.sin(oFPData2.get_T()));this._Y = (oFPData2.get_Y() + oFPData2.get_K() * x * Math.sin(oFPData2.get_T()) + oFPData2.get_K() * y * Math.cos(oFPData2.get_T()));if (!this._GP.AntiProjection(this._X, this._Y, 114).booleanValue()){this._ErrMsg = ("高斯投影失敗：" + this._GP.get_ErrMsg());bRtn = false;}this._L = (114.0D + this._GP.get_L());bRtn = true;}this._B = this._GP.get_Lat();return bRtn;}public double get_X(){return this._X;}public double get_Y(){return this._Y;}public double get_B(){return this._B;}public double get_L(){return this._L;}public String get_ErrMsg(){return this._ErrMsg;} }

GaussProjection.java

package PkgCover;public class GaussProjection {double f = 0.0D;double e2 = 0.0D;double e1 = 0.0D;double b;double a = 0.0D;private String _ErrMsg;private int _ZoneWide;private double _l;private int _Maridian;private double _Lng;private double _Lat;private double _Y;private double _X;private double _A8;private double _A6;private double _A4;private double _A2;private double _A0;private double _M8;private double _M6;private double _M4;private double _M2;private double _M0;public GaussProjection(enumEllipseSys ellipse, int zoneWide){this._ZoneWide = zoneWide;switch (ellipse){case WGS84: this.a = 6378245.0D;this.b = 6356863.0187730473D;this.e1 = (Math.sqrt(Math.pow(this.a, 2.0D) - Math.pow(this.b, 2.0D)) / this.a);this.e2 = (Math.sqrt(Math.pow(this.a, 2.0D) - Math.pow(this.b, 2.0D)) / this.b);break;case BeiJing54: this.a = 6378137.0D;this.b = 6356752.3141999999D;this.e1 = (Math.sqrt(Math.pow(this.a, 2.0D) - Math.pow(this.b, 2.0D)) / this.a);this.e2 = (Math.sqrt(Math.pow(this.a, 2.0D) - Math.pow(this.b, 2.0D)) / this.b);break;case XiAn80: this.a = 6378140.0D;this.b = 6356755.2881575283D;this.e1 = (Math.sqrt(Math.pow(this.a, 2.0D) - Math.pow(this.b, 2.0D)) / this.a);this.e2 = (Math.sqrt(Math.pow(this.a, 2.0D) - Math.pow(this.b, 2.0D)) / this.b);break;}this.f = ((this.a - this.b) / this.a);CalMidPara();}public Boolean Projection(double B, double L, double H){Boolean bRtn = Boolean.valueOf(true);try{int n = 0;if (this._ZoneWide == 6){n = (int)L / 6;if (L % 6.0D > 0.0D) {n++;}this._Maridian = (6 * n - 3);this._l = (L - this._Maridian);}if (3 == this._ZoneWide){n = (int)(L - 1.5D) / 3;if ((L - 1.5D) % 3.0D > 0.0D) {n++;}this._Maridian = (3 * n);this._l = (L - this._Maridian);}B = B * 3.141592653589793D / 180.0D;L = L * 3.141592653589793D / 180.0D;this._l = (this._l * 3.141592653589793D / 180.0D);double X = this._A0 * B - this._A2 * Math.sin(2.0D * B) / 2.0D + this._A4 * Math.sin(4.0D * B) / 4.0D - this._A6 * Math.sin(6.0D * B) / 6.0D + this._A8 * Math.sin(8.0D * B) / 8.0D;double T = Math.tan(B);double ita = this.e2 * Math.cos(B);double SinB = Math.sin(B);double CosB = Math.cos(B);double N = this.a * Math.pow(1.0D - Math.pow(this.e1 * SinB, 2.0D), -0.5D);this._X = (X + N * T * Math.pow(CosB, 2.0D) * Math.pow(this._l, 2.0D) / 2.0D + N * T * (5.0D - Math.pow(T, 2.0D) + 9.0D * Math.pow(ita, 2.0D) + 4.0D * Math.pow(ita, 2.0D)) * Math.pow(CosB * this._l, 4.0D) / 24.0D + N * T * (61.0D - 58.0D * Math.pow(T, 2.0D) + Math.pow(T, 4.0D)) * Math.pow(CosB * this._l, 6.0D) / 720.0D);this._Y = (N * CosB * this._l + N * (1.0D - Math.pow(T, 2.0D) + Math.pow(ita, 2.0D)) * Math.pow(CosB * this._l, 3.0D) / 6.0D + N * (5.0D - 18.0D * Math.pow(T, 2.0D) + Math.pow(T, 4.0D) + 14.0D * Math.pow(ita, 2.0D) - 58.0D * Math.pow(ita * T, 2.0D)) * Math.pow(CosB * this._l, 5.0D) / 120.0D);}catch (Exception ex){this._ErrMsg = ex.getMessage();bRtn = Boolean.valueOf(false);}return bRtn;}public Boolean AntiProjection(double x, double y, int maridian){this._Maridian = maridian;boolean bRtn = true;try{double Bf0 = x / this._A0;double Bf = GetBf(Bf0, x);double Nf = this.a * Math.pow(1.0D - Math.pow(this.e1 * Math.sin(Bf), 2.0D), -0.5D);double Z = y / (Nf * Math.cos(Bf));double Mf = this.a * (1.0D - this.e1 * this.e1) * Math.pow(1.0D - Math.pow(this.e1 * Math.sin(Bf), 2.0D), -1.5D);double Tf = Math.tan(Bf);double ITAf2 = Math.pow(this.e2 * Math.cos(Bf), 2.0D);double B2 = Tf * Math.pow(y, 2.0D) / (2.0D * Mf * Nf);double B4 = Tf * Math.pow(y, 4.0D) * (5.0D + 3.0D * Math.pow(Tf, 2.0D) + ITAf2 - 9.0D * Math.pow(Tf, 2.0D) * ITAf2) / (24.0D * Mf * Math.pow(Nf, 3.0D));double B6 = Tf * Math.pow(y, 6.0D) * (61.0D + 90.0D * Math.pow(Tf, 2.0D) + 45.0D * Math.pow(Tf, 4.0D)) / (720.0D * Mf * Math.pow(Nf, 5.0D));double B3 = (1.0D + 2.0D * Math.pow(Tf, 2.0D) + ITAf2) * Math.pow(y, 3.0D) / (6.0D * Math.pow(Nf, 3.0D) * Math.cos(Bf));double B5 = (5.0D + 28.0D * Math.pow(Tf, 2.0D) + 24.0D * Math.pow(Tf, 4.0D) + 6.0D * ITAf2 + 8.0D * ITAf2) * Math.pow(y, 5.0D) / (120.0D * Math.pow(Nf, 5.0D) * Math.cos(Bf));this._Lat = ((Bf - B2 + B4 - B6) * 180.0D / 3.141592653589793D);this._l = ((Z - B3 + B5) * 180.0D / 3.141592653589793D);}catch (Exception ex){this._ErrMsg = (this._ErrMsg + "高斯投影反算出錯！" + ex.getMessage());bRtn = false;}return Boolean.valueOf(bRtn);}private double GetBf(double bf, double x){double M = -this._A2 * Math.sin(2.0D * bf) / 2.0D + this._A4 * Math.sin(4.0D * bf) / 4.0D - this._A6 * Math.sin(6.0D * bf) / 6.0D + this._A8 * Math.sin(8.0D * this.f) / 8.0D;double Bf = (x - M) / this._A0;if (Math.abs(bf - Bf) > 3.141592653589793D * Math.pow(10.0D, -8.0D) / 648.0D) {Bf = GetBf(Bf, x);}return Bf;}private void CalMidPara(){this._M0 = (this.a * (1.0D - Math.pow(this.e1, 2.0D)));this._M2 = (3.0D * Math.pow(this.e1, 2.0D) * this._M0 / 2.0D);this._M4 = (5.0D * Math.pow(this.e1, 2.0D) * this._M2 / 4.0D);this._M6 = (7.0D * Math.pow(this.e1, 2.0D) * this._M4 / 6.0D);this._M8 = (9.0D * Math.pow(this.e1, 2.0D) * this._M6 / 8.0D);this._A0 = (this._M0 + this._M2 / 2.0D + 3.0D * this._M4 / 8.0D + 5.0D * this._M6 / 16.0D + 35.0D * this._M8 / 128.0D);this._A2 = (this._M2 / 2.0D + this._M4 / 2.0D + 15.0D * this._M6 / 32.0D + 7.0D * this._M8 / 16.0D);this._A4 = (this._M4 / 8.0D + this._M6 * 3.0D / 16.0D + this._M8 * 7.0D / 32.0D);this._A6 = (this._M6 / 32.0D + this._M8 / 16.0D);this._A8 = (this._M8 / 128.0D);}public double get_X(){return this._X;}public double get_Y(){return this._Y;}public double get_Lat(){return this._Lat;}public double get_Lng(){return this._Lng;}public String get_ErrMsg(){return this._ErrMsg;}public double get_L(){return this._l;}public int get_Maridian(){return this._Maridian;} }

案例

首先經緯度轉換Util，把相應的城建坐標"499549.528", "317657.996"轉換為經緯度

import java.io.BufferedReader; import java.io.IOException; import java.io.InputStreamReader; import java.net.URL; import java.net.URLConnection;import PkgCover.CoodinateCover; import PkgCover.Coordinate; import PkgCover.CoordinateSys; import PkgCover.enumEllipseSys; import PkgCover.enumLocalCoordinateSys;public class Utils {public static Double[] getCoordinate(String coordx, String coordy) {Double[] returnCoordinate = new Double[2];if (coordx == null || coordy == null) {return null;}CoodinateCover cc = new CoodinateCover(3, enumEllipseSys.WGS84, enumLocalCoordinateSys.BeiJingCJ);if (!cc.FunCover(Double.valueOf(coordy).doubleValue(), Double.valueOf(coordx).doubleValue())) {try {throw new RuntimeException(cc.get_ErrMsg());} catch (Exception e) {System.out.println("捕獲轉換異常：" + e);}}returnCoordinate[0] = cc.get_Y();returnCoordinate[1] = cc.get_X();return returnCoordinate;}public static Double[] getAinCoordinate(String coordx, String coordy) {Double[] returnCoordinate = new Double[2];if (coordx == null || coordy == null) {return null;}CoodinateCover cc = new CoodinateCover(3, enumEllipseSys.WGS84, enumLocalCoordinateSys.BeiJingCJ);if (!cc.FunAntiCover(Double.valueOf(coordy).doubleValue(), Double.valueOf(coordx).doubleValue())) {try {throw new RuntimeException(cc.get_ErrMsg());} catch (Exception e) {System.out.println("捕獲轉換異常：" + e);}}if (cc.get_L() != 0.0D && cc.get_B() != 0.0D) {Coordinate _Coordinate = new Coordinate(cc.get_B(), cc.get_L(), CoordinateSys.WGS84);if (_Coordinate.WGS84ToGCJ02()) {returnCoordinate[0] = _Coordinate.get_GCJLat();returnCoordinate[1] = _Coordinate.get_GCJLng();} else {returnCoordinate[0] = cc.get_L();returnCoordinate[1] = cc.get_B();}} else {returnCoordinate[0] = cc.get_L();returnCoordinate[1] = cc.get_B();}return returnCoordinate;}public static void main(String[] args) {Double[] dd = Utils.getAinCoordinate("499549.528", "317657.996");System.out.println(dd[0] + "," + dd[1]);} }

總結

以上是生活随笔為你收集整理的城建坐标与经纬度转换工具的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。

上一篇： android 代码布局设置wrap_c
下一篇： python3函数调用时间_Python