MATLAB CGCS2000高斯投影坐标转经纬度坐标

坐标系转换这边需要用到mapping toolbox

首先根据原始（x,y）坐标对应的投影坐标系查询EPSG编号

例如这边CGCS2000 / 3-degree Gauss-Kruger CM 123E的编号就是4450

对应的编号可以https://blog.csdn.net/qq_41441896/article/details/104525296在这篇博客中查到

如果是其他投影坐标系也可以在https://spatialreference.org/ref/epsg/ Spatial Reference List中查询任意的投影坐标系编号

随后定义好投影坐标系后即可用projinv函数转换坐标

% 空间坐标系EPSG编号 code = 4550;% CGCS2000 / 3-degree Gauss-Kruger CM 123E % --- 定义 CGCS2000 投影坐标系 --- cgcs2000_proj = projcrs(code); % --- 坐标转化 --- [lat,lon] = projinv(cgcs2000_proj,x,y);

反过来，lon,lat转为CGCS2000的代码为

[x,y] = projfwd(cgcs2000_proj,lat,lon);