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绘制人人网好友关系图示例
Apr 01 Python
使用url_helper简化Python中Django框架的url配置教程
May 30 Python
Python读取键盘输入的2种方法
Jun 16 Python
python操作MySQL 模拟简单银行转账操作
Sep 27 Python
python基础教程项目五之虚拟茶话会
Apr 02 Python
Python socket套接字实现C/S模式远程命令执行功能案例
Jul 06 Python
Django Python 获取请求头信息Content-Range的方法
Aug 06 Python
win10环境下配置vscode python开发环境的教程详解
Oct 16 Python
Pandas+Matplotlib 箱式图异常值分析示例
Dec 09 Python
python同义词替换的实现(jieba分词)
Jan 21 Python
Python中import导入不同目录的模块方法详解
Feb 18 Python
python实现时间序列自相关图(acf)、偏自相关图(pacf)教程
Jun 03 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 运行效率总结(提示程序速度)
2009/11/26 PHP
PHP学习笔记之一
2011/01/17 PHP
phpadmin如何导入导出大数据文件及php.ini参数修改
2013/02/18 PHP
一个简单的PHP验证码实现代码
2014/05/10 PHP
PHP的邮件群发系统phplist配置方法详细总结
2016/03/30 PHP
用正则表达式 动态创建/增加css style script 兼容IE firefox
2009/03/10 Javascript
用nodejs写的一个简单项目打包工具
2013/05/11 NodeJs
JS 实现点击a标签的时候让其背景更换
2013/10/15 Javascript
转换字符串为json对象的方法详解
2013/11/29 Javascript
快速解决Canvas.toDataURL 图片跨域的问题
2016/05/10 Javascript
laydate日历控件使用方法详解
2017/11/20 Javascript
简单的Vue SSR的示例代码
2018/01/12 Javascript
webpack里使用jquery.mCustomScrollbar插件的方法
2018/05/30 jQuery
AngularJS 事件发布机制
2018/08/28 Javascript
深入理解Vue 的钩子函数
2018/09/05 Javascript
node.js 模块和其下载资源的镜像设置的方法
2018/09/06 Javascript
详解async/await 异步应用的常用场景
2019/05/13 Javascript
JS中的算法与数据结构之列表(List)实例详解
2019/08/16 Javascript
在Vue中使用this.$store或者是$route一直报错的解决
2019/11/08 Javascript
node.js使用zlib模块进行数据压缩和解压操作示例
2020/02/12 Javascript
python使用PyFetion来发送短信的例子
2014/04/22 Python
Pycharm学习教程(5) Python快捷键相关设置
2017/05/03 Python
django中的HTML控件及参数传递方法
2018/03/20 Python
Python模拟登入的N种方式(建议收藏)
2020/05/31 Python
基于Keras中Conv1D和Conv2D的区别说明
2020/06/19 Python
澳大利亚宠物食品和用品商店:PETstock
2020/01/02 全球购物
美国最大的购物网站:Amazon.com(亚马逊美国)
2020/05/23 全球购物
简单叙述一下MYSQL的优化
2016/05/09 面试题
《我的第一本书》教学反思
2014/02/15 职场文书
社区节水倡议书
2015/04/29 职场文书
小学教研工作总结2015
2015/05/13 职场文书
离婚起诉书怎么写
2015/05/19 职场文书
2019年朋友圈经典励志语录50条
2019/07/05 职场文书
利用ajax+php实现商品价格计算
2021/03/31 PHP
Apache Calcite 实现方言转换的代码
2021/04/24 Servers
Win11开始菜单添加休眠选项
2022/04/19 数码科技