[R-sig-Geo] doubt with readGDAL

[Barry Rowlingson]

2009/10/4 Roger Bivand <Roger.Bivand at nhh.no>:

> Thanks for a clear report - the y range and the cell sizes were wrong when
> all three of offset=, region.dim= and output.dim= were used together. With
> current rgdal 0.6-18, using the Arc ASCII file (not GTiff), I see for the
> undecimated subscene:
>
>> gridparameters(catb)
>
>  cellcentre.offset     cellsize cells.dim
> x      4.839851e-08 0.0008333333      5000
> y      4.041750e+01 0.0008333333      2600
>>
>> bbox(catb)
>
>            min      max
> x -0.0004166183  4.16625
> y 40.4170835754 42.58375
>
> but for the decimated subscene:
>
>> gridparameters(cata)
>
>  cellcentre.offset   cellsize cells.dim
> x       0.004584215 0.01000167       500
> y     -15.769325078 0.01923397       260
>>
>> bbox(cata)
>
>            min        max
> x -4.166183e-04   5.000417
> y -1.577894e+01 -10.778109
>
>
> In forthcoming 0.6-19 (committed to the sourceforge rgdal repository) and
> the undecimated subscene:

 I'm using: Rgdal=0.6-18 with R=2.9.2, sp=0.9-37and I don't see this
problem (see my previous post).

Barry