开发者社区

文档建议反馈控制台

最新优惠活动

文章/答案/技术大牛

发布

GDAL ReadAsArray()一直返回None

GDAL是地理空间数据抽象库（Geospatial Data Abstraction Library）的简称，它是一个开源的库，用于读取、写入和处理各种地理空间数据格式。GDAL提供了一系列的API和工具，可以在不同的编程语言中使用，包括C/C++、Python等。

在GDAL中，ReadAsArray()是一个函数，用于从栅格数据集中读取数据，并返回一个NumPy数组。然而，当ReadAsArray()一直返回None时，可能有以下几个原因：

数据集路径错误：首先，需要确保提供给ReadAsArray()函数的数据集路径是正确的。可以检查路径是否包含正确的文件名、文件格式以及文件的存储位置。
数据集打开失败：如果数据集无法成功打开，ReadAsArray()函数将返回None。这可能是由于数据集文件损坏、文件格式不受支持或者权限问题导致的。可以尝试使用其他工具或方法打开数据集，以确定是否存在问题。
数据集中没有有效的数据：如果数据集中没有有效的数据，ReadAsArray()函数也会返回None。可以使用其他工具或方法检查数据集中是否存在有效的数据。
数据集的波段索引错误：ReadAsArray()函数可以指定要读取的波段索引，如果指定的波段索引超出了数据集的范围，函数将返回None。可以检查波段索引是否正确，并确保在读取数据之前了解数据集中的波段信息。

总结起来，当GDAL的ReadAsArray()函数一直返回None时，需要检查数据集路径、数据集打开情况、数据集中是否存在有效数据以及波段索引是否正确。如果问题仍然存在，可以考虑使用其他方法或工具来读取和处理地理空间数据。

腾讯云相关产品和产品介绍链接地址：

腾讯云地理空间计算服务：https://cloud.tencent.com/product/gis
腾讯云云服务器（CVM）：https://cloud.tencent.com/product/cvm
腾讯云对象存储（COS）：https://cloud.tencent.com/product/cos
腾讯云人工智能（AI）：https://cloud.tencent.com/product/ai
腾讯云物联网（IoT）：https://cloud.tencent.com/product/iot
腾讯云区块链（BCS）：https://cloud.tencent.com/product/bcs
腾讯云云原生应用引擎（TKE）：https://cloud.tencent.com/product/tke

请注意，以上链接仅供参考，具体产品选择应根据实际需求和情况进行评估。

页面内容是否对你有帮助？

有帮助

没帮助

相关·内容

打开栅格数据的正确方式

(band.DataType)}') # DataType属性返回的是数字 print(f'NoData值：{band.GetNoDataValue()}') # 很多影像都是NoData，...我们在做数据处理时要特别对待 print(f'统计值（最大值最小值）：{band.ComputeRasterMinMax()}') # 有些数据本身就存储了统计信息，有些数据没有需要计算 # 关闭数据集 ds = None...函数直接进行（相当于第一种转换） import gdal # 打开栅格数据集 ds = gdal.Open('example.tif') # 在数据集层面转换 image = ds.ReadAsArray...() print(f'数据的尺寸：{image.shape}') # 输出结果为：数据的尺寸：(3, 4800, 4800) # 这说明ReadAsArray方法将每个波段都转换为了一个二维数组 #...ds = None 输出结果：数据的尺寸：(3, 4800, 4800) 波段大小：(4800, 4800) 波段大小：(4800, 4800) 波段大小：(4800, 4800) 使用gdal_array

1.1K2 0

一堆的遥感影像里求平均（逐个像元），求几个平均。。。

cols = ds.RasterXSize rows = ds.RasterYSize band = ds.GetRasterBand(1) data = band.ReadAsArray...cols = ds.RasterXSize rows = ds.RasterYSize band = ds.GetRasterBand(1) data = band.ReadAsArray...def Readxy(RasterFile): #读取每个图像的信息 ds = gdal.Open(RasterFile,GA_ReadOnly) if ds is None: print ('Cannot...RasterFile sys.exit(1) cols = ds.RasterXSize rows = ds.RasterYSize band = ds.GetRasterBand(1) data = band.ReadAsArray...rows = ds.RasterYSize band = ds.GetRasterBand(1) data = band.ReadAsArray(0,0,cols,rows) return data

1.4K3 0

对Python3+gdal 读取tiff格式数据的实例讲解

