怎么用python计算遥感影像中的NDVI
时间: 2023-12-22 07:04:32 浏览: 171
要用Python计算遥感影像中的NDVI,可以使用以下步骤:
1. 导入需要用到的Python库,如gdal、numpy等。
```python
import gdal
import numpy as np
```
2. 打开需要计算NDVI的遥感影像。
```python
ds = gdal.Open('image.tif')
```
3. 读取红光波段和近红外波段数据,并将其转换为numpy数组格式。
```python
red_band = ds.GetRasterBand(3).ReadAsArray().astype(float)
nir_band = ds.GetRasterBand(4).ReadAsArray().astype(float)
```
4. 计算NDVI值。
```python
ndvi = (nir_band - red_band) / (nir_band + red_band)
```
5. 将计算出的NDVI值保存为遥感影像。
```python
driver = ds.GetDriver()
out_ds = driver.CreateCopy('ndvi.tif', ds) # 创建一个新的遥感影像,用于保存计算出的NDVI值
out_band = out_ds.GetRasterBand(1)
out_band.WriteArray(ndvi)
out_band.SetNoDataValue(-999) # 设置无效值
out_ds = None # 关闭遥感影像
```
以上就是用Python计算遥感影像中的NDVI的步骤。
阅读全文