GDAL实现标准差拉伸渲染影像

2024-04-12 09:18

本文主要是介绍GDAL实现标准差拉伸渲染影像,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

import os
from sys import path
from osgeo import gdal
from osgeo import osr
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.colors as mpc
import argparsedef stdStretch(imgFile, colorList, flag, outPath):# imgFile 影像# colorList 色带参数# flag=True 将数据投影到3857坐标系进行显示# flag=False 将数据按原坐标系进行显示gdal.AllRegister()ds = gdal.Open(imgFile)if flag:dsrs = osr.SpatialReference()dsrs.ImportFromEPSG(3857)vrt_ds = gdal.AutoCreateWarpedVRT(ds, None, dsrs.ExportToWkt(), gdal.GRA_Bilinear)band = gdal.Dataset.GetRasterBand(vrt_ds, 1)else:band = gdal.Dataset.GetRasterBand(ds, 1)array = gdal.Band.ReadAsArray(band)avg = np.nanmean(array)std = np.nanstd(array)# 默认按2.5倍标准差进行拉伸v_max = avg+2.5*stdv_min = avg-2.5*stdnorm = plt.Normalize(v_min, v_max)colors = colorListcm = mpc.LinearSegmentedColormap.from_list('colormap', colors, 512)plt.imshow(array, cm, norm)plt.axis('off')plt.subplots_adjust(bottom=0, top=1, left=0, right=1, hspace=0, wspace=0)plt.colorbar()outname = os.path.basename(imgFile)outname = outname.replace(outname[outname.rindex('.'):], '.png')outname = os.path.join(outPath, outname)plt.savefig(outname, transparent=True, dpi=100)plt.close()# 封装为命令行工具 批量处理目录下的影像
if __name__ == '__main__':parser = argparse.ArgumentParser('标准差拉伸渲染影像')parser.add_argument('inputPath', type=str, help='The directory  stores the images or a image file to process!')parser.add_argument('colorList', type=str,help='Color map used in render!')parser.add_argument('outPath', type=str,help='The directory stores the result!')parser.add_argument('flag', type=bool, help='A flag which indicates whether the image should be prjected to epsg3857!')parser.add_argument('extension',type=str,help='The extension of image file to be processed!',default=None)args=parser.parse_args()inputPath=args.inputPathcolors=args.colorList.split(',')outPath=args.outPathflag=args.flagextension=args.extension# extension=None 则处理目录下的所有影像if os.path.isfile(inputPath):# 处理单张影像stdStretch(inputPath,colors,flag,outPath)else:# 处理整个目录imgFiles=os.listdir(inputPath)for imgFile in imgFiles:if not(extension) or imgFile.endswith(extension):stdStretch(os.path.join(inputPath,imgFile),colors,flag,outPath)# colors = '#0000FF, #3B7FFF, #00FFFF,#B4FF91, #FFFF00, #FF9100, #FF0000'print('finished!')
调用方式
user@DESKTOP-IFIB8V0 MINGW64 /i/毕业设计整理/LST
$ python color.py --help
usage: 标准差拉伸渲染影像 [-h] inputPath colorList outPath flag extensionpositional arguments:inputPath   The directory stores the images or a image file to process!colorList   Color map used in render!outPath     The directory stores the result!flag        A flag which indicates whether the image should be prjected to epsg3857!extension   The extension of image file to be processed!optional arguments:-h, --help  show this help message and exit
(base) 
渲染结果

在这里插入图片描述
在这里插入图片描述

这篇关于GDAL实现标准差拉伸渲染影像的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/896720

相关文章

mybatis执行insert返回id实现详解

《mybatis执行insert返回id实现详解》MyBatis插入操作默认返回受影响行数,需通过useGeneratedKeys+keyProperty或selectKey获取主键ID,确保主键为自... 目录 两种方式获取自增 ID:1. ​​useGeneratedKeys+keyProperty(推

Spring Boot集成Druid实现数据源管理与监控的详细步骤

《SpringBoot集成Druid实现数据源管理与监控的详细步骤》本文介绍如何在SpringBoot项目中集成Druid数据库连接池,包括环境搭建、Maven依赖配置、SpringBoot配置文件... 目录1. 引言1.1 环境准备1.2 Druid介绍2. 配置Druid连接池3. 查看Druid监控

Linux在线解压jar包的实现方式

《Linux在线解压jar包的实现方式》:本文主要介绍Linux在线解压jar包的实现方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录linux在线解压jar包解压 jar包的步骤总结Linux在线解压jar包在 Centos 中解压 jar 包可以使用 u

c++ 类成员变量默认初始值的实现

《c++类成员变量默认初始值的实现》本文主要介绍了c++类成员变量默认初始值,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 目录C++类成员变量初始化c++类的变量的初始化在C++中,如果使用类成员变量时未给定其初始值,那么它将被

Qt使用QSqlDatabase连接MySQL实现增删改查功能

《Qt使用QSqlDatabase连接MySQL实现增删改查功能》这篇文章主要为大家详细介绍了Qt如何使用QSqlDatabase连接MySQL实现增删改查功能,文中的示例代码讲解详细,感兴趣的小伙伴... 目录一、创建数据表二、连接mysql数据库三、封装成一个完整的轻量级 ORM 风格类3.1 表结构

基于Python实现一个图片拆分工具

《基于Python实现一个图片拆分工具》这篇文章主要为大家详细介绍了如何基于Python实现一个图片拆分工具,可以根据需要的行数和列数进行拆分,感兴趣的小伙伴可以跟随小编一起学习一下... 简单介绍先自己选择输入的图片,默认是输出到项目文件夹中,可以自己选择其他的文件夹,选择需要拆分的行数和列数,可以通过

Python中将嵌套列表扁平化的多种实现方法

《Python中将嵌套列表扁平化的多种实现方法》在Python编程中,我们常常会遇到需要将嵌套列表(即列表中包含列表)转换为一个一维的扁平列表的需求,本文将给大家介绍了多种实现这一目标的方法,需要的朋... 目录python中将嵌套列表扁平化的方法技术背景实现步骤1. 使用嵌套列表推导式2. 使用itert

Python使用pip工具实现包自动更新的多种方法

《Python使用pip工具实现包自动更新的多种方法》本文深入探讨了使用Python的pip工具实现包自动更新的各种方法和技术,我们将从基础概念开始,逐步介绍手动更新方法、自动化脚本编写、结合CI/C... 目录1. 背景介绍1.1 目的和范围1.2 预期读者1.3 文档结构概述1.4 术语表1.4.1 核

在Linux中改变echo输出颜色的实现方法

《在Linux中改变echo输出颜色的实现方法》在Linux系统的命令行环境下,为了使输出信息更加清晰、突出,便于用户快速识别和区分不同类型的信息,常常需要改变echo命令的输出颜色,所以本文给大家介... 目python录在linux中改变echo输出颜色的方法技术背景实现步骤使用ANSI转义码使用tpu

Python使用python-can实现合并BLF文件

《Python使用python-can实现合并BLF文件》python-can库是Python生态中专注于CAN总线通信与数据处理的强大工具,本文将使用python-can为BLF文件合并提供高效灵活... 目录一、python-can 库:CAN 数据处理的利器二、BLF 文件合并核心代码解析1. 基础合