1、遇到的问题：numpy版本 im_data = dataset.ReadAsArray(0,0,im_width,im_height)#获取数据这句报错升级numpy：pip install -...U numpy 但是提示已经是最新版本解决：卸载numpy 重新安装 2.直接从压缩包中读取tiff图像参考：http://gdal.org/gdal_virtual_file_systems.html...#gdal_virtual_file_systems_vsizip 当前情况是2层压缩： /’/vsitar/C:/Users/summer/Desktop/a_PAN1.tiff’ 3.读tiff def...= gdal.Open(fileName) if dataset == None: print(fileName+ "掩膜失败，文件无法打开") return im_width = dataset.RasterXSize...您可能感兴趣的文章: 利用python GDAL库读写geotiff格式的遥感影像方法在python中利用GDAL对tif文件进行读写的方法浅谈python下tiff图像的读取和保存方法用python

2.1K2 0

栅格数据创建与保存

思路与方法使用Python进行栅格数据处理，很多时候，我们会将GDAL的Dataset对象转化为NumPy的ndarray对象，这样我们可以使用很多通用的Python库对数据进行处理，然后再借助GDAL..., xsize=None, ysize=None, transform=None, projection=None,...: dtype = gdal.GDT_Float32 if xsize is None: xsize = np_array.shape[-...ds.RasterXSize # 图像列数 y_size = ds.RasterYSize # 图像行数 proj = ds.GetProjection() # 返回的是WKT格式的字符串 trans...= ds.GetGeoTransform() # 返回的是六个参数的tuple # 在数据集层面ReadAsArray方法将每个波段都转换为了一个二维数组 image = ds.ReadAsArray

1.6K1 1

【Python】GDAL基本操作遥感大图显示

/~gohlke/pythonlibs/#gdal 有些老版本gdal的引入方式是直接import： import gdal 新版本的gdal引入方式如下： from osgeo import gdal...：(0,1,0,0,0,1) 空间参照系统信息 ·GetProjection()方法返回栅格数据的坐标转换参数，示例： from osgeo import gdal data = gdal.Open(...()函数返回(rows, columns) 对于多波段栅格数据，ReadAsArray()函数返回(bands, rows, columns) 按块读取栅格 ReadAsArray同样支持按块读取栅格信息...，即读取部分区域图像信息，示例： from osgeo import gdal data = gdal.Open("xdu.tif") data_array = data.ReadAsArray(xoff...class Tiles: def __init__(self, dataset, size=1000, image_type=None): """ 初始化影像瓦片

2.5K3 1

Python自动计算大量遥感数据的NDVI

band_red = dataset.GetRasterBand(3) data_red = band_red.ReadAsArray() data_red...= data_red.astype(float) band_nir = dataset.GetRasterBand(4) data_nir = band_nir.ReadAsArray...output_dataset = None print(filename, "finished!") ...dataset.GetRasterBand()用以获取指定的栅格波段，而band.ReadAsArray()则将波段数据读取为数组；同时，我这里还用了astype()转换数组的格式，避免原本遥感影像的数据格式带来的问题...同时，需要清理和关闭数据集，将数据集和输出数据集设置为None以释放资源。还可以打印文件名和finished!，表示当前文件处理完成。

1321 0

Python | GDAL处理影像

支持的数据类型 # 导入gdal，注意导入的名称 import os from osgeo import gdal #或者直接用import gdal from matplotlib import pyplot...offsets = gdal.ApplyGeoTransform(inv_gt, 465200, 5296000) xoff, yoff = map(int, offsets) value = band.ReadAsArray...) # print(success, inv_gt) # GDAL 2.x: You get the geotransform or None inv_gt = gdal.InvGeoTransform...[行, 列]，与GDAL的恰恰相反，GDAL为[列，行]！...[1] else: raise RuntimeError('Inverse geotransform failed') elif inv_gt is None: raise

4.4K3 1

Pymysql cur.fetchall() 返回 None

该函数返回一个包含查询结果的元组列表。如果 cur.fetchall() 返回 None，可能是由于以下多种问题导致的。...预期在下一个循环中，cur.fetchall() 方法应该返回一个包含所有表名的元组，其中应该包括新创建的表。但实际上，cur.fetchall() 方法却返回了 None。...version sql += "INSERT INTO %s VALUES ('HFRHWHARDWAREVERSION',结论检查 SQL 语句是否成功执行：确保查询是正确的，并且确实返回了结果...通过这些步骤，我们可以排查 pymysql 中 cur.fetchall() 返回 None 的问题。

1081 0

SVM检测影像的云

img_ds.GetRasterBand(1).DataType)) for b in range(img.shape[2]): img[:, :, b] = img_ds.GetRasterBand(b + 1).ReadAsArray...num_leaves=31, max_depth=-1, learning_rate=0.1, n_estimators=100, subsample_for_bin=200000, objective=None..., class_weight=None, min_split_gain=0.0, min_child_weight=0.001, min_child_samples=20, subsample=1.0,...subsample_freq=0, colsample_bytree=1.0, reg_alpha=0.0, reg_lambda=0.0, random_state=None, n_jobs=-1,.../cloubtest.tif') band = ds.GetRasterBand(2) arr = band.ReadAsArray() x0, dx, dxdy, y0, dydx, dy = ds.GetGeoTransform

