Python&&GDAL實(shí)現(xiàn)NDVI的計算方式
更新時間：2020年01月09日 09:04:12 作者：蘭小莫
今天小編就為大家分享一篇Python&&GDAL實(shí)現(xiàn)NDVI的計算方式，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧
很短的幾句代碼，可是我卻花了很長的時間才寫出來，因?yàn)閍rray那里的除法運(yùn)算結(jié)果老是不對，正常在-1-1之間的。從別的資料摘來處理NDVI計算的array代碼處，出現(xiàn)了很多問題，可能它用了什么優(yōu)化計算的函數(shù)，但是結(jié)果不對，果斷放棄了。
直接硬算，影像波段是整數(shù)，轉(zhuǎn)成浮點(diǎn)型數(shù)字就行，然后再參與運(yùn)算得出了正確的結(jié)果范圍。
這個效率還是不行，用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實(shí)現(xiàn)NDVI的計算方式就是小編分享給大家的全部內(nèi)容了，希望能給大家一個參考，也希望大家多多支持腳本之家。
您可能感興趣的文章:
相關(guān)文章
這篇文章主要介紹了python?turtle庫畫圣誕樹詳細(xì)代碼教程,圣誕節(jié)快到了，下面小編就來利用python?turtle庫畫一顆圣誕樹,?主要成分有圣誕樹的本體、大小蝴蝶結(jié)、星星、圣誕帽和襪子，需要的朋友可以參考一下2021-12-12
Python的Tkinter點(diǎn)擊按鈕觸發(fā)事件的例子
今天小編就為大家分享一篇Python的Tkinter點(diǎn)擊按鈕觸發(fā)事件的例子，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2019-07-07
從matplotlib的角度來說,繪圖風(fēng)格也算是圖像類型的一部分,所以這篇文章小編想帶大家了解一下Python中matplotlib的繪圖風(fēng)格,有需要的可以參考下2023-09-09
這篇文章主要給大家介紹了利用Python判斷變量是否為Json格式的字符串的相關(guān)資料，文中給出了詳細(xì)的示例代碼供大家參考學(xué)習(xí)，需要的朋友們下面來一起看看吧。2017-05-05
Python定時任務(wù)框架APScheduler原理及常用代碼
這篇文章主要介紹了Python定時任務(wù)框架APScheduler原理及常用代碼,文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友可以參考下2020-10-10
Pytorch實(shí)現(xiàn)LSTM和GRU示例
今天小編就為大家分享一篇Pytorch實(shí)現(xiàn)LSTM和GRU示例，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2020-01-01
本文給大家分享的是走著的一些強(qiáng)制殺掉python線程經(jīng)驗(yàn)教訓(xùn)，如果你使用強(qiáng)制手段干掉線程，那么很大幾率出現(xiàn)意想不到的bug。 請記住一點(diǎn)，鎖資源不會因?yàn)榫€程退出而釋放鎖資源 ！2017-02-02
本文主要介紹了pycharm安裝opencv-python報錯的解決，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2023-07-07
Pycharm配置Qt Designer及Pyuic的實(shí)現(xiàn)方法
本文介紹了如何安裝Qt designer和Pyuic以及他們的基本用法，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2023-07-07
這篇文章主要介紹了python根據(jù)文件大小打log日志的方法,可實(shí)現(xiàn)循環(huán)打日志,即第一個文件達(dá)到maxBytes大小后就寫入第二個文件,以此類推,是非常實(shí)用的技巧,需要的朋友可以參考下2014-10-10