如何为 python cartopy 几何设置偏移量？

Question

所以问题是当我将 .shp 文件中的几何图形添加到 cartopy 图时，有一个偏移量，我不知道如何设置偏移量。

我是 python 新手，因此感谢任何帮助。

picture here

import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
#from cartopy.feature import GSHHSFeature 

from cartopy.io.shapereader import Reader


canada_east = -63
canada_west = -123
canada_north = 75
canada_south = 37

standard_parallels = (49, 77)
central_longitude = -(91 + 52 / 60)

data = Reader('icitw_wgs84')

fig = plt.figure()
ax = fig.add_subplot(1, 1, 1,
                     projection=ccrs.LambertConformal(central_longitude=central_longitude,
                                                      standard_parallels=standard_parallels))
ax.set_extent([-79.65, -79.1, 43.57, 43.87])


ax.add_feature(cfeature.LAKES.with_scale('10m'))
ax.add_feature(cfeature.LAND.with_scale('10m'))
ax.add_feature(cfeature.RIVERS.with_scale('10m'))


ax.add_geometries(data.geometries(), crs=ccrs.Geodetic(), edgecolor='k', facecolor='none')

Answer 1

我认为您所看到的是低分辨率的陆地/湖泊数据集。对于这种比例的地图，您最好使用地图图块而不是 NaturalEarth 土地要素。 cartopy 中已有多种选择，Stamen Terrain 或 Open Street Map 可能是不错的选择：

import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs

from cartopy.io.shapereader import Reader
from cartopy.io.img_tiles import StamenTerrain, OSM

standard_parallels = (49, 77)
central_longitude = -(91 + 52 / 60)

data = Reader('citygcs')

fig = plt.figure()
ax = fig.add_subplot(1, 1, 1,                  
                     projection=ccrs.LambertConformal(central_longitude=central_longitude,                                     
                                                      standard_parallels=standard_parallels))
ax.set_extent([-79.65, -79.1, 43.57, 43.87])

tiler = OSM()
ax.add_image(tiler, 10)

ax.add_geometries(data.geometries(), crs=ccrs.Geodetic(), edgecolor='k', 
facecolor='none')
plt.show()

或者使用StamenTerrain:

关于引用椭圆可能还有其他问题（我注意到 shapefile 名称中有 WGS84），这里有一个很好的参考：https://scitools.org.uk/cartopy/docs/v0.16/gallery/effects_of_the_ellipse.html。

如果您的代码示例最少并且所有数据都可用（我必须自己去找一个类似的 shapefile 来重现），这将有所帮助，请参阅此处获取指南：https://stackoverflow.com/help/mcve。