通过经纬度转换为XY值
输入的参数为经度,纬度
得到的结果是[X,Y]
代码内部longitude0代表中央子午线的值,需要根据实际内容配置,也可以自己提到参数中
public static double[] GaussToBLToGauss(double longitude, double latitude)//参数 经度,纬度
{
double[] output = new double[2];
double longitude1,latitude1, longitude0, X0,Y0, xval,yval;
//NN曲率半径,测量学里面用N表示
//M为子午线弧长,测量学里用大X表示
//fai为底点纬度,由子午弧长反算公式得到,测量学里用Bf表示
//R为底点所对的曲率半径,测量学里用Nf表示
double a,f, e2,ee, NN, T,C,A, M, iPI;
iPI = 0.0174532925199433; //3.1415926535898/180.0;
a=6378137.0; f=1/298.257222101; //CGCS2000坐标系参数
//a=6378137.0; f=1/298.2572236; //wgs84坐标系参数
longitude0 = 117.0;//中央子午线 根据实际进行配置
longitude0 = longitude0 * iPI ;//中央子午线转换为弧度
longitude1 = longitude * iPI ; //经度转换为弧度
latitude1 =