地理坐标平面化
首先是地理坐标的平面化转化。在一个小范围内(例如是方圆几公里内),我们可以假设地面是平的,而不是弯的。如果经纬度都用弧度表示,那么1纬度对应的长度是:
1lat_len=R*lat
其中R是地球半径。
而相同经度间的距离会随着纬度的增加而减少,在lat这一纬度下,1经度对应的长度是:
1lon_len=Rloncos(lat)
那么,(lon,lat)这个坐标平面化后的坐标就是:(Rloncos(lat),R*lat)
加上弧度,(lon,lat)取值示例(116,30)
对应的x,y坐标分别为
double x = R * lon * Math.cos(lat/180 *Math.PI)/180 *Math.PI;
double y = R * lat/180 * Math.PI;
其中R可取值 R = 6378137; // 赤道半径 6378137
WGS84坐标转平面坐标代码
public class Gx {
static double R = 6378137; // 赤道半径 6378137
static double R1 = 6356752;
public static GetX x(double lon, double lat) {
double x = R * lon * Math.cos(lat/180 *Math.PI)/180 *Math.PI;
double y = R * lat/180 * Math.PI;
GetX xy = new GetX(x,y);
System.out.println(x+ "," + y);
return xy ;
}
}
平面坐标转为 WGS84坐标代码
public class Gp {
static double R = 6378137;
public static GetX p(double x, double y) {
double lat = y/R *180 /Math.PI;
double lon = x/(R * Math.cos(lat/180 *Math.PI))*180 /Math.PI;
GetX xy = new GetX(lat,lon);
System.out.println(x+ "," + y);
return xy ;
}
}