|
Size: 692
Comment:
|
Size: 2518
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 5: | Line 5: |
| Line 7: | Line 6: |
| {{http://www.scipy.org/Cookbook/Matplotlib/Maps?action=AttachFile&do=get&target=basemap3c.png||height="399px",width="532px"}} | |
| Line 10: | Line 10: |
| == How to do it in Python == | == Map and projection basics == |
| Line 12: | Line 12: |
| '''References:''' * [[http://www.scipy.org/Cookbook/Matplotlib/Maps | Cookbook / Matplotlib / Maps]] * [[http://matplotlib.sourceforge.net/basemap/doc/html/ | Matplotlib Basemap Toolkit documentation]] |
geoid ???? take slides from cartography exercise some words Tissot's Indicatrix Geodetic datum and coordinate transformations |
| Line 17: | Line 19: |
| === Projections === {{attachment:tmerc.png||height="370px",width="494px"}} {{attachment:nplaea.png||height="370px"}} |
|
| Line 18: | Line 22: |
| {{attachment:vandg.png||height="370px"}} {{attachment:sinu.png||height="370px"}} | |
| Line 19: | Line 24: |
| {{attachment:tmerc_tissot.png||height="370px",width="494px"}} Example conformal map | |
| Line 20: | Line 26: |
| {{attachment:nplaea_tissot.png||height="370px"}} Example equal area map == 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. === Define projection === 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. ??? generate all maps with Tissots indicatrix ??? {{http://matplotlib.sourceforge.net/basemap/doc/html/_images/merc.png}} {{http://matplotlib.sourceforge.net/basemap/doc/html/_images/sinu.png}} === Load some vector data === '''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 |
Visualising spatial data
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
geoid ???? take slides from cartography exercise
some words Tissot's Indicatrix
Geodetic datum and coordinate transformations
Projections
Example conformal map
Example equal area map
How to do it in Python
There are different ways to generate maps with geographic information within Python. These are either based on the PyNGL (Gallery) or 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.
Define projection
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.
??? generate all maps with Tissots indicatrix ???
Load some vector data
References:
- Snyder: Map projections - a reference manual