Revision: 5757
http://matplotlib.svn.sourceforge.net/matplotlib/?rev=5757&view=rev
Author: jswhit
Date: 2008-07-13 05:59:45 -0700 (Sun, 13 Jul 2008)
Log Message:
-----------
add tissot's indicatrix to examples to indicate distortion.
Modified Paths:
--------------
trunk/toolkits/basemap/doc/users/figures/aea.py
trunk/toolkits/basemap/doc/users/figures/eqdc.py
trunk/toolkits/basemap/doc/users/figures/laea.py
trunk/toolkits/basemap/doc/users/figures/lcc.py
trunk/toolkits/basemap/doc/users/figures/stere.py
Modified: trunk/toolkits/basemap/doc/users/figures/aea.py
===================================================================
--- trunk/toolkits/basemap/doc/users/figures/aea.py 2008-07-13 12:31:53 UTC
(rev 5756)
+++ trunk/toolkits/basemap/doc/users/figures/aea.py 2008-07-13 12:59:45 UTC
(rev 5757)
@@ -1,18 +1,28 @@
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
+from matplotlib.patches import Polygon
# setup albers equal area conic basemap
# lat_1 is first standard parallel.
# lat_2 is second standard parallel.
# lon_0,lat_0 is central point.
-m = Basemap(width=12000000,height=9000000,
+m = Basemap(width=8000000,height=6000000,
resolution='l',projection='aea',\
- lat_1=45.,lat_2=55,lat_0=50,lon_0=-107.)
+ lat_1=40.,lat_2=60,lon_0=35,lat_0=50)
m.drawcoastlines()
+m.drawcountries()
m.fillcontinents(color='coral',lake_color='aqua')
# draw parallels and meridians.
m.drawparallels(np.arange(-80.,81.,20.))
m.drawmeridians(np.arange(-180.,181.,20.))
m.drawmapboundary(fill_color='aqua')
+# draw tissot's indicatrix to show distortion.
+ax = plt.gca()
+for y in np.linspace(m.ymax/20,19*m.ymax/20,9):
+ for x in np.linspace(m.xmax/20,19*m.xmax/20,12):
+ lon, lat = m(x,y,inverse=True)
+ seg = m.tissot(lon,lat,1.,100)
+ poly = Polygon(seg,facecolor='green',zorder=10,alpha=0.5)
+ ax.add_patch(poly)
plt.title("Albers Equal Area Projection")
plt.savefig('aea.png')
Modified: trunk/toolkits/basemap/doc/users/figures/eqdc.py
===================================================================
--- trunk/toolkits/basemap/doc/users/figures/eqdc.py 2008-07-13 12:31:53 UTC
(rev 5756)
+++ trunk/toolkits/basemap/doc/users/figures/eqdc.py 2008-07-13 12:59:45 UTC
(rev 5757)
@@ -1,6 +1,7 @@
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
+from matplotlib.patches import Polygon
# setup equidistant conic basemap.
# lat_1 is first standard parallel.
# lat_2 is second standard parallel.
@@ -15,5 +16,12 @@
m.drawparallels(np.arange(-80.,81.,20.))
m.drawmeridians(np.arange(-180.,181.,20.))
m.drawmapboundary(fill_color='aqua')
+ax = plt.gca()
+for y in np.linspace(m.ymax/20,19*m.ymax/20,9):
+ for x in np.linspace(m.xmax/20,19*m.xmax/20,12):
+ lon, lat = m(x,y,inverse=True)
+ seg = m.tissot(lon,lat,1.5,100)
+ poly = Polygon(seg,facecolor='green',zorder=10,alpha=0.5)
+ ax.add_patch(poly)
plt.title("Equidistant Conic Projection")
plt.savefig('eqdc.png')
Modified: trunk/toolkits/basemap/doc/users/figures/laea.py
===================================================================
--- trunk/toolkits/basemap/doc/users/figures/laea.py 2008-07-13 12:31:53 UTC
(rev 5756)
+++ trunk/toolkits/basemap/doc/users/figures/laea.py 2008-07-13 12:59:45 UTC
(rev 5757)
@@ -1,6 +1,7 @@
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
+from matplotlib.patches import Polygon
# setup lambert azimuthal equal area basemap.
# lat_ts is latitude of true scale.
# lon_0,lat_0 is central point.
@@ -13,5 +14,13 @@
m.drawparallels(np.arange(-80.,81.,20.))
m.drawmeridians(np.arange(-180.,181.,20.))
m.drawmapboundary(fill_color='aqua')
+# draw tissot's indicatrix to show distortion.
+ax = plt.gca()
+for y in np.linspace(m.ymax/20,19*m.ymax/20,9):
+ for x in np.linspace(m.xmax/20,19*m.xmax/20,12):
+ lon, lat = m(x,y,inverse=True)
+ seg = m.tissot(lon,lat,1.5,100)
+ poly = Polygon(seg,facecolor='green',zorder=10,alpha=0.5)
+ ax.add_patch(poly)
plt.title("Lambert Azimuthal Equal Area Projection")
plt.savefig('laea.png')
Modified: trunk/toolkits/basemap/doc/users/figures/lcc.py
===================================================================
--- trunk/toolkits/basemap/doc/users/figures/lcc.py 2008-07-13 12:31:53 UTC
(rev 5756)
+++ trunk/toolkits/basemap/doc/users/figures/lcc.py 2008-07-13 12:59:45 UTC
(rev 5757)
@@ -1,6 +1,7 @@
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
+from matplotlib.patches import Polygon
# setup lambert conformal basemap.
# lat_1 is first standard parallel.
# lat_2 is second standard parallel (defaults to lat_1).
@@ -18,5 +19,13 @@
m.drawparallels(np.arange(-80.,81.,20.))
m.drawmeridians(np.arange(-180.,181.,20.))
m.drawmapboundary(fill_color='aqua')
+# draw tissot's indicatrix to show distortion.
+ax = plt.gca()
+for y in np.linspace(m.ymax/20,19*m.ymax/20,9):
+ for x in np.linspace(m.xmax/20,19*m.xmax/20,12):
+ lon, lat = m(x,y,inverse=True)
+ seg = m.tissot(lon,lat,1.5,100)
+ poly = Polygon(seg,facecolor='green',zorder=10,alpha=0.5)
+ ax.add_patch(poly)
plt.title("Lambert Conformal Projection")
plt.savefig('lcc.png')
Modified: trunk/toolkits/basemap/doc/users/figures/stere.py
===================================================================
--- trunk/toolkits/basemap/doc/users/figures/stere.py 2008-07-13 12:31:53 UTC
(rev 5756)
+++ trunk/toolkits/basemap/doc/users/figures/stere.py 2008-07-13 12:59:45 UTC
(rev 5757)
@@ -1,6 +1,7 @@
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
+from matplotlib.patches import Polygon
# setup stereographic basemap.
# lat_ts is latitude of true scale.
# lon_0,lat_0 is central point.
@@ -13,5 +14,13 @@
m.drawparallels(np.arange(-80.,81.,20.))
m.drawmeridians(np.arange(-180.,181.,20.))
m.drawmapboundary(fill_color='aqua')
+# draw tissot's indicatrix to show distortion.
+ax = plt.gca()
+for y in np.linspace(m.ymax/20,19*m.ymax/20,9):
+ for x in np.linspace(m.xmax/20,19*m.xmax/20,12):
+ lon, lat = m(x,y,inverse=True)
+ seg = m.tissot(lon,lat,1.5,100)
+ poly = Polygon(seg,facecolor='green',zorder=10,alpha=0.5)
+ ax.add_patch(poly)
plt.title("Stereographic Projection")
plt.savefig('stere.png')
This was sent by the SourceForge.net collaborative development platform, the
world's largest Open Source development site.
-------------------------------------------------------------------------
Sponsored by: SourceForge.net Community Choice Awards: VOTE NOW!
Studies have shown that voting for your favorite open source project,
along with a healthy diet, reduces your potential for chronic lameness
and boredom. Vote Now at http://www.sourceforge.net/community/cca08
_______________________________________________
Matplotlib-checkins mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/matplotlib-checkins
