1

我正在尝试将 Coastlines() 地图的投影与 shapefile 的投影合成，其 .prj 文件显示：

GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",
SPHEROID["WGS_1984",6378137.0,298.257223563]],
PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]

我的尝试是：

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
from cartopy.io import shapereader

# set up a map with coastlines around Auckland:
plt.figure(figsize=(10, 10))
platecarree = ccrs.PlateCarree(globe=ccrs.Globe(datum='WGS84'))

ax = plt.axes(projection=platecarree)
extent = [174.25, 175.25, -37.5, -36.5]
ax.set_extent(extent)
ax.coastlines('10m',color='red')

# read in shapefile and plot the polygons:
shp2 = shapereader.Reader('auckland_geology_wgs84gcs.shp')
formations = shp2.records()

for formation in formations:
     # plot water blue, and all other rocks yellow
     if formation.attributes['MAIN_ROCK'] == b'                                ':
         ax.add_geometries(formation.geometry, ccrs.PlateCarree(),facecolor='blue',alpha=.1)
     else:
         ax.add_geometries(formation.geometry, ccrs.PlateCarree(), facecolor='yellow',alpha=.1)
plt.show()

我尝试在我的 platecarree 定义中为地球参数提供半径和 prj 文件的反向展平，但如果我设置甚至更改这些数字，我没有看到输出有任何变化。

此外，使用定义的“platecarree”投影（使用 WGS84 调用地球）作为 add_geometries 调用中的 crs，我的输出为空白。

照原样，结果在我看来就像投影不匹配

4

1 回答 1

1

我尝试使用 QGIS 和从 Natural Earth（10m 海岸线）和 GADM（NZ adm0 级别）下载的数据来重现您的问题。看起来 NE10m 海岸线是罪魁祸首！GADM 与您的地质层完美对齐，而 NE10m 关闭（并变形）。带有地质图和海岸线的 QGIS 屏幕截图

于 2016-07-04T10:01:45.510 回答