C++使用GDAL库实现Tiff文件的读取
目录
- 库编译
- 项目配置
- 测试
- 测试图片
- 项目结构
- 测试代码
- 运行结果
库编译
首先需要编译对应版本的GDAL库,官方链接:https://gdal.org/ ;
我有GDAL源码、编译好的win32和win64,百度网盘链接:
链接:https://pan.baidu.com/s/1trMAX-Ar6HFPxiy3ghejhg
提取码:prra
解压密码: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库依赖的
测试
测试图片
采用清华大学宫鹏教授学科组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*)GDALO编程客栈pen(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->GetRasterCo开发者_Go入门unt() << 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*1ProjectionisGEOGCS["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,python2.000009)PixelSize = (0.000090,-0.000090)
解析成功。下一篇将会介绍如何读取指定像元值
到此这篇关于C++使用GDAL库实现Tiff文件的读取的文章就介绍到这了,更多相关C++ GDAL Tiff文件内容请搜索我们以前的文章或继续浏览编程下面的相关文章希望大家以后多多支持我们!
精彩评论