C++语言GDAL批量裁剪多波段栅格图像：基于像元个数裁剪

news2025/1/11 14:27:34

  本文介绍基于C++ 语言的GDAL模块，按照给定的像元行数与列数，批量裁剪大量多波段栅格遥感影像文件，并将所得到的裁剪后新的多波段遥感影像文件保存在指定路径中的方法。

  在之前的文章中，我们多次介绍了在不同平台，或基于不同代码语言，对栅格遥感影像加以裁剪、批量裁剪的方法，主要包括Python中ArcPy基于矢量范围批量裁剪大量栅格遥感影像（https://blog.csdn.net/zhebushibiaoshifu/article/details/128307676）、Google Earth Engine谷歌地球引擎GEE矢量数据裁剪栅格数据（https://blog.csdn.net/zhebushibiaoshifu/article/details/117390431）、基于ENVI实现栅格遥感影像按图层行列号与像元数量划定矩形研究区域并裁剪（https://blog.csdn.net/zhebushibiaoshifu/article/details/118978851）等；而本文，我们就介绍一下基于**C++**语言的GDAL模块，实现批量裁剪需求的方法。

  首先，来看一下本文的需求。现在有一个文件夹，如下图所示，其中具有多个.tiff格式的多波段遥感影像文件（为了方便，我们这里文件夹内就只有2个文件，但实际上一般我们批量处理的需求肯定远远大于这个数量）。

  我们希望实现的，就是基于这个文件夹内每一景遥感影像，将其左上角100 * 100像元的这一部分给裁剪下来（如下图所示），并分别保存为新的遥感影像文件（其中，新的文件名称就在原有文件名称后加一个_C后缀即可），并存放在另一个指定的结果文件夹中。我们希望裁剪后的遥感影像，和原有的遥感影像对比起来，呈现如下图所示的情况。

  本文所用代码如下。

#include <iostream>
#include <string>
#include <gdal/gdal.h>
#include <gdal/gdal_priv.h>
using namespace std;

