from mpl_toolkits.basemap import Basemap
from pylab import title, show, arange, figure

# map with continents drawn and filled.
fig = figure()
lon_0=57
m1 = Basemap(projection='geos',lon_0=lon_0,rsphere=(6378137.00,6356752.3142),resolution=None)
ax = fig.add_axes([0.1,0.1,0.8,0.8],axisbg='k')
# plot just upper right quadrant.
m = Basemap(projection='geos',lon_0=lon_0,rsphere=(6378137.00,6356752.3142),resolution='l',llcrnrx=0.,llcrnry=0.,urcrnrx=m1.urcrnrx/2.,urcrnry=m1.urcrnry/2.)
print m.projparams
m.drawcoastlines()
m.drawmapboundary(fill_color='aqua')
m.fillcontinents(color='coral',lake_color='aqua')
m.drawcountries()
# draw parallels and meridians.
m.drawparallels(arange(-90.,120.,30.))
m.drawmeridians(arange(0.,420.,60.))
m.drawmapboundary()
title('Geostationary Map Centered on Lon=%s' % lon_0)

fig = figure()
m1 = Basemap(projection='ortho',lon_0=lon_0,lat_0=10,resolution=None)
ax = fig.add_axes([0.1,0.1,0.8,0.8],axisbg='k')
# plot just upper right quadrant.
m = Basemap(projection='ortho',lon_0=lon_0,lat_0=10,resolution='l',llcrnrx=0.,llcrnry=0.,urcrnrx=m1.urcrnrx/2.,urcrnry=m1.urcrnry/2.)
print m.projparams
m.drawcoastlines()
m.drawmapboundary(fill_color='aqua')
m.fillcontinents(color='coral',lake_color='aqua')
m.drawcountries()
# draw parallels and meridians.
m.drawparallels(arange(-90.,120.,30.))
m.drawmeridians(arange(0.,420.,60.))
m.drawmapboundary()
title('Orthographic Map Centered on Lon=%s' % lon_0)

show()

@@ -85,48 +85,8 @@
 re%0A%0A
-# map with continents drawn and filled.%0A
 fig 
@@ -105,16 +105,34 @@
 on_0=57%0A
+# global geos map%0A
 m1 = Bas
@@ -289,32 +289,73 @@
 r right quadrant
+ (coordinates determined from global map)
 .%0Am = Basemap(pr
@@ -498,35 +498,16 @@
 nry/2.)%0A
-print m.projparams%0A
 m.drawco
@@ -820,16 +820,45 @@
 igure()%0A
+# global ortho map%0Alat_0=10.%0A
 m1 = Bas
@@ -891,33 +891,36 @@
 n_0=lon_0,lat_0=
-1
+lat_
 0,resolution=Non
@@ -1002,16 +1002,53 @@
 quadrant
+ (corners determined from global map)
 .%0Am = Ba
@@ -1090,17 +1090,20 @@
 0,lat_0=
-1
+lat_
 0,resolu
@@ -1182,27 +1182,8 @@
 2.)%0A
-print m.projparams%0A
 m.dr
@@ -1472,25 +1472,41 @@
 n Lon=%25s
+, Lat=%25s
 ' %25 
+(
 lon_0
+,lat_0)
 )%0A%0Ashow(

add some comments