【已解决】使用cartopy读取shp，显示空白

gaishibo60

大神们，求助，我使用cartopy读取shp并画图但是缺显示空白，请问这是为什么呢？附件是使用的shp地图（百度链接：https://pan.baidu.com/s/1IwBI8u66Ek3lspDxZBOMrQ    提取码：zqty），使用arcgis可正常打开显示。

我所使用的代码如下，谢谢各位大神！


import cartopy.crs as ccrs
import cartopy.io.shapereader as shpreader
from cartopy.feature import ShapelyFeature
import matplotlib.pyplot as plt
# define map extent
lllat=0
urlat=60
lllon=70
urlon=140
# Set up cartopy instance
m = ccrs.Mercator()
# m = ccrs.PlateCarree()
ax = plt.axes(projection=m)
plt.gcf().set_size_inches(20, 10)
extent=[lllon,urlon,lllat,urlat]
ax.set_extent(extent)
reader = shpreader.Reader(r'C:\Users\G\Desktop\shp\国界线.shp')
# reader = shpreader.Reader(r'C:\Users\G\Desktop\shp\province.shp')
shape_feature = ShapelyFeature(reader.geometries(), ccrs.PlateCarree(), facecolor="w",
                               edgecolor='black', lw=1)
ax.add_feature(shape_feature)
plt.savefig("图片.png")

1099221723

set_extent 没换投影

1099221723

地图也不是经纬度的

rqq

把坐标系转换为WGS84

gaishibo60

1099221723 发表于 2022-6-6 09:00
地图也不是经纬度的

谢谢大佬，已解决

gaishibo60

rqq 发表于 2022-6-6 09:37
把坐标系转换为WGS84

谢谢大佬，已解决

gaishibo60

1099221723 发表于 2022-6-6 08:55
set_extent 没换投影

谢谢大佬，已解决

Leonf

感谢大佬！！遇到了相同的问题，转换为WGS1984后就解决了

聆念

请问怎么解决的呀？方便不方便贴代码出来学习一下