重庆分公司,新征程启航
为企业提供网站建设、域名注册、服务器等服务
本文介绍了java 地心坐标系(ECEF)和WGS-84坐标系(WGS84)互转的实现,分享给大家,具体如下:
成都创新互联公司服务项目包括尚志网站建设、尚志网站制作、尚志网页制作以及尚志网络营销策划等。多年来,我们专注于互联网行业,利用自身积累的技术优势、行业经验、深度合作伙伴关系等,向广大中小型企业、政府机构等提供互联网行业的解决方案,尚志网站推广取得了明显的社会效益与经济效益。目前,我们服务的客户以成都为中心已经辐射到尚志省份的部分城市,未来相信会继续扩大服务区域并继续获得客户的支持与信任!
public static String WGS84toECEF(double latitude, double longitude, double height) { double X; double Y; double Z; double a = 6378137; double b = 6356752.314245; double E = (a * a - b * b) / (a * a); double COSLAT = Math.cos(latitude * Math.PI / 180); double SINLAT = Math.sin(latitude * Math.PI / 180); double COSLONG = Math.cos(longitude * Math.PI / 180); double SINLONG = Math.sin(longitude * Math.PI / 180); double N = a / (Math.sqrt(1 - E * SINLAT * SINLAT)); double NH = N + height; X = NH * COSLAT * COSLONG; Y = NH * COSLAT * SINLONG; Z = (b * b * N / (a * a) + height) * SINLAT; return X + "," + Y + "," + Z; } public static String ECEFtoWGS84(double x, double y, double z) { double a, b, c, d; double Longitude;//经度 double Latitude;//纬度 double Altitude;//海拔高度 double p, q; double N; a = 6378137.0; b = 6356752.31424518; c = Math.sqrt(((a * a) - (b * b)) / (a * a)); d = Math.sqrt(((a * a) - (b * b)) / (b * b)); p = Math.sqrt((x * x) + (y * y)); q = Math.atan2((z * a), (p * b)); Longitude = Math.atan2(y, x); Latitude = Math.atan2((z + (d * d) * b * Math.pow(Math.sin(q), 3)), (p - (c * c) * a * Math.pow(Math.cos(q), 3))); N = a / Math.sqrt(1 - ((c * c) * Math.pow(Math.sin(Latitude), 2))); Altitude = (p / Math.cos(Latitude)) - N; Longitude = Longitude * 180.0 / Math.PI; Latitude = Latitude * 180.0 / Math.PI; return Longitude + "," + Latitude + "," + Altitude; }
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持创新互联。