使用RasterIO快速显示图像

         这篇算是RasterIO的一个演示程序吧，源代码可以去资源中心去下载。在之前有一篇博客中（地址为：http://blog.csdn.net/liminlu0314/article/details/7072224），GDALRasterIO函数的使用中，有很多人留言，之后又在文章后面更新了一些例子代码。但是还是有很多人不是很清楚，就写个显示的小程序进行说明下吧。

        首先简单的介绍下，我的界面用的QT，VS2008 SP1的IDE环境。程序只能显示8bit的数据，其他的数据可以显示，估计会出问题。下面就对核心的代码进行说明，其他关于界面的编写不做任何说明，有问题请自行Google。先贴代码，后面进行说明：

void CImageDisplayDlg::ShowRaster()
{
	int iStartCol = spinStartCol->value();
	int iStartRow = spinStartRow->value();

	int dataWidth = spinEndCol->value() - iStartCol;
	int dataHeight = spinEndRow->value() - iStartRow;
	int	dataBands =  m_pDataset->GetRasterCount();

	int	band_list[3] = {1,2,3};

	m_dScale = dataHeight > dataWidth ? dataHeight : dataWidth;
	int iViewHeight = 541;
	m_dScale = iViewHeight/m_dScale;

	int iSize = GDALGetDataTypeSize(GDT_Byte) / 8;
	int iScaleWidth = static_cast(dataWidth*m_dScale+0.5);
	int iScaleHeight = static_cast(dataHeight*m_dScale+0.5);

	iScaleWidth = (iScaleWidth*8+31)/32*4;

	unsigned char* pBuffer = new unsigned char[iScaleWidth*iScaleHeight*dataBands];
	CPLErr err = m_pDataset->RasterIO(GF_Read, 0, 0, dataWidth, dataHeight, pBuffer, iScaleWidth, iScaleHeight,
		GDT_Byte, dataBands, band_list, iSize*dataBands, iSize*iScaleWidth*dataBands, iSize);	//读取所有波段的数据

	unsigned char* pDataBuffer = NULL;
	if (dataBands >=3 )
	{
		pDataBuffer = pBuffer;
	}
	else
	{
		pDataBuffer = new unsigned char[iScaleWidth*iScaleHeight*3];
		for (int i=0; isetPixmap(pixmap);
}

        好了，最核心的就是这一个函数，叫ShowRaster。首先获取图像的大小和波段数，我的显示控件的宽度是541个像素，所以根据图像的宽高按照比例缩放到541个像元以内。然后计算缩放后的图像高和宽，就是上面函数中的iScaleWidth和iScaleHeight两个。知道图像缩放后的高度和宽度，这个高度是好使的，但是宽度需要按照公式iScaleWidth = (iScaleWidth*8+31)/32*4再计算一下，要不然图像显示出来可能不正确。关于这个公式可以搜关于位图的一些知识就知道了。

        如果波段数大于3个，就只去前三个进行读取，按照RGBRGB ...RGB的方式进行排列，如果波段数小于3个，那么就读取第一个波段的数据，然后将这个波段的数据复制两次，也就是RGB的值都是同一个。

        这样读取到图像的像元值，然后使用QImage类进行构造一个QImage，再用QImage构造一个QPixmap，最后将这个QPixmap设置给显示的控件就可以了。

        程序代码下载地址：http://download.csdn.net/detail/liminlu0314/4337917。有问题请留言。

        程序截图一：全图显示

程序运行截图二：指定范围显示