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分析apache访问日志脚本分享
Feb 26 Python
python实现爬虫统计学校BBS男女比例(一)
Dec 31 Python
Java Web开发过程中登陆模块的验证码的实现方式总结
May 25 Python
利用Tkinter(python3.6)实现一个简单计算器
Dec 21 Python
python实现自动发送报警监控邮件
Jun 21 Python
详解关于Django中ORM数据库迁移的配置
Oct 08 Python
Python基于滑动平均思想实现缺失数据填充的方法
Feb 21 Python
Python中的相关分析correlation analysis的实现
Aug 29 Python
python之array赋值技巧分享
Nov 28 Python
python实现从wind导入数据
Dec 03 Python
Python模块_PyLibTiff读取tif文件的实例
Jan 13 Python
Python列表list操作相关知识小结
Jan 29 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
叶罗丽:为什么大家对颜冰这对CP非常关心,却对金茉两人十分冷漠
2020/03/17 国漫
phpmyadmin操作流程
2006/10/09 PHP
php 中的str_replace 函数总结
2007/04/27 PHP
php数组函数序列之array_unshift() 在数组开头插入一个或多个元素
2011/11/07 PHP
Android App中DrawerLayout抽屉效果的菜单编写实例
2016/03/21 PHP
PHP实现求连续子数组最大和问题2种解决方法
2017/12/26 PHP
javascript 二分法(数组array)
2010/04/24 Javascript
推荐30个新鲜出炉的精美 jQuery 效果
2012/03/26 Javascript
js 获取时间间隔实现代码
2014/05/12 Javascript
js简单实现标签云效果实例
2015/08/06 Javascript
jquery实现鼠标经过显示下划线的渐变下拉菜单效果代码
2015/08/24 Javascript
JS锚点的设置与使用方法
2016/09/05 Javascript
AngularJS 自定义过滤器详解及实例代码
2016/09/14 Javascript
利用n 升级工具升级Node.js版本及在mac环境下的坑
2017/02/15 Javascript
bootstrap的常用组件和栅格式布局详解
2017/05/02 Javascript
jQuery中可见性过滤器简单用法示例
2018/03/31 jQuery
解决vue router组件状态刷新消失的问题
2018/08/01 Javascript
Vuex 单状态库与多模块状态库详解
2018/12/11 Javascript
JS实现鼠标拖拽盒子移动及右键点击盒子消失效果示例
2019/01/29 Javascript
图解JS原型和原型链实现原理
2020/09/15 Javascript
python多线程方式执行多个bat代码
2016/06/07 Python
Python 中Pickle库的使用详解
2018/02/24 Python
详解Python self 参数
2019/08/30 Python
详解Python 中sys.stdin.readline()的用法
2019/09/12 Python
python滑块验证码的破解实现
2019/11/10 Python
Python制作数据预测集成工具(值得收藏)
2020/08/21 Python
python基于pexpect库自动获取日志信息
2021/02/01 Python
英国家庭家具、照明和花园家具购物网站:Furniture123
2018/12/31 全球购物
青年创业培训欢迎词
2014/01/10 职场文书
房屋买卖协议书
2014/04/10 职场文书
企业人事任命书
2014/06/05 职场文书
2014乡镇党政班子四风问题思想汇报
2014/09/14 职场文书
2015年毕业生实习评语
2015/03/25 职场文书
Python OpenCV快速入门教程
2021/04/17 Python
简单谈谈Python面向对象的相关知识
2021/06/28 Python
深入理解go slice结构
2021/09/15 Golang