[gdal-dev] "Syncing" two rasters together...

[Jonathan Greenberg]

GDALers:

What is the most efficient way, given a "reference raster", and an
arbitrary raster (we'll call it "unsynced") synced together to allow
them to be stacked: the output of this should be the unsynced raster
with the same number of rows, columns, pixel size, upper left
coordinates, and projection as the reference raster.

Here are the assumptions:
- both rasters already have the same pixel size and projection, but
the offset of the upper left coordinate may not be "perfect"

- for each pixel location defined by the reference raster, the output
should have the closest (nearest neighbor) pixel from the unsynced
raster.  It should have some NA value in regions where the reference
and unsynced do not overlap.

Thoughts?  The application of this would be e.g. taking a set of
Landsat images from different time periods (each image has slightly
different numbers of rows and columns), and stacking them together to
perform time series analysis.

Cheers!

--j

-- 
Jonathan A. Greenberg, PhD
Assistant Professor
Global Environmental Analysis and Remote Sensing (GEARS) Laboratory
Department of Geography and Geographic Information Science
University of Illinois at Urbana-Champaign
259 Computing Applications Building, MC-150
605 East Springfield Avenue
Champaign, IL  61820-6371
Phone: 217-300-1924
http://www.geog.illinois.edu/~jgrn/
AIM: jgrn307, MSN: jgrn307 at hotmail.com, Gchat: jgrn307, Skype: jgrn3007