用GDAL和OpenCV编写RGB 16位图像

write a RGB 16 bits image with GDAL and OpenCV

本文关键字：RGB 16位 图像 编写 OpenCV GDAL      更新时间：2023-10-16

我需要用R,G,B组件(OpenCV的Mat类型)编写RGB图像。我有一个与GDT_Byte工作的代码，但如果我改为GDT_UInt16不再工作。

metadata_info=(GDALDataset*)GDALOpen(fname.c_str(),GA_ReadOnly);
metadata_info->GetGeoTransform(adfGeoTransform);
target = poDriver->Create(fname_t.c_str(),sX,sY,3,GDT_UInt16,NULL);
target->SetGeoTransform(GT);
target->SetProjection( metadata_info->GetProjectionRef() );
band = target->GetRasterBand(1);
band->RasterIO(GF_Write,0,0,sX,sY,(void *)R.data,sX,sY,GDT_UInt16,0,0);
band = target->GetRasterBand(2);
band->RasterIO(GF_Write,0,0,sX,sY,(void *)G.data,sX,sY,GDT_UInt16,0,0);
band = target->GetRasterBand(3);
band->RasterIO(GF_Write,0,0,sX,sY,(void *)B.data,sX,sY,GDT_UInt16,0,0);

如果我使用8位图像并将GDT_UInt16更改为GDT_Byte和GUInt16更改为GByte，则我没有问题。

*注意:问题不是图像的负载(我认为)，因为我写了一个先验的opencv图像RGB。我阅读图像的组件:(我在架构x64中工作)

img=imread(fname_B.c_str(),CV_LOAD_IMAGE_UNCHANGED|CV_LOAD_IMAGE_ANYDEPTH);
B = Mat(img);

我也试过这个代码来读取R, G和B:

poDataset=(GDALDataset*) GDALOpen(fname.c_str(),GA_ReadOnly);
poBand=poDataset->GetRasterBand(1);
sX=poBand->GetXSize(); sY= poBand->GetYSize();
B=Mat(nYSize, nXSize, CV_16UC1);
poBand->RasterIO(GF_Read,0,0,sX,sY,(void*)B.data,sX,sY,GDT_UInt16,0,0);

同样，GDALINFO的结果是:

Driver: GTiff/GeoTIFF
Files: test_gdal.tif
Size is 7611, 7761
...
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
...
Band 1 Block=7611x1 Type=UInt16, ColorInterp=Gray
Band 2 Block=7611x1 Type=UInt16, ColorInterp=Undefined
Band 3 Block=7611x1 Type=UInt16, ColorInterp=Undefined

我注意到问题只是解释，因为，当我打开创建的RGB文件时，我可以分离频带。

值得注意的是，OpenCV 3.0 beta支持使用GDAL绑定直接用imread加载图像。

https://github.com/Itseez/opencv/blob/master/samples/cpp/tutorial_code/HighGUI/GDAL_IO/gdal-image.cpp 

/// Load generic image as 8-bit color
cv::Mat image = cv::imread(argv[1], cv::IMREAD_LOAD_GDAL | cv::IMREAD_COLOR ); 
/// load a 16-bit image grayscale
cv::Mat dem = cv::imread(argv[2], cv::IMREAD_LOAD_GDAL | cv::IMREAD_ANYDEPTH );
相关文章：