load_from_gdal

load_from_gdal是gempy.core.grid_modules.topography.Topography包中的一个函数，用于从GDAL数据源加载地形数据。它将地形数据集加载为numpy数组后再进行处理。

语法

load_from_gdal(filename, crs=None, resolution=None, extent=None, mask=None, **kwargs)

参数

• filename：文件名或URL，支持多种GDAL数据源格式，如GeoTiff，HDF5等。
• crs：坐标系统，默认为None。
• resolution：分辨率，默认为None。
• extent：范围，默认为None。
• mask：指定的遮罩文件名或遮罩numpy数组，用于掩盖地形数据。

返回值

返回一个numpy.ndarray，其中包含从GDAL数据源加载的地形数据。

描述

load_from_gdal函数通过GDAL库打开并读取文件或链接中的地形数据，然后将其转换为numpy数组。如果提供了范围或分辨率，则仅读取数据的子集。GeoTiff格式的文件可以自动检测其坐标系统和范围，但其他格式的文件必须明确指定。

示例

import gempy as gp

# 访问GeoTiff数据并加载
geo_data = gp.create_data('Tutorial_ch7')
gp.topography.load_from_gdal(geo_data, filename='my_elevation.tiff')

# 更改坐标系统和范围
gp.topography.load_from_gdal(geo_data, filename='my_elevation.tiff', crs='EPSG:4326', extent=[-125, -116, 32, 42])

# 应用掩码来遮挡数据
gp.topography.load_from_gdal(geo_data, filename='my_elevation.tiff', mask='my_mask.tiff')