(一)GDAL开发&OGRSpatialReference投影

2 篇文章 1 订阅

这几天在使用GDAL来读取遥感数据,其中涉及到投影的转换。一开始我是在自己的电脑上可以运行得成功的。之后在我们同事的电脑上运行就出现了问题。具体问题是创建创建投影的时候返回了空。我们知道用GDAL定义投影是OGRSpatialReference类来创建,其中如果需要转换,还需要proj.dll库,否则我们的运行程序是无法运行的,或者返回空的。现在使用下面的三种方式来创建了投影参考。

   char    *EPSG4326_WKT = NULL;
	OGRSpatialReference   *myespEPSG4326osr = new OGRSpatialReference("");
	myespEPSG4326osr->SetWellKnownGeogCS("EPSG:4326");

	myespEPSG4326osr->exportToWkt(&EPSG4326_WKT);


	char    *WGS84_WKT = NULL;
	OGRSpatialReference   myespWGS84osr;
	myespWGS84osr.SetWellKnownGeogCS("WGS84");
	myespWGS84osr.exportToWkt(&WGS84_WKT);

	char    *EPSGID_WKT = NULL;
	OGRSpatialReference   myespEPSGIDosr;
	myespEPSGIDosr.importFromEPSG((int)4326);

	myespEPSGIDosr.exportToWkt(&EPSGID_WKT);

使用exportToWkt函数后第一和第三种方式返回的wkt的字符都是为空的。而只有第二种方式不为空,按照了网络相关资料,需要设置环境变量,经过测试后,我的电脑也设置了环境变量但是还是没有成功。其中成功的值如下所示。

GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4326"]]

 

评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

yGIS

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值