पाइथन बेसमैप

Question

में देशों को भरें हाय मैं कुछ निश्चित देशों में भरे कुछ देशों के साथ अजगर बेसमैप का उपयोग करके एक नक्शा प्लॉट करने की कोशिश कर रहा हूं।

क्या वहां कोई त्वरित और आसान समाधान है ??

Answer 1

जैसा कि @unutbu द्वारा पहले ही कहा जा चुका है, थॉमस 'here ठीक है जो आप के बाद हैं।

import cartopy.crs as ccrs 
import matplotlib.pyplot as plt 
import cartopy.io.shapereader as shpreader 
import itertools 
import numpy as np 

shapename = 'admin_0_countries' 
countries_shp = shpreader.natural_earth(resolution='110m', 
             category='cultural', name=shapename) 

# some nice "earthy" colors 
earth_colors = np.array([(199, 233, 192), 
           (161, 217, 155), 
           (116, 196, 118), 
           (65, 171, 93), 
           (35, 139, 69), 
           ])/255. 
earth_colors = itertools.cycle(earth_colors) 



ax = plt.axes(projection=ccrs.PlateCarree()) 
for country in shpreader.Reader(countries_shp).records(): 
    print country.attributes['name_long'], earth_colors.next() 
    ax.add_geometries(country.geometry, ccrs.PlateCarree(), 
         facecolor=earth_colors.next(), 
         label=country.attributes['name_long']) 

plt.show() 

Answer 2

pelson से जवाब से प्रेरित होकर, मैं पोस्ट:

आप Cartopy, इसी कोड (v0.7 में) के साथ ऐसा करने http://scitools.org.uk/cartopy/docs/latest/tutorials/using_the_shapereader.html थोड़ा से अनुकूलित किया जा सकता इच्छुक हों मेरे पास समाधान है। मैं इसे आपके पास छोड़ दूंगा जो सबसे अच्छा काम करता है, इसलिए मैं इस समय कोई जवाब स्वीकार नहीं करूंगा।);

#! /usr/bin/env python 

import sys 
import os 
from pylab import * 
from mpl_toolkits.basemap import Basemap 
import matplotlib as mp 

from shapelib import ShapeFile 
import dbflib 
from matplotlib.collections import LineCollection 
from matplotlib import cm 

def get_shapeData(shp,dbf): 
    for npoly in range(shp.info()[0]): 
    shpsegs = [] 
    shpinfo = [] 

    shp_object = shp.read_object(npoly) 
    verts = shp_object.vertices() 
    rings = len(verts) 
    for ring in range(rings): 
     if ring == 0: 
      shapedict = dbf.read_record(npoly) 
     name = shapedict["name_long"] 
     continent = shapedict["continent"] 
     lons, lats = zip(*verts[ring]) 
     if max(lons) > 721. or min(lons) < -721. or max(lats) > 91. or min(lats) < -91: 
      raise ValueError,msg 
     x, y = m(lons, lats) 
     shpsegs.append(zip(x,y)) 
     shapedict['RINGNUM'] = ring+1 
     shapedict['SHAPENUM'] = npoly+1 
     shpinfo.append(shapedict) 

    lines = LineCollection(shpsegs,antialiaseds=(1,)) 
    lines.set_facecolors(cm.jet(np.random.rand(1))) 
    lines.set_edgecolors('k') 
    lines.set_linewidth(0.3) 
    ax.add_collection(lines) 


if __name__=='__main__': 

    f=figure(figsize=(10,10)) 
    ax = plt.subplot(111) 
    m = Basemap(projection='merc',llcrnrlat=30,urcrnrlat=72,\ 
      llcrnrlon=-40,urcrnrlon=50,resolution='c') 
    m.drawcountries(linewidth=0.1,color='w') 

    sfile = 'ne_10m_admin_0_countries' 

    shp = ShapeFile(sfile) 
    dbf = dbflib.open(sfile) 
    get_shapeData(shp,dbf) 

    show() 
    sys.exit(0) 

यह यहाँ मेरी उदाहरण है कि कैसे अल्बानिया में सही रंग (बहुत ही सुंदर नहीं मुझे पता है) को भरने के लिए परिणाम

है।

#HACK for Albania 
    shpsegs = [] 
    shpinfo = [] 

    shp_object = shp.read_object(9) 
    verts = shp_object.vertices() 
    rings = len(verts) 
    for ring in range(rings): 
     if ring == 0: 
      shapedict = dbf.read_record(9) 
     name = shapedict["name_long"] 
     continent = shapedict["continent"] 
     lons, lats = zip(*verts[ring]) 
     if max(lons) > 721. or min(lons) < -721. or max(lats) > 91. or min(lats) < -91: 
      raise ValueError,msg 
     x, y = m(lons, lats) 
     shpsegs.append(zip(x,y)) 
     shapedict['RINGNUM'] = ring+1 
     shapedict['SHAPENUM'] = npoly+1 
     shpinfo.append(shapedict) 
    lines = LineCollection(shpsegs,antialiaseds=(1,)) 
    if name == 'Albania': 
    lines.set_facecolors('w') 
    lines.set_edgecolors('k') 
    lines.set_linewidth(0.3) 
    ax.add_collection(lines) 

यह महत्वपूर्ण है कि आप अन्य सभी आकारों के बाद ऐसा करते हैं। शायद आप इस कोड के कुछ हिस्सों से छुटकारा पा सकते हैं, लेकिन जैसा कि मैंने कहा था कि यह मेरे लिए पर्याप्त था।

नाम या महाद्वीप द्वारा अपने आवेदन मैं रंग contries के लिए, इसलिए इन पंक्तियों:

name = shapedict["name_long"] 
    continent = shapedict["continent"] 

डेटा का उपयोग किया है कि मैं इस वेबसाइट से मिले: http://www.naturalearthdata.com/