gdal（1）GF-1影像提取NDVI

GF-1号多光谱影像 band3为红波段 band4为近红外波段

话不多说，直接上代码

# -*- coding:utf_8 -*-
"""
作者：ZHOU
日期：2021年09月21日
计算NDVI数值
openfile函数为读遥感影像
ndvi函数为计算NDVI
write函数为写tif数据
"""
try:import gdal
except:from osgeo import gdal
from osgeo import ogr
from osgeo.gdalconst import *
import numpy as npdef openfile(filepath):dataset = gdal.Open(filepath)dataset_xsize = dataset.RasterXSizedataset_ysize = dataset.RasterYSizedataset_data = dataset.ReadAsArray(0,0,dataset_xsize,dataset_ysize)print(dataset_data.shape)return dataset,dataset_datadef ndvi(red,nir,dataset_data):red_band = dataset_data[red]nir_abnd = dataset_data[nir]red_band,nir_abnd = red_band.astype(np.float64),nir_abnd.astype(np.float64)add = red_band+nir_abndsub = nir_abnd-red_bandndvi_data = np.divide(sub,add,out=np.zeros_like(sub),where=add!=0)print(ndvi_data.shape)return ndvi_datadef writefile(dataset,in_data,path):in_prj= dataset.GetProjection()in_geo = dataset.GetGeoTransform()if len(in_data) == 3:in_bands,in_ysize,in_xsize = in_data.shapeelse:in_xsize , in_ysize = in_data.shapein_bands = 1driver = gdal.GetDriverByName('GTiff')new_data = driver.Create(path,in_xsize,in_ysize,in_bands,GDT_CFloat64)print(in_bands)if new_data != None:new_data.SetProjection(in_prj)new_data.SetGeoTransform(in_geo)if in_bands ==1:new_data.GetRasterBand(1).WriteArray(in_data)else:for i in range(in_bands):new_data.GetRasterBand(i+1).WriteArray(in_data[i])del datasetif __name__ == '__main__':dataset,dataset_data = openfile(r'')ndvi_data = ndvi(2,3,dataset_data)writefile(dataset,ndvi_data,r'')

gdal（1）GF-1影像提取NDVI相关推荐

TM影像提取NDVI值小记
昨天下载了Landset8的TM影像,云很多,并做预处理,然后融合,提取NDVI值,竟然达到了几百多,很吓人,之前处理GF1WFV的影像并没有出现这种情况,和同事讨论了一下,觉得是融合的时候出了问题, ...
ENVI中基于Modis影像提取NDVI处理步骤
1.数据下载数据为MODIS-15天叶面积指数数据,可在地理空间数据云中进行下载,行列号可以参照下图. 2.提取NDVI数据在Envi classic或者4.8,安装"MODIS Con ...
ENVI软件对Landsat-8数据进行辐射定标、大气校正、提取NDVI、估算植被覆盖度等操作【图说GIS】
目录一.前言二.数据下载第一步访问网站第二步登录账户第三步点击高级检索, 下载数据三.软件下载及安装四.辐射定标第一步打开数据第二步辐射定标五.计算区域平均高程第一步 ...
ENVI提取NDVI与植被覆盖度估算
公众号原文链接,欢迎关注 gENVI提取NDVI与植被覆盖度估算ENVI中NDVI计算与植被覆盖度提取https://mp.weixin.qq.com/s/Uc_14ObiYKKafQrtKqFjvg ...
【ArcGIS遇上Python】ArcGIS批量处理栅格影像（NDVI）归一化完整案例代码
文章<ArcGIS栅格影像(NDVI)归一化处理的两种方法精解>详细介绍了ArcGIS中进行NDVI影像归一化处理的两种方法,但是无法批量进行,本文就放大招用Python实现批量NDVI归 ...
基于矢量成果从影像提取中深度学习样本库
大数据之:影像提取中深度学习样本库获取的思考话说,虾神一直是做空间统计和数据分析的,对于深度学习这个热门学科,一直以来也就停留在"了解"阶段,虽然这个平展开来,里面比较核心的技术 ...
Python遥感图像处理应用篇(二十二)：Python+GDAL 批量等距离裁剪影像-续
之前写过一篇按照指定行列号数量来进行影像等距离裁剪的博客,链接如下: Python遥感图像处理应用篇(二十二):Python+GDAL 批量等距离裁剪影像_空中旋转篮球的博客-CSDN博客_pytho ...
【ArcGIS风暴】ArcGIS栅格影像（NDVI）归一化处理的两种方法精解
[ArcGIS遇上Python]ArcGIS批量处理栅格影像(NDVI)归一化完整案例代码图像归一化是指对图像进行了一系列标准的处理变换,使之变换为一固定标准形式的过程,该标准图像称作归一化图像. ...
8影像计算ndvi landsat_使用GDAL读取遥感影像的信息
读取影像数据集的元数据 GDAL已经提供了足够方便的函数,可以读取影像的一些元数据信息, 从而方便对数据进行处理.GDAL一般是以字典的形式对元数据进行组织的, 但是对于不同的栅格数据类型,元数据的类 ...

gdal（1）GF-1影像提取NDVI

gdal（1）GF-1影像提取NDVI相关推荐

最新文章

热门文章