ESRI Ascii Grid
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]