使用光谱 python 读取 .img envi 文件

Posted 2023-03-12

【中文标题】使用光谱 python 读取 .img envi 文件

【英文标题】：reading .img envi file with spectral python

【发布时间】：2017-12-24 00:23:28

【问题描述】：

我正在尝试打开一个env .img 文件并且有一个同名的.hdr 文件。在 .img 文件中有两个图像，我可以使用以下代码读取它们。

from spectral import *
img = open_image('LC08_L1TP_029029_20130330_20170310_01_T1_sensor_B05.img')

img(BSQ文件)的属性如下图

In[352] img

Out[352]: 

Data Source:   '.\LC08_L1TP_029029_20130330_20170310_01_T1_sensor_B05.img'
# Rows:           7311
# Samples:        7371
# Bands:             2
Interleave:        BSQ
Quantization:  16 bits
Data format:     int16

我想从 img 中提取这两个图像。但是当我尝试使用时

img[:,:,1]

它给了我一个大小为 (7311,7371,1) 的数组，但数组中的所有值都为零，但我知道它们应该是非零值。

我的问题是如何从 BSQ 文件中提取这两个图像？

【参考方案1】：

你可以试试这个变种：

from spectral import * 

img = open_image('LC08_L1TP_029029_20130330_20170310_01_T1_sensor_B05.img') 
img_open = img.open_memmap(writeable = True) 
img_band = img_open[:,:,1] 
envi.save_image('new_image.bsq', ext='bsq', interleave = 'BSQ')

或

这个变种需要通过 hdr 文件打开图像。但它应该像以前的变体一样工作。

from spectral import *
img = envi.open('<name of hdr file>') 
img_open = img.open_memmap(writeable = True) 
img_band = img_open[:,:,1] 
envi.save_image('new_image.bsq', ext='bsq', interleave = 'BSQ')

【参考方案2】：

您可以从同名的 hdr 文件中读取 envi 图像。

import numpy as np
from spectral import*    
img1=open_image("<path to file.hdr>").read_band(0)
img2=open_image("<path to file.hdr>").read_band(1)

现在您已经提取了两个波段到 img1 和 img2，您可以根据自己的判断保存或显示它们。