GDAL通过RasterIO读写图像文件，数据集的创建，图像相关信息的获取

最新推荐文章于 2024-07-04 08:59:40 发布

mjyywd

最新推荐文章于 2024-07-04 08:59:40 发布

阅读量4.8k

点赞数 3

分类专栏： GIS 文章标签：数据 gdal 地图

本文链接：https://blog.csdn.net/mjyywd/article/details/20371739

版权

GIS 专栏收录该内容

3 篇文章 1 订阅

订阅专栏

首先创建数据集

GDALDataset *poDataset; //创建整个图像数据集
GDALRasterBand *poRasterBand;  //创建单个波段数据集
GDALAllRegister(); //注册驱动
poDataset=(GDALDataset *)GDALOpen(  //打开文件
"e:\\124036zuhe.tif"         //文件路径
,GA_ReadOnly);              //打开方式  包括GA_ReadOnly 和  GA_Update
if(poDataset!=NULL)
{
double adfGeoTransform[6];
cout<<"RasterXSize:"<<poDataset->GetRasterXSize()<<endl;
cout<<"RasterYSize:"<<poDataset->GetRasterYSize()<<endl;
cout<<"RasterBand:"<<poDataset->GetRasterBand(1)<<endl;
cout<<"RasterCount:"<<poDataset->GetRasterCount()<<endl;
cout<<"Projetion information:"<<poDataset->GetProjectionRef()<<endl; 
cout<<"仿射:"<<poDataset->GetGeoTransform(adfGeoTransform)<<endl;
cout<<"GcpCount:"<<poDataset->GetGCPCount()<<endl;
cout<<"GcpProjection:"<<poDataset->GetGCPProjection()<<endl;
cout<<"Gcp:"<<poDataset->GetGCPs()<<endl;
poRasterBand=poDataset->GetRasterBand(1);
cout<<"BandXSize:"<<poRasterBand->GetXSize()<<endl;
cout<<"BandYSize:"<<poRasterBand->GetYSize()<<endl;
cout<<"BandDataType:"<<GDALGetDataTypeName(poRasterBand->GetRasterDataType())<<endl;
cout<<"ColorType:"<<GDALGetColorInterpretationName(poRasterBand->GetColorInterpretation())<<endl;
//获取相关信息
float  *buff=new float[1000*1000*3]; //创建读写缓冲区
int num[3]={1,2,3};  //读取的波段顺序
poDataset->RasterIO(GF_Read   //设置读写方式GF_Read和GF_Write
,3000,3000,                   //读取图像的起始位置X，Y
1000,1000,                    //要读取区域的大小，分别为X，Y方向
buff,                         //读入的缓冲区指针
1000,1000,                    //指定缓冲区大小
GDT_Float32,                  //读取图像的数据类型
3,                           //读取的波段数
num                          //读取的波段顺序指针
,0,0,0);                     //先不管
//读取图像到缓冲区
const char *pszFormat = "GTiff";
char *SavePath="D:\\3.tif";
        poDriver = GetGDALDriverManager()->GetDriverByName(pszFormat);  //获取指定驱动类型
        poDatasetCreate=poDriver->Create(                              //创建待写入文件的数据集
        SavePath,                                                      //文件路径
        1000,1000,                                                     //图像文件的大小
        3,                                                             //图像的波段数
        GDT_Float32                                                    //图像数据类型
        ,NULL);                                                        //图像储存方式
        
        double  GeoTransform[6]={ 444720, 30, 0, 3751320, 0, -30 };
poDatasetCreate->RasterIO(GF_Write,0,0,1000,1000,buff,
1000,1000,GDT_Float32,3,num,0,0,0);
//设置图像坐标

poDatasetCreate->RasterIO(                                     //缓冲区内容写入创建的数据集
GF_Write,                                                  //数据集打开方式
0,0,                                                       //写入的起始位置X,Y
1000,1000,                                                 //写入的数据大小X,Y
buff,                                                      //缓冲区指针
1000,1000,                                                 //缓冲区大小
GDT_Float32,                                               //图像数据类型
3,                                                         //波段数目
num,                                                       //波段顺序
0,0,0);