Debian12虚拟环境下Python-GDAL无法导入gdal_array的问题求助

环境背景 #

• 系统：Debian 12
• 已安装系统级GDAL：

  sudo apt install gdal-bin
  sudo apt-get install libgdal-dev
  

• 系统GDAL版本：GDAL 3.6.2, released 2023/01/02（通过ogrinfo --version验证）
• Python环境：Python 3.11.2虚拟环境，已安装numpy和GDAL

出错场景 #

运行以下Python代码处理BIL文件时，在data = band.ReadAsArray()行报错：

from osgeo import gdal
import matplotlib.pyplot as plt
import numpy as np
    
# 注册EHdr驱动（处理头文件）
gdal.GetDriverByName('EHdr').Register()
    
# 打开BIL文件
filename = 'my_file.bil'
dataset = gdal.Open(filename)

# 读取栅格波段（假设单波段）
band = dataset.GetRasterBand(1)
data = band.ReadAsArray()
# 显示数据
plt.imshow(data, cmap='gray')  # 可调整cmap
plt.colorbar()
plt.show()

错误堆栈：

Traceback (most recent call last):
File "/test.py", line 16, in <module>
data = band.ReadAsArray()
       ^^^^^^^^^^^^^^^^^^
File "/venv/lib/python3.11/site-packages/osgeo/gdal.py", line 4148, in ReadAsArray
from osgeo import gdal_array
File "/venv/lib/python3.11/site-packages/osgeo/gdal_array.py", line 13, in <module>
from . import _gdal_array
ImportError: cannot import name '_gdal_array' from 'osgeo' (venv/lib/python3.11/site-packages/osgeo/__init__.py)

已尝试的无效方法 #

• 使用Conda安装（未成功，转而用venv/pip）
• 先安装numpy再指定版本安装GDAL：pip install gdal==3.6.2
• 通过wheels安装GDAL
• 匹配系统全局numpy版本（卸载虚拟环境中numpy 2.25，安装numpy==1.24.2）

注：系统默认Python（不启用虚拟环境）下代码可正常运行。

解决方法 #

放弃直接使用GDAL，改用rasterio（依赖GDAL）：

1. 安装rasterio：

pip install rasterio

2. 运行以下代码处理BIL文件：

import rasterio

# 打开BIL文件
with rasterio.open('image.bil') as src:
    # 读取栅格数据为NumPy数组
    data = src.read()

    # 获取元数据
    profile = src.profile

    # 处理数据（如显示、分析等）
    print(data)
    print(profile)

内容的提问来源于stack exchange，提问作者melev