C语言GPS经纬度距离高德地图距离计算.pdfVIP

C语⾔GPS经纬度距离⾼德地图距离计算 #include math.h #include stdio.h double gcj_ lng, gcj_ lat ; double pi = 3. 4 5926; double a = 6378245.0; double es = 0 62296594323 ; double transformlat(double lng, double lat) { double ret ; ret = - 00.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0. * lng * lat + 0.2 * sqrt(fabs(lng)) ; ret += (20.0 * sin(6.0 * lng * pi) + 20.0 * sin(2.0 * lng * pi)) * 2.0 / 3.0; ret += (20.0 * sin(lat * pi) + 40.0 * sin(lat / 3.0 * pi)) * 2.0 / 3.0; ret += ( 60.0 * sin(lat / 2.0 * pi) + 320 * sin(lat * pi / 30.0)) * 2.0 / 3.0; return (ret) ; } double transformlng(double lng, double lat) { double ret ; ret = 300.0 + lng + 2.0 * lat + 0. * lng * lng + 0. * lng * lat + 0. * sqrt(fabs(lng)) ; ret += (20.0 * sin(6.0 * lng * pi) + 20.0 * sin(2.0 * lng * pi)) * 2.0 / 3.0; ret += (20.0 * sin(lng * pi) + 40.0 * sin(lng / 3.0 * pi)) * 2.0 / 3.0; ret += ( 50.0 * sin(lng / 2.0 * pi) + 300.0 * sin(lng / 30.0 * pi)) * 2.0 / 3.0; return (ret) ; } void WGS84_to_GCJ02(double lng, double lat) { / * * 实现WGS84坐标系向GCJ02坐标系的转换 * :param lng : WGS84坐标系下的经度 * :param lat: WGS84坐标系下的纬度 * :return: 转换后的GCJ02下经纬度 */ double dlat, dlng, radlat, magic, sqrtmagic; dlat = transformlat(lng - 05.0, lat - 35.0) ; dlng = transformlng(lng - 05.0, lat - 35.0) ; radlat = lat / 80.0 * pi; magic = sin(radlat) ; magic = - es * magic * magic; sqrtmagic = sqrt(magic) ; dlat = (dlat * 80.0) / ((a * ( - es)) / (magic * sqrtmagic) * pi) ; dlng = (dlng * 80.0) / (a / sqrtmagic * cos(radlat) * pi) ; gcj_ lng = lat + dlat ; gcj_ lat = lng + dlng ; } / * p oint_ a ,p oint_ b是经纬度,格式为[lng,lat] */ double st_distance(double point_alon, double point_alat, double point_ blon, double point_ blat) { double radia