GDAL输出影像

6 篇文章 0 订阅
5 篇文章 0 订阅

//输出影像
/**
* @brief	ReadImageToBuff.
* 把数据从内存中输出到文件中,默认保存为tif格式。
* @created:	2012/03/25
* @author:	Zhengwen.Fu
* @param 	szPath		输出路径
* @param 	img			数据块内存块指针
* @param 	nSamples	数据块内存块代表的影像的列数
* @param 	nLines		数据块内存块代表的影像的行数
* @param 	nBandCount	数据块内存块代表的影像的波段数
* @param 	pszFormat	GDAL文件格式标识串,默认是"GTiff",表示存储为tiff格式,其它格式还有:HDF,ENVI等,具体可参考GDAL说明文档中的:GDAL Raster Formats专题,在index.html页面的Supported Formats里;
* @exception 	
* @return 	true	成功
* @return 	false	失败
* @note 	注意szPath指定的文件后缀名要和pszFormat指定的文件格式要求的后缀名一致;
* @remarks 	
* @history:	1.Created by Zhengwen.Fu on 3 25th,2012.      	
*/
//template<typename T>
bool SaveImageToFile(char *szPath,float* img,int nSamples,int nLines,int nBandCount,const char *pszFormat ="GTiff")
{
	GDALAllRegister();         //利用GDAL读取图片,先要进行注册
	CPLSetConfigOption("GDAL_FILENAME_IS_UTF8","NO");   //设置支持中文路径

	//准备保存图片
	GDALDriver *poDriver=GetGDALDriverManager()->GetDriverByName(pszFormat);
	if( poDriver == NULL ) 
	{
		return false;
	}
	char **papszOptions = NULL; 
	GDALDataset *WriteDataSet = poDriver->Create(szPath, nSamples,nLines,nBandCount,GDT_Float32,papszOptions);

	if(WriteDataSet->RasterIO(GF_Write,0,0,nSamples,nLines,img,nSamples,nLines,GDT_Float32,nBandCount,NULL,0,0,0)==CE_Failure )
	{
		return false;
	}
	delete WriteDataSet;WriteDataSet=NULL;
	return true;
}
用法示例:
	float data[10*10*1];
	memset(data,0,100*sizeof(float));

	//输出为tiff格式
	SaveImageToFile("F:\\Work\\数据\\tif.tif",data,10,10,1,"GTiff");

	//HDF4
	SaveImageToFile("F:\\Work\\数据\\HDF4.hdf4",data,10,10,1,"HDF4");
	//不支持

	//HDF5
	SaveImageToFile("F:\\Work\\数据\\HDF5.hdf5",data,10,10,1,"HDF5");
	//不支持

	//ENVI标准影像
	SaveImageToFile("F:\\Work\\ENVI.img",data,10,10,1,"ENVI");
	//ENVI .hdr Labelled Raster采用头文件加主文件的方式保存,头文件以.hdr做后缀名,主文件的后缀可以任意,但通常都使用.img,这正是造成与Erdas Imagine混淆的原因。
	//文件名不能用hdr做后缀,因为写img时会写同名的hdr文件,如果用hdr则造成头文件和主文件是同一文件,造成读写异常

	SaveImageToFile("F:\\Work\\数据\\ENVI2.envi",data,10,10,1,"ENVI");
	//主文件可以使用除hdr之外的任何后缀做文件后缀名

	//Erdas Image文件
	SaveImageToFile("F:\\Work\\数据\\ErdasImagine.img",data,10,10,1,"HFA");
	//输出的文件仅有一个.img文件

	SaveImageToFile("F:\\Work\\数据\\ErdasImagine.eds",data,10,10,1,"HFA");
	//可以不使用.img做后缀,但是你必须知道你输出的这个文件时Erdas Image文件,否则可能用Envi等是无法打开的。


  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
输出栅格影像数据通常需要使用 GDAL 库,下面是一些步骤: 1. 打开需要输出的栅格影像数据文件,可以使用 GDAL 的 Open 函数打开。 2. 创建输出文件,可以使用 GDAL 的 Create 函数创建一个新的栅格影像数据文件,也可以使用 GDAL 的 CreateCopy 函数从现有的栅格影像数据文件创建一个副本。 3. 设置输出文件的地理参考信息和投影信息,可以使用 GDAL 的 SetProjection 和 SetGeoTransform 函数设置。 4. 写入数据到输出文件中,可以使用 GDALWriteRaster 函数写入数据。 5. 关闭输出文件,可以使用 GDAL 的 Close 函数关闭输出文件,释放资源。 下面是一个简单的示例代码: ```python from osgeo import gdal # 打开需要输出的栅格影像数据文件 input_filename = 'input.tif' input_dataset = gdal.Open(input_filename) # 创建输出文件 output_filename = 'output.tif' driver = gdal.GetDriverByName('GTiff') output_dataset = driver.Create(output_filename, input_dataset.RasterXSize, input_dataset.RasterYSize, input_dataset.RasterCount, input_dataset.GetRasterBand(1).DataType) # 设置输出文件的地理参考信息和投影信息 output_dataset.SetProjection(input_dataset.GetProjection()) output_dataset.SetGeoTransform(input_dataset.GetGeoTransform()) # 写入数据到输出文件中 for i in range(input_dataset.RasterCount): band = input_dataset.GetRasterBand(i+1) output_band = output_dataset.GetRasterBand(i+1) output_band.WriteRaster(0, 0, input_dataset.RasterXSize, input_dataset.RasterYSize, band.ReadRaster(0, 0, input_dataset.RasterXSize, input_dataset.RasterYSize)) # 关闭输出文件 output_dataset = None ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

深蓝静音

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

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

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

打赏作者

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

抵扣说明:

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

余额充值