rasterio介绍以及代码

简介

官网链接：https://rasterio.readthedocs.io/en/latest/intro.html

1、原理

与地理空间抽象数据库（GDAL）绑定的，用于访问不同类型栅格数据的包同时为C语言服务，而C语言存在指针问题，不精通C语言的人在使用与gdal绑定的包访问栅格数据时，很容易发生指针悬空问题。
Rasterio的目标是为python量身打造地理空间抽象数据库，——即使用更少的非惯用扩展类和更多的惯用python类型来表达gdal数据模型，同时提高python代码的效率。高性能，较低的认知负担，更清晰，更透明的代码。这就是Rasterio的目的。

2、安装

windows安装输入以下代码即可（二选一）

conda install rasterio

pip install rasterio

3、代码示例

用rasterio打开并显示栅格数据,结果如下

import os 
import gdal
import rasterio
import matplotlib.pyplot as plt

path = r'G:\裁剪\GZQ_clip\20190312.tif'     #文件路径
dataset = rasterio.open(path)               #打开文件
fig = plt.figure(figsize=(16,9))            #画布大小vmax=1,vmin=0,
show((dataset,1),cmap='RdYlBu')

相对于imread而言，rasterio可取读取空间数据，而imread只能读取非空间数据。另附imread读取图形代码

import matplotlib.pyplot as plt
from matplotlib.pyplot import imread
img = imread(r'E:\picture\0G2A0425.jpg')
plt.imshow(img)
plt.show