如何将遥感图像读为Emgucv的Image<gray,float>形式

基于VS2010的遥感图像处理，我们通常将遥感图像读为Emgucv中的Image<gray,float>形式。故处理遥感常用的库有Emgucv（基于.NET的opencv）和GDAL（GDAL提供对多种栅格数据的支持）库。

1、GDAL读取遥感数据

GDAL库是一个跨平台的栅格地理数据格式库，包括读取、写入、转换、处理各种栅格数据格式（有些特定的格式对一些操作如写入等不支持）。它使用了一个单一的抽象数据模型就支持了大多数的栅格数据。

GDAL常用的读取遥感数据的函数：

OSGeo.GDAL.Gdal.AllRegister();//对所有支持的文件格式进行注册

采用GDAL读取图像成功后，可以获取图像的一些基本信息，如下：

描述信息：const char* GDALDataset::GetDriver()->GetDescription()，通常是图像的格式

图像大小： 图像宽度 int GDALDataset::GetRasterXSize()

图像高度 int GDALDataset::GetRasterYSize()

波段数：int GDALDataset::GetRasterCount() 波段数即图像每个像素点所含的颜色种类，物理中的光学中学过颜色就是某频率的光波。波段少则一个，多则很多个，在遥感影象中波段通常有多个。

投影信息：GDALDataset::GetProjectionRef() 有的图像没有投影信息，不如一般的JPG、BMP格式图像。

地理坐标信息：double adfGeoTransform[6]

GDALDataset::GetGeoTransform(adfGeoTransform)

地理坐标信息是一个含6个double型数据的数组，adfGeoTransform[1]和adfGeoTransform[5]表示东西和南北方向一个像素对应的距离，adfGeoTransform[0]和adfGeoTransform[3]表示左上角的坐标。

波段信息：数据集中重要的信息，有波段尺寸、数据类型、颜色信息等。

获取波段的方法： GDALRasterBand *poBand;

poBand = poDataset->GetRasterBand(i) poBand为指向第i个波段的指针

波段尺寸：int poBand->GetXSize()

int poBand->GetYSize()

数据类型：const char* GDALGetDataTypeName(poBand->GetRasterDataType())

颜色信息：const char* GDALGetColorInterpretationName(poBand->GetColorInterpretation())

2、读取遥感数据到Image<gray,float>中

//从原图像数据中读取需要的数据到托管数组中
float[] buffer = new float[roi.Width * roi.Height];
srcBand.ReadRaster(
roi.X, roi.Y,
roi.Width, roi.Height,
buffer,
roi.Width, roi.Height,
0, 0);
//分配非托管内存
IntPtr myPtr = System.Runtime.InteropServices.Marshal.AllocHGlobal(
System.Runtime.InteropServices.Marshal.SizeOf(buffer[0]) * buffer.Length);
//将托管数组中的数据赋值到非托管内存中
System.Runtime.InteropServices.Marshal.Copy(buffer, 0, myPtr, buffer.Length);
//利用非托管内存创建EmguCV的Image类的对象，并将数据复制到托管内存中
Image<Gray, float> srcImg = new Image<Gray, float>(
roi.Width, roi.Height,
System.Runtime.InteropServices.Marshal.SizeOf(buffer[0]) * roi.Width,
myPtr)
.Copy();
//释放从进程的非托管内存中分配的内存
System.Runtime.InteropServices.Marshal.FreeHGlobal(myPtr);