#!/usr/bin/env python

import matplotlib.pyplot as plt           # pyplot module import
from mpl_toolkits.basemap import Basemap  # basemap import
import numpy as np                        # Numpy import

# Lambert Conformal map of USA lower 48 states
m = Basemap(llcrnrlon=-119, 
            llcrnrlat=22, 
            urcrnrlon=-64,
            urcrnrlat=49, 
            projection='lcc', 
            lat_1=33, 
            lat_2=45,
            lon_0=-95, 
            resolution='h', 
            area_thresh=10000)

# display blue marble image (from NASA) as map background
m.bluemarble()

# draw the coastlines of continental area
m.drawcoastlines()
# draw country boundaries
m.drawcountries(linewidth=2)
# draw states boundaries (America only)
m.drawstates()

plt.savefig('fig_USmapSat.png')
plt.show()