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使用bar绘制堆积/带误差棒柱形图的实现

    Python使用bar绘制堆积/带误差棒柱形图的实现

    本文先讲解bar参数如何使用,然后分别演示堆积柱形图和带误差柱形图画法。具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2021-09-09
  • Python爬虫框架Scrapy简介

    Python爬虫框架Scrapy简介

    这篇文章介绍了Python爬虫框架Scrapy,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2022-06-06
  • Python中np.linalg.norm()用法实例总结

    Python中np.linalg.norm()用法实例总结

    在线性代数中一个向量通过矩阵转换成另一个向量时,原有向量的大小就是向量的范数,这个变化过程的大小就是矩阵的范数,下面这篇文章主要给大家介绍了关于Python中np.linalg.norm()用法的相关资料,需要的朋友可以参考下
    2022-07-07
  • Python实现获取视频时长功能

    Python实现获取视频时长功能

    这篇文章主要介绍了Python如何实现获取视频时长功能,可以精确到毫秒。文中的示例代码简洁易懂,对我们的学习有一定的帮助,感兴趣的可以了解一下
    2021-12-12
  • 详解python中DRF框架的数据校验方式

    详解python中DRF框架的数据校验方式

    这篇文章主要为大家详细介绍了python中DRF框架的数据校验方式,文中的示例代码讲解详细,具有一定的学习价值,感兴趣的小伙伴可以跟随小编一起了解一下
    2023-10-10
  • 基于Tensorflow一维卷积用法详解

    基于Tensorflow一维卷积用法详解

    这篇文章主要介绍了基于Tensorflow一维卷积用法详解,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2020-05-05
  • django实现用户登陆功能详解

    django实现用户登陆功能详解

    这篇文章主要介绍了django实现用户登陆功能详解,具有一定借鉴价值,需要的朋友可以参考下。
    2017-12-12
  • Python OrderedDict字典排序方法详解

    Python OrderedDict字典排序方法详解

    这篇文章主要介绍了Python OrderedDict字典使用方法详解,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
    2020-05-05
  • 详解Python如何获取视频文件的大小和时长

    详解Python如何获取视频文件的大小和时长

    这篇文章主要为大家详细介绍了Python如何实现获取视频文件的大小和时长,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起了解一下
    2023-03-03
  • Python创建xml文件示例

    Python创建xml文件示例

    这篇文章主要介绍了Python创建xml文件的方法,结合实例形式分析了Python针对xml格式数据及文件读写相关操作技巧,需要的朋友可以参考下
    2017-03-03

最新评论