int main()
{
    GDALAllRegister();

    string inputFolder = "/home/cppGDAL/TIF/WFV1_2021_STB/test";
    string outputFolder = "/home/cppGDAL/TIF/WFV1_2021_STB_C";

    CPLStringList fileList;
    fileList = VSIReadDir(inputFolder.c_str());

    for (int i = 0; i < fileList.size(); i++)
    {
        string inputImagePath = fileList[i];
        if (!EQUAL(CPLGetExtension(inputImagePath.c_str()), "tiff"))
        {
            continue;
        }

        string full_path = inputFolder + "/" + inputImagePath;
        GDALDataset *poDataset = (GDALDataset *)GDALOpen(full_path.c_str(), GA_ReadOnly);

        int width = poDataset->GetRasterXSize();
        int height = poDataset->GetRasterYSize();

        int xOffset = 0;
        int yOffset = 0;
        int xSize = 100;
        int ySize = 100;
        cout << full_path;
        string outputImageName = (outputFolder + "/" + inputImagePath.substr(0, inputImagePath.size() - 5) + "_C.tiff");
        cout << outputImageName;
        GDALDriver *poDriver = GetGDALDriverManager()->GetDriverByName("GTiff");
        GDALDataset *poOutputDataset = poDriver->Create(outputImageName.c_str(), xSize, ySize, 4, GDT_Float32, NULL);

        poOutputDataset->SetProjection(poDataset->GetProjectionRef());
        double adfGeoTransform[6];
        poDataset->GetGeoTransform(adfGeoTransform);
        // adfGeoTransform[1] = adfGeoTransform[1] / (width / xSize);
        // adfGeoTransform[5] = adfGeoTransform[5] / (height / ySize);
        poOutputDataset->SetGeoTransform(adfGeoTransform);

        for (int bandIndex = 1; bandIndex <= 4; bandIndex++)
        {
            float *buffer = new float[xSize * ySize];
            GDALRasterBand *poBand = poDataset->GetRasterBand(bandIndex);
            GDALRasterBand *poOutputBand = poOutputDataset->GetRasterBand(bandIndex);
            poBand->RasterIO(GF_Read, xOffset, yOffset, xSize, ySize, buffer, xSize, ySize, GDT_Float32, 0, 0, NULL);
            poOutputBand->RasterIO(GF_Write, 0, 0, xSize, ySize, buffer, xSize, ySize, GDT_Float32, 0, 0, NULL);
            delete[] buffer;
        }

        GDALClose(poOutputDataset);
        GDALClose(poDataset);
    }

    GDALDestroyDriverManager();

    return 0;
}

  我们来介绍一下上述代码的具体内容。

  首先，我们需要通过GDALAllRegister();，来注册所有的GDAL驱动器。同时，我们定义了输入和输出文件夹路径——inputFolder就是存储输入遥感影像（待裁剪的遥感影像）的文件夹路径，outputFolder则是存储结果遥感影像的文件夹路径。

  其次，我们通过CPLStringList fileList;定义一个字符串列表，用于存储文件夹中的文件列表；并使用VSIReadDir函数读取输入文件夹中的所有文件，并将结果存储在fileList中。

  接下来，我们使用循环迭代处理每个文件。首先，通过string inputImagePath = fileList[i];获取当前文件的路径；如果文件的扩展名不是tiff，则跳过该文件。接下来，对于文件的扩展名是tiff的，我们构建完整的输入文件路径，并使用GDALOpen函数打开输入文件，返回一个GDALDataset对象，存储在poDataset中。

  接下来，我们即可获取输入文件的宽度和高度，并定义裁剪区域的偏移量（左上角像元的位置）、宽度和高度。前面提到了，我这里就是需要在原本遥感影像的最左上角（所以偏移量均为0），裁剪下来100 * 100像元的这一部分。其次，构建输出文件的路径，并使用GetGDALDriverManager()->GetDriverByName函数获取GTiff驱动器对象，存储在poDriver中。随后，我们使用poDriver->Create函数创建输出文件，返回一个GDALDataset对象，存储在poOutputDataset中。

  接下来这个部分需要稍微注意一下。首先，我们使用poOutputDataset->SetProjection设置输出文件的投影信息，即与输入文件相同的投影；其次，使用poDataset->GetGeoTransform获取输入文件的地理变换参数，存储在adfGeoTransform数组中。由于在我这里，裁剪后遥感影像的像元大小（即单个像元的长度与宽度）没有改变，且裁剪前后栅格遥感影像的左上角像元没有发生变化，所以新的栅格遥感影像的地理变换参数和老的栅格遥感影像比起来，无需有任何改变；但是如果大家的裁剪需求不是这样的话（比如像元大小发生变化了，或者是裁剪并不是从左上角像元开始的），那么就需要调整这6个地理变换参数——至于怎么变，这就比较复杂了，我也还没完全搞清楚，大家就结合自己的实际需求，到GDAL官网查阅即可。最后，我们使用poOutputDataset->SetGeoTransform，设置输出文件的地理变换参数，在我这里就是与输入文件完全相同的地理变换参数。

  代码最后，我们使用循环迭代处理每个波段（我这里每一个遥感影像都是共4个波段）。首先，创建一个大小为xSize * ySize的浮点型缓冲区，并使用poBand->RasterIO从输入文件中读取对应波段的像元数据到缓冲区；接下来，使用poOutputBand->RasterIO将缓冲区中的数据写入到输出文件对应波段中。随后，即可释放缓冲区内存，并关闭输出文件和输入文件。

  运行上述代码，我们即可在结果文件夹中看到已经裁剪好的遥感影像文件，且新的文件的文件名称也符合我们的要求；如下图所示。

  至此，大功告成。

欢迎关注：疯狂学习GIS

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.coloradmin.cn/o/1058659.html

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈，一经查实，立即删除！

相关文章

TouchGFX之后端通信

TouchGFX之后端通信

在大多数应用中&#xff0c;UI需以某种方式连接到系统的其余部分&#xff0c;并发送和接收数据。 它可能会与硬件外设&#xff08;传感器数据、模数转换和串行通信等&#xff09;或其他软件模块进行交互通讯。 Model类​ 所有TouchGFX应用都有Model类&#xff0c;Model类除了存…
阅读更多...
色彩一致性自动处理方法在遥感图像中的应用

色彩一致性自动处理方法在遥感图像中的应用

前言 在获取卫星遥感影像时&#xff0c;由于受不均匀的光照、不同的大气条件和不同的传感器设备等因素的影响&#xff0c;遥感影像中会存在局部亮度和色彩分布不均匀的现象&#xff0c;下面是在BigMap地图下载器中收集的几幅谷歌卫星影像&#xff0c;像下面这种都是拼接好的影像…
阅读更多...
从零开始的C++（四）

