Python&&GDAL实现NDVI的计算方式

 更新时间:2020年01月09日 09:04:12   作者:兰小莫  
今天小编就为大家分享一篇Python&&GDAL实现NDVI的计算方式,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧

很短的几句代码,可是我却花了很长的时间才写出来,因为array那里的除法运算结果老是不对,正常在-1-1之间的。从别的资料摘来处理NDVI计算的array代码处,出现了很多问题,可能它用了什么优化计算的函数,但是结果不对,果断放弃了。

直接硬算,影像波段是整数,转成浮点型数字就行,然后再参与运算得出了正确的结果范围。

这个效率还是不行,用c++&&GDAL计算速率快得多了。

from osgeo import gdal_array as ga
import gdal, ogr, os, osr
import numpy as np
b3=r'C:\Users\suns\Desktop\b4.TIF'
b4=r'C:\Users\suns\Desktop\b5.TIF'
arr=ga.LoadFile(b3)
arr1=ga.LoadFile(b4)
ga.numpy.seterr(all="ignore")
ndvi=((arr1-arr)*1.0)/((arr1+arr)*1.0)
ndvi1=ga.numpy.nan_to_num(ndvi)
target=r'C:\Users\suns\Desktop\ndvi1.tif'
out=ga.SaveArray(ndvi1,target,format = "GTiff",prototype = b4)
out=None

以上这篇Python&&GDAL实现NDVI的计算方式就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持脚本之家。

相关文章

  • python turtle库画圣诞树详细代码教程

    python turtle库画圣诞树详细代码教程

    这篇文章主要介绍了python turtle库画圣诞树详细代码教程,圣诞节快到了,下面小编就来利用python turtle库画一颗圣诞树, 主要成分有圣诞树的本体、大小蝴蝶结、星星、圣诞帽和袜子,需要的朋友可以参考一下
    2021-12-12
  • Python的Tkinter点击按钮触发事件的例子

    Python的Tkinter点击按钮触发事件的例子

    今天小编就为大家分享一篇Python的Tkinter点击按钮触发事件的例子,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2019-07-07
  • Python matplotlib绘图风格详解

    Python matplotlib绘图风格详解

    从matplotlib的角度来说,绘图风格也算是图像类型的一部分,所以这篇文章小编想带大家了解一下Python中matplotlib的绘图风格,有需要的可以参考下
    2023-09-09
  • Python判断变量是否为Json格式的字符串示例

    Python判断变量是否为Json格式的字符串示例

    这篇文章主要给大家介绍了利用Python判断变量是否为Json格式的字符串的相关资料,文中给出了详细的示例代码供大家参考学习,需要的朋友们下面来一起看看吧。
    2017-05-05
  • Python定时任务框架APScheduler原理及常用代码

    Python定时任务框架APScheduler原理及常用代码

    这篇文章主要介绍了Python定时任务框架APScheduler原理及常用代码,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
    2020-10-10
  • Pytorch实现LSTM和GRU示例

    Pytorch实现LSTM和GRU示例

    今天小编就为大家分享一篇Pytorch实现LSTM和GRU示例,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2020-01-01
  • 不要用强制方法杀掉python线程

    不要用强制方法杀掉python线程

    本文给大家分享的是走着的一些强制杀掉python线程经验教训,如果你使用强制手段干掉线程,那么很大几率出现意想不到的bug。 请记住一点,锁资源不会因为线程退出而释放锁资源 !
    2017-02-02
  • pycharm安装opencv-python报错的解决

    pycharm安装opencv-python报错的解决

    本文主要介绍了pycharm安装opencv-python报错的解决,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2023-07-07
  • Pycharm配置Qt Designer及Pyuic的实现方法

    Pycharm配置Qt Designer及Pyuic的实现方法

    本文介绍了如何安装Qt designer和Pyuic以及他们的基本用法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2023-07-07
  • python根据文件大小打log日志

    python根据文件大小打log日志

    这篇文章主要介绍了python根据文件大小打log日志的方法,可实现循环打日志,即第一个文件达到maxBytes大小后就写入第二个文件,以此类推,是非常实用的技巧,需要的朋友可以参考下
    2014-10-10

最新评论