Java实现基于GDAL将单波段影像转为三波段影像-唯一值渲染

在处理遥感影像的渲染时,经常需要处理单波段影像。单波段影像没有任何颜色,只有一个波段的值。渲染时只能采用色带拉伸、离散颜色、唯一值渲染这几种方式。直接将单波段影像转成三波段的影像,并将三个波段转为颜色对应的rgb值,这样可以加速渲染、切片的过程。这里我有一张单波段影像,需要按照唯一值的方式,进行渲染,这里记录一下实现过程。

Java依赖

        <dependency><groupId>org.gdal</groupId><artifactId>gdal</artifactId><version>3.5.0</version></dependency>

注意我本地部署的是gdal3以上的版本,不同版本的api会有不同

实现代码

package map.tile.server.tool;import org.gdal.gdal.Band;
import org.gdal.gdal.Dataset;
import org.gdal.gdal.Driver;
import org.gdal.gdal.gdal;
import org.gdal.gdalconst.gdalconstConstants;import java.awt.*;import java.util.HashMap;
import java.util.Map;import static org.gdal.gdalconst.gdalconstConstants.GDT_Byte;public class RasterTool {public static void main(String[] args) throws Exception {// 注册所有支持的格式gdal.AllRegister();String inputFile = "E:\\栅格\\2023_20230430.tif";Dataset dataset = gdal.Open(inputFile, gdalconstConstants.GA_ReadOnly);if (dataset == null) {System.out.println("无法打开输入图像");return;}int width = dataset.GetRasterXSize();int height = dataset.GetRasterYSize();Band band = dataset.GetRasterBand(1);System.out.println(band.GetRasterDataType());//设立每类值对应样式,样色Color colorLevel1 = new Color(0, 92, 230);Color colorLevel2 = new Color(56, 168, 0);Color colorLevel3 = new Color(85, 255, 0);Color colorLevel4 = new Color(255, 170, 0);Map<Short, Color> styles = new HashMap<>(4);styles.put((short) 1, colorLevel1);styles.put((short) 2, colorLevel2);styles.put((short) 3, colorLevel3);styles.put((short) 4, colorLevel4);// 创建新的RGB图像对象Double[] noDataValue = new Double[1];band.GetNoDataValue(noDataValue);Driver driver = gdal.GetDriverByName("GTiff");Dataset rgbDataset = driver.Create("E:\\栅格\\output_rgb_image.tif", width, height, 3, GDT_Byte);rgbDataset.SetProjection(dataset.GetProjection());rgbDataset.SetGeoTransform(dataset.GetGeoTransform());for (int i = 0; i < 3; ++i) {Band outputBand = rgbDataset.GetRasterBand(i + 1);outputBand.SetNoDataValue(noDataValue[0]);short[] buffer = new short[width];//            一次读取一行,防止数据溢出for (int j = 0; j < height; j++) {band.ReadRaster(0, j, width, 1, buffer);// 转换颜色processData(buffer, i + 1, styles);outputBand.WriteRaster(0, j, width, 1, buffer);}outputBand.FlushCache();outputBand.delete();}rgbDataset.FlushCache();rgbDataset.delete();dataset.delete();}private static void processData(short[] buffer, int bandIndex, Map<Short, Color> styles) {if (bandIndex == 1) {for (int i = 0; i < buffer.length; i++) {Color color = styles.get(buffer[i]);if (color != null) {buffer[i] = (short) color.getRed();}}} else if (bandIndex == 2) {for (int i = 0; i < buffer.length; i++) {Color color = styles.get(buffer[i]);if (color != null) {buffer[i] = (short) color.getGreen();}}} else if (bandIndex == 3) {for (int i = 0; i < buffer.length; i++) {Color color = styles.get(buffer[i]);if (color != null) {buffer[i] = (short) color.getBlue();}}}}
}

效果对比

将转换后的文件,直接拖入拖入QGIS展示效果如下:

  • 转换前

转换前

  • 转换后
    转换后

注意,如果使用ArcGIS查看效果的话,需要将拉伸类型选择“无”,将应用Gamma拉伸取消勾选。

ArcGIS配置

感想

  1. GDAL的Java交互真的不太友好
  2. 太大的影像居然没有办法直接读取,整张影像读取,因为可能或超过数组的长度限制。

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

C++ 中关键字 Static

1、什么是static? static 是C中很常用的修饰符&#xff0c;它被用来控制变量的存储方式和可见性。由关键字static修饰类中成员&#xff0c;成为类的静态成员。类的静态成员为其所有对象共享&#xff0c;不管有多少对象&#xff0c;静态成员只有一份存于公用内存中。静态…

2023高级人工智能期末总结

1、人工智能概念的一般描述 人工智能是那些与人的思维相关的活动&#xff0c;诸如决策、问题求解和学习等的自动化&#xff1b; 人工智能是一种计算机能够思维&#xff0c;使机器具有智力的激动人心的新尝试&#xff1b; 人工智能是研究如何让计算机做现阶段只有人才能做得好的…

接口自动化测试要做什么?

作者&#xff1a;不辣的皮皮 链接&#xff1a;https://www.zhihu.com/question/384727359/answer/1124441469 来源&#xff1a;知乎 著作权归作者所有。商业转载请联系作者获得授权&#xff0c;非商业转载请注明出处。 可以分为四个步骤/阶段。 原理 》 业务逻辑》 工具》 …

【现代密码学】笔记3.1-3.3 --规约证明、伪随机性《introduction to modern cryphtography》

【现代密码学】笔记3.1-3.3 --规约证明、伪随机性《introduction to modern cryphtography》 写在最前面私钥加密与伪随机性 第一部分密码学的计算方法论计算安全加密的定义&#xff1a;对称加密算法 伪随机性伪随机生成器&#xff08;PRG&#xff09; 规约法规约证明 构造安全…

Python算法例34 寻找丢失的数

1. 问题描述 给一个由1~n的整数随机组成的一个字符串序列&#xff0c;其中丢失了一个整数&#xff0c;本例将找到它。 2. 问题示例 给出n20&#xff0c;str19201234567891011121314151618&#xff0c;丢失的数是17。 3. 代码实现 def find_missing_number(n, string):nums…

正则表达式手册

文章目录 一、表达式全集二、常见的正则表达式1. 校验数字的表达式2. 校验字符的表达式3. 特殊需求表达式 一、表达式全集 字符描述\将下一个字符标记为一个特殊字符、或一个原义字符、或一个向后引用、或一个八进制转义符。例如&#xff0c;“n”匹配字符“n”。“\n”匹配一…

【代码随想录】刷题笔记Day47

前言 又过了个愉快的周末~大组会终于不用开了&#xff0c;理论上已经可以回家了&#xff01;但是我多留学校几天吧&#xff0c;回家实在太无聊了&#xff0c;也没太多学习的氛围 198. 打家劫舍 - 力扣&#xff08;LeetCode&#xff09; dp[i]含义 考虑下标i&#xff08;包括…

【leetcode】力扣算法之两数相加【中等难度】

题目描述 给你两个 非空 的链表&#xff0c;表示两个非负的整数。它们每位数字都是按照 逆序 的方式存储的&#xff0c;并且每个节点只能存储 一位 数字。 请你将两个数相加&#xff0c;并以相同形式返回一个表示和的链表。 你可以假设除了数字 0 之外&#xff0c;这两个数都…

Qt5.14.2实现将html文件转换为pdf文件

文章目录 简介源码widget.cppwidget.uihtml文件演示效果简介 QPdfWriter是Qt框架中用于创建和写入PDF文件的类。它允许您在您的Qt应用程序中动态生成并输出PDF文档,以便进行打印、保存或导出。 QPdfWriter类提供了以下一些常用的函数和方法,可以让您创建和定制PDF文件: 构…

ASM磁盘管理:从初始化参数到自动化管理的全面解析

文章目录 一、引言二、ASM初始化参数三、ASM三大系统权限四、ASM实例的启停1.Oracle ASM的启停可以通过两种方式进行2.查看集群中的资源状态3.配置 ASM资源随着系统启动而启动4.配置数据库实例随着ASM启动而启动 五、数据库实例与ASM的交互六、 启动策略详解七、 ASM后台进程八…

Hugging face库

1. Transformers 是NLP, CV, audio&#xff0c;speech processing 任务的库。 它包含了超过10000个预训练的模型&#xff0c;涵盖了各种NLP任务&#xff0c;如文本分类&#xff0c;问答&#xff0c;文本生成&#xff0c;情感分析等。Transformers库支持多种深度学习框架&#…

MySQL决战:MySQL数据导入导出

目录 前言 一.navact数据导入导出&#xff08;第三方工具&#xff09; 1.导入数据 2.数据导出 二. mysqldump命令导入导出数据 1.mysqldump介绍 2.数据导出 3.数据导入 三.load data file进行数据导入导出&#xff08;只限于单表&#xff09; 1.数据导出 增加导出权…

python算法每日一练:连续子数组的最大和

这是一道关于动态规划的算法题&#xff1a; 题目描述&#xff1a; 给定一个整数数组 nums&#xff0c;请找出该数组中连续子数组的最大和&#xff0c;并返回这个最大和。 示例&#xff1a; 输入&#xff1a;[-2, 1, -3, 4, -1, 2, 1, -5, 4] 输出&#xff1a;6 解释&#xff…

linux系统基础知识-基础IO

IO 概念引入位图的概念IO的系统调用函数openwriteread()close简单使用样例&#xff1a; 文件描述符fd默认文件流stdin/stdout/stderr文件描述符的分配规则 重定向的概念输出重定向输入重定向追加重定向dup2()系统调用总结 文件缓冲区深入理解缓冲区的概念输出缓冲区部分代码解释…

WEB 3D技术 three.js 线框几何体

本文 我们说一下 线框几何体 想将一个几何体 以线框形式展现 threeJS中 有两种类可以实现 第一种 WireframeGeometry 这种几何体 其实就类似于 将材质中的 wireframe 开启 这种方法 之前我们也用过 还有一种 就是 EdgesGeometry 边缘几何体 我们先将代码写成这样 import .…

【深度学习每日小知识】Data Augmentation 数据增强

数据增强是通过对原始数据进行各种转换和修改来人工生成附加数据的过程&#xff0c;旨在增加机器学习模型中训练数据的大小和多样性。这对于计算机视觉领域尤为重要&#xff0c;因为图像经常被用作输入数据。 计算机视觉中的数据增强 数据增强的主要目标是解决过拟合问题&…

使用迭代优化递归程

王有志&#xff0c;一个分享硬核Java技术的互金摸鱼侠加入Java人的提桶跑路群&#xff1a;共同富裕的Java人 今天我们将会分析上篇文章中递归算法存在的问题&#xff0c;并通过迭代去优化。 递归存在的问题 上一篇中&#xff0c;我们计算了序号10以内的斐波那契数。今天为了清…

【Leetcode】236.二叉树的最近公共祖先

一、题目 1、题目描述 给定一个二叉树, 找到该树中两个指定节点的最近公共祖先。 百度百科中最近公共祖先的定义为:“对于有根树 T 的两个节点 p、q,最近公共祖先表示为一个节点 x,满足 x 是 p、q 的祖先且 x 的深度尽可能大(一个节点也可以是它自己的祖先)。” 示例1…

商品秒杀总结

秒杀模块编写思维总结 分为两种模式一个是限时购买&#xff0c;一个是限裤购买。 我们这里使用的是指定时间段下面的限库购买 单独使用一个库来存储数据&#xff0c;下面有两张表&#xff0c;一张表是具体的商品&#xff0c;一张表是订单。用户下的订单数据都要放到redis 中…

【Linux】journalctl和dmesg日志的区别

journalctl 和 dmesg 是两个在 Linux 系统中查看日志信息的工具&#xff0c;它们提供了不同层次和角度的系统日志。 journalctl&#xff1a; journalctl 是 systemd-journald 服务的一个前端&#xff0c;用于检查和查询系统日志。它使用 systemd 的日志系统&#xff0c;将日志信…