#acl AdminGroup:read,write,delete,revert  EditorGroup:read,write All:read
#format wiki
#language en
#pragma section-numbers off
= Visualising spatial data =
{{http://www.scipy.org/Cookbook/Matplotlib/Maps?action=AttachFile&do=get&target=basemap3c.png||height="399px",width="532px"}}

High level programming languages like e.g. Python, Matlab, R, IDL allow for the easy generation of maps by using specific modules. Alternatives to produce maps is special commerical or non-commercial software like e.g. ArcGIS, Generic mapping tool (GMT), ...

== Map and projection basics ==

[[EarthShape | The Earth shape and it's approximation ]]



== How to do it in Python ==
There are different ways to generate maps with geographic information within Python. These are either based on the [[http://www.pyngl.ucar.edu/|PyNGL]] ([[http://www.pyngl.ucar.edu/Examples/gallery.shtml|Gallery]]) or [[http://matplotlib.sourceforge.net/basemap/doc/html/|Basemap]] packages. PyNGL is very flexible to handle and allows for the generation of high quality maps. However, it takes a while to generate a nice looking map in PyNGL. If you just want to visualise your data, Basemap might be a very good alternative in many cases.

=== The Basemap way ===
The first step is the generation of an Basemap object which includes also the definition of the projection. Which projection to choose depends on your application. You might want to have a map with specific properties like e.g. equal area or equal angular and optimized for a specific region like e.g. the polar regions.

You can then draw pre-defined data or use your own data.

{{{#!python
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np

#/// define a map projection ///
map = Basemap(projection='ortho', lat_0 = 50, lon_0 = -100,
              resolution = 'l', area_thresh = 1000.)

#to see a list of potential projections, type help(Basemap)
#be aware that each projection has different types of parameters

#... you like coastlines?
map.drawcoastlines()

#... and countries
map.drawcountries()

#... fill continents
map.fillcontinents(color = 'coral')

#... draw coordinates
map.drawmeridians(np.arange(0, 360, 30))
map.drawparallels(np.arange(-90, 90, 30))

}}}




You want to plot some point data?
{{{#!python

clf()
map.bluemarble()

# lat/lon coordinates of five cities.
lats = [40.02, 32.73, 38.55, 48.25, 17.29]
lons = [-105.16, -117.16, -77.00, -114.21, -88.10]
cities=['Boulder, CO','San Diego, CA',
        'Washington, DC','Whitefish, MT','Belize City, Belize']
# compute the native map projection coordinates for cities.
x,y = map(lons,lats)
# plot filled circles at the locations of the cities.
map.plot(x,y,'ro')
# plot the names of those five cities.
for name,xpt,ypt in zip(cities,x,y):
    plt.text(xpt+50000,ypt+50000,name,color='white')


}}}



=== References: ===
 * [[http://www.scipy.org/Cookbook/Matplotlib/Maps|Cookbook / Matplotlib / Maps]]
 * [[http://matplotlib.sourceforge.net/basemap/doc/html/|Matplotlib Basemap Toolkit documentation]]
 * Snyder: Map projections -  a reference manual