Albers-Flächentreue Projektion¶
Eine flächentreue Projektion. Die auf der Karte gezeichneten grünen Formen sind flächentreue Kreise auf der Erdoberfläche. Bekannt als "Tissot'sche Indikatrix", können sie verwendet werden, um die Winkel- und Flächenverzerrung einer Kartenprojektion zu zeigen. 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.
Die Verzerrung ist bei dieser Projektion in der Nähe der Pole sehr groß.
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
# 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=8000000,height=7000000,
resolution='l',projection='aea',\
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,10):
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.25,100,\
facecolor='green',zorder=10,alpha=0.5)
plt.title("Albers Equal Area Projection")
plt.show()