从零开始的C++（四）

上篇链接&#xff1a;http://t.csdnimg.cn/3nyT9 1.拷贝构造函数&#xff1a; 上篇中介绍了析构函数&#xff0c;即在对象销毁时自动调用的函数&#xff0c;常用于含有malloc、fopen等成员变量的对象。然而&#xff0c;在将对象做函数实参进行值传递的时候&#xff0c;可能会…
阅读更多...
Unity宣布自2024年起将根据游戏安装量收费，你对此有何看法？

Unity宣布自2024年起将根据游戏安装量收费，你对此有何看法？

文章目录 每日一句正能量前言Unity的来历Unity的应用对于收费的看法个人角度&#xff1a;公司角度&#xff1a; 后记 每日一句正能量 水与水之间有距离&#xff0c;但地心下直相牵&#xff0c;人与人之间有距离&#xff0c;但心里时刻挂念&#xff0c;发条短信道声晚安&#xf…
阅读更多...
【强化学习】05 —— 基于无模型的强化学习（Prediction）

【强化学习】05 —— 基于无模型的强化学习（Prediction）

文章目录 简介蒙特卡洛算法时序差分方法Example1 MC和TD的对比偏差&#xff08;Bias&#xff09;/方差&#xff08;Variance&#xff09;的权衡Example2 Random WalkExample3 AB 反向传播(backup)Monte-Carlo BackupTemporal-Difference BackupDynamic Programming Backup Boot…
阅读更多...
深入浅出，SpringBoot整合Quartz实现定时任务与Redis健康检测（一）

深入浅出，SpringBoot整合Quartz实现定时任务与Redis健康检测（一）

