GDAL坐标系转换最好用设置文件名的方式

最新推荐文章于 2023-08-02 19:14:30 发布

directx3d_beginner

最新推荐文章于 2023-08-02 19:14:30 发布

阅读量169

点赞数

CC 4.0 BY-SA版权

分类专栏：遇到的技术问题

本文链接：https://blog.youkuaiyun.com/directx3d_beginner/article/details/116157493

遇到的技术问题专栏收录该内容

104 篇文章

订阅专栏

Arcgis的坐标参考文件一般是以.prj为后缀

即

std::string strPrj1 = "ESRI::" + strProjFileName1;
OGRSpatialReference oSourceSRS;
oSourceSRS.SetFromUserInput(strPrj1.c_str());

   std::string strPrj2 = "ESRI::" + strProjFileName2;
   OGRSpatialReference oTargetSRS;
   oTargetSRS.SetFromUserInput(strPrj2.c_str());

//坐标系转换

   double srcX = 120.235842;
   double srcY = 30.222458;
   double x = srcX;
   double y = srcY;
   OGRCoordinateTransformation *poCT = OGRCreateCoordinateTransformation(&oSourceSRS,
       &oTargetSRS);

   if (poCT == NULL || !poCT->Transform(1, &x, &y))
       printf("Transformation failed.\n");
   else
   {
       printf("(%f,%f) ->(%f,%f)\n", srcX, srcY, x, y);
   }

因为如果解析字符串，容易出错