I'm working with NOAA's PRECipitation REConstruction over Land dataset, and am looking for an efficient way to sample the data by pixels on NASA's MODIS grid, using open source tools. The PRECL data uses 0.5 degree blocks, so I'll have quite a few repeated values.
I'm currently using a mapping generated by ArcGIS that specifies the array index within the PRECL data for every MODIS pixel. This works (and it's fast, as I'm just doing a series of array lookups), but it'd be nice to get away from that extra data set, especially if there were some way to do the transformation in code, without Arc.
Thanks in advance, all!