栅格:访问地理空间栅格数据¶
地理信息系统使用geotiff和其他格式来组织和存储栅格栅格数据集,如卫星图像和地形模型。Rasterio读取和写入这些格式,并提供基于numpy n维数组和geojson的python API。
下面是一个示例程序,它提取栅格有效数据足迹的geojson形状。
import rasterio
import rasterio.features
import rasterio.warp
with rasterio.open('example.tif') as dataset:
# Read the dataset's valid data mask as a ndarray.
mask = dataset.dataset_mask()
# Extract feature shapes and values from the array.
for geom, val in rasterio.features.shapes(
mask, transform=dataset.transform):
# Transform shapes from the dataset's own coordinate
# reference system to CRS84 (EPSG:4326).
geom = rasterio.warp.transform_geom(
dataset.crs, 'EPSG:4326', geom, precision=6)
# Print GeoJSON shapes to stdout.
print(geom)
程序输出:
{'type': 'Polygon', 'coordinates': [[(-77.730817, 25.282335), ...]]}
Rasterio支持2.7和3.3或更高版本的python。