WGS84_UTM.zip

#include "coorconv.h" #include <iostream> #include <string> using namespace std; CCoorConv::CCoorConv() : m_nUTMZone(-10000), m_bNorth(false), m_pj_utm(nullptr) { const char *wgs84 = "+proj=longlat +datum=WGS84 +no_defs ";//GPS所用坐标系,EPSG:4326 m_pj_wgs84 = pj_init_plus(wgs84); } CCoorConv::~CCoorConv() { if (m_pj_utm != nullptr) { pj_free(m_pj_utm); m_pj_utm = nullptr; } if (m_pj_wgs84 != nullptr) { pj_free(m_pj_wgs84); m_pj_wgs84 = nullptr; } } bool CCoorConv::LatLonToUTMXY(const WGS84Corr wgs, UTMCorr& utm) { int nUTMZone = (int)((wgs.log + 186.0) / 6.0); bool bNorth = wgs.lat > 0 ? true : false; if (m_nUTMZone != nUTMZone || m_bNorth != bNorth) { //"+proj=utm +zone=50 +datum=WGS84 +units=m +no_defs"; // "+proj=utm +zone=48 +south +datum=WGS84 +units=m +no_defs" std::string qstrUTM = "+proj=utm +zone=" + std::to_string(nUTMZone); if (!bNorth) { qstrUTM += " +south "; } qstrUTM += " +datum=WGS84 +units=m +no_defs"; if (m_pj_utm != nullptr) { pj_free(m_pj_utm); m_pj_utm = NULL; } string strUTM = qstrUTM; const char *pUTM = strUTM.c_str(); m_pj_utm = pj_init_plus(pUTM); m_nUTMZone = nUTMZone; m_bNorth = bNorth; } utm.x = wgs.log * DEG_TO_RAD; utm.y = wgs.lat * DEG_TO_RAD; pj_transform(m_pj_wgs84, m_pj_utm, 1, 1, &utm.x, &utm.y, nullptr); return true; } bool CCoorConv::UTMXYToLatLon(const UTMCorr utm, WGS84Corr& wgs) { if (m_pj_utm == nullptr) { return false; } wgs.log = utm.x; wgs.lat = utm.y; pj_transform(m_pj_utm, m_pj_wgs84, 1, 1, &wgs.log, &wgs.lat, nullptr); wgs.log = wgs.log * RAD_TO_DEG; wgs.lat = wgs.lat * RAD_TO_DEG; return true; }