[gdal-dev] Re: gdal_rasterize -tr and -te

[Even Rouault]

Le mercredi 06 octobre 2010 18:42:14, Eli Adam a écrit :
> The results of -tap are what I expect, pixels aligned with 0,0.  Although I
> have a separate question at the end.

ok, I've applied the patch in trunk in r20778

> 
> Prior to the patch I ran this:
>  gdal_rasterize -tr 100 100 -a Elevation -l town town.shp townrasterb.tif
> and got this:
> gdalinfo townrasterb.tif
> ...
> Origin = (7267911.066051597706974,512434.991065477312077)
> Pixel Size = (100.000000000000000,-100.000000000000000)
> Metadata:
>   AREA_OR_POINT=Area
> Image Structure Metadata:
>   INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  ( 7267911.066,  512434.991) (124d 6'49.77"W, 45d 0'53.75"N)
> Lower Left  ( 7267911.066,  290634.991) (124d 4'33.15"W, 44d24'25.68"N)
> Upper Right ( 7392811.066,  512434.991) (123d37'52.26"W, 45d 1'45.22"N)
> Lower Right ( 7392811.066,  290634.991) (123d35'53.84"W, 44d25'16.61"N)
> Center      ( 7330361.066,  401534.991) (123d51'17.02"W, 44d43' 6.24"N)
> Band 1 Block=1249x1 Type=Float64, ColorInterp=Gray
> 
> After the patch I ran this:
> gdal_rasterize -tr 100 100 -tap -a Elevation -l town town.shp
> townrasterd.tif and got this:
> gdalinfo townrasterd.tif
> ...
> Origin = (7267900.000000000000000,512500.000000000000000)
> Pixel Size = (100.000000000000000,-100.000000000000000)
> Metadata:
>   AREA_OR_POINT=Area
> Image Structure Metadata:
>   INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  ( 7267900.000,  512500.000) (124d 6'49.97"W, 45d 0'54.39"N)
> Lower Left  ( 7267900.000,  290600.000) (124d 4'33.28"W, 44d24'25.33"N)
> Upper Right ( 7392900.000,  512500.000) (123d37'51.05"W, 45d 1'45.90"N)
> Lower Right ( 7392900.000,  290600.000) (123d35'52.60"W, 44d25'16.29"N)
> Center      ( 7330400.000,  401550.000) (123d51'16.49"W, 44d43' 6.41"N)
> Band 1 Block=1250x1 Type=Float64, ColorInterp=Gray
> 
> town.shp is a point layer with 12 points, 11 of those points have values in
> -a Elevation (1 is null), and in the output tif, 9 pixels have values. 
> The far south and far east points do not have corresponding pixels with
> values.  If this needs review, I can open a ticket and attach a very small
> example.
> 
> I include this info in this thread since that is the case before -tap and
> after -tap I get the same thing, except the far east point has a
> corresponding pixel with -tap.  This probably doesn't have anything (much)
> to do with -tap.
> 

Yes, I'm pretty confident this has nothing to do with -tap. More likely an 
issue with the rasterize code itself. You can probably file a ticket with data 
and command on how to reproduce to keep track of the issue.