1

次のコードを機能させるのに途方に暮れています。何らかの理由で、GeoPandas *.plot() は機能しませんが、いくつかの単純なプロットには Pandas と GeoPandas の両方を使用したいと考えています。

GeoPandas から Shapely オブジェクトを取得して、ベースマップにプロットしようとしています。問題は、ポリゴンがプロットされないことです。GeoPandas.geometry からそれらを反復処理し、それらを軸コレクションに追加してから、plot() を使用しますが、役に立ちません。ベースマップは正常に動作しているようで、コードはエラーを出していませんが、ポリゴン (郡) は表示されません...

お手伝いありがとう！

import geopandas as gpd
from descartes import PolygonPatch
import matplotlib as mpl
import mpl_toolkits.basemap as base
import matplotlib.pyplot as plt

counties_file = r'C:\Users\...\UScounties\UScounties.shp'
counties = gpd.read_file(counties_file)

#new plot
fig = plt.figure(figsize=(5,5),dpi=300)
#ax = fig.add_subplot(111)
ax = ax = plt.gca()

minx, miny, maxx, maxy = counties.total_bounds

#map
m = base.Basemap(llcrnrlon=minx, llcrnrlat=miny,
             urcrnrlon=maxx, urcrnrlat=maxy,
             resolution='h', area_thresh=100000,
             projection='merc')

patches = []

#add polygons
for poly in counties.geometry:
    #deal with single polygons and multipolygons
    if poly.geom_type == 'Polygon':
        p = PolygonPatch(poly, facecolor='blue', alpha=1)
        #plt.gca().add_patch(p)
        #ax.add_patch(p)
        patches.append(p)

    elif poly.geom_type == 'MultiPolygon':
        for single in poly:
            q = PolygonPatch(single,facecolor='red', alpha=1)
            #ax.add_patch(p)
            patches.append(q)

m.drawcoastlines(linewidth=.1)
m.fillcontinents()
m.drawcountries(linewidth=.25,linestyle='solid')
m.drawstates(linewidth=.25,linestyle='dotted')
m.drawmapboundary(fill_color='white')

ax.add_collection(mpl.collections.PatchCollection(patches, match_original=True))
ax.plot()

plt.show()
4

1 に答える 1

3

シェープファイルが正しい投影システムにあるかどうかを確認してください。ベースマップは現在、メルカトル図法に設定されています。その後、それは私のために働いた。

于 2014-04-23T20:30:41.957 に答える