[R-sig-Geo] ESRI Ascii Grid

[Dylan Beaudette]

On Wed, Jan 14, 2009 at 1:14 AM, Roger Bivand <Roger.Bivand at nhh.no> wrote:
> On Tue, 13 Jan 2009, Dylan Beaudette wrote:
>
>> On Tuesday 13 January 2009, Edzer Pebesma wrote:
>>>
>>> This sounds very familiar; I ran into this about 6 years ago or so.
>>
>> Could it be related to this ticket:
>>
>> http://trac.osgeo.org/gdal/ticket/2369
>
> Yes, exactly! I have checked with readGDAL() with GDAL 1.6.0, and the
> problem is resolved (in fact it defaults to int32), while the problem
> remains in the 1.5.* GDAL series (including 1.5.4 released last week) for
> the sample file in the ticket. May I add the sample file to rgdal to give an
> indication of this?

Yes, please do include the sample file. This particular bug bit me
while processing a bunch of climate data about 4 months ago.

Cheers,

Dylan

>
> The next Windows binary packages will be GDAL 1.5.4, not 1.6.0 - I'm working
> on getting 1.6.0 to check through rgdal.
>
> If Brian is working on something other than Windows with the CRAN binary
> package, upgrading GDAL will fix the problem, because rgdal will use the
> installed GDAL.
>
> Roger
>
> PS. GDAL < 1.6 can be tricked into doing the right thing by using a
> NODATA_value in the file that is either less than the minimum int16 or
> greater than the maximum int16, so with a text editor this can be a
> "solution". The GDAL driver is more flexible in letting you choose
> subscenes, etc, so this may be prefered to readAsciiGrid().
>
>>
>> ?

[snip]