Create map with Primary Productivity

[mriedijk]

Data from 2002-2016 can be downloaded here in hdf or xyz format: http://www.science.oregonstate.edu/ocean.productivity/standard.product.php

This would result in a separate map, but the same timeslider.

[roblabs]

HDF notes

• http://stackoverflow.com/questions/22703498/read-open-a-modis-aqua-hdf-file-and-display-plot-the-output-in-gdal-and-matplot
• https://trac.osgeo.org/gdal/wiki/HDF
• Sample from the Oregon State web site.
• HDF appears gray scale

$ gdalinfo vgpm.2002335.hdf
Driver: HDF4Image/HDF4 Dataset
Files: vgpm.2002335.hdf
Size is 2160, 1080
Coordinate System is `'
Metadata:
  Central_Azimuth=0
  Hole Value=-9999
  Intercept=0
  Isotropic=0
  Latitude Center=0
  Limit=-90, -180, 90, 180
  Longitude Center=0
  Position=0, 0, 1, 1
  Projection Category=IDL
  Projection ID=8
  Projection Name=Cylindrical
  Rotation=0
  Scaling=linear
  Slope=1
  Source=not specified
  Start Time String=12/01/2002 00:00:00
  Start Time=1038729600
  Stop Time String=12/31/2002 23:59:59
  Stop Time=1041407999
  Units=mgC m-2 day-1
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 1080.0)
Upper Right ( 2160.0,    0.0)
Lower Right ( 2160.0, 1080.0)
Center      ( 1080.0,  540.0)
Band 1 Block=2160x462 Type=Float32, ColorInterp=Gray

xyz data example

lon lat value (-9999 = no data)
-179.916672 89.916664 -9999
-179.750000 89.916664 -9999
-10.249995 64.583336 1741.54
-10.083323 64.583336 1583.7
-9.916667 64.583336 1727.24
-9.749995 64.583336 1742.46
-9.583323 64.583336 1772.91
-9.416667 64.583336 1659.97
-9.249995 64.583336 1648.34
-9.083323 64.583336 1707.72