QGIS二次开发：数值表示的墨卡托坐标和经纬度表示的WSG84坐标相互转化

本文链接：https://blog.csdn.net/danshiming/article/details/127393692

本博文用到的QGIS 版本说明：

QGIS 3.26.2。

经常在QGIS中看到类似如下很大数值表示的坐标，如：

auto mapCenterPt = QgsPointXY(11472426, 4197692);

而大多数情况下，我们对经纬度表示的坐标比较直观和熟悉，如：

103.05855621730184, 35.248420842233159

前者是基于墨卡托坐标下的坐标表示，后者是基于WSG84坐标系下的坐标表示。如下代码可以将墨卡托坐标下的坐标转为WSG84坐标系下的坐标：

/* 功能：将墨卡托坐标系下的点坐标转为WGS84坐标系下的坐标
   curPoint:墨卡托坐标下的点。输入值
   newPoint：WGS84坐标下的点。输出值
*/
void CMyGIS::convertMercatorToWGS84(const QgsPoint& curPoint, QgsPoint& newPoint)
{
    QgsCoordinateReferenceSystem srcCrs("EPSG:3857"); // 墨卡托坐标
	QgsCoordinateReferenceSystem dstCrs("EPSG:4326"); // WSG84坐标
	const QgsCoordinateTransformContext context;
	QgsCoordinateTransform crsTras(srcCrs, dstCrs, context);

	auto tmpPt = crsTras.transform(curPoint);
	newPoint.setX(tmpPt.x());
	newPoint.setY(tmpPt.y());
}

如下代码可以将WSG84坐标下的坐标转为墨卡托坐标系下的坐标：

/* 功能：将WGS84坐标系下的点坐标转为墨卡托坐标系下的坐标
   curPoint:WGS84坐标下的点。输入值
   newPoint：墨卡托坐标下的点。输出值
*/
void CMyGIS::convertWGS84ToMercator(const QgsPoint& curPoint, QgsPoint& newPoint)
{
	QgsCoordinateReferenceSystem srcCrs("EPSG:4326"); // WSG84坐标
	QgsCoordinateReferenceSystem dstCrs("EPSG:3857"); // 墨卡托坐标
	const QgsCoordinateTransformContext context;
	QgsCoordinateTransform crsTras(srcCrs, dstCrs, context);

	auto tmpPt = crsTras.transform(curPoint);
	newPoint.setX(tmpPt.x());
	newPoint.setY(tmpPt.y());
}

说明：上述代码用到的QGIS版本为QGIS 3.26.2。

关于墨卡托坐标系和WSG84坐标系之间的理论原理讲解，请参见：墨卡托坐标以及墨卡托坐标转经纬度_荆楚闲人的博客-CSDN博客