目录 前言 环境配置 Quartz 什么是Quartz&#xff1f; 应用场景 核心组件 Job JobDetail Trigger CronTrigger SimpleTrigger Scheduler 任务存储 RAM JDBC 导入依赖 定时任务 销量统计 Redis检测 使用 ​编辑 注意事项 前言 在悦享校园1.0中引入了Quart…
阅读更多...
番外3：下载+安装VMware(前期准备）

番外3：下载+安装VMware(前期准备）

step1: 查看自己笔记本电脑配置&#xff1b; step2: 下载并安装VMware&#xff08;下载地址www..kkx.net/soft/16841.html&#xff09;这里选择本地普通下载&#xff1b; step3: 安装VMware过程中需要填写密钥&#xff08;本人用的最后一个&#xff09;; #UU54R-FVD91-488PP-7N…
阅读更多...
国庆day4

国庆day4

运算符重载代码 #include <iostream> using namespace std; class Num { private:int num1; //实部int num2; //虚部 public:Num(){}; //无参构造Num(int n1,int n2):num1(n1),num2(n2){}; //有参构造~Num(){}; //析构函数const Num operator(const Num &other)cons…
阅读更多...
微服务技术栈-Ribbon负载均衡和Nacos注册中心

微服务技术栈-Ribbon负载均衡和Nacos注册中心

文章目录 前言一、Ribbon负载均衡1.LoadBalancerInterceptor&#xff08;负载均衡拦截器&#xff09;2.负载均衡策略IRule 二、Nacos注册中心1.Nacos简介2.搭建Nacos注册中心3.服务分级存储模型4.环境隔离5.Nacos与Eureka的区别 总结 前言 在上面那个文章中介绍了微服务架构的…
阅读更多...
【C语言】函数的定义、传参与调用（一）

【C语言】函数的定义、传参与调用（一）

目录 导读&#xff1a; 1. 为什么要用函数 2. C语言中函数的分类 2.1 库函数 2.1.1 什么是库函数 2.1.2 C语言常用的库函数 2.2 自定义函数 2.2.1 什么是自定义函数 2.2.2 定义函数的方法 2.2.3 举例 3. 函数的参数 3.1 传参不同的对比 3.2 形式参数&#xff08;形…
阅读更多...
以太网基础学习（一）——以太网概述

以太网基础学习（一）——以太网概述

一、以太网概述 以太网(Ethernet)指的是由 Xerox公司创建并由Xerox、Intel和 DEC公司联合开发的基带局域网规范&#xff0c;通用的以太网标准于1980年9月30日出台&#xff0c;是当今现有局域网采用的最通用的通信协议标准&#xff08;是局域网的一种&#xff09;。 以太网是一种…
阅读更多...
libevent源码学习笔记

libevent源码学习笔记

libevent源码学习笔记 libevent安装libevent源码解析&#xff08;1&#xff09;事件对象&#xff08;2&#xff09;事件操作&#xff08;3&#xff09;事件循环&#xff08;4&#xff09;事件处理 常用指令问题记录问题一&#xff1a;长连接的管理问题二&#xff1a;连接关闭问…
阅读更多...
WebSocket实战之三遇上PAC

WebSocket实战之三遇上PAC

一、前言 前两天销售数据实时刷新功能开发测试完成&#xff0c;开开心心部署到生产环境&#xff0c;然后直接懵逼傻眼了&#xff0c;竟然连接不上WebSocket服务端&#xff0c;浏览器端请求头报 Provisional headers are shown 信息&#xff0c;然后采用一系列操作排查问题。 …
阅读更多...
DS线性表之链表

DS线性表之链表

前言 我们上一期介绍了顺序表&#xff0c;它的底层就是数组&#xff0c;我们也分别对顺序表的动态版本和静态版本进行了实现&#xff01;并且分析了顺序表的优缺点&#xff0c;优点是&#xff1a;尾插、尾删效率很高&#xff0c;其时间复杂度是O(1)&#xff1b;缺点是&#xff…
阅读更多...
用于数据增强的十个Python库

用于数据增强的十个Python库

数据增强是人工智能和机器学习领域的一项关键技术。它涉及到创建现有数据集的变体&#xff0c;提高模型性能和泛化。Python是一种流行的AI和ML语言&#xff0c;它提供了几个强大的数据增强库。在本文中&#xff0c;我们将介绍数据增强的十个Python库&#xff0c;并为每个库提供…
阅读更多...
（二）激光线扫描-相机标定

（二）激光线扫描-相机标定

1. 何为相机标定？ 当相机拍摄照片时，我们看到的图像通常与我们实际看到的不完全相同。这是由相机镜头引起的，而且发生的频率比我们想象的要高。 这种图像的改变就是我们所说的畸变。一般来说，畸变是指直线在图像中出现弯曲或弯曲。 这种畸变我们可以通过相机标定来进行解…
阅读更多...
轮询与中断

轮询与中断

中断控制器 #include"exynos_4412.h"int main() {/*产生一个中断信号*//*1.属于外设层次&#xff0c;让外部的硬件控制器能产生一个中断信号并发送给中断控制器*//*将GPX1_1设置成中断功能*/GPX1.CON GPX1.CON |( 0xF << 4);/*设置GPX1_1中断的触发方式---下降…
阅读更多...
栈的应用场景(一)

栈的应用场景(一)

逆波兰表达式 1.题目2.思路3.代码 1.题目 2.思路 3.代码 class Solution {public int evalRPN(String[] tokens) {//创建一个栈Stack<Integer> stack new Stack<>();//对字符串数组进行遍历for(String x : tokens){//分数操作符和操作数两种情况,我们得判断//一下…
阅读更多...
JavaScript系列从入门到精通系列第十二篇：JavaScript中对象的简介和对象的基本操作以及JavaScript中的属性值和属性名

JavaScript系列从入门到精通系列第十二篇：JavaScript中对象的简介和对象的基本操作以及JavaScript中的属性值和属性名

文章目录 前言 一&#xff1a;对象分类 1&#xff1a;内建对象 2&#xff1a;宿主对象 3&#xff1a;自建对象 二&#xff1a;对象的基本操作 1&#xff1a;创建对象 2&#xff1a;向对象中添加属性 3&#xff1a;读取对象中的属性 4&#xff1a;修改对象中的属性 三…
阅读更多...
使用晶体管做布尔逻辑和逻辑门

使用晶体管做布尔逻辑和逻辑门

目录 二进制&#xff0c;三进制&#xff0c;五进制 true&#xff0c;false表示0&#xff0c;1 早期计算机采用进制 布尔逻辑 三个基本操作&#xff1a;NOT,AND,OR 基础“真值表” NOT 如何实现&#xff1f; AND如何实现&#xff1f; OR如何实现&#xff1f; 图标表示…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023