Jonathan,

You should set an header in the first file. Dunne_soil.dat is recognized as
an Ascii Grid by GDAL because it has an header :
ncols         720
nrows         360
xllcorner     -180
yllcorner     -90
cellsize      0.5
NODATA_value  -9999

soilwhc.dat has no header. Also NA is coded with no space (-99.0-99.0).
Maybe if you search and replace -99.0 by -9999[space] this would make it
more standard. The lines are also rather short (36), but very long (7200).
I guess each row is indeed 10 lines. It would share the same header than
dunne_soil.dat. So you have some reformatting to get it working.

Etienne


2013/5/10 Jonathan Greenberg <jgrn@illinois.edu>

> r-sig-geo'ers:
>
> I'm trying some random hacks, but was wondering if there is a more
> efficient way to import the following file:
> ftp://ftp.daac.ornl.gov/data/global_soil/Dunne/data/soilwhc.dat
>
> as a raster().  The file is a fixed width file (width=5) of floating point
> values, but had the Arc header stripped off.  A reference file that DOES
> work with raster is:
>
> ftp://ftp.daac.ornl.gov/data/global_soil/Dunne/data/dunne_soil.dat
> reference_raster <- raster(dunne_soil.dat)
> # This has the same number of rows and columns and resolution as
> # the file above.
>
> Thoughts on the most efficient way to convert the first file to a workable
> raster?
>
> --j
>
> --
> Jonathan A. Greenberg, PhD
> Assistant Professor
> Global Environmental Analysis and Remote Sensing (GEARS) Laboratory
> Department of Geography and Geographic Information Science
> University of Illinois at Urbana-Champaign
> 607 South Mathews Avenue, MC 150
> Urbana, IL 61801
> Phone: 217-300-1924
> http://www.geog.illinois.edu/~jgrn/
> AIM: jgrn307, MSN: jgrn307@hotmail.com, Gchat: jgrn307, Skype: jgrn3007
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

	[[alternative HTML version deleted]]

