卫星影像地理坐标系的转换

最新推荐文章于 2024-06-27 21:48:06 发布

三月微暖寻春笋

最新推荐文章于 2024-06-27 21:48:06 发布

阅读量4.3k

点赞数

分类专栏： GDAL 文章标签： gdal 地理坐标系

本文链接：https://blog.csdn.net/u013230291/article/details/60961943

版权

GDAL 专栏收录该内容

9 篇文章 1 订阅

订阅专栏

影像地理坐标获取以及坐标系之间的坐标转换

    GDALAllRegister();
	CPLSetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");//支持中文路径！
	///读入输入影像
	const char * pszFile = imgfile.c_str();
	GDALDataset *poDataset = (GDALDataset*)GDALOpen(pszFile, GA_ReadOnly);//使用只读方式打开图像 
	if (poDataset == NULL)
	{
		printf("File: %s不能打开！\n", pszFile);
	}
	GDALRasterBand *pBand_temp = NULL;
	pBand_temp = poDataset->GetRasterBand(1);
	///获取四角点经纬度
	double adfGeoTransform[6];
	if (poDataset->GetGeoTransform(adfGeoTransform) == CE_Failure)
	{
		printf("获取仿射变换参数失败");
	}
	std::cout << "仿射变换参数:" << endl << adfGeoTransform[0] << endl << adfGeoTransform[1] << endl << adfGeoTransform[2] << endl << adfGeoTransform[3] << endl << adfGeoTransform[4] << endl << adfGeoTransform[5] << endl;
	double X1 = adfGeoTransform[0] ;//经度
	double Y1 = adfGeoTransform[3];//纬度
	double X2 = adfGeoTransform[0] + iw * adfGeoTransform[1];
	double Y2 = adfGeoTransform[3];
	double X3 = adfGeoTransform[0] + iw  * adfGeoTransform[1];
	double Y3 = adfGeoTransform[3] + ih  * adfGeoTransform[5];
	double X4 = adfGeoTransform[0];
	double Y4 = adfGeoTransform[3] + ih * adfGeoTransform[5];
	/输出影像四个角点的地理坐标
	std::cout << Y1 << endl << X1 << endl;
	std::cout << Y2 << endl << X2 << endl;
	std::cout << Y3 << endl << X3 << endl;
	std::cout << Y4 << endl << X4 << endl;
	double X[4] = { X1, X2, X3 ,X4};
	double Y[4] = { Y1, Y2, Y3, Y4};
	int* flag = NULL;
	///若坐标系为非GWS84,则进行空间坐标系转化
	OGRSpatialReference obj_OSRS;
	obj_OSRS.SetWellKnownGeogCS("WGS84");

	const char* spr=poDataset->GetProjectionRef();
	OGRSpatialReference OSRS(spr);
	if (obj_OSRS.IsSame(&OSRS) == TRUE)
	{
		std::cout << "数据为GWS84地理坐标系统" << std::endl;
		///地理边界
		dataxml.topleft_latitude = Y1;
		dataxml.topleft_longitude = X1;
		dataxml.topright_latitude = Y2;
		dataxml.topright_longitude = X2;
		dataxml.bottomright_latitude = Y3;
		dataxml.bottomright_longitude = X3;
		dataxml.bottomleft_latitude = Y4;
		dataxml.bottomleft_longitude = X4;
	}
	else
	{
		OGRCoordinateTransformation *transf1=OGRCreateCoordinateTransformation(&OSRS, &obj_OSRS);
		if (transf1 == NULL)
		{
			std::cout << "数据不支持转化" << std::endl;
		}
		else
		{
			std::cout << "支持坐标系统之间转化" << std::endl;
		}
		transf1->TransformEx(4, X, Y, NULL, flag);
		///地理边界
		dataxml.topleft_latitude = Y[0];
		dataxml.topleft_longitude = X[0];
		dataxml.topright_latitude = Y[1];
		dataxml.topright_longitude = X[1];
		dataxml.bottomright_latitude = Y[2];
		dataxml.bottomright_longitude = X[2];
		dataxml.bottomleft_latitude = Y[3];
		dataxml.bottomleft_longitude = X[3];
		OGRCoordinateTransformation::DestroyCT(transf1);
	}
	
	/*if (OSRS.IsGeographic() == TRUE)
	{
		cout << "镶嵌数据为地理坐标系统" << endl;
	}
	if (OSRS.IsProjected()==TRUE)
	{

		cout << "镶嵌数据为投影坐标系统" << endl;
	}*/
	
	/*std::string srs_txt = srcdata + "\\" + dataname + ".txt";
	ofstream out_sr(srs_txt);
	out_sr << spr;
	out_sr.close();*/
	
	/*OGRSpatialReference::DestroySpatialReference(&obj_OSRS);
	OGRSpatialReference::DestroySpatialReference(&OSRS);*/

	/输出影像四个角点的地理坐标
	std::cout << Y[0] << endl << X[0] << endl;
	std::cout << Y[1] << endl << X[1] << endl;
	std::cout << Y[2] << endl << X[2] << endl;
	std::cout << Y[3] << endl << X[3] << endl;