Python：如何将 shapefile 放在一个绘图中的光栅文件顶部，然后将绘图保存为 Jpeg 文件格式

Posted 2023-02-18

【中文标题】Python：如何将 shapefile 放在一个绘图中的光栅文件顶部，然后将绘图保存为 Jpeg 文件格式

【英文标题】：Python : How to place a shapefile on top of raster file in one plot, and then save the plot in a Jpeg file format

【发布时间】：2020-09-10 18:32:09

【问题描述】：

我在网上搜索了三天后发布了这个问题，但没有成功。希望能在这里得到答案。请不要删除帖子，因为我在这里也没有找到答案。谢谢。

我有 2 个文件：

光栅图像文件（即气温 2020-01-01.tif） 世界国家边界 shapefile（（即 World_Countries_base_map.shp）

目标：我想在光栅文件的顶部绘制 shapefile，然后将绘图保存为 Jpeg 文件格式，最终得到这样的结果：

我是 Python 的新手，使用 Spyder 编写了这个简单的代码：

# Import needed packages
import os
import rasterio
import matplotlib.pyplot as plt
import geopandas as gpd
import earthpy as et
from matplotlib import pyplot

## list all raster images in tiff format in the folder:
list_files = [f for f in 
       os.listdir('C:/Users/Desktop/Question/Raster_Air_temp') 
       if '.tif' in f]
print(list_files[1])  # checking the 1st file in the list

## reading the first tiff file:    
raster_image = rasterio.open(list_files[1])

## plot it
draft_output = pyplot.imshow(raster_image.read(1), cmap='jet')

## importing world shapefile
World_map = gpd.read_file('C:/Users/Desktop/Question/World_shapefile/World_Countries_base_map.shp')

# plot World shapefile
fig, ax = plt.subplots(figsize = (30,30))  # image size and quality can be controled by figsize
ax.set_title('The Glob Map', fontsize=50); 
World_map.plot(ax=ax, color='white', edgecolor='black')     # colors note at  https://matplotlib.org/tutorials/colors/colormaps.html
plt.show()

## Plot both World shapefile and raster image in one graph:

????  

但是，如上所示，此代码仅在控制台中为我生成 2 个分隔图。

问题：如何在???中输入正确的代码？实现我的目标的代码部分（如上所述）？ 感谢所有 cmets 和帮助。

在这里，我分享这两个文件是为了方便那些需要帮助的人。 Download the files from my Dropbox

.

【参考方案1】：

由于我无法访问您的数据，因此我使用来自 geopandas 的一些样本数据和一个随机 numpy ndarray 作为 tiff 代理来展示原理。

关键是用 rasterios rasterplot 显示 tiff，不要忘记设置 DEM 的范围！

import rasterio
import numpy as np
from rasterio import plot as rasterplot
import geopandas as gpd
from matplotlib import pyplot as plt


# this is how you'd open the raster dataset if you have one
#tiff = rasterio.open('example.tif')
#tiff_extent = [tiff.bounds[0], tiff.bounds[2], tiff.bounds[1], tiff.bounds[3]]

# i am making this array up
tiff_band_1 = np.random.randint(0, 10, size=(65, 64))
tiff_extent = [4159200.0, 4808100.0, 2828000.0, 3482600.0]

shapefile = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
shapefile = shapefile.to_crs('epsg:3035')
shapefile = shapefile[shapefile.name == 'Germany']

f, ax = plt.subplots()

# plot DEM
rasterplot.show(
    tiff_band_1,  # use tiff.read(1) with your data
    extent=tiff_extent,
    ax=ax,

)
# plot shapefiles
shapefile.plot(ax=ax, facecolor='w', edgecolor='k')
plt.show()