[gdal-dev] gdalwarp coordinate alignment with -tap option

Even Rouault even.rouault at mines-paris.org
Thu Feb 13 10:43:08 PST 2014


Le jeudi 13 février 2014 17:03:11, Gregory, Matthew a écrit :
> Hi Eli,
> 
> Eli Adam wrote:
> > <matt.gregory at oregonstate.edu> wrote:
> > > Hi all,
> > > 
> > > We have been using -tap as a command-line argument to gdalwarp to
> > > force pixels to be aligned with the origin (0,0) in the output
> > 
> > coordinate system.
> > 
> > > We now have a situation where we want to specify a different 'snap'
> > > coordinate for pixel alignment (e.g. (15,15)).  Is there a simple
> > > mechanism for specifying this within gdalwarp that I'm overlooking?
> > 
> > Not exactly sure what you mean.  But you might try using -tr as well.
> > I'm not exactly sure of how -tr and -tap relate from the docs,
> > 
> > -tr xres yres:set output file resolution (in target georeferenced units)
> > -tap:(GDAL >= 1.8.0) (target aligned pixels) align the coordinates of
> > the extent of the output file to the values of the -tr, such that the
> > aligned extent includes the minimum extent.
> 
> Thanks for your help.  I left out some key information from my initial post
> which was that I am also passing -tr 30 30 as an argument to gdalwarp. 
> When using -tap along with this, the pixels are aligned with the origin. 
> What I really want is for the 30m pixels to have a corner at (15, 15) in
> the output projection.
> 
> The only way I've figured out how to do this is to first find the projected
> extents of the output image (via osr in python), modify these extents to
> snap to my desired snap window and then use the -te command line option to
> explicitly set the extent.
> 
> I was just looking for an obvious way to specify this just using gdalwarp
> command line options.

Your method is the best method I can think of currently.

> 
> thanks, matt
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
Geospatial professional services
http://even.rouault.free.fr/services.html


More information about the gdal-dev mailing list