Contour Filled SurfaceΒΆ

../../_images/earth_proj1.png
from matplotlib import pyplot as plt
from matplotlib import cm
import s3dlib.surface as s3d

#.. Earth Elevation Contour Projection

# 2. Setup and map surfaces .........................................
rez, lrez  = 5, 5

earth = s3d.SphericalSurface(rez,color='k')
earth.map_geom_from_image('data/elevation.png',1)

elev = earth.contourLines(1.1,coor=2)

surface = elev.get_filled_surface(coor='s',dist=0.05,lrez=lrez)
surface.map_cmap_from_op( lambda c : s3d.SphericalSurface.coor_convert(c,False)[0] , 'jet')

# 3. Construct figure, add surfaces, and plot ......................

fig = plt.figure(figsize=plt.figaspect(1), facecolor='black' )
fig.text(0.02,0.02,str(surface), ha='left', va='bottom',
        fontsize='smaller', color='white')
ax = plt.axes(projection='3d')
ax.set_title('Elevation Contour Projection',color='w')
minmax = (-0.85,0.85)
ax.set(xlim=minmax, ylim=minmax, zlim=minmax)
ax.set_facecolor('black')
ax.set_axis_off()
ax.view_init(25,-150)

ax.add_collection3d(surface)

fig.tight_layout(pad=0)
plt.show()