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


Posted in Python onJanuary 09, 2020

很短的几句代码,可是我却花了很长的时间才写出来,因为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&&GDAL实现NDVI的计算方式就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python批量提交沙箱问题实例
Oct 08 Python
Python六大开源框架对比
Oct 19 Python
Python基本语法经典教程
Mar 11 Python
python简单实例训练(21~30)
Nov 15 Python
Tensorflow环境搭建的方法步骤
Feb 07 Python
PyQt 实现使窗口中的元素跟随窗口大小的变化而变化
Jun 18 Python
使用python实现对元素的长截图功能
Nov 14 Python
Python timer定时器两种常用方法解析
Jan 20 Python
Python文件夹批处理操作代码实例
Jul 21 Python
详解python中GPU版本的opencv常用方法介绍
Jul 24 Python
python3 os进行嵌套操作的实例讲解
Nov 19 Python
将不规则的Python多维数组拉平到一维的方法实现
Jan 11 Python
Python-numpy实现灰度图像的分块和合并方式
Jan 09 #Python
Pytorch 搭建分类回归神经网络并用GPU进行加速的例子
Jan 09 #Python
使用PyTorch训练一个图像分类器实例
Jan 08 #Python
pytorch 实现将自己的图片数据处理成可以训练的图片类型
Jan 08 #Python
pytorch下大型数据集(大型图片)的导入方式
Jan 08 #Python
Python 实现训练集、测试集随机划分
Jan 08 #Python
Pyecharts绘制全球流向图的示例代码
Jan 08 #Python
You might like
php替换超长文本中的特殊字符的函数代码
2012/05/22 PHP
PHP整数取余返回负数的相关解决方法
2014/05/15 PHP
javascript模拟地球旋转效果代码实例
2013/12/02 Javascript
解决JS中乘法的浮点错误的方法
2014/01/03 Javascript
用js模拟struts2的多action调用示例
2014/05/19 Javascript
JS与jQ读取xml文件的方法
2015/12/08 Javascript
如何使用jquery实现文字上下滚动效果
2016/10/12 Javascript
Vue.js 递归组件实现树形菜单(实例分享)
2016/12/21 Javascript
webpack打包单页面如何引用的js
2017/06/07 Javascript
详解Vue组件实现tips的总结
2017/11/01 Javascript
js循环map 获取所有的key和value的实现代码(json)
2018/05/09 Javascript
微信小程序实现单列下拉菜单效果
2019/04/25 Javascript
Nodejs监控事件循环异常示例详解
2019/09/22 NodeJs
使用Vant完成通知栏Notify的提示操作
2020/11/11 Javascript
python操作MongoDB基础知识
2013/11/01 Python
Python递归函数定义与用法示例
2017/06/02 Python
python定时利用QQ邮件发送天气预报的实例
2017/11/17 Python
django使用html模板减少代码代码解析
2017/12/12 Python
如何在python中使用selenium的示例
2017/12/26 Python
python实现从本地摄像头和网络摄像头截取图片功能
2019/07/11 Python
python函数参数(必须参数、可变参数、关键字参数)
2019/08/16 Python
python读取ini配置的类封装代码实例
2020/01/08 Python
使用python采集Excel表中某一格数据
2020/05/14 Python
ITK 实现多张图像转成单个nii.gz或mha文件案例
2020/07/01 Python
django 将自带的数据库sqlite3改成mysql实例
2020/07/09 Python
自我评价怎么写好呢?
2013/12/05 职场文书
超市业务员岗位职责
2013/12/05 职场文书
运动会100米解说词
2014/01/23 职场文书
师生聚会感言
2014/01/26 职场文书
计算机科学系职业生涯规划书
2014/03/08 职场文书
李敖北大演讲稿
2014/05/24 职场文书
大学拉赞助协议书范文
2014/09/26 职场文书
毕业设计论文评语
2014/12/31 职场文书
初中开学典礼新闻稿
2015/07/17 职场文书
公司年会主持词范文!
2019/05/07 职场文书
oracle DGMGRL ORA-16603报错的解决方法(DG Broker)
2021/04/06 Oracle