[GRASS-user] r.to.rast3
bks at centrum.cz
bks at centrum.cz
Fri Jul 27 11:31:16 EDT 2007
Thanks for the answer.
The reason I am trying to do that is really simple. Perhaps, can you give me some recommandations...
I would like to get a simple visualization of only one soil horizon depth in few transects. 2d is good enough...
I have his depth raster (depth from surface) and surface elevation.
Thank you in advance.
Stanislav
______________________________________________________________
> Od: soerengebbert at googlemail.com
> Komu: "bks at centrum.cz" <bks at centrum.cz>
> Datum: 27.07.2007 12:28
> Předmět: Re: [GRASS-user] r.to.rast3
>
>Hi,
>i guess the number of open file handler is
>the limiting factor in this case. Looks like the limit is 1024?
>I don't know if r.to.rast3 has any influence on this.
>And i don't know if you can increase the number of
>open file handler's on your system.
>
>I don't think that i am able to provide a patch
>or a workaround in the next weeks or months.
>
>What you can do is to try a lower z resolution, maybe 0.034
>or something to avoid 3100 depths. Because each map is
>opened for each depth.
>
>Soeren
>
>2007/7/26, bks at centrum.cz <bks at centrum.cz>:
>>
>> Hi,
>> I have installed 6.2.2 but it didn't solved the problem "r.to.rast3". I
do
>> not know, what I am doing badly....
>>
>> Here are details:
>>
>> GRASS 6.2.2cvs (lokvrst102):~ > r.to.rast3 input=elev_bas,elev_surf
>> output=rast3_epipedon --overwrite
>>
>> Open raster map elev_surf - one time for each depth (1/3100)
>> .....
>> Open raster map elev_surf - one time for each depth (1017/3100)
>> Open raster map elev_surf - one time for each depth (1018/3100)
>> Open raster map elev_surf - one time for each depth (1019/3100)
>> Open raster map elev_surf - one time for each depth (1020/3100)
>> WARNING: Unable to open
>> /home/bekas/lokvrst102/PERMANENT/cell_misc/elev_surf/f_format
>> WARNING: quantization file [elev_surf] in mapset [PERMANENT] missing
>> WARNING: Can't open header file for [MASK in PERMANENT]
>> WARNING: Unable to open
>> /home/bekas/lokvrst102/PERMANENT/cell_misc/elev_surf/f_format
>> ERROR: Input maps have to be from the same type. CELL, FCELL or DCELL!
>> [Raster MASK present]
>>
>> GRASS 6.2.2cvs (lokvrst102):~ > g.region -p3
>> projection: 0 (x,y)
>> zone: 0
>> north: -1188163.21777146
>> south: -1188428.35974525
>> west: -581335.49616366
>> east: -581170.41708042
>> top: 246.00000000
>> bottom: 215.00000000
>> nsres: 1.00053575
>> nsres3: 1.00053575
>> ewres: 1.00047929
>> ewres3: 1.00047929
>> tbres: 0.01
>> rows: 265
>> rows3: 265
>> cols: 165
>> cols3: 165
>> depths: 3100
>> cells: 43725
>> 3dcells: 135547500
>> [Raster MASK present]
>>
>>
>> Layer: elev_bas Date: Thu Jul 26 19:27:20
>> 2007 |
>> | Mapset: PERMANENT Login of Creator:
>> bekas |
>> | Location:
>> lokvrst102 |
>> | DataBase:
>> /home/bekas |
>> | Title: ( elev_bas
>> ) |
>> | timestamp:
>> none |
>>
>>
|----------------------------------------------------------------------------|
>>
>> |
|
>> | Type of Map: raster Number of Categories:
>> 255 |
>> | Data
>> Type: DCELL |
>> | Rows:
>> 265 |
>> |
>> Columns: 165
|
>> | Total
>> Cells: 43725 |
>> | Projection: x,y (zone
>> 0) |
>> | N: -1188163.21777146 S: -1188428.35974525 Res:
>> 1.00053575 |
>> | E: -581170.41708042 W: -581335.49616366 Res:
1.00047929
>> |
>> | Range of data: min = 216.996098 max = 243.984776
>> |
>>
>> |
|
>> | Data
>> Source: |
>>
>> |
|
>>
>> |
|
>>
>> |
|
>> | Data
>> Description: |
>> | generated by r.mapcalc
>> |
>>
>> |
|
>> |
>> Comments:
|
>> | elev_d20 - epiped.cm
>> |
>> |
>>
>>
+----------------------------------------------------------------------------+
>> | Layer: elev_surf Date: Thu Jul 26 19:28:20
>> 2007 |
>> | Mapset: PERMANENT Login of Creator:
>> bekas |
>> | Location:
>> lokvrst102 |
>> | DataBase:
>> /home/bekas |
>> | Title: ( elev_surf
>> ) |
>> | timestamp:
>> none |
>>
>>
|----------------------------------------------------------------------------|
>>
>> |
|
>> | Type of Map: raster Number of Categories:
>> 255 |
>> | Data
>> Type: DCELL |
>> | Rows:
>> 265 |
>> |
>> Columns: 165
|
>> | Total
>> Cells: 43725 |
>> | Projection: x,y (zone
>> 0) |
>> | N: -1188163.21777146 S: -1188428.35974525 Res:
>> 1.00053575 |
>> | E: -581170.41708042 W: -581335.49616366 Res:
1.00047929
>> |
>> | Range of data: min = 217.848999 max = 244.296890
>> |
>>
>> |
|
>> | Data
>> Source: |
>>
>> |
|
>>
>> |
|
>>
>> |
|
>> | Data
>> Description: |
>> | generated by r.mapcalc
>> |
>>
>> |
|
>> |
>> Comments:
|
>> | elev_bas + epiped.cm
>> |
>>
>> |
|
>>
>>
+----------------------------------------------------------------------------+
>>
>>
>>
>>
>>
>> ______________________________________________________________
>> > Od: soerengebbert at gmx.de
>> > Komu: grassuser at grass.itc.it
>> > CC: bks at centrum.cz
>> > Datum: 23.07.2007 20:22
>> > Předmět: Re: [GRASS-user] r.to.rast3
>> >
>> >Hi,
>> >i'm not able to reproduce this error with the latest cvs version of
>> grass6.
>> >It seems to work with integer, float and double raster maps.
>> >
>> >Can you please provide more information?
>> >
>> >1.) What version of grass are you using
>> >2.) The command line you have used
>> >3.) The projection and size of your mapset
>> >
>> >Best regards
>> >Soeren
>> >
>> >Am Freitag, 20. Juli 2007 18:26 schrieb bks at centrum.cz:
>> >> Hi,
>> >> I have a problem in runing r.to.rast3.
>> >> I cannot go over the following error message: "Input maps have to be
>> from
>> >> the same type. CELL, FCELL, DCELL!"
>> >>
>> >> The input 2D slice maps were made by r.mapcalc.
>> >> To fix the problem I have tried to: r.mapcalc "map1=float(map1)"
and
>> the
>> >> with map2.
>> >> r.info .... from both of them show: data type: FCELL
>> >> It doesn't help.
>> >>
>> >> Wenn I work with integer (CELL), I receive this error message:
"Could
>> not
>> >> open map ... (allway the second one
>> >>
>> >> >from the two input ones)"
>> >>
>> >> Thanks in advance for your help!
>> >> Stanislav
>> >>
>> >> _______________________________________________
>> >> grassuser mailing list
>> >> grassuser at grass.itc.it
>> >> http://grass.itc.it/mailman/listinfo/grassuser
>> >
>>
>
More information about the grass-user
mailing list