4483 0

Basemap系列教程：读取WRF模式数据

+'":XLONG') ds_lat = gdal.Open('NETCDF:"'+wrf_out_file+'":XLAT') ds_t2 = gdal.Open('NETCDF:"'+wrf_out_file...map.contourf(ds_lon.ReadAsArray()[1], ds_lat.ReadAsArray()[1], ds_t2.ReadAsArray()[1]) map.drawcoastlines...注意 GDAL 是如何读取 NetCDF文件的。...x, y = map(ds_lon.ReadAsArray()[1], ds_lat.ReadAsArray()[1]) map.contourf(x, y, ds_t2.ReadAsArray()[...x, y = map(ds_lon.ReadAsArray()[1], ds_lat.ReadAsArray()[1]) u = ds_u.ReadAsArray()[1] v = ds_v.ReadAsArray

1.9K2 0

MOD043km的hdf数据转为tif

MOD043km的hdf数据转为tifc尝试了两种办法，1是靠 gdal.open(r'D:/Thesis/ML/modis3km/MOD04_3K.A2018001.0320.061.2018003202214....hdf') 代码如下： import gdal, osr import numpy as np import os import rasterio from rasterio.warp import...import convertModisGDAL ds = gdal.Open(r'D:/Thesis/ML/modis3km/MOD04_3K.A2018001.0320.061.2018003202214....hdf') sub=ds.GetSubDatasets() aod=gdal.Open(sub[10][0]).ReadAsArray()##这步有个问题aod=gdal.Open(sub[10])....ReadAsArray()的时候立刻报错，还是要[10][0]才可以，后续的问题不知道是不是这步有毛病 band = ds.GetRasterBand(1)#2 #arr = band.ReadAsArray

1.9K4 0

Basemap系列教程：绘图

ds = gdal.Open(".....为了看到结果，linewidth 不能是 None linewidths 六边形的线宽，默认为 None，即不绘制边界 alpha 表示 layer的透明度注意：旧版本的库不支持 hexbin...ds = gdal.Open(".....ds = gdal.Open(".....text(x, y, s, fontdict=None, withdash=False, **kwargs) text 方法不属于 Basemap，而是属于 matplotlib。

4.3K1 0

GDAL对缺失投影定义的AIG文件根据经纬度坐标提取像元值

文件直接在上述教程进行测试发现能够顺利读取AIG，但是根据正确坐标返回的坐标为像素值为空（或者在行列计算时就不存在），思考该问题应该是投影系统出现了问题。...gcs = None shape = None # gdal.AllRegister() dataset = gdal.Open(in_file_path)...): ''' 根据单个图像坐标，或者依据GDAL的六参数模型将给定的投影、地理坐标转为影像图上坐标后，返回对应像元的像素值 :param file_path: 图像的文件路径...): ''' 根据多个图像坐标，或者依据GDAL的六参数模型将给定的投影、地理坐标转为影像图上坐标后，返回对应像元的像素值 :param file_path: 图像的文件路径..., gcs, pcs, extend, shape = get_file_info(file_path, prj_config) img = dataset.GetRasterBand(1).ReadAsArray

1.8K0 0

GDAL计算平均值的图（tif）

ds.GetRasterBand(1).SetNoDataValue(noDataValue) ds.GetRasterBand(1).WriteArray(data) ds = None...(RasterFile,GA_ReadOnly) if ds is None: print ('Cannot open '),RasterFile sys.exit...cols = ds.RasterXSize rows = ds.RasterYSize band = ds.GetRasterBand(1) data = band.ReadAsArray...(RasterFile,GA_ReadOnly) if ds is None: print ('Cannot open '),RasterFile sys.exit...cols = ds.RasterXSize rows = ds.RasterYSize band = ds.GetRasterBand(1) data = band.ReadAsArray

1.9K2 0

利用Python对Landsat进行去云

下面我们看一下具体代码： import numpy as np from osgeo import gdal out_path='D:\DATA\cloud_mask.tif' QA_band =...gdal.Open(r'D:\DATA\LC09_L2SP_055016_20220703_20220705_02_T1\LC09_L2SP_055016_20220703_20220705_02_T1...in_band = red_band.GetRasterBand(1) x_size = in_band.XSize y_size = in_band.YSize QA_array = QA_band.ReadAsArray...().astype(np.float) red_array = red_band.ReadAsArray().astype(np.float) #识别云及阴影位置 cloud_mask=(np.bitwise_and...()) outBand = outimg.GetRasterBand(1) outBand.WriteArray(out_red_array) outBand.FlushCache outimg = None

1.5K2 0

一个是idw插值（idw源码），另一个是补缺插值~其实还是一个东西

=None, leafsize=10): if not X is None: self.tree = cKDTree(X, leafsize=leafsize ) if not z is None: self.z...= np.array(z) def fit(self, X=None, z=None, leafsize=10): return self....(r'D:\Thesis\samrast1.tif') bandg = ds.GetRasterBand(1) elevationg = bandg.ReadAsArray() [cols, rows]...= elevationg.shape format = "GTiff"#5 driver = gdal.GetDriverByName(format)#6 pmkrigrid=[] pmdata=np.array...kriging.predict(test_data=testpoint, l=.5, sigma=.2) #pre=list(predict) #pre.reverse()#这里面有个顺序问题，我一直没搞懂

6243 0

Python辐射校正遥感图像并以一列的形式导出Excel

本文介绍基于Python语言中的gdal模块，读取一景.tif格式的栅格遥感影像文件，提取其中每一个像元的像素数值，对像素值加以计算（辐射定标）后，再以一列数据的形式将计算后的各像元像素数据保存在一个...(file_path, gdal.GA_ReadOnly) band = dataset.GetRasterBand(1) data = band.ReadAsArray() dataset = None...dataset = gdal.Open(file_path, gdal.GA_ReadOnly)意味着我们以只读方式打开遥感影像文件，并将返回的Dataset对象赋值给变量dataset；随后，获取第一个波段的像元值...，这可以通过band = dataset.GetRasterBand(1)来完成（需要注意，这里波段编号的索引是从1开始的）；随后，data = band.ReadAsArray()意思是将波段的像元值读取为一个二维数组...随后，我们需要关闭遥感影像文件，通过将dataset变量设为None，释放对遥感影像文件的引用，从而关闭文件。其次，我们对像元值进行处理。

1481 0

python实现3D地图可视化

库我们使用了多个库： 1.gdal；主要是用于读取地图信息，这个库在GIS中很常用，使用C++代码编写的，如果安装不了需要在pypi里面找一下对应的资源。...from osgeo import gdal import cv2 gdal.UseExceptions() ds = gdal.Open('E:/Pythoncode/读取地理信息/无标题.tif'...) bandg = ds.GetRasterBand(1) elevationg = bandg.ReadAsArray() bandr = ds.GetRasterBand(2) elevationr...= bandr.ReadAsArray() bandb = ds.GetRasterBand(3) elevationb = bandb.ReadAsArray() import matplotlib.pyplot

1.8K1 0

为什么我的递归函数返回None

Type "a" or "b": a got input: None 我不明白为什么 get_input() 函数返回的是 None，因为它本应只返回 my_var。这个 None 是从哪里来的？...答：它返回 None 是因为当你递归调用它时: if my_var != "a" and my_var !...没有返回那个值。因此，尽管递归确实发生了，但返回值却被丢弃了，然后你会从函数末尾退出。...在函数末尾退出意味着 Python 会隐式地返回 None，就像下面这样： >>> def f(x): ......Python3 documentation 因此，除了在 if 语句中调用 get_input() 之外，还需要返回递归调用返回的内容。

1341 0

利用gdal、rasterio将modis文件进行格式转换、投影转换

("D:\Thesis\ML\maiac_02.tif")#1 band = ds2.GetRasterBand(1)#2 arr = band.ReadAsArray()#3， [cols, rows...() aod55=gdal.Open(sub[1][0]).ReadAsArray() if aod55.ndim == 4: a0=aod55[0,:]#提取第一个波段...("D:/minxinan/data/temp/ndvitemp.tif")#1 band = ds2.GetRasterBand(1)#2 arr = band.ReadAsArray()#3， [cols...ds.GetSubDatasets() ndvi=gdal.Open(sub[0][0]).ReadAsArray() outDataRaster = driver.Create(outp...('D:/minxinan/data/NDVI/ndvihdf/'+i) sub=ds.GetSubDatasets() aod=gdal.Open(sub[0][0]).ReadAsArray

2K3 0

点击加载更多

扫码

添加站长进交流群

领取专属 10元无门槛券

手把手带您无忧上云

扫码加入开发者社群

相关资讯

热门标签

活动推荐

运营活动

活动名称

广告关闭