I have DEMs and 1m Aerial imagery. I plan on using gdal_calc to shade the imagery. gdal_calc requires each image to be the same size and no projection checking is performed. so I'd like to retile each dataset using gdal_retile but the output tiles must have the same boundaries.
I was hoping to tile the data in their native resolutions and end up with tiles that have the same boundaries. Any ideas?