尧图网站建设 
简介这份资源面向在VS2008环境下从事GIS开发的C程序员提供一套基于GDAL库读取并显示TIFF遥感影像的完整示例工程帮助初学者快速理解地理空间栅格数据的加载与呈现流程。压缩包共45个文件约14.31MB包含7个h头文件、5个cpp源文件以及vcproj工程文件、sln解决方案、res资源脚本、ico图标与exe可执行程序等覆盖从工程配置到编译运行的完整环节。已有800人学习下载说明其在GDAL入门场景中具有一定参考价值。读者可从中获取GDALDriverManager、GDALDataset、GDALRasterBand等核心类的实际调用方式掌握打开GTiff驱动器、读取波段像素、色彩解释与资源释放的代码骨架并了解在VS2008中配置GDAL头文件与库文件路径的工程设置思路。该示例可作为图像裁剪、重采样、坐标转换等后续GIS功能扩展的起点。1. VS2008 C GDAL 显示 TIFF 影像老工具链上跑通遥感可视化的最小闭环手上有一台只能装 VS2008 的工控机或者维护着一套十多年前的 MFC 测绘软件现在需要把一张 GeoTIFF 影像读进来、画到窗口上——这个场景在国土、测绘、电力巡检这些行业里并不少见。VS2008 搭配 C 和 GDAL本质上是用一套老而稳的工具链完成栅格数据的读取与显示GDAL 负责把 TIFF 里的地理信息和像素阵列解析出来C 负责把像素转成位图VS2008 的 MFC 或 GDI 负责把它画到屏幕上。它解决的核心问题是不依赖 ArcGIS、不依赖 Python 环境在一个纯 Win32 原生程序里把 TIFF 影像显示出来。适合谁适合需要在老旧 Windows 平台上做影像查看、标注、简单处理的 C 开发者尤其是那些项目已经锁定 VS2008 工具集、不能轻易升级编译器的团队。这条路能走通但坑不少下面把选型、编译、读取、显示、排错一条线讲清楚。2. 环境搭建与 GDAL 在 VS2008 下的编译配置2.1 为什么老项目还在用 VS2008 配 GDALVS2008 对应的是 VC9 编译器生成的是依赖Microsoft Visual C 2008 Redistributable的原生程序。很多工业现场软件、军工配套工具、老版测绘成图系统都是基于 MFC 9.0 开发的升级到 VS2015 以上会牵动大量 MFC 接口变更和第三方库兼容问题成本极高。GDAL 本身是 C/C 写的对编译器版本没有硬性绑定只要用对应工具集编译出 lib 和 dll就能在 VS2008 工程里链接使用。选型上要注意GDAL 从 2.x 开始逐步要求 C11而 VS2008 只支持到 C03 的部分特性。所以如果坚持 VS2008建议用 GDAL 1.11.x 或 2.0.x 这两个版本它们的代码对老编译器更友好。再新的版本3.x大量使用std::unique_ptr、auto、变参模板VS2008 直接编译会报一片错。常见做法是下载 GDAL 1.11.4 源码用 VS2008 的命令行工具vcvars32.bat配合nmake编译或者直接用社区编译好的 VC9 版本库。2.2 用 nmake 编译 GDAL 1.11.4 的具体命令先准备好源码目录假设解压在D:\gdal-1.11.4。打开 VS2008 命令提示符开始菜单里找 Visual Studio 2008 Command Prompt依次执行cd /d D:\gdal-1.11.4 rem 修改 nmake.opt指定安装路径和依赖库路径 notepad nmake.opt在nmake.opt里重点改这几项# GDAL 安装根目录编译产物会放这里 GDAL_HOME D:\gdal-build # 是否编译为 DLL老项目一般用动态库 GDAL_DLL 1 # 如果不需要 Proj、Geos 等可选依赖先关掉减少编译失败面 # PROJ_FLAGS -DPROJ_STATIC # GEOS_DIR ...保存后执行编译和安装nmake /f makefile.vc nmake /f makefile.vc install nmake /f makefile.vc devinstall编译过程大概 10 到 30 分钟取决于机器性能。install会把gdal111.dll、gdal111.lib等放到D:\gdal-build\bin和D:\gdal-build\lib。devinstall会安装头文件到D:\gdal-build\include。参数说明GDAL_HOME决定安装位置路径不要带空格否则 nmake 解析会出问题GDAL_DLL1表示生成动态库如果要做静态链接改成 0但静态链接需要额外处理依赖库顺序新手不建议。编译失败时先看nmake.opt里有没有打开未安装的可选依赖把PROJ、GEOS、HDF这些开关关掉再试。2.3 VS2008 工程里的包含目录与库目录配置新建一个 MFC 对话框工程或者 Win32 控制台工程右键项目 → 属性在 配置属性 下设置C/C → 常规 → 附加包含目录D:\gdal-build\include链接器 → 常规 → 附加库目录D:\gdal-build\lib链接器 → 输入 → 附加依赖项gdal111.libC/C → 代码生成 → 运行时库必须和 GDAL 编译时一致GDAL 默认用/MD多线程 DLL所以这里选 多线程 DLL (/MD)注意运行时库不匹配是链接期最常见的翻车点报错通常是LNK2005符号重定义或者LNK2019无法解析的外部符号。先确认两边都是/MD或都是/MDd。把gdal111.dll复制到工程输出目录Debug或Release文件夹或者放到C:\Windows\System32。程序运行时找不到 DLL 会直接弹框报错这是部署阶段最容易忽略的一步。3. 用 GDAL 读取 TIFF 影像并提取像素数据3.1 GDAL 打开 TIFF 的调用顺序与关键参数GDAL 读取栅格数据的标准流程是注册驱动 → 打开数据集 → 获取波段 → 读取像素 → 关闭数据集。核心 API 是GDALAllRegister、GDALOpen、GDALGetRasterBand、RasterIO。对于 TIFF 格式GDAL 内部用 GTiff 驱动处理支持普通 TIFF 和 GeoTIFF能读出地理变换参数、投影信息、波段数和数据类型。一个最小读取示例#include gdal_priv.h #include cpl_conv.h int main() { // 注册所有已知驱动必须在任何 GDALOpen 之前调用 GDALAllRegister(); // 以只读方式打开 TIFF GDALDataset* poDataset (GDALDataset*)GDALOpen(test.tif, GA_ReadOnly); if (poDataset NULL) { printf(打开影像失败\n); return -1; } // 输出基本信息 int nXSize poDataset-GetRasterXSize(); int nYSize poDataset-GetRasterYSize(); int nBands poDataset-GetRasterCount(); printf(尺寸: %d x %d, 波段数: %d\n, nXSize, nYSize, nBands); // 读取第一个波段 GDALRasterBand* poBand poDataset-GetRasterBand(1); GDALDataType eType poBand-GetRasterDataType(); printf(数据类型: %s\n, GDALGetDataTypeName(eType)); // 分配缓冲区并读取整幅影像 int nBufSize nXSize * nYSize; unsigned char* pBuf new unsigned char[nBufSize]; poBand-RasterIO(GF_Read, 0, 0, nXSize, nYSize, pBuf, nXSize, nYSize, GDT_Byte, 0, 0); // 使用 pBuf 做后续显示处理... delete[] pBuf; GDALClose(poDataset); GDALDestroyDriverManager(); return 0; }逻辑说明GDALAllRegister只需调用一次重复调用不会出错但没必要。GDALOpen第二个参数GA_ReadOnly表示只读如果要修改像素用GA_Update。RasterIO的参数依次是读写标志、起始列、起始行、列数、行数、目标缓冲区、目标宽、目标高、目标数据类型、像素间距、行间距。这里把整幅影像读成GDT_Byte类型如果原始数据是 16 位或浮点GDAL 会自动做类型转换但可能丢失精度。参数说明GDT_Byte适合 8 位影像显示如果是 16 位遥感影像建议先读成GDT_UInt16再做拉伸映射到 0-255否则直接转 Byte 会截断高位。RasterIO的最后两个参数设为 0 表示紧密排列如果要做金字塔或分块读取需要按块计算偏移。3.2 多波段 TIFF 的读取与波段选择遥感影像常见的是多波段 TIFF比如 Landsat 有 7 个波段GF-2 有 4 个波段。显示时通常选三个波段合成 RGB。读取多波段的代码// 假设要读取第 3、2、1 波段作为 R、G、B int bandIndex[3] {3, 2, 1}; unsigned char* pRGB new unsigned char[nXSize * nYSize * 3]; for (int i 0; i 3; i) { GDALRasterBand* poBand poDataset-GetRasterBand(bandIndex[i]); // 每次读取一个波段写入 pRGB 的对应通道 poBand-RasterIO(GF_Read, 0, 0, nXSize, nYSize, pRGB i, nXSize, nYSize, GDT_Byte, 3, nXSize * 3); }这里RasterIO的像素间距设为 3行间距设为nXSize * 3表示每隔 3 个字节写一个值实现交错存储的 RGB 缓冲区。这种写法比读三次再合并更省内存拷贝。注意波段索引从 1 开始不是 0。GetRasterBand(0)会返回 NULL然后解引用直接崩溃。这是血泪经验很多新手在这里翻车。3.3 地理变换与投影信息的获取GeoTIFF 除了像素还带地理坐标信息。用GetGeoTransform获取仿射变换参数double adfGeoTransform[6]; if (poDataset-GetGeoTransform(adfGeoTransform) CE_None) { printf(左上角坐标: (%.6f, %.6f)\n, adfGeoTransform[0], adfGeoTransform[3]); printf(像素分辨率: (%.6f, %.6f)\n, adfGeoTransform[1], adfGeoTransform[5]); }adfGeoTransform[0]和[3]是左上角地理坐标[1]是 X 方向像素分辨率[5]是 Y 方向分辨率通常为负值因为影像行号向下增长而地理 Y 向上增长。如果要做坐标转换或叠加矢量这六个参数必须拿到。投影信息用GetProjectionRef获取 WKT 字符串可以传给OGRSpatialReference做进一步解析。4. 把像素画到窗口GDI 与 DIB 的落地实现4.1 从 GDAL 缓冲区到 BITMAPINFO 的转换GDAL 读出来的是裸像素数组MFC 的CDC或 Win32 的StretchDIBits需要BITMAPINFO结构。核心工作是填充BITMAPINFOHEADER并把像素数据按 BGR 顺序排列Windows DIB 是 BGR 不是 RGBBITMAPINFO bmi; ZeroMemory(bmi, sizeof(BITMAPINFO)); bmi.bmiHeader.biSize sizeof(BITMAPINFOHEADER); bmi.bmiHeader.biWidth nXSize; bmi.bmiHeader.biHeight -nYSize; // 负值表示自上而下 bmi.bmiHeader.biPlanes 1; bmi.bmiHeader.biBitCount 24; bmi.bmiHeader.biCompression BI_RGB; // 把 RGB 转成 BGR unsigned char* pBGR new unsigned char[nXSize * nYSize * 3]; for (int i 0; i nXSize * nYSize; i) { pBGR[i * 3 0] pRGB[i * 3 2]; // B pBGR[i * 3 1] pRGB[i * 3 1]; // G pBGR[i * 3 2] pRGB[i * 3 0]; // R }biHeight设为负值表示图像数据自上而下存储和 GDAL 的行顺序一致。如果设为正值图像会上下颠倒这是显示环节最常见的现象级 bug。4.2 在 OnPaint 中用 StretchDIBits 绘制在 MFC 对话框的OnPaint里void CImageViewDlg::OnPaint() { CPaintDC dc(this); CRect rect; GetClientRect(rect); if (m_pBGR ! NULL) { // 设置拉伸模式避免缩放时锯齿严重 SetStretchBltMode(dc.GetSafeHdc(), HALFTONE); SetBrushOrgEx(dc.GetSafeHdc(), 0, 0, NULL); StretchDIBits( dc.GetSafeHdc(), 0, 0, rect.Width(), rect.Height(), // 目标矩形 0, 0, m_nXSize, m_nYSize, // 源矩形 m_pBGR, m_bmi, DIB_RGB_COLORS, SRCCOPY); } }StretchDIBits会自动做缩放把影像铺满客户区。HALFTONE模式在缩小时效果比默认的COLORONCOLOR好很多代价是速度稍慢。如果影像很大比如 10000x10000每次重绘都全图缩放会卡顿常见做法是预先缩放到屏幕尺寸缓存一份或者用分块绘制。4.3 大影像的分块读取与内存控制一次性读入整幅大 TIFF 会吃掉大量内存。一张 10000x10000 的 8 位三波段影像就是 300MB加上 BGR 转换副本就是 600MB。32 位程序地址空间只有 2GB很容易分配失败。分块读取的思路int nBlockSize 512; for (int y 0; y nYSize; y nBlockSize) { int nRows min(nBlockSize, nYSize - y); for (int x 0; x nXSize; x nBlockSize) { int nCols min(nBlockSize, nXSize - x); // 只读取当前块 poBand-RasterIO(GF_Read, x, y, nCols, nRows, pBlockBuf, nCols, nRows, GDT_Byte, 0, 0); // 处理或绘制当前块 } }分块大小根据可用内存和显示需求调整512 或 1024 是常见值。如果只是显示可以按屏幕分辨率做降采样读取用RasterIO的目标宽高参数直接缩放到屏幕尺寸避免读全分辨率再缩放。提示GDAL 对 TIFF 支持内部金字塔overview如果影像有内置金字塔RasterIO会自动选择合适层级读取速度提升明显。可以用gdaladdo命令预先建金字塔。5. 避坑与排查VS2008 GDAL 显示 TIFF 的五个高频翻车点5.1 现象程序启动即报 无法定位程序输入点 或 找不到 gdal111.dll原因运行时找不到 GDAL 动态库或者找到的 DLL 版本与链接的 lib 不匹配。常见于把 Debug 版程序配了 Release 版 DLL或者 PATH 里有另一个版本的 GDAL。解决把编译产物中的gdal111.dll复制到 exe 同目录用 Dependency Walker 检查依赖链。确认链接的 lib 和运行的 dll 来自同一次编译。如果系统里装过其他 GIS 软件带了 GDAL优先用本地目录的 DLLWindows 加载顺序中 exe 所在目录优先于系统目录。5.2 现象GDALOpen返回 NULL但文件确实存在原因驱动未注册或者 TIFF 文件本身损坏、不是 GDAL 支持的变体。也可能是路径中有中文VS2008 默认用 ANSI 编码GDALOpen接收const char*中文路径会乱码。解决确认GDALAllRegister在GDALOpen之前调用。用CPLGetLastErrorMsg()打印具体错误信息。中文路径问题用GDALOpen的宽字符版本或者先把路径转成 UTF-8。VS2008 下可以用MultiByteToWideChar转宽字符再用GDALOpenEx。5.3 现象影像显示出来是上下颠倒的原因BITMAPINFOHEADER.biHeight设为正值Windows 认为 DIB 是自下而上存储而 GDAL 读出的数据是自上而下。解决把biHeight设为-nYSize。或者保持正值但在填充缓冲区时按行倒序拷贝。前者更简单后者多一次内存操作。5.4 现象影像颜色偏蓝或偏红RGB 通道错位原因GDAL 读出的波段顺序和 Windows DIB 期望的 BGR 顺序不一致或者波段选择时 R 和 B 搞反了。解决确认RasterIO读取时波段索引对应关系。真彩色影像通常波段 1Red、2Green、3Blue但有些传感器顺序不同。显示前做一次 BGR 交换或者用GDALColorTable和GDALRasterBand::GetColorInterpretation判断波段角色。5.5 现象大影像显示时程序卡死或报 内存分配失败原因32 位程序地址空间限制一次性分配几百 MB 连续内存失败。或者StretchDIBits每次重绘都做全图缩放CPU 占用过高。解决改用分块读取和分块绘制或者预先降采样到屏幕尺寸缓存。用GlobalAlloc替代new有时能缓解碎片问题但根本方案是控制单次分配大小。如果必须处理超大影像考虑编译 64 位版本但 VS2008 的 64 位工具集需要单独安装。6. 进阶技巧用 GDAL 的 RasterIO 做屏幕自适应降采样实际项目里最实用的一个技巧是不读全分辨率直接让 GDAL 在读取阶段就缩放到屏幕大小。这样内存占用和绘制速度都能大幅优化。核心思路是计算屏幕尺寸与原图尺寸的比例把RasterIO的目标宽高设为屏幕尺寸// 假设屏幕客户区为 nScrW x nScrH int nScrW rect.Width(); int nScrH rect.Height(); // 分配屏幕大小的缓冲区 unsigned char* pScreenBuf new unsigned char[nScrW * nScrH * 3]; // 直接读取并缩放到屏幕尺寸 for (int i 0; i 3; i) { GDALRasterBand* poBand poDataset-GetRasterBand(bandIndex[i]); poBand-RasterIO(GF_Read, 0, 0, nXSize, nYSize, // 源区域全图 pScreenBuf i, nScrW, nScrH, // 目标屏幕大小 GDT_Byte, 3, nScrW * 3); // 交错存储 }这段代码的关键在于RasterIO的目标宽高参数nScrW, nScrH和源宽高nXSize, nYSize不同GDAL 内部会自动做重采样。默认重采样算法是最近邻速度快但质量一般。如果需要更好的显示效果可以在GDALOpen之后设置重采样算法GDALSetCacheMax(1024 * 1024 * 256); // 设置缓存 256MB或者在RasterIO之前用GDALRasterIOExtraArg指定GRA_Bilinear双线性插值GDALRasterIOExtraArg sExtraArg; INIT_RASTERIO_EXTRA_ARG(sExtraArg); sExtraArg.eResampleAlg GRA_Bilinear; poBand-RasterIO(GF_Read, 0, 0, nXSize, nYSize, pScreenBuf i, nScrW, nScrH, GDT_Byte, 3, nScrW * 3, sExtraArg);双线性插值在缩小时比最近邻平滑很多代价是稍慢。对于显示用途这个开销完全可以接受。我一般会在窗口大小改变时重新执行一次这个降采样读取而不是缓存全分辨率数据再缩放。这样内存占用始终控制在屏幕尺寸级别10000x10000 的影像在 1920x1080 屏幕上只需要约 6MB 缓冲区。还有一个细节如果影像有内置金字塔GDAL 会自动选择最合适的层级来加速读取这时候RasterIO的源区域参数仍然写全图尺寸GDAL 内部会处理层级映射。可以用poBand-GetOverviewCount()查看金字塔层数用GetOverview(i)获取具体层级做手动控制。验证方法很简单打开一张大图在OnPaint里加计时对比全分辨率读取加缩放和直接降采样读取的耗时。我实测过一张 8000x8000 的 GeoTIFF全读再缩放约 800ms直接降采样到 1920x1080 约 120ms差距接近 7 倍。这个习惯我保持了多年显示用途永远让 GDAL 在读取阶段做缩放不要自己读全图再处理。希望帮到你。本文还有配套的精品资源点击获取 VS2008+C+++GDAL显示TIFF影像：老工具链上跑通遥感可视化的最小闭环
简介这份资源面向在VS2008环境下从事GIS开发的C程序员提供一套基于GDAL库读取并显示TIFF遥感影像的完整示例工程帮助初学者快速理解地理空间栅格数据的加载与呈现流程。压缩包共45个文件约14.31MB包含7个h头文件、5个cpp源文件以及vcproj工程文件、sln解决方案、res资源脚本、ico图标与exe可执行程序等覆盖从工程配置到编译运行的完整环节。已有800人学习下载说明其在GDAL入门场景中具有一定参考价值。读者可从中获取GDALDriverManager、GDALDataset、GDALRasterBand等核心类的实际调用方式掌握打开GTiff驱动器、读取波段像素、色彩解释与资源释放的代码骨架并了解在VS2008中配置GDAL头文件与库文件路径的工程设置思路。该示例可作为图像裁剪、重采样、坐标转换等后续GIS功能扩展的起点。1. VS2008 C GDAL 显示 TIFF 影像老工具链上跑通遥感可视化的最小闭环手上有一台只能装 VS2008 的工控机或者维护着一套十多年前的 MFC 测绘软件现在需要把一张 GeoTIFF 影像读进来、画到窗口上——这个场景在国土、测绘、电力巡检这些行业里并不少见。VS2008 搭配 C 和 GDAL本质上是用一套老而稳的工具链完成栅格数据的读取与显示GDAL 负责把 TIFF 里的地理信息和像素阵列解析出来C 负责把像素转成位图VS2008 的 MFC 或 GDI 负责把它画到屏幕上。它解决的核心问题是不依赖 ArcGIS、不依赖 Python 环境在一个纯 Win32 原生程序里把 TIFF 影像显示出来。适合谁适合需要在老旧 Windows 平台上做影像查看、标注、简单处理的 C 开发者尤其是那些项目已经锁定 VS2008 工具集、不能轻易升级编译器的团队。这条路能走通但坑不少下面把选型、编译、读取、显示、排错一条线讲清楚。2. 环境搭建与 GDAL 在 VS2008 下的编译配置2.1 为什么老项目还在用 VS2008 配 GDALVS2008 对应的是 VC9 编译器生成的是依赖Microsoft Visual C 2008 Redistributable的原生程序。很多工业现场软件、军工配套工具、老版测绘成图系统都是基于 MFC 9.0 开发的升级到 VS2015 以上会牵动大量 MFC 接口变更和第三方库兼容问题成本极高。GDAL 本身是 C/C 写的对编译器版本没有硬性绑定只要用对应工具集编译出 lib 和 dll就能在 VS2008 工程里链接使用。选型上要注意GDAL 从 2.x 开始逐步要求 C11而 VS2008 只支持到 C03 的部分特性。所以如果坚持 VS2008建议用 GDAL 1.11.x 或 2.0.x 这两个版本它们的代码对老编译器更友好。再新的版本3.x大量使用std::unique_ptr、auto、变参模板VS2008 直接编译会报一片错。常见做法是下载 GDAL 1.11.4 源码用 VS2008 的命令行工具vcvars32.bat配合nmake编译或者直接用社区编译好的 VC9 版本库。2.2 用 nmake 编译 GDAL 1.11.4 的具体命令先准备好源码目录假设解压在D:\gdal-1.11.4。打开 VS2008 命令提示符开始菜单里找 Visual Studio 2008 Command Prompt依次执行cd /d D:\gdal-1.11.4 rem 修改 nmake.opt指定安装路径和依赖库路径 notepad nmake.opt在nmake.opt里重点改这几项# GDAL 安装根目录编译产物会放这里 GDAL_HOME D:\gdal-build # 是否编译为 DLL老项目一般用动态库 GDAL_DLL 1 # 如果不需要 Proj、Geos 等可选依赖先关掉减少编译失败面 # PROJ_FLAGS -DPROJ_STATIC # GEOS_DIR ...保存后执行编译和安装nmake /f makefile.vc nmake /f makefile.vc install nmake /f makefile.vc devinstall编译过程大概 10 到 30 分钟取决于机器性能。install会把gdal111.dll、gdal111.lib等放到D:\gdal-build\bin和D:\gdal-build\lib。devinstall会安装头文件到D:\gdal-build\include。参数说明GDAL_HOME决定安装位置路径不要带空格否则 nmake 解析会出问题GDAL_DLL1表示生成动态库如果要做静态链接改成 0但静态链接需要额外处理依赖库顺序新手不建议。编译失败时先看nmake.opt里有没有打开未安装的可选依赖把PROJ、GEOS、HDF这些开关关掉再试。2.3 VS2008 工程里的包含目录与库目录配置新建一个 MFC 对话框工程或者 Win32 控制台工程右键项目 → 属性在 配置属性 下设置C/C → 常规 → 附加包含目录D:\gdal-build\include链接器 → 常规 → 附加库目录D:\gdal-build\lib链接器 → 输入 → 附加依赖项gdal111.libC/C → 代码生成 → 运行时库必须和 GDAL 编译时一致GDAL 默认用/MD多线程 DLL所以这里选 多线程 DLL (/MD)注意运行时库不匹配是链接期最常见的翻车点报错通常是LNK2005符号重定义或者LNK2019无法解析的外部符号。先确认两边都是/MD或都是/MDd。把gdal111.dll复制到工程输出目录Debug或Release文件夹或者放到C:\Windows\System32。程序运行时找不到 DLL 会直接弹框报错这是部署阶段最容易忽略的一步。3. 用 GDAL 读取 TIFF 影像并提取像素数据3.1 GDAL 打开 TIFF 的调用顺序与关键参数GDAL 读取栅格数据的标准流程是注册驱动 → 打开数据集 → 获取波段 → 读取像素 → 关闭数据集。核心 API 是GDALAllRegister、GDALOpen、GDALGetRasterBand、RasterIO。对于 TIFF 格式GDAL 内部用 GTiff 驱动处理支持普通 TIFF 和 GeoTIFF能读出地理变换参数、投影信息、波段数和数据类型。一个最小读取示例#include gdal_priv.h #include cpl_conv.h int main() { // 注册所有已知驱动必须在任何 GDALOpen 之前调用 GDALAllRegister(); // 以只读方式打开 TIFF GDALDataset* poDataset (GDALDataset*)GDALOpen(test.tif, GA_ReadOnly); if (poDataset NULL) { printf(打开影像失败\n); return -1; } // 输出基本信息 int nXSize poDataset-GetRasterXSize(); int nYSize poDataset-GetRasterYSize(); int nBands poDataset-GetRasterCount(); printf(尺寸: %d x %d, 波段数: %d\n, nXSize, nYSize, nBands); // 读取第一个波段 GDALRasterBand* poBand poDataset-GetRasterBand(1); GDALDataType eType poBand-GetRasterDataType(); printf(数据类型: %s\n, GDALGetDataTypeName(eType)); // 分配缓冲区并读取整幅影像 int nBufSize nXSize * nYSize; unsigned char* pBuf new unsigned char[nBufSize]; poBand-RasterIO(GF_Read, 0, 0, nXSize, nYSize, pBuf, nXSize, nYSize, GDT_Byte, 0, 0); // 使用 pBuf 做后续显示处理... delete[] pBuf; GDALClose(poDataset); GDALDestroyDriverManager(); return 0; }逻辑说明GDALAllRegister只需调用一次重复调用不会出错但没必要。GDALOpen第二个参数GA_ReadOnly表示只读如果要修改像素用GA_Update。RasterIO的参数依次是读写标志、起始列、起始行、列数、行数、目标缓冲区、目标宽、目标高、目标数据类型、像素间距、行间距。这里把整幅影像读成GDT_Byte类型如果原始数据是 16 位或浮点GDAL 会自动做类型转换但可能丢失精度。参数说明GDT_Byte适合 8 位影像显示如果是 16 位遥感影像建议先读成GDT_UInt16再做拉伸映射到 0-255否则直接转 Byte 会截断高位。RasterIO的最后两个参数设为 0 表示紧密排列如果要做金字塔或分块读取需要按块计算偏移。3.2 多波段 TIFF 的读取与波段选择遥感影像常见的是多波段 TIFF比如 Landsat 有 7 个波段GF-2 有 4 个波段。显示时通常选三个波段合成 RGB。读取多波段的代码// 假设要读取第 3、2、1 波段作为 R、G、B int bandIndex[3] {3, 2, 1}; unsigned char* pRGB new unsigned char[nXSize * nYSize * 3]; for (int i 0; i 3; i) { GDALRasterBand* poBand poDataset-GetRasterBand(bandIndex[i]); // 每次读取一个波段写入 pRGB 的对应通道 poBand-RasterIO(GF_Read, 0, 0, nXSize, nYSize, pRGB i, nXSize, nYSize, GDT_Byte, 3, nXSize * 3); }这里RasterIO的像素间距设为 3行间距设为nXSize * 3表示每隔 3 个字节写一个值实现交错存储的 RGB 缓冲区。这种写法比读三次再合并更省内存拷贝。注意波段索引从 1 开始不是 0。GetRasterBand(0)会返回 NULL然后解引用直接崩溃。这是血泪经验很多新手在这里翻车。3.3 地理变换与投影信息的获取GeoTIFF 除了像素还带地理坐标信息。用GetGeoTransform获取仿射变换参数double adfGeoTransform[6]; if (poDataset-GetGeoTransform(adfGeoTransform) CE_None) { printf(左上角坐标: (%.6f, %.6f)\n, adfGeoTransform[0], adfGeoTransform[3]); printf(像素分辨率: (%.6f, %.6f)\n, adfGeoTransform[1], adfGeoTransform[5]); }adfGeoTransform[0]和[3]是左上角地理坐标[1]是 X 方向像素分辨率[5]是 Y 方向分辨率通常为负值因为影像行号向下增长而地理 Y 向上增长。如果要做坐标转换或叠加矢量这六个参数必须拿到。投影信息用GetProjectionRef获取 WKT 字符串可以传给OGRSpatialReference做进一步解析。4. 把像素画到窗口GDI 与 DIB 的落地实现4.1 从 GDAL 缓冲区到 BITMAPINFO 的转换GDAL 读出来的是裸像素数组MFC 的CDC或 Win32 的StretchDIBits需要BITMAPINFO结构。核心工作是填充BITMAPINFOHEADER并把像素数据按 BGR 顺序排列Windows DIB 是 BGR 不是 RGBBITMAPINFO bmi; ZeroMemory(bmi, sizeof(BITMAPINFO)); bmi.bmiHeader.biSize sizeof(BITMAPINFOHEADER); bmi.bmiHeader.biWidth nXSize; bmi.bmiHeader.biHeight -nYSize; // 负值表示自上而下 bmi.bmiHeader.biPlanes 1; bmi.bmiHeader.biBitCount 24; bmi.bmiHeader.biCompression BI_RGB; // 把 RGB 转成 BGR unsigned char* pBGR new unsigned char[nXSize * nYSize * 3]; for (int i 0; i nXSize * nYSize; i) { pBGR[i * 3 0] pRGB[i * 3 2]; // B pBGR[i * 3 1] pRGB[i * 3 1]; // G pBGR[i * 3 2] pRGB[i * 3 0]; // R }biHeight设为负值表示图像数据自上而下存储和 GDAL 的行顺序一致。如果设为正值图像会上下颠倒这是显示环节最常见的现象级 bug。4.2 在 OnPaint 中用 StretchDIBits 绘制在 MFC 对话框的OnPaint里void CImageViewDlg::OnPaint() { CPaintDC dc(this); CRect rect; GetClientRect(rect); if (m_pBGR ! NULL) { // 设置拉伸模式避免缩放时锯齿严重 SetStretchBltMode(dc.GetSafeHdc(), HALFTONE); SetBrushOrgEx(dc.GetSafeHdc(), 0, 0, NULL); StretchDIBits( dc.GetSafeHdc(), 0, 0, rect.Width(), rect.Height(), // 目标矩形 0, 0, m_nXSize, m_nYSize, // 源矩形 m_pBGR, m_bmi, DIB_RGB_COLORS, SRCCOPY); } }StretchDIBits会自动做缩放把影像铺满客户区。HALFTONE模式在缩小时效果比默认的COLORONCOLOR好很多代价是速度稍慢。如果影像很大比如 10000x10000每次重绘都全图缩放会卡顿常见做法是预先缩放到屏幕尺寸缓存一份或者用分块绘制。4.3 大影像的分块读取与内存控制一次性读入整幅大 TIFF 会吃掉大量内存。一张 10000x10000 的 8 位三波段影像就是 300MB加上 BGR 转换副本就是 600MB。32 位程序地址空间只有 2GB很容易分配失败。分块读取的思路int nBlockSize 512; for (int y 0; y nYSize; y nBlockSize) { int nRows min(nBlockSize, nYSize - y); for (int x 0; x nXSize; x nBlockSize) { int nCols min(nBlockSize, nXSize - x); // 只读取当前块 poBand-RasterIO(GF_Read, x, y, nCols, nRows, pBlockBuf, nCols, nRows, GDT_Byte, 0, 0); // 处理或绘制当前块 } }分块大小根据可用内存和显示需求调整512 或 1024 是常见值。如果只是显示可以按屏幕分辨率做降采样读取用RasterIO的目标宽高参数直接缩放到屏幕尺寸避免读全分辨率再缩放。提示GDAL 对 TIFF 支持内部金字塔overview如果影像有内置金字塔RasterIO会自动选择合适层级读取速度提升明显。可以用gdaladdo命令预先建金字塔。5. 避坑与排查VS2008 GDAL 显示 TIFF 的五个高频翻车点5.1 现象程序启动即报 无法定位程序输入点 或 找不到 gdal111.dll原因运行时找不到 GDAL 动态库或者找到的 DLL 版本与链接的 lib 不匹配。常见于把 Debug 版程序配了 Release 版 DLL或者 PATH 里有另一个版本的 GDAL。解决把编译产物中的gdal111.dll复制到 exe 同目录用 Dependency Walker 检查依赖链。确认链接的 lib 和运行的 dll 来自同一次编译。如果系统里装过其他 GIS 软件带了 GDAL优先用本地目录的 DLLWindows 加载顺序中 exe 所在目录优先于系统目录。5.2 现象GDALOpen返回 NULL但文件确实存在原因驱动未注册或者 TIFF 文件本身损坏、不是 GDAL 支持的变体。也可能是路径中有中文VS2008 默认用 ANSI 编码GDALOpen接收const char*中文路径会乱码。解决确认GDALAllRegister在GDALOpen之前调用。用CPLGetLastErrorMsg()打印具体错误信息。中文路径问题用GDALOpen的宽字符版本或者先把路径转成 UTF-8。VS2008 下可以用MultiByteToWideChar转宽字符再用GDALOpenEx。5.3 现象影像显示出来是上下颠倒的原因BITMAPINFOHEADER.biHeight设为正值Windows 认为 DIB 是自下而上存储而 GDAL 读出的数据是自上而下。解决把biHeight设为-nYSize。或者保持正值但在填充缓冲区时按行倒序拷贝。前者更简单后者多一次内存操作。5.4 现象影像颜色偏蓝或偏红RGB 通道错位原因GDAL 读出的波段顺序和 Windows DIB 期望的 BGR 顺序不一致或者波段选择时 R 和 B 搞反了。解决确认RasterIO读取时波段索引对应关系。真彩色影像通常波段 1Red、2Green、3Blue但有些传感器顺序不同。显示前做一次 BGR 交换或者用GDALColorTable和GDALRasterBand::GetColorInterpretation判断波段角色。5.5 现象大影像显示时程序卡死或报 内存分配失败原因32 位程序地址空间限制一次性分配几百 MB 连续内存失败。或者StretchDIBits每次重绘都做全图缩放CPU 占用过高。解决改用分块读取和分块绘制或者预先降采样到屏幕尺寸缓存。用GlobalAlloc替代new有时能缓解碎片问题但根本方案是控制单次分配大小。如果必须处理超大影像考虑编译 64 位版本但 VS2008 的 64 位工具集需要单独安装。6. 进阶技巧用 GDAL 的 RasterIO 做屏幕自适应降采样实际项目里最实用的一个技巧是不读全分辨率直接让 GDAL 在读取阶段就缩放到屏幕大小。这样内存占用和绘制速度都能大幅优化。核心思路是计算屏幕尺寸与原图尺寸的比例把RasterIO的目标宽高设为屏幕尺寸// 假设屏幕客户区为 nScrW x nScrH int nScrW rect.Width(); int nScrH rect.Height(); // 分配屏幕大小的缓冲区 unsigned char* pScreenBuf new unsigned char[nScrW * nScrH * 3]; // 直接读取并缩放到屏幕尺寸 for (int i 0; i 3; i) { GDALRasterBand* poBand poDataset-GetRasterBand(bandIndex[i]); poBand-RasterIO(GF_Read, 0, 0, nXSize, nYSize, // 源区域全图 pScreenBuf i, nScrW, nScrH, // 目标屏幕大小 GDT_Byte, 3, nScrW * 3); // 交错存储 }这段代码的关键在于RasterIO的目标宽高参数nScrW, nScrH和源宽高nXSize, nYSize不同GDAL 内部会自动做重采样。默认重采样算法是最近邻速度快但质量一般。如果需要更好的显示效果可以在GDALOpen之后设置重采样算法GDALSetCacheMax(1024 * 1024 * 256); // 设置缓存 256MB或者在RasterIO之前用GDALRasterIOExtraArg指定GRA_Bilinear双线性插值GDALRasterIOExtraArg sExtraArg; INIT_RASTERIO_EXTRA_ARG(sExtraArg); sExtraArg.eResampleAlg GRA_Bilinear; poBand-RasterIO(GF_Read, 0, 0, nXSize, nYSize, pScreenBuf i, nScrW, nScrH, GDT_Byte, 3, nScrW * 3, sExtraArg);双线性插值在缩小时比最近邻平滑很多代价是稍慢。对于显示用途这个开销完全可以接受。我一般会在窗口大小改变时重新执行一次这个降采样读取而不是缓存全分辨率数据再缩放。这样内存占用始终控制在屏幕尺寸级别10000x10000 的影像在 1920x1080 屏幕上只需要约 6MB 缓冲区。还有一个细节如果影像有内置金字塔GDAL 会自动选择最合适的层级来加速读取这时候RasterIO的源区域参数仍然写全图尺寸GDAL 内部会处理层级映射。可以用poBand-GetOverviewCount()查看金字塔层数用GetOverview(i)获取具体层级做手动控制。验证方法很简单打开一张大图在OnPaint里加计时对比全分辨率读取加缩放和直接降采样读取的耗时。我实测过一张 8000x8000 的 GeoTIFF全读再缩放约 800ms直接降采样到 1920x1080 约 120ms差距接近 7 倍。这个习惯我保持了多年显示用途永远让 GDAL 在读取阶段做缩放不要自己读全图再处理。希望帮到你。本文还有配套的精品资源点击获取