Re: [Matplotlib-users] Draw paths on map using matplotlib-basemap

[Thomas L. <thl...@ms...>]

hi,

You could just add a scatter point at the first x[0],y[0] of each group, maybe annotated with a name label ?

note, for efficiency purposes, you should label the axes outside of your loop (+ set the title too), you only need to do this once!

Cheers,

Thomas

**********************
Dr Thomas Lecocq 
Geologist
Royal Observatory of Belgium

- Seismology -
**********************


From: gl...@co...
To: mat...@li...
Date: Tue, 12 Feb 2013 19:47:19 +0000
Subject: [Matplotlib-users] Draw paths on map using matplotlib-basemap
















Thank you for your ideas, I leave the correct code to plot trajectories of any object, in my case I have drawn the trajectories of convective storms.
# --- Construimos el mapa ---
import numpy as npfrom mpl_toolkits.basemap import Basemapimport matplotlib.pyplot as pltfrom PIL import *fig = plt.figure(figsize=(12,12))
ax = fig.add_axes([0.1,0.1,0.8,0.8])
m = Basemap(projection='cyl', llcrnrlat=12, urcrnrlat=35,llcrnrlon=-120, urcrnrlon=-80, resolution='c', area_thresh=1000.)
m.bluemarble()m.drawcoastlines(linewidth=0.5)m.drawcountries(linewidth=0.5)m.drawstates(linewidth=0.5)
# --- Dibujamos paralelos y meridianos ---
m.drawparallels(np.arange(10.,35.,5.),labels=[1,0,0,1])m.drawmeridians(np.arange(-120.,-80.,5.),labels=[1,0,0,1])m.drawmapboundary(fill_color='aqua')
# --- Abrimos el archivo que contiene los datos ---
import pandas as pd
df = pd.read_csv('scm-2004.csv')for evento, group in df.groupby(['evento']):    latitude = group.lat.values    longitude = group.lon.values    x,y = m(longitude, latitude)    plt.plot(x,y,'y-',linewidth=2 )    plt.xlabel('Longitud')    plt.ylabel('Latitud')    plt.title('Trayectorias de Sistemas Convectivos 2004')


plt.savefig('track-2004.jpg', dpi=100)
With the above code, I get
the desired figure. 60
paths drawn on 
the map of México.


I have only one last question: how could indicate the start of each of the storms, someone has an idea how I can do this?











------------------------------------------------------------------------------
Free Next-Gen Firewall Hardware Offer
Buy your Sophos next-gen firewall before the end March 2013 
and get the hardware for free! Learn more.
http://p.sf.net/sfu/sophos-d2d-feb
_______________________________________________
Matplotlib-users mailing list
Mat...@li...
https://lists.sourceforge.net/lists/listinfo/matplotlib-users