QT GDAL加载显示tif图像

已于 2024-05-11 22:57:17 修改
阅读量115 收藏 1
点赞数 1
文章标签： qt 图像图像
于 2024-05-11 22:49:30 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u012156872/article/details/138732597
版权

1、GDAL库初始化

    // 初始化
    GDALAllRegister();
    // 设置支持中文路径和文件名
    CPLSetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");

在此之前记得引用头文件

#include <gdal_priv.h> 

2、加载tiff图像文件 

void gdalLoadTiffImage()
{
    QString filePath = QFileDialog::getOpenFileName(this,
                                                    QString::fromLocal8Bit("选择图像"),
                                                    "",
                                                    QString::fromLocal8Bit("图像文件(*.tif *.tiff)"));
    if (filePath.isEmpty())
        return;

    GDALDataset *dataset = (GDALDataset *)GDALOpen(filePath.toLocal8Bit().constData(), GA_ReadOnly);
    if (!dataset) {
        // 错误处理
        return;
    }

    int width = dataset->GetRasterXSize();
    int height = dataset->GetRasterYSize();
    int bands = dataset->GetRasterCount();

    // 分配QImage的内存
    QImage::Format format = QImage::Format_RGB888;
    QImage image = QImage(width, height, format);

    // 假设我们总是处理RGB图像，并且波段顺序是红、绿、蓝
    int redBand = 1, greenBand = 2, blueBand = 3;
    if (bands == 1) {
        redBand = greenBand = blueBand = 1;
    }

    // 读取波段数据并转换为QImage
    for (int y = 0; y < height; ++y) {
        for (int x = 0; x < width; ++x) {
            uint8_t red, green, blue;

            GDALRasterBand *bandRed = dataset->GetRasterBand(redBand);
            GDALRasterBand *bandGreen = dataset->GetRasterBand(greenBand);
            GDALRasterBand *bandBlue = dataset->GetRasterBand(blueBand);

            // 像素值范围
            int success = 0;
            double minimum = bandRed->GetMinimum(&success);
            double maximum = bandRed->GetMaximum(&success);
            double noData = bandRed->GetNoDataValue(&success);

            GDALDataType dataType = bandRed->GetRasterDataType();
            switch (dataType) {
                case GDT_Byte:{
                    bandRed->RasterIO(GF_Read, x, y, 1, 1, &red, 1, 1, GDT_Byte, 0, 0);
                    bandGreen->RasterIO(GF_Read, x, y, 1, 1, &green, 1, 1, GDT_Byte, 0, 0);
                    bandBlue->RasterIO(GF_Read, x, y, 1, 1, &blue, 1, 1, GDT_Byte, 0, 0);

                    // 设置QImage的像素值
                    QRgb pixel = qRgb(red, green, blue);
                    image.setPixel(x, y, pixel);
                    break;
                } break;
                case GDT_UInt16: {
                    uint16_t red16, green16, blue16;
                    bandRed->RasterIO(GF_Read, x, y, 1, 1, &red16, 1, 1, GDT_UInt16, 0, 0);
                    bandGreen->RasterIO(GF_Read, x, y, 1, 1, &green16, 1, 1, GDT_UInt16, 0, 0);
                    bandBlue->RasterIO(GF_Read, x, y, 1, 1, &blue16, 1, 1, GDT_UInt16, 0, 0);

                    // 假设值在0-65535之间，缩放到0-255
                    red = static_cast<uint8_t>(red16 * 255.0 / 65535.0);//65535.0 -> (maximum - minimum));
                    green = static_cast<uint8_t>(green16 * 255.0 / 65535.0);
                    blue = static_cast<uint8_t>(blue16 * 255.0 / 65535.0);

                    // 设置QImage的像素值
                    QRgb pixel = qRgb(red, green, blue);
                    image.setPixel(x, y, pixel);
                    break;
                }
                case GDT_Int16: {
                    int16_t red16, green16, blue16;
                    bandRed->RasterIO(GF_Read, x, y, 1, 1, &red16, 1, 1, GDT_Int16, 0, 0);
                    bandGreen->RasterIO(GF_Read, x, y, 1, 1, &green16, 1, 1, GDT_Int16, 0, 0);
                    bandBlue->RasterIO(GF_Read, x, y, 1, 1, &blue16, 1, 1, GDT_Int16, 0, 0);

                    // 假设值在-32768到32767之间，先转换为0-65535再缩放到0-255
                    // 注意这里需要处理负数，并可能需要进行偏移和缩放
                    uint16_t redUnsigned = static_cast<uint16_t>(red16 + 32768);
                    uint16_t greenUnsigned = static_cast<uint16_t>(green16 + 32768);
                    uint16_t blueUnsigned = static_cast<uint16_t>(blue16 + 32768);

                    red = static_cast<uint8_t>(redUnsigned * 255.0 / 65535.0);
                    green = static_cast<uint8_t>(greenUnsigned * 255.0 / 65535.0);
                    blue = static_cast<uint8_t>(blueUnsigned * 255.0 / 65535.0);

                    // 设置QImage的像素值
                    QRgb pixel = qRgb(red, green, blue);
                    image.setPixel(x, y, pixel);
                    break;
                }
                case GDT_UInt32:{
                    uint32_t red32, green32, blue32;
                    bandRed->RasterIO(GF_Read, x, y, 1, 1, &red32, 1, 1, GDT_UInt32, 0, 0);
                    bandGreen->RasterIO(GF_Read, x, y, 1, 1, &green32, 1, 1, GDT_UInt32, 0, 0);
                    bandBlue->RasterIO(GF_Read, x, y, 1, 1, &blue32, 1, 1, GDT_UInt32, 0, 0);

                    // 假设值在0-4294967295之间，缩放到0-255
                    red = static_cast<uint8_t>(red32 * 255.0 / 4294967295.0);
                    green = static_cast<uint8_t>(green32 * 255.0 / 4294967295.0);
                    blue = static_cast<uint8_t>(blue32 * 255.0 / 4294967295.0);

                    // 设置QImage的像素值
                    QRgb pixel = qRgb(red, green, blue);
                    image.setPixel(x, y, pixel);
                    break;
                }
                case GDT_Int32:{
                    int32_t red32, green32, blue32;
                    bandRed->RasterIO(GF_Read, x, y, 1, 1, &red32, 1, 1, GDT_Int32, 0, 0);
                    bandGreen->RasterIO(GF_Read, x, y, 1, 1, &green32, 1, 1, GDT_Int32, 0, 0);
                    bandBlue->RasterIO(GF_Read, x, y, 1, 1, &blue32, 1, 1, GDT_Int32, 0, 0);

                    // 假设值在-2147483648到2147483647之间，先转换为0-4294967295再缩放到0-255
                    uint32_t redUnsigned = static_cast<uint32_t>(red32 + 2147483648);
                    uint32_t greenUnsigned = static_cast<uint32_t>(green32 + 2147483648);
                    uint32_t blueUnsigned = static_cast<uint32_t>(blue32 + 2147483648);

                    red = static_cast<uint8_t>(redUnsigned * 255.0 / 4294967295.0);
                    green = static_cast<uint8_t>(greenUnsigned * 255.0 / 4294967295.0);
                    blue = static_cast<uint8_t>(blueUnsigned * 255.0 / 4294967295.0);

                    // 设置QImage的像素值
                    QRgb pixel = qRgb(red, green, blue);
                    image.setPixel(x, y, pixel);
                    break;
                }
                case GDT_Float32:{
                    float redFloat, greenFloat, blueFloat;
                    bandRed->RasterIO(GF_Read, x, y, 1, 1, &redFloat, 1, 1, GDT_Float32, 0, 0);
                    bandGreen->RasterIO(GF_Read, x, y, 1, 1, &greenFloat, 1, 1, GDT_Float32, 0, 0);
                    bandBlue->RasterIO(GF_Read, x, y, 1, 1, &blueFloat, 1, 1, GDT_Float32, 0, 0);

                    // 假设值在0.0到1.0之间，直接乘以255
                    // 如果不是0.0到1.0，你需要知道数据的具体范围并进行适当的缩放
                    red = static_cast<uint8_t>(redFloat * 255.0);
                    green = static_cast<uint8_t>(greenFloat * 255.0);
                    blue = static_cast<uint8_t>(blueFloat * 255.0);

                    // 如果需要处理可能的超出范围的值
                    red = qMin(qMax((uint8_t)red, (uint8_t)0), (uint8_t)255);
                    green = std::min(std::max(green, (uint8_t)0), (uint8_t)255);
                    blue = std::min(std::max(blue, (uint8_t)0), (uint8_t)255);

                    // 设置QImage的像素值
                    QRgb pixel = qRgb(red, green, blue);
                    image.setPixel(x, y, pixel);
                    break;
                }
                case GDT_Float64:{
                    double redDouble, greenDouble, blueDouble;
                    bandRed->RasterIO(GF_Read, x, y, 1, 1, &redDouble, 1, 1, GDT_Float64, 0, 0);
                    bandGreen->RasterIO(GF_Read, x, y, 1, 1, &greenDouble, 1, 1, GDT_Float64, 0, 0);
                    bandBlue->RasterIO(GF_Read, x, y, 1, 1, &blueDouble, 1, 1, GDT_Float64, 0, 0);

                    // 假设值在0.0到1.0之间，直接乘以255
                    // 如果不是0.0到1.0，你需要知道数据的具体范围并进行适当的缩放
                    red = static_cast<uint8_t>(redDouble * 255.0);
                    green = static_cast<uint8_t>(greenDouble * 255.0);
                    blue = static_cast<uint8_t>(blueDouble * 255.0);

                    // 如果需要处理可能的超出范围的值
                    red = std::min(std::max(red, (uint8_t)0), (uint8_t)255);
                    green = std::min(std::max(green, (uint8_t)0), (uint8_t)255);
                    blue = std::min(std::max(blue, (uint8_t)0), (uint8_t)255);

                    // 设置QImage的像素值
                    QRgb pixel = qRgb(red, green, blue);
                    image.setPixel(x, y, pixel);
                    break;
                }

                // 对于复数类型（GDT_CInt16, GDT_CInt32, GDT_CFloat32, GDT_CFloat64），
                // GDAL通常不直接用于存储图像数据（如RGB），而是用于科学计算。
                case GDT_CInt16:{
                    int16_t redComplexValue[2], greenComplexValue[2], blueComplexValue[2]; // 数组中的第一个元素是实部，第二个元素是虚部
                    bandRed->RasterIO(GF_Read, x, y, 1, 1, redComplexValue, 1, 1, GDT_CInt16, 0, 0);
                    bandGreen->RasterIO(GF_Read, x, y, 1, 1, greenComplexValue, 1, 1, GDT_CInt16, 0, 0);
                    bandBlue->RasterIO(GF_Read, x, y, 1, 1, blueComplexValue, 1, 1, GDT_CInt16, 0, 0);

                    // 只使用实部，并像GDT_Int16那样处理
                    int16_t red16 = redComplexValue[0];
                    int16_t green16 = greenComplexValue[0];
                    int16_t blue16 = blueComplexValue[0];

                    // 假设值在-32768到32767之间，先转换为0-65535再缩放到0-255
                    // 注意这里需要处理负数，并可能需要进行偏移和缩放
                    uint16_t redUnsigned = static_cast<uint16_t>(red16 + 32768);
                    uint16_t greenUnsigned = static_cast<uint16_t>(green16 + 32768);
                    uint16_t blueUnsigned = static_cast<uint16_t>(blue16 + 32768);

                    red = static_cast<uint8_t>(redUnsigned * 255.0 / 65535.0);
                    green = static_cast<uint8_t>(greenUnsigned * 255.0 / 65535.0);
                    blue = static_cast<uint8_t>(blueUnsigned * 255.0 / 65535.0);

                    // 设置QImage的像素值
                    QRgb pixel = qRgb(red, green, blue);
                    image.setPixel(x, y, pixel);

                    break;
                }
                case GDT_CInt32:{
                    // 参考GDT_CInt16
                    break;
                }
                case GDT_CFloat32:{
                    // 参考GDT_CInt16
                    break;
                }
                case GDT_CFloat64:{
                    // 参考GDT_CInt16
                    break;
                }
            }
        }
    }

    GDALClose(dataset);
}

3、GDAL库资源释放

GDALDestroyDriverManager();

在QT安装路径下有个tif文件目录可供测试使用：C:\Qt\Qt5.15.2\5.15.2\Src\qtimageformats\tests\shared\images\tiff

按照常理来说，如果波段为1，数据类型是GDT_Byte，数据范围正常来说应该是0~255，最终显示出来的图像是灰度图，但是如下图片的数据范围却为0~65535，用自己代码显示当然也是灰度图像，但是用其它工具软件显示确实彩色图像（如PS），究竟是哪里出了问题，有待日后去研究了！！！

老狼主
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
博客
TightVNC H264编解码（一）
12-31 1488
经过尽二个月的不断摸索，TightVNC终于支持H264编解码了，前期真正编写H264编解码器只花了一周左右时间，但是测试发现效果并不是太理想（帧率太低，延迟大），这大概是所有远程控制软件要处理的通病吧。作为一位经验丰富的程序员，造成这种问题的原因，有以下几点：1、抓屏？2、H264编码？3、TCP网络发送？4、H264解码？5、视频渲染？为了排除到底哪个环节出了问题，看来不摸透底层事务流程是行不通的，不得已摸透代码画出了下面这个简陋的流程图。Server端事务处理流程图.
博客
QT HttpServer
07-31 823
资源从网上收集，包含httpserver、sslserver、http-parser三个模块，解决了编译过程中由于路径引用错误问题。下载链接：https://download.csdn.net/download/u012156872/20664799，下载编译即可应用到自己项目中，简单快速，不用造车。...
博客
QT 界面设计篇（水波纹进度条QProgressBarWater）
11-30 1637
效果图：头文件声明（QProgressBarWater.h）：#ifndef QPROGRESSBARWATER_H#define QPROGRESSBARWATER_H#include <QWidget>#include <QPainter>#include <QPaintEvent>class QProgressBarWater : public QWidget{ Q_OBJECTpublic: explic.
博客
EasyPR编译指南
08-02 3188
1、下载源码https://github.com/liuruoze/EasyPRDownloadZIP或者git（https://github.com/liuruoze/EasyPR.git）将代码解压到硬盘目录下。例如F:\车牌识别\EasyPR（车牌识别）\EasyPR2、双击解决方案文件EasyPR.sln，用VSIDE（这里使用VS2019）打开，提示点击确定。3、选择libeasypr项目，编译出现如下错误这是因为没有配置OpenCV环境所致。4、...
博客
QT MySQL驱动
05-14 111
在某些场合QT自带的QSQLDatabase驱动有效，但是open的时候，程序直接崩溃了，原因未知。不得已将上次使用过的。根据QT接口进行改动，进行无缝衔接。
博客
QT vlc音视频开发踩坑
05-11 175
为此我检查错误耗时将近1整天时间，甚至于都一度怀疑人生（之前在VC++下开发一切顺利）。SDK初始化一切顺利，最终在调用如下代码的时候始终失败，返回NULL值。问题解决，原来如此简单啊！
博客
QT 在线安装失败
05-07 115
QT在线安装使用的国外服务器，所以导致安装极其缓慢或者失败，唯一能想到的解决方法就是切换到国内的镜像源（除非翻墙）。
博客
QT 无标题栏窗口的拖动和拉伸
12-14 470
如下2种方案，2种方案都不完美（拉伸窗口都会抖动），但是如果一定要使用，强烈推荐第1种方案（可减小拉伸时窗口抖动）
博客
C++ libVNC开源库
12-11 520
VNC具体用途，不用过多介绍了。如果需要实时更新服务器端画面，while循环就可以了。
博客
QT QCefView自编译
12-09 511
6、回到QCefView目录，执行批处理generate-win-x86_64.bat即可，出现如下信息表示生成成功。5、设置QT版本，修改配置文件QCefView/cmake/QtConfig.cmake，如下。4、请根据对应的版本修改编译文件CefViewCore/CMakeLists.txt，如下。2、下载QCefViewCore，放到QCefViewCore目录下，地址。请注意最新版本兼容性可能存在问题，我下载的版本是。1、下载QCefView，地址。3、下载Cef库文件，地址。
博客
Ubuntu 配置共享文件夹
10-26 257
安装虚拟机经常会遇到已经共享了文件夹，但是就是找不到目录文件，搜集网上资料，汇总一下指令，留作以后备用吧：vmware-hgfsclient在文件最后一行添加如下内容：
博客
C++ OleDB数据库驱动（MySql驱动姊妹篇）
10-11 221
这里只介绍OelDB调用数据的方法。OleDB因为兼容ODBC和ADO，所以受大众追捧。大家也都非常熟悉了，这里就直接上代码吧。
博客
Visual Assist X插件安装失败总结
09-27 932
2、安装成功以后，却显示插件无法启用，或者是插件启动了，但是实则没有启动，看不到插件选项。1、如果是运行VA_X_SETUP.exe无反应，无法启动，请右键属性，修改兼容模式。3、此时重启应该就正常了。
博客
常用时间日期格式符号定义
09-27 1773
编程过程中经常会用到时间日期格式化，这里是一些常见的日期时间格式字符定义的列表：年份：yyyyyy月份：MMMMMM天数：ddd小时：hhh分钟：mmm秒数：sss毫秒：zzzz星期几：dddddddAP时区：t。
博客
C++ linux下使用X11实现屏幕截图（续）
09-18 241
使用libpng库将XGetImage获取到的内容保存为PNG格式的文件。
博客
C++ MFC中如何根据不同文件编码格式读取和显示文件内容
09-17 741
这里假设VC++的字符集编码格式是Unicode，要显示的文本内容是strResult。上一篇，介绍了写文件，这一篇我要介绍的如何读取文件内容并显示出来。3、读取unicode-le编码格式的文本文件。4、读取unicode-be编码格式的文本文件。2、读取utf-8编码格式的文本文件。1、读取ANSI编码格式的文本文件。
博客
C++ MFC中如何根据不同文件编码格式创建和写文本内容
09-17 499
这里假设VC++的字符集编码格式是Unicode，要写入的文本内容是m_strTEXT。3、保存unicode-le编码格式的文本文件。4、保存unicode-be编码格式的文本文件。2、保存utf-8编码格式的文本文件。1、保存ANSI编码格式的文本文件。
博客
C++ COM接口调用
09-14 333
那如果这个dll更新特频繁，每次都需要反注册，再重新注册，似乎有那么一点点的烦人哦。既然COM组件也是一种库文件，那有没有一种调用方法，像调用动态链接库那样先加载，再寻址，最后调用呢？代码虽然复杂很多，但是目的我们是达到了，庆祝一下吧！前提条件：COM是VSC++编写的，先引用头文件。
博客
C++ OSD水印叠加
09-06 440
视频相关领域（监控、LED大屏）很多场合可能会涉及到在一幅画面上进行文字或时间的叠加，常规做法都是在后端实现，即先渲染图像，然后叠加OSD文字信息，这种方法简单且高效。但是有些场合必须要求在前端进行叠加，比如监控领域中视频名称信息、时间信息等这些都需要在前端图像编码的时候就已经叠加到图像上，以防止用户修改和验伪。示例中采用jpeg图像作为底图，所以需要借助libjpeg将图像解码还原RGB原始图像数据，叠加OSD文字信息需借助freetype生成点阵RGB图像数据，然后对2个图像进行像素合并或重写。
博客
C++ MySql驱动
09-06 270
5. 扩展性：MySQL支持水平扩展和垂直扩展，可以方便地增加或减少服务器的数量、内存和存储容量，以满足不断增长的需求。8. 社区支持：MySQL拥有庞大的开源社区，用户可以通过社区获取帮助、分享经验和参与开发，获得及时的技术支持。6. 兼容性：MySQL遵循SQL标准，与其他常见的数据库系统兼容，能够无缝迁移和集成现有的数据和应用程序。4. 高性能：MySQL采用了多种优化策略和技术，如索引、缓存、分区等，以提供快速的读写操作和响应时间。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值