|
马上注册,结交更多好友,享用更多功能^_^
您需要 登录 才可以下载或查看,没有账号?立即注册
x
想请问一下报错是不是因为numpy和gdal版本不兼容,如果是的话,我numpy版本应该安装哪一个版本,才能对应GDAL.如果不是因为版本,又因为什么?谢谢
numpy==1.14.5,GDAL-3.4.2
代码:
from tifffile import imread
import numpy as np
from osgeo import gdal, osr
# 读取TIFF格式无人机影像数据
image_path = 'F:\\duneline\\dune\\dune.tif'
output_path = 'overlay_image.tif'
image = imread(image_path)
# 将输入图像转换为灰度图像
gray_image = np.mean(image, axis=2)
# 自适应阈值二值化
binary = np.zeros_like(gray_image, dtype=np.uint8)
window_size = 11 # 窗口大小
k = 0.1 # 控制阈值的参数
for i in range(window_size//2, gray_image.shape[0]-window_size//2):
for j in range(window_size//2, gray_image.shape[1]-window_size//2):
window = gray_image[i-window_size//2:i+window_size//2+1, j-window_size//2:j+window_size//2+1]
threshold = np.mean(window) - k * np.std(window)
if gray_image[i, j] > threshold:
binary[i, j] = 255
# 配置输出的空间参考信息
in_ds = gdal.Open(image_path)
projection = in_ds.GetProjection()
geotransform = in_ds.GetGeoTransform()
# 创建新的输出图像
driver = gdal.GetDriverByName('GTiff')
out_ds = driver.Create(output_path, image.shape[1], image.shape[0], 1, gdal.GDT_Byte)
# 设置输出图像的空间参考信息
out_ds.SetProjection(projection)
out_ds.SetGeoTransform(geotransform)
# 将沙脊线部分写入输出图像
out_band = out_ds.GetRasterBand(1)
out_band.WriteArray(binary)
# 释放资源
out_band.FlushCache()
out_ds = None
in_ds = None
print("保存成功!")
报错:
ModuleNotFoundError: No module named 'numpy.core._multiarray_umath'
Traceback (most recent call last):
File "F:/duneline/duneline.py", line 143, in <module>
out_band.WriteArray(binary)
File "C:\安装软件\Python\lib\site-packages\osgeo\gdal.py", line 3668, in WriteArray
from osgeo import gdal_array
File "C:\安装软件\Python\lib\site-packages\osgeo\gdal_array.py", line 13, in <module>
from . import _gdal_array
ImportError: numpy.core.multiarray failed to import
Process finished with exit code 1
报错信息提示了一个模块导入错误:ModuleNotFoundError: No module named 'numpy.core._multiarray_umath'。这个错误是由于numpy的版本问题导致的。在你的代码中,你使用的numpy版本是1.14.5,而GDAL的版本是3.4.2。这两个版本之间可能存在不兼容的情况。
解决这个问题的方法是将numpy降级到与GDAL兼容的版本。根据GDAL的官方文档,GDAL 3.4.2与numpy 1.20.0兼容。所以你可以尝试将numpy升级到1.20.0版本。
你可以使用以下命令来升级numpy:
- pip install numpy==1.20.0
复制代码
安装完成后,重新运行你的代码,应该就能解决这个问题了。如果还有其他问题,请随时提问。
|
|