如果我有标准NMEA格式的纬度或经度读数有一个简单的方法/公式将该读数转换为米,然后我可以用Java(J9)实现?
编辑:好吧似乎我想做的事情不容易,但我真正想做的是:
假设我有一个lat和一个长途的点和一个lat和很长的用户有一个简单的方法来比较它们来决定何时告诉用户他们在一个合理的距离点的方式点?我觉得合理是主题,但是这很容易做到还是仍然过于数学?
这是一个javascript函数:
function measure(lat1, lon1, lat2, lon2){ // generally used geo measurement function
var R = 6378.137; // Radius of earth in KM
var dLat = lat2 * Math.PI / 180 - lat1 * Math.PI / 180;
var dLon = lon2 * Math.PI / 180 - lon1 * Math.PI / 180;
var a = Math.sin(dLat/2) * Math.sin(dLat/2) +
Math.cos(lat1 * Math.PI / 180) * Math.cos(lat2 * Math.PI / 180) *
Math.sin(dLon/2) * Math.sin(dLon/2);
var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));
var d = R * c;
return d * 1000; // meters
}
说明:https://en.wikipedia.org/wiki/Haversine_formula
根据其经度和纬度,半正公式确定球体上两点之间的大圆距离.
鉴于你正在寻找一个简单的公式,这可能是最简单的方法,假设地球是一个40075公里的周长球.
纬度1°的长度(以米为单位)=始终为111.32 km
以1°经度为单位的长度= 40075 km*cos(纬度)/ 360
为了近似两个坐标之间的短距离,我使用了http://en.wikipedia.org/wiki/Lat-lon中的公式 :
m_per_deg_lat = 111132.954 - 559.822 * cos( 2 * latMid ) + 1.175 * cos( 4 * latMid); m_per_deg_lon = 111132.954 * cos ( latMid );
.
在下面的代码中,我留下了原始数字,以显示它们与维基百科公式的关系.
double latMid, m_per_deg_lat, m_per_deg_lon, deltaLat, deltaLon,dist_m; latMid = (Lat1+Lat2 )/2.0; // or just use Lat1 for slightly less accurate estimate m_per_deg_lat = 111132.954 - 559.822 * cos( 2.0 * latMid ) + 1.175 * cos( 4.0 * latMid); m_per_deg_lon = (3.14159265359/180 ) * 6367449 * cos ( latMid ); deltaLat = fabs(Lat1 - Lat2); deltaLon = fabs(Lon1 - Lon2); dist_m = sqrt ( pow( deltaLat * m_per_deg_lat,2) + pow( deltaLon * m_per_deg_lon , 2) );
维基百科条目指出,纵向100km的距离计算在0.6m范围内,纬度100km处的距离计算在1cm范围内,但我没有证实这一点,因为我的使用精度很高.
纬度和经度指定点,而不是距离,所以你的问题有点荒谬.如果你问的是两个(纬度,经度)点之间的最短距离,请参阅这篇关于大圆距离的维基百科文章.
地球是一个令人讨厌的不规则表面,所以没有简单的公式可以做到这一点.你必须使用地球的近似模型,并将坐标投射到它上面.我通常看到用于此的模型是WGS 84.这就是GPS设备通常用来解决完全相同的问题.
NOAA有一些软件可以在他们的网站上下载以帮助解决这个问题.
有许多工具可以使这很容易.有关涉及的内容的详细信息,请参阅monjardin的答案.
但是,这样做并不一定困难.听起来你正在使用Java,所以我建议你研究类似GDAL的东西.它为它们的例程提供了java包装器,它们具有从Lat/Lon(地理坐标)转换为UTM(投影坐标系)或其他一些合理的地图投影所需的所有工具.
UTM很好,因为它是米,所以很容易使用.但是,您需要获得适当的UTM区域才能做得很好.通过谷歌搜索可以获得一些简单的代码,以便为纬度/经度对找到合适的区域.
这是bh函数的R版本,以防万一:
measure <- function(lon1,lat1,lon2,lat2) {
R <- 6378.137 # radius of earth in Km
dLat <- (lat2-lat1)*pi/180
dLon <- (lon2-lon1)*pi/180
a <- sin((dLat/2))^2 + cos(lat1*pi/180)*cos(lat2*pi/180)*(sin(dLon/2))^2
c <- 2 * atan2(sqrt(a), sqrt(1-a))
d <- R * c
return (d * 1000) # distance in meters
}