(一)GDAL开发&OGRSpatialReference投影

这几天在使用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"]]

 

你可能感兴趣的:(C++,GDAL)