Lambert Azimuthal Equal Area Projektion¶
Eine flächentreue Projektion. Die auf der Karte gezeichneten grünen Formen sind flächentreue Kreise auf der Erdoberfläche. Bekannt als „Tissots Indikatrix“ können sie zur Darstellung der Winkel- und Flächenverzerrung einer Kartenprojektion verwendet werden. Bei einer winkeltreuen Projektion bleibt die Form der Kreise erhalten, aber nicht die Fläche. Bei einer flächentreuen Projektion bleibt die Fläche erhalten, aber nicht die Form.
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
# setup lambert azimuthal equal area basemap.
# lat_ts is latitude of true scale.
# lon_0,lat_0 is central point.
m = Basemap(width=12000000,height=8000000,
resolution='l',projection='laea',\
lat_ts=50,lat_0=50,lon_0=-107.)
m.drawcoastlines()
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)
poly = m.tissot(lon,lat,1.5,100,\
facecolor='green',zorder=10,alpha=0.5)
plt.title("Lambert Azimuthal Equal Area Projection")
plt.show()
