Hi Roger
first of all thanks a lot for your prompt reply. I tried it with two
other raster layers, and one worked.
I played around, and the problem seems top be the length of the output
variable:
GRASS 6.2.2 (grass):~/Documents/Projects/AlienSpread/R > r.out.bin -b
input=Hakea_dc_2007
output=/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee
null=-1
Creating BIL support files ...
Header File = eeeeee.wld <<<<<<------HERE IT IS NOT
World File =
/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee.wld
Exporting Raster as integer values (bytes=2)
Using the Current Region settings ...
north=6247000.000000
south=6195000.000000
east=270200.000000
west=251000.000000
r=520
c=192
Percent complete: 100%
r.out.bin complete.
[Raster MASK present]
GRASS 6.2.2 (grass):~/Documents/Projects/AlienSpread/R > ll /home/rkrug/eee*
-rw-r--r-- 1 rkrug users 199680 2007-10-16 12:18
/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee
-rw-r--r-- 1 rkrug users 60 2007-10-16 12:18
/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee.wld
[Raster MASK present]
GRASS 6.2.2 (grass):~/Documents/Projects/AlienSpread/R > r.out.bin -b
input=Hakea_dc_2007
output=/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee
null=-1
Creating BIL support files ...
Header File =
/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee.hdr
<<<<<<------HERE IT IS
World File =
/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee.wld
Exporting Raster as integer values (bytes=2)
Using the Current Region settings ...
north=6247000.000000
south=6195000.000000
east=270200.000000
west=251000.000000
r=520
c=192
Percent complete: 100%
r.out.bin complete.
[Raster MASK present]
GRASS 6.2.2 (grass):~/Documents/Projects/AlienSpread/R > ll /home/rkrug/eee*
-rw-r--r-- 1 rkrug users 199680 2007-10-16 12:18
/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee
-rw-r--r-- 1 rkrug users 199680 2007-10-16 12:18
/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee
-rw-r--r-- 1 rkrug users 60 2007-10-16 12:18
/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee.wld
-rw-r--r-- 1 rkrug users 83 2007-10-16 12:18
/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee.hdr
-rw-r--r-- 1 rkrug users 60 2007-10-16 12:18
/home/rkrug/eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee.wld
[Raster MASK present]
GRASS 6.2.2 (grass):~/Documents/Projects/AlienSpread/R >
So it is a problem in grass itself.
I should be able to find a way around, but I should file a bug report.
Thanks a million,
Rainer
Roger Bivand wrote:
> On Tue, 16 Oct 2007, Rainer M. Krug wrote:
>
>> Hi
>>
>> I just upgraded to the new version, and I get an error message when
>> importing an existing raster. The output is below.
>
> This is not related to the sp update. Please do debug(readRAST6) and
> re-run, stepping through until you pass the the setting of the command
> string for r.out.bin. Please then print the cmd variable to screen and
debug: tull <- ifelse(.Platform$OS.type == "windows", system(cmd),
system(cmd,
ignore.stderr = ignore.stderr))
Browse[1]> cmd
[1] "r.out.bin -b input=Hakea_seeds_2007
output=/home/rkrug/Documents/Projects/AlienSpread/R/../grass/simulation/.tmp/rkrug-linux/Hakea_seeds_2007
null=-1"
Browse[1]>
> copy&paste it. It should have a -b flag. Please also give the output of
> sessionInfo() for your system, and the output of system("g.version") too.
> sessionInfo()
R version 2.6.0 (2007-10-03)
i686-pc-linux-gnu
locale:
LC_CTYPE=en_GB.UTF-8;LC_NUMERIC=C;LC_TIME=en_GB.UTF-8;LC_COLLATE=en_GB.UTF-8;LC_MONETARY=en_GB.UTF-8;LC_MESSAGES=en_GB.UTF-8;LC_PAPER=en_GB.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_GB.UTF-8;LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] gstat_0.9-40 spgrass6_0.4-3 maptools_0.6-17 rgdal_0.5-15
[5] sp_0.9-16 foreign_0.8-23
loaded via a namespace (and not attached):
[1] grid_2.6.0 lattice_0.16-5
> system("g.version")
GRASS 6.2.2 (2007)
>
>
> The *.wld file got put in the correct place, so the *.hdr ought to have
> been there too.
>
> If you do debug(readBinGrid), you'll see where it fails now, about 4
> lines in. You could also do list.files(dirname(fname)) just before it
> fails in readBinGrid(), to check whether there ia anything else odd
> going on. Have you checked to see whether it has been written anywhere
> else?
>
> Roger
>
>>
>>
>>> readRAST6("Hakea_seeds_2007")
>> Creating BIL support files ...
>> Header File = ld
>> World File =
>> /home/rkrug/Documents/Projects/AlienSpread/R/../grass/simulation/.tmp/rkrug-linux/Hakea_seeds_2007.wld
>>
>> Exporting Raster as double values(bytes=8)
>> Using the Current Region settings ...
>> north=6247000.000000
>> south=6195000.000000
>> east=270200.000000
>> west=251000.000000
>> r=520
>> c=192
>>
>> Percent complete:
>> r.out.bin complete.
>> Error in readBinGrid(rtmpfl11, colname = vname[i], proj4string = p4,
>> integer = to_int) :
>> no such file:
>> /home/rkrug/Documents/Projects/AlienSpread/R/../grass/simulation/.tmp/rkrug-linux/Hakea_seeds_2007.hdr
>>
>>>
>>
>> The contents of the folder are as follow:
>>
>> (grass):~/Documents/Projects/AlienSpread/grass/simulation/.tmp/rkrug-linux
>>
>>> ll Hakea*
>> -rw-r--r-- 1 rkrug users 798720 2007-10-16 11:10 Hakea_seeds_2007
>> -rw-r--r-- 1 rkrug users 60 2007-10-16 11:10 Hakea_seeds_2007.wld
>> -rw-r--r-- 1 rkrug users 798720 2007-10-16 10:52
>> Hakea_seedsDispersed_2007
>> -rw-r--r-- 1 rkrug users 60 2007-10-16 10:52
>> Hakea_seedsDispersed_2007.wld
>> [Raster MASK present]
>> GRASS 6.2.2
>> (grass):~/Documents/Projects/AlienSpread/grass/simulation/.tmp/rkrug-linux
>>
>>>
>>
>> Is there something wrong in my paths or os there a bug?
>>
>> Thanks a lot,
>>
>>
>> Rainer
>>
>>
>> Roger Bivand wrote:
>>> A new version of the foundation sp package is reaching CRAN mirrors,
>>> and
>>> Windows and OSX binaries should follow within a day or so. There is a
>>> significant change in how SpatialGrid objects are created, which we
>>> hope
>>> will not affect users negatively, and which is giving a two orders of
>>> magnitude speedup for reading larger rasters (single band 900 by
>>> 3600 down
>>> from 12 minutes to 6 seconds on an older 1.5GHz machine).
>>>
>>> If anyone sees any negative impacts, please contact the authors
>>> immediately directly.
>>>
>>> Roger
>>>
>>
>