Stereografische Projektion¶
Eine winkeltreue Projektion. Die auf der Karte gezeichneten grünen Formen sind flächentreue Kreise auf der Erdoberfläche. Bekannt als „Tissots Indikatrix“, können sie verwendet werden, um die Winkel- und Flächenverzerrung einer Kartenprojektion zu zeigen. Bei einer winkeltreuen Projektion wird die Form der Kreise beibehalten, aber nicht die Fläche. Bei einer flächentreuen Projektion wird die Fläche beibehalten, aber nicht die Form.
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
# setup stereographic basemap.
# lat_ts is latitude of true scale.
# lon_0,lat_0 is central point.
m = Basemap(width=12000000,height=8000000,
resolution='l',projection='stere',\
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("Stereographic Projection")
plt.show()
