关于gdal numpy计算NDVI等指数的问题

最新推荐文章于 2023-12-16 22:51:18 发布
原创 最新推荐文章于 2023-12-16 22:51:18 发布 · 303 阅读 · 0 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
文章标签：

#numpy

本文讲述了在使用Python的GDAL和NumPy计算NDVI时遇到的NaN和Inf问题，作者提供了如何在QGIS外部处理这些异常值的方法，包括用nan_to_num函数替换NaN，将Inf设置为0，并对结果进行二值化处理。

在利用python gdal numpy计算NDVI等运算栅格数组的指数时出现nan和inf的问题。

_band3_np_data = np.array(_band3_ds_data)
_band6_np_data = np.array(_band6_ds_data)

以上得到np多维数组
直接进行运算

 _band_ndsi = (_band3_ds_data - _band6_ds_data) / (_band3_ds_data + _band6_ds_data)

会出现nan和inf，在QGIS中加载出现问

最低0.47元/天 解锁文章
【ArcPy】 Python批量提取多张遥感影像的NDVI值
陈书予
09-05 685
在处理遥感影像时，理解影像的波段结构和坐标系统非常重要。一般来说，遥感影像会包含多个波段，每个波段代表特定的光谱范围。Red波段：通常对应于第三波段（在Landsat影像中为Band 3）。NIR波段：通常对应于第四波段（在Landsat影像中为Band 4）。通过本文的介绍，你已经学会了如何使用Python批量提取多张遥感影像的NDVI值。这一技能在遥感影像分析中具有广泛的应用前景，可以帮助你更高效地分析植被覆盖情况。
Python遥感图像处理应用篇(八)：批量计算landsat8中NDVI指数数据置信区间并保存结果
soderayer的博客
05-02 3875
""" @soderayer date:20220428 to read NDVI files testfile Landsat8 collection2 level2 NDVI data calculate the confidence interval to csv data info table： Attribute Value Long Name Normalized Difference Vegetation Index Short Name LC8NDVI, LE7ND...
arcgis批量求ndvi_arcgis批量ndvi_arcgis_批量ndvi_
09-30
利用编写的Python代码，批量求navi。
Python批量计算遥感图像NDVI（归一化植被指数）
热门推荐
qq_43177210的博客
07-11 1万+
归一化植被指数：遥感影像中，近红外波段的反射值与红光波段的反射值之差比上两者之和。 计算公式：NDVI = (NIR-R)/(NIR+R) 注：NIR为近红外波段的反射值 R为红光波段的反射值。 由于学习需要，从网上寻找方法试图批量计算该指数，防止以后忘记。附代码如下： （主要使用gdal库） import os from PIL import Image import numpy as np from osgeo import gdal import glob import ...
【Python】GDAL库计算NDVI指数简易教程
weixin_47262549的博客
03-10 746
使用gdal库计算NDVI的简易版教程
gdal+python计算NDVI
Harry233的博客
10-23 1260
使用GDAL+python写NDVI
Python&&GDAL实现NDVI的计算方式
09-18
在地理信息系统和遥感领域，归一化植被指数（NDVI）是一个衡量植被生长状况的重要指标，它通过比较近红外波段（NIR）和红光波段（RED）的反射率差异来计算。Python作为一种流行的数据分析编程语言，结合GDAL...
Python遥感图像处理应用篇(二十一)：Python+GDAL 批量计算遥感图像NDVI指数
soderayer的博客
07-03 2351
之前也写过NDVI批量计算的实现方法，采用的是Acrpy python2.7实现的，这里我们采用GDAL实现NDVI指数的批量计算。使用数据，多波段遥感图像数据。Landsat05 Colection2 Level2数据。NDVI=(NIR-RED)/(NIR+RED) 就算之前需要搞清楚多光谱数据中的波段对应关系。影像数据地表反射率修正系数：0.0000275+（-0.2）运行环境：运行环境：windows10 pycharm python3.7.7 GDAL-3.2.3-cp37-cp37m-w
基于Python的gdal库读取遥感影像TIFF，Python批量对NDVI植被指数计算（源代码）
qq_41240287的博客
03-17 2238
一.单次针对波段分开的TIFF影像的NDVI植被指数的计算 1.代码如下：（代码中已经将NDVI指数异常值进行了剔除，取值范围最终在[-1, 1]区间内） （为了方便程序运行的进度查看，在函数内增添处理过程进度条显示） from osgeo import gdal, gdalconst import os import numpy as np def NDVI_calculte(Red_band_file, NIR_band_file, NDVI_output_file): """ ..
python环境下利用单波段影像计算NDVI并可视化
qq_26809799的博客
01-14 1074
from osgeo import gdal_array import gdal, ogr, os, osr import numpy as np import cv2 #文件位置 b1=r'C:\Users\Administrator\Desktop\python_RandomText\maqu_redband_3.TIF' b4=r'C:\Users\Administrator\Desktop\python_RandomText\maqu_Nearredband_3.TIF' #将单波段影像转换成数组
利用python计算NDVI
qq_26809799的博客
01-15 1785
一、代码如下 import osgeo import gdal from matplotlib import pyplot as plt import numpy as np import cv2 from PIL import Image import cv2 #--1)读取遥感影像：尺寸（宽度和长度）、波段数、投影信息、数组 def NDVI_tif(filename): dataset = gdal.Open(filename) # 打开文件 im_width = datase
使用python实现ndvi归一化植被指数计算
木易GIS的博客
02-07 1422
使用python实现ndvi归一化植被指数计算
python+gdal+numpy计算ndvi
llc的博客
09-05 2134
from osgeo import gdal import numpy as np ds: gdal.Dataset = gdal.Open('gf.tif') # 红波段 red = ds.GetRasterBand(3).ReadAsArray() * 0.0001 # 近红外波段 nir = ds.GetRasterBand(4).ReadAsArray() * 0.0001 np.seterr(all="ignore") ndvi = (nir - red) / (nir + red) #
Python遥感图像处理指南(2)-卫星影像分析MNDWI提取
gis收藏家的博客
12-16 1470
在第一篇文章中，我们使用了我们导入的软件包中已定义的函数，如 rasterio 中的 open 和 read 以及 matplotlib 中的 imshow。现在，我们可以编写自己的函数了。函数是一段代码，它接收一些参数作为输入并返回输出（经过一些计算），可以在代码中多次重复使用，也可以被其他代码重用。在 Python 中定义一个新函数非常简单。return v在此语法示例中，arg1和arg2强制参数，应该传递给函数，否则会引发错误，而arg3是一个可选参数，如果是，则假定默认值 1从函数调用中省略。
计算ndvi值需要的数据_大气校正前后山地流域NDVI变化特征对比
weixin_29164185的博客
12-31 8696
摘要大气校正是遥感影像处理的重要环节,探究大气校正前后归一化植被指数(normalized difference vegetation index,NDVI)特征及其地形梯度特征具有一定理论意义。以贵州省境内包溪河流域为例,基于Landsat8遥感影像和FLAASH大气校正,探究NDVI分等级和分地形计算误差。得出以下结论: ①未进行大气校正所得NDVI在空间上总体被低估; ②未经大气校...
【ArcGIS遇上Python】ArcGIS Python实现Modis NDVI批量求年最大值
「 刘一哥与GIS的故事」
08-29 1万+
一年中的12个月份的月最大合成（mvc）数据放在“F:\\Vegetation Change\\Data\\GIMMS Data\\1MVC\\"，数据名称格式为mvc_198801,mvc_198802........mvc_198812。处理年份为1981-2006，代码为： import arcpy arcpy.CheckOutExtension("spatial") for i in r
Arcgis快速计算NDVI
所有专家都曾始于新手
10-02 7551
如何使用Arcgis像ENVI一样波段计算NDVI的值，事实上，Arcgis更快速一些。
(结果不对)python调用ArcGIS工具批量mvc最大化合成NDVI每月数据
qq_37254633的博客
01-24 1963
数据：近十年MOD13Q1的NDVI数据，每隔16天一期，每个月有两期 数据命名方式： 使用ArcGIS的 ArcToolbox—Spatial Analys工具—局部分析—像元统计数据 窗口进行最大化合成得到有效值正确结果如图： 使用python代码批量处理有效值错误： import arcpy from arcpy import env from arcpy.sa import * arcpy.CheckOutExtension("spatial") env.workspace = "C:/MOD13
python处理遥感数据（NVDI计算、辐射校正）
baidu_34785556的博客
12-27 9872
前言：作为遥感、GIS从业者，会用python将能提高你的核心竞争力，同时也证明你可以摆脱专业软件的束缚，迈入自主研发的新阶段。ps:其实ENVI、Arcgis等相关软件已经非常完善了，想要突破比较困难，但是这一小步的迈出也决定了你已跨入自主研发的道路。 目录 1.辐射校正 1.1传感器端辐射校正 1.2大气校正 1.3地表辐射校正 2.代码 3.结束 1.辐射校正 辐射校正的目...
python gdal ndvi计算
09-05
要使用Python gdal来计算NDVI，首先需要导入gdal库和numpy库。然后，使用gdal打开需要计算NDVI的遥感图像。 接下来，可以使用gdal库的GetRasterBand方法来读取图像的红光波段和近红外波段数据。可以使用gdal库的...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值