使用GDAL库读取Tiff文件

2023-06-20 192

版权

本文内容由阿里云实名注册用户自发贡献，版权归原作者所有，阿里云开发者社区不拥有其著作权，亦不承担相应法律责任。具体规则请查看《阿里云开发者社区用户服务协议》和《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容，填写侵权投诉表单进行举报，一经查实，本社区将立刻删除涉嫌侵权内容。

简介： 使用GDAL库读取Tiff文件

库编译

首先需要编译对应版本的GDAL库，官方链接：https://gdal.org/ ；

我有GDAL源码、编译好的win32和win64,百度网盘链接：

链接：https://pan.baidu.com/s/1-7NE56Uue-Xq7DBuqcRukw?pwd=6666

提取码：6666

解压密码：LinXi

项目配置

依赖库

将上面解压好的文件放在一个固定的路径，注意不要有中文。

我的存放路径如下：

项目配置

1、在项目配置中的配置属性-》常规-》输出目录：.\bin_D\

2、在项目配置中的配置属性-》C/C++》常规-》附加包含目录：这里填写你刚刚存放对应平台的include文件。例如：32位你就选择gdal_vs(86);64位你就选择gdal_vs。

3、在项目配置中的配置属性-》连接器》常规-》附加库目录：这里填写你刚刚存放对应平台的lib文件。例如：32位你就选择gdal_vs(86);64位你就选择gdal_vs。

4、在项目配置中的配置属性-》连接器》输入》附加依赖项：这里不与平台有关，直接写gdal_i_lib即可。

5、然后将gdal*.dll动态库拷贝到工程的输入路径下，proj动态库是gdal库依赖的，具体编译可以看我这篇博客：PROJ 9.1.1源码下载编译（Win10+VS2022）

测试

测试图片

采用清华大学宫鹏教授学科组10m土地覆盖数据，具体可以在官网下载。

项目结构

测试代码

#include "gdal_priv.h"
#include<iostream>
#include <iostream>
#include <list>
#include "gdal_priv.h"   
#include "cpl_conv.h" //for CPLMalloc()   
using namespace std;
void ReadTiffFile(double lon, double lat)
{
  //注册文件格式   
  GDALAllRegister();
  const char* pszFile = "./fromglc10v01_0_100.tif"; // 读取tiff文件
  //const char* pszFile = "./01C_20210101-20220101.tif";//"C:\\Test.img";  //1.jpg
  //使用只读方式打开图像 
  GDALDataset* poDataset = (GDALDataset*)GDALOpen(pszFile, GA_ReadOnly);
  if (poDataset == NULL)
  {
    printf("File: %s不能打开！\n", pszFile);
    return;
  }
  //输出图像的格式信息   
  std::cout << "Driver:" << poDataset->GetDriver()->GetDescription() << "/"
    << poDataset->GetDriver()->GetMetadataItem(GDAL_DMD_LONGNAME) << std::endl;
  //输出图像的大小和波段个数   
  std::cout << "Size is" << poDataset->GetRasterXSize() << "*"
    << poDataset->GetRasterYSize() << "*" << poDataset->GetRasterCount() << std::endl;
  //输出图像的投影信息   
  if (poDataset->GetProjectionRef() != NULL)
  {
    std::cout << "Projectionis" << poDataset->GetProjectionRef() << std::endl;
  }
  GDALRasterBand* poBand;
  int            nBlockXSize, nBlockYSize;
  int            bGotMin, bGotMax;
  double         adfMinMax[2];
  //读取第一个波段   
  poBand = poDataset->GetRasterBand(1);
  int   nXSize = poBand->GetXSize();
  int nYSize = poBand->GetYSize();
  //输出图像的坐标和分辨率信息   
  double adfGeoTransform[6];
  if (poDataset->GetGeoTransform(adfGeoTransform) == CE_None)
  {
    printf("Origin =(%.6f,%.6f)\n",
      adfGeoTransform[0], adfGeoTransform[3]);
    printf("PixelSize = (%.6f,%.6f)\n",
      adfGeoTransform[1], adfGeoTransform[5]);
  }
  CPLFree(poDataset);
  //关闭文件   
  GDALClose((GDALDatasetH)poDataset);
}
int main()
{
  double lon = 101.8615835729;
  double lat = 1.7389243989;
  ReadTiffFile(lon, lat);
  return 0;
}

运行结果

Driver:GTiff/GeoTIFF
Size is22265*22265*1
ProjectionisGEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]
Origin =(99.999918,2.000009)
PixelSize = (0.000090,-0.000090)