GDAL直接操纵内存中的影像数据的办法

 

GDAL直接操纵内存中的影像数据的办法

转载自：http://cryolite.iteye.com/blog/358929

 

没有一个简单现成的函数可以做到这个，Frank Warmerdam提示了一种办法，使用“内存文件系统”，但不是通常linux下的那种内存文件系统，而是GDAL自定义的。

首先需要安装一个内存文件处理器
这个特殊的文件处理器可以将一段内存块作为文件处理，这些文件位于"/vsimem/" 目录下（不用找了，在shell下是看不到的），也就是说该目录下的所有文件都被GDAL的内存文件处理器处理。不过目前还不支持目录操作，所以vsimem目录是平的。另外不同的进程有各自的/vsimem/目录，而同一进程下的所有线程共享同一/vsimem/目录。

可以用那些VSI*L函数创建和销毁内存里的数据，就像是操作文件一样。还有些方法可以高效的创建内存文件：不必拷贝原始的数据，这样同样的数据就不会重复出现在多个地方了；可以将某段内存块关联成“GDAL内存文件系统"中的一个文件。

A. 将内存中的一块二进制影像数据buffer读入到GDALDataset，然后进行各种GDAL操作。
一个例子如下：
1. 影像二进制数据在内存中的位置为pabyInData，数据的长度为nInDataLength，我们为它创建一个内存影像文件work.dat
VSIFCloseL( VSIFileFromMemBuffer( "/vsimem/work.dat", pabyInData,
                                       nInDataLength, FALSE ) );
2. 打开这个内存影像文件，得到一个GDALDataset
GDALDatasetH hDS = GDALOpen( "/vsimem/work.dat", GA_ReadOnly );

用完后记得要调用下列函数删除内存中的这个文件：
VSIUnlink( "/vsimem/work.dat" );

B.将GDALDataset转换成某种格式（比如GTiff)的二进制影像数据流。一个例子如下
hDS为想输出为二进制数据的GDALDataset
1. 转换成目标数据格式的GDALDataset，这个GDALDataset是建立在“内存影像文件系统”上的
GDALDriverH hDriver = GDALGetDriverByName( "GTiff" );
GDALDatasetH hOutDS= GDALCreateCopy( hDriver, "/vsimem/out.tif", hDS, TRUE, NULL,
                              NULL, NULL );
hOutDS是一个指针，用完后记得要关闭，不然可能会发生内存泄漏
GDALClose(hOutDS);

2. 调用VSIGetMemFileBuffer函数从内存影像文件中读出二进制数据流

C代码  

1. vsi_l_offset outDataLength;  // 读出的二进制数据流的长度放在这个变量里  
2. int bUnlinkAndSeize = TRUE;   
3. GByte * binData = VSIGetMemFileBuffer("/vsimem/out.tif", &outDataLength, bUnlinkAndSeize);  

vsi_l_offset outDataLength;  // 读出的二进制数据流的长度放在这个变量里
int bUnlinkAndSeize = TRUE;
GByte * binData = VSIGetMemFileBuffer("/vsimem/out.tif", &outDataLength, bUnlinkAndSeize);

binData为转换后的Gtif格式的二进制影像数据，outDataLength为影像数据的长度
bUnlinkAndSeize为TRUE的意思是，内存影像文件将会被自动删除，所以不必再调用VSIUnlink("/vsimem/out.tif")删除内存影像文件了；如果为FALSE，则会保留内存影像文件。

VSIGetMemFileBuffer函数分配的内存要通过CPLFree函数释放：

C代码  

1. CPLFree(binData);  

CPLFree(binData);

VSI*函数内部会调用安装函数VSIInstallMemFileHandler();这个安装函数一般不需要直接调用，反复调用这个安装函数也没什么害处。

虽然还是原来的文件操作，但实际上是对内存直接操作的，只不过披上文件的外衣而已。
http://www.gdal.org/cpl__vsi_8h.html#66e2e6f093fd42f8a941b962d4c8a19e