Re: [R-sig-Geo] A new ASTER Global DEM data set

2009-11-20 Thread Barry Rowlingson
2009/11/20 Yong Li yong...@unimelb.edu.au: Dear all folks, I was informed in the 6th Digital Earth Conference that there is a better place to acquire high resolution of global DEM developed by ASTER, called GDEM with 30 m resolution, and fantastically free of charge. I tried some here

[R-sig-Geo] Error, exporting ascii grid generated by DTM grid

2009-11-20 Thread Tobin Cara
Hello, I am receiving a bizarre error recently: Erreur dans write.asciigrid(data_int_ked[1], name1) : Asciigrid does not support grids with non-square cells I do not understand this as I am putting interpolated data back into the cells of an ascii DTM grid with square cells. Information on

Re: [R-sig-Geo] A new ASTER Global DEM data set

2009-11-20 Thread Ashton Shortridge
On Friday 20 November 2009 07:09:35 Barry Rowlingson wrote: 2009/11/20 Yong Li yong...@unimelb.edu.au: Dear all folks, I was informed in the 6th Digital Earth Conference that there is a better place to acquire high resolution of global DEM developed by ASTER, called GDEM with 30 m

Re: [R-sig-Geo] A new ASTER Global DEM data set

2009-11-20 Thread Tomislav Hengl
Hi all, FYI: I've run a small comparison between ASTER DEM and LIDAR DEMs (say 'true' topography): http://geomorphometry.org/content/gdem-quick-assessment The 4 case studies can be downloaded from here: http://geomorphometry.org/system/files/GDEM_assessment.zip I'm not too happy with what

Re: [R-sig-Geo] A new ASTER Global DEM data set

2009-11-20 Thread Karl Ove Hufthammer
On Fri, 20 Nov 2009 16:53:15 +0100 Tomislav Hengl he...@spatial- analyst.net wrote: The 4 case studies can be downloaded from here: http://geomorphometry.org/system/files/GDEM_assessment.zip GDEM is a frankestein - once you zoom in and add some shading, The maps with shading look quite

Re: [R-sig-Geo] Error, exporting ascii grid generated by DTM grid

2009-11-20 Thread Robert J. Hijmans
Shouldn't this: write.asciigrid(data_int_ked[1], name1) be: write.asciigrid(data_int_ked, name1) I think you are writing a single cell, which in this case could be interpreted as rectangular (if it covers the full grid) On Fri, Nov 20, 2009 at 6:06 AM, Tobin Cara cara.to...@epfl.ch wrote:

Re: [R-sig-Geo] changing the data type of a gdal dataset

2009-11-20 Thread Robert J. Hijmans
Alan, You can have a look at the code in the raster package on R-forge (see the writeGDAL function). Or perhaps use this package: one of the primary reasons for developing it was to automate (and hide) block by block raster processing --- and quite successfully, it is being used to process raster

[R-sig-Geo] Writing an ASCII grid with nodata values

2009-11-20 Thread Ned Horning
I am trying to export a SpatialPixelsDataFrame object with floating point values to an ASCII grid file and want to specify a nodata value. I tried write.asciigrid from the sp package and I couldn't figure out how to output float values. The pixel values exported as integers. For now I am