;Used to generate a color map using the standard lambert projection. Requires a map of the earth, BlueMarble.png, 5400x2700pix.
;Generates 3 images, color_lambert_[0-2].png which must be combined in the GIMP to create the color image.
pro color_lambert,_extra=extra
  Filename="BlueMarble.png"
  cityloc=[[-147.870  ,    64.8200],[-105.200   ,   40.0000],[-76.3300   ,   37.0200],[-111.820     , 41.7600],[0.000000  ,   0.000000]]
  cityname=["Fairbanks, AK","Boulder, CO","Hampton, VA","Logan UT"]
  for channel=0,2 do begin
    image=reform((read_png(Filename))[channel,0:5399,1800:2699])
    lat=float(indgen(900))*60./900.+30.
    lat=rebin(transpose(lat),5400,900)
    lon=float(indgen(5400))*360./5400.-180.
    lon=rebin(lon,5400,900)
    limage=map_patch_lambert(Image, Lon, Lat, bbox=bbox,_extra=extra,missing=255)
    erase
    tv,255-limage
    draw_lambert,bbox=[0,0],_extra=extra,linestyle=0,cityloc=cityloc,cityname=cityname,citycolor=0,linecolor=0
    saveimage,string(channel,format='(%"color_lambert_%d.png")'),/png,/invert
  end
end