下載了好幾個(gè)源程序都不能用,在他們的基礎(chǔ)上經(jīng)過我與公式的核對(duì).現(xiàn)在將代碼改成了java

以下是完整代碼:絕對(duì)能用,我現(xiàn)在的項(xiàng)目中就用的這個(gè).

package tms.base.systemlib;

@SuppressWarnings("unchecked")

public class GaussXYDeal {

//? 由高斯投影坐標(biāo)反算成經(jīng)緯度

public static double[] GaussToBL(double X, double Y)//, double *longitude, double *latitude)

{

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年北京坐標(biāo)系參數(shù)

a=6378140.0; f=1/298.257; //80年西安坐標(biāo)系參數(shù)

ZoneWide = 6; 6度帶寬

ProjNo = (int)(X/1000000L) ; //查找?guī)?hào)

longitude0 = (ProjNo-1) * ZoneWide + ZoneWide / 2;

longitude0 = longitude0 * iPI ; //中央經(jīng)線

X0 = ProjNo*1000000L+500000L;

Y0 = 0;

xval = X-X0; yval = Y-Y0; //帶內(nèi)大地坐標(biāo)

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;

//計(jì)算經(jīng)度(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);

//轉(zhuǎn)換為度 DD

output[0] = longitude1 / iPI;

output[1] = latitude1 / iPI;

return output;

//*longitude = longitude1 / iPI;

//*latitude = latitude1 / iPI;

}

? 由經(jīng)緯度反算成高斯投影坐標(biāo)

public void GaussToBLToGauss(double longitude, double latitude)

{

int ProjNo=0; int ZoneWide; 帶寬

double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval;

double a,f, e2,ee, NN, T,C,A, M, iPI;

iPI = 0.0174532925199433; 3.1415926535898/180.0;

ZoneWide = 6; 6度帶寬

a=6378245.0; f=1.0/298.3; //54年北京坐標(biāo)系參數(shù)

a=6378140.0; f=1/298.257; //80年西安坐標(biāo)系參數(shù)

ProjNo = (int)(longitude / ZoneWide) ;

longitude0 = ProjNo * ZoneWide + ZoneWide / 2;

longitude0 = longitude0 * iPI ;

latitude0 = 0;

System.out.println(latitude0);

longitude1 = longitude * iPI ; //經(jīng)度轉(zhuǎn)換為弧度

latitude1 = latitude * iPI ; //緯度轉(zhuǎn)換為弧度

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;

//*X = xval;

//*Y = yval;

System.out.println("x："+xval);

System.out.println("y："+yval);

}

}

“七”樂無窮，盡在新浪新版博客，快來體驗(yàn)啊~~~請(qǐng)點(diǎn)擊進(jìn)入~

總結(jié)

以上是生活随笔為你收集整理的java 高斯投影_贡献高斯投影坐标与经纬度转换的java源码的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。