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 过滤字符串的技巧,map与itertools.imap
Sep 06 Python
web.py在模板中输出美元符号的方法
Aug 26 Python
Python描述器descriptor详解
Feb 03 Python
用Python创建声明性迷你语言的教程
Apr 13 Python
Django中URL视图函数的一些高级概念介绍
Jul 20 Python
python基于隐马尔可夫模型实现中文拼音输入
Apr 01 Python
Python实现的简单计算器功能详解
Aug 25 Python
python实现LBP方法提取图像纹理特征实现分类的步骤
Jul 11 Python
python自带tkinter库实现棋盘覆盖图形界面
Jul 17 Python
pandas 对日期类型数据的处理方法详解
Aug 08 Python
Python 获取项目根路径的代码
Sep 27 Python
python中uuid模块实例浅析
Dec 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
论建造顺序的重要性
2020/03/04 星际争霸
4.与数据库的连接
2006/10/09 PHP
基于PHP读取TXT文件向数据库导入海量数据的方法
2013/04/23 PHP
Codeigniter框架的更新事务(transaction)BUG及解决方法
2014/07/25 PHP
PHP往XML中添加节点的方法
2015/03/12 PHP
Javascript异步表单提交,图片上传,兼容异步模拟ajax技术
2010/05/10 Javascript
Textbox控件注册回车事件及触发按钮提交事件具体实现
2013/03/04 Javascript
使用JavaScript的ActiveXObject对象检测应用程序是否安装的方法
2014/04/15 Javascript
javascript与有限状态机详解
2014/05/08 Javascript
基于jQuery的JavaScript模版引擎JsRender使用指南
2014/12/29 Javascript
javascript常见数据验证插件大全
2015/08/03 Javascript
基于jQuery和CSS3制作数字时钟附源码下载(jquery篇)
2015/11/24 Javascript
浅谈js的url解析函数封装
2016/06/28 Javascript
详解JSON1:使用TSQL查询数据和更新JSON数据
2016/11/21 Javascript
深入理解基于vue-cli的vuex配置
2017/07/24 Javascript
基于js 各种排序方法和sort方法的区别(详解)
2018/01/03 Javascript
在antd Table中插入可编辑的单元格实例
2020/10/28 Javascript
py2exe 编译ico图标的代码
2013/03/08 Python
Python编程之列表操作实例详解【创建、使用、更新、删除】
2017/07/22 Python
Python数据分析之获取双色球历史信息的方法示例
2018/02/03 Python
Python中的单行、多行、中文注释方法
2018/07/19 Python
python中的常量和变量代码详解
2018/07/25 Python
Django添加bootstrap框架时无法加载静态文件的解决方式
2020/03/27 Python
解决import tensorflow as tf 出错的原因
2020/04/16 Python
Python中SQLite如何使用
2020/05/27 Python
CSS3中currentColor关键字的妙用
2016/02/27 HTML / CSS
优秀体育委员自荐书
2014/01/31 职场文书
中班中秋节活动反思
2014/02/18 职场文书
护士上岗前培训自我鉴定
2014/04/20 职场文书
国际语言毕业生求职信
2014/07/08 职场文书
地道战观后感2000字
2015/06/04 职场文书
物资采购管理制度
2015/08/06 职场文书
2016年师德学习心得体会
2016/01/12 职场文书
三好学生竞选稿范文
2019/08/21 职场文书
Vue监视数据的原理详解
2022/02/24 Vue.js
Go 内联优化让程序员爱不释手
2022/06/21 Golang