[GRASS-user] 3D raster color coding

Francois Chartier fra.chartier at gmail.com
Sat Jun 2 11:17:42 PDT 2018


and here is g.region -p3
projection: 1 (UTM)
zone:       17
datum:      nad83
ellipsoid:  grs80
north:      4870795
south:      4825715
west:       609204
east:       651669
top:        352.00000000
bottom:     45.00000000
nsres:      500.88888889
nsres3:     500.88888889
ewres:      499.58823529
ewres3:     499.58823529
tbres:      1
rows:       90
rows3:      90
cols:       85
cols3:      85
depths:     307
cells:      7650
cells3:     2348550

Le sam. 2 juin 2018 à 14:12, Francois Chartier <fra.chartier at gmail.com> a
écrit :

> see export below of color tab.
> (Sat Jun 02 14:11:39 2018)
>
> r3.colors.out map=BHALL_May24v2 at Toronto
>
> -8.99099 0:0:255
> 5.30354 0:255:255
> 19.5981 255:255:0
> 33.8926 255:0:0
> nv 255:255:255
> default 255:255:255
> (Sat Jun 02 14:11:39 2018) Command finished (0 sec)
>
> Le sam. 2 juin 2018 à 14:05, Francois Chartier <fra.chartier at gmail.com> a
> écrit :
>
>> C:\Users\Francois Chartier>g.region -p
>> projection: 1 (UTM)
>> zone:       17
>> datum:      nad83
>> ellipsoid:  grs80
>> north:      4870795
>> south:      4825715
>> west:       609204
>> east:       651669
>> nsres:      500.88888889
>> ewres:      499.58823529
>> rows:       90
>> cols:       85
>> cells:      7650
>>
>> metadata from 3d raster
>> Title:    BHALL_May24v2
>>   |
>>  | Units:    none
>>      |
>>  | Vertical unit: units
>>      |
>>  | Timestamp: none
>>     |
>>
>>  |----------------------------------------------------------------------------|
>>  |
>>     |
>>  |   Type of Map:  raster_3d            Number of Categories: 0
>>      |
>>  |   Data Type:    FCELL
>>     |
>>  |   Rows:         90
>>      |
>>  |   Columns:      85
>>      |
>>  |   Depths:       307
>>     |
>>  |   Total Cells:  2348550
>>     |
>>  |   Total size:           3180868 Bytes
>>     |
>>  |   Number of tiles:      336
>>     |
>>  |   Mean tile size:       9466 Bytes
>>      |
>>  |   Tile size in memory:  30420 Bytes
>>     |
>>  |   Number of tiles in x, y and  z:   6, 7, 8
>>     |
>>  |   Dimension of a tile in x, y, z:   15, 13, 39
>>      |
>>  |
>>     |
>>  |        Projection: UTM (zone 17)
>>      |
>>  |            N:    4870795    S:    4825715   Res: 500.88888889
>>     |
>>  |            E:     651669    W:     609204   Res: 499.58823529
>>     |
>>  |            T:        352    B:         45   Res:     1
>>      |
>>  |   Range of data:   min = -8.99099445 max = 33.89262009
>>      |
>>  |
>>     |
>>  |   Data Source:
>>      |
>>  |
>>     |
>>  |
>>     |
>>  |
>>     |
>>  |   Data Description:
>>     |
>>  |    generated by v.vol.rst
>>     |
>>  |
>>     |
>>  |   Comments:
>>     |
>>  |    v.vol.rst input="BHMat1All_May12 at Toronto" wcolumn="cat"
>> tension=40. \   |
>>  |    smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=2 wscale=1.0
>> zscale=1.\   |
>>  |    0 elevation="BHALL_May24v2"
>>      |
>>  |
>>     |
>>
>>  +----------------------------------------------------------------------------+
>>
>> Le sam. 2 juin 2018 à 14:01, Veronica Andreo <veroandreo at gmail.com> a
>> écrit :
>>
>>> g.region manual page:
>>> https://grass.osgeo.org/grass74/manuals/g.region.html
>>>
>>> # print current region info
>>> g.region -p
>>>
>>> # set to a raster 3d
>>> g.region raster_3d=nameofyourmap
>>>
>>> El sáb., 2 jun. 2018 19:58, Francois Chartier <fra.chartier at gmail.com>
>>> escribió:
>>>
>>>> how do you ensure that region is selected to the extent of the raster.
>>>> I think I already did this when making v.vol.rst.
>>>>
>>>> Le sam. 2 juin 2018 à 13:56, Veronica Andreo <veroandreo at gmail.com> a
>>>> écrit :
>>>>
>>>>> What about the region? Is it set to the extent of your raster?
>>>>>
>>>>> El sáb., 2 jun. 2018 19:52, Francois Chartier <fra.chartier at gmail.com>
>>>>> escribió:
>>>>>
>>>>>> i removed the 'r' and used elevation and bcyr as the color tables,
>>>>>> see below.  Display window is blank.
>>>>>>
>>>>>> to recap the steps:
>>>>>> 1 - i selected the Mapset  and also tried Map that corresponds to the
>>>>>> 3d raster
>>>>>> 2 - then in Define tab, i only choose the Name of Color Table: bcyr,
>>>>>> i have not selected anything else,
>>>>>> 3 - then i just Run.
>>>>>>
>>>>>> (Sat Jun 02 13:45:02 2018)
>>>>>>
>>>>>> r3.colors map=BHALL_May24v2 at Toronto color=elevation
>>>>>>
>>>>>> Color table for raster3d map <BHALL_May24v2 at Toronto> set to
>>>>>> 'elevation'
>>>>>> (Sat Jun 02 13:45:02 2018) Command finished (0 sec)
>>>>>>
>>>>>> (Sat Jun 02 13:46:14 2018)
>>>>>>
>>>>>> r3.colors map=BHALL_May24v2 at Toronto color=bcyr
>>>>>>
>>>>>> Color table for raster3d map <BHALL_May24v2 at Toronto> set to 'bcyr'
>>>>>> (Sat Jun 02 13:46:15 2018) Command finished (0 sec)
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> Le sam. 2 juin 2018 à 09:01, Francois Chartier <
>>>>>> fra.chartier at gmail.com> a écrit :
>>>>>>
>>>>>>> Is there any way to set the range for the color table?
>>>>>>>
>>>>>>> On Sat, Jun 2, 2018, 02:33 Veronica Andreo, <veroandreo at gmail.com>
>>>>>>> wrote:
>>>>>>>
>>>>>>>> Try removing -r from your command (unchecking the box in the GUI).
>>>>>>>>
>>>>>>>> That flag is used to remove a color palette.
>>>>>>>>
>>>>>>>> Best,
>>>>>>>> Vero
>>>>>>>>
>>>>>>>> El sáb., 2 jun. 2018 3:48, Francois Chartier <
>>>>>>>> fra.chartier at gmail.com> escribió:
>>>>>>>>
>>>>>>>>>
>>>>>>>>> Hi Veronica,
>>>>>>>>>
>>>>>>>>> something is odd.  here are the steps - Display window is still
>>>>>>>>> blank and the 3D raster has data ranging from -8 to 34.
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Step 2 - set as viridis - according to manual this
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> *result in command window*
>>>>>>>>> (Fri Jun 01 21:44:05 2018)
>>>>>>>>>
>>>>>>>>> r3.colors -r map=BHALL_May24v2 at Toronto color=viridis
>>>>>>>>>
>>>>>>>>> (Fri Jun 01 21:44:05 2018) Command finished (0 sec)
>>>>>>>>>
>>>>>>>>> *Display is still blank.   *
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> *Data is there*
>>>>>>>>>
>>>>>>>>>           |
>>>>>>>>>  |   Type of Map:  raster_3d            Number of Categories: 0
>>>>>>>>>            |
>>>>>>>>>  |   Data Type:    FCELL
>>>>>>>>>             |
>>>>>>>>>  |   Rows:         90
>>>>>>>>>            |
>>>>>>>>>  |   Columns:      85
>>>>>>>>>            |
>>>>>>>>>  |   Depths:       307
>>>>>>>>>             |
>>>>>>>>>  |   Total Cells:  2348550
>>>>>>>>>             |
>>>>>>>>>  |   Total size:           3180868 Bytes
>>>>>>>>>             |
>>>>>>>>>  |   Number of tiles:      336
>>>>>>>>>             |
>>>>>>>>>  |   Mean tile size:       9466 Bytes
>>>>>>>>>            |
>>>>>>>>>  |   Tile size in memory:  30420 Bytes
>>>>>>>>>             |
>>>>>>>>>  |   Number of tiles in x, y and  z:   6, 7, 8
>>>>>>>>>             |
>>>>>>>>>  |   Dimension of a tile in x, y, z:   15, 13, 39
>>>>>>>>>            |
>>>>>>>>>  |
>>>>>>>>>             |
>>>>>>>>>  |        Projection: UTM (zone 17)
>>>>>>>>>            |
>>>>>>>>>  |            N:    4870795    S:    4825715   Res: 500.88888889
>>>>>>>>>             |
>>>>>>>>>  |            E:     651669    W:     609204   Res: 499.58823529
>>>>>>>>>             |
>>>>>>>>>  |            T:        352    B:         45   Res:     1
>>>>>>>>>            |
>>>>>>>>>  |   Range of data:   min = -8.99099445 max = 33.89262009
>>>>>>>>>            |
>>>>>>>>>  |
>>>>>>>>>             |
>>>>>>>>>  |   Data Source:
>>>>>>>>>            |
>>>>>>>>>  |
>>>>>>>>>             |
>>>>>>>>>  |
>>>>>>>>>             |
>>>>>>>>>  |
>>>>>>>>>             |
>>>>>>>>>  |   Data Description:
>>>>>>>>>             |
>>>>>>>>>  |    generated by v.vol.rst
>>>>>>>>>             |
>>>>>>>>>  |
>>>>>>>>>             |
>>>>>>>>>  |   Comments:
>>>>>>>>>             |
>>>>>>>>>  |    v.vol.rst input="BHMat1All_May12 at Toronto" wcolumn="cat"
>>>>>>>>> tension=40. \   |
>>>>>>>>>  |    smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=2 wscale=1.0
>>>>>>>>> zscale=1.\   |
>>>>>>>>>  |    0 elevation="BHALL_May24v2"
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Le ven. 1 juin 2018 à 12:54, Francois Chartier <
>>>>>>>>> fra.chartier at gmail.com> a écrit :
>>>>>>>>>
>>>>>>>>>> Thanks, I'll check it out.
>>>>>>>>>>
>>>>>>>>>> On Fri, Jun 1, 2018, 03:05 Veronica Andreo, <veroandreo at gmail.com>
>>>>>>>>>> wrote:
>>>>>>>>>>
>>>>>>>>>>> Hello,
>>>>>>>>>>>
>>>>>>>>>>> That's probably because that palette is for integer values
>>>>>>>>>>> ranging from 0 to 360 [0]. Try with a relative one such as viridis or bcyr.
>>>>>>>>>>>
>>>>>>>>>>> HTH,
>>>>>>>>>>> Vero
>>>>>>>>>>>
>>>>>>>>>>> [0] https://grass.osgeo.org/grass74/manuals/r.colors.html
>>>>>>>>>>>
>>>>>>>>>>> El jue., 31 may. 2018 a las 3:18, Francois Chartier (<
>>>>>>>>>>> fra.chartier at gmail.com>) escribió:
>>>>>>>>>>>
>>>>>>>>>>>> Hi,
>>>>>>>>>>>>
>>>>>>>>>>>> I have created a 3D raster from soil properties with about 70K
>>>>>>>>>>>> points, with a* lot of help from the mailing list i must say*
>>>>>>>>>>>> .
>>>>>>>>>>>> from the copy of the metadata of the 3d raster (below) it would
>>>>>>>>>>>> confirm that the 3D raster worked and the grid is populated with data.
>>>>>>>>>>>> I have copied below the r3.colors and assigned  'aspectcolr'
>>>>>>>>>>>> as a trial.  Unfortunately in the Display window nothing shows up, whether
>>>>>>>>>>>> in 2D or in 3D view.  I do not know what i am doing incorrectly.
>>>>>>>>>>>>
>>>>>>>>>>>> Once I am done and able to develop the 3D raster of the soil
>>>>>>>>>>>> properties, I will prepare a simple step by step for developing such a
>>>>>>>>>>>> raster which may be useful for others.
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> +----------------------------------------------------------------------------+
>>>>>>>>>>>>  | Layer:    BHALL_May24v2 at Toronto          Date: Sun May 27
>>>>>>>>>>>> 16:43:46 2018    |
>>>>>>>>>>>>  | Mapset:   Toronto                        Login of Creator:
>>>>>>>>>>>> Francois Charti |
>>>>>>>>>>>>  | Location: Toronto
>>>>>>>>>>>>               |
>>>>>>>>>>>>  | DataBase: C:\Users\Francois Chartier\Documents\grassdata
>>>>>>>>>>>>                |
>>>>>>>>>>>>  | Title:    BHALL_May24v2
>>>>>>>>>>>>               |
>>>>>>>>>>>>  | Units:    none
>>>>>>>>>>>>                |
>>>>>>>>>>>>  | Vertical unit: units
>>>>>>>>>>>>                |
>>>>>>>>>>>>  | Timestamp: none
>>>>>>>>>>>>               |
>>>>>>>>>>>>
>>>>>>>>>>>>  |----------------------------------------------------------------------------|
>>>>>>>>>>>>  |
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Type of Map:  raster_3d            Number of Categories:
>>>>>>>>>>>> 0               |
>>>>>>>>>>>>  |   Data Type:    FCELL
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Rows:         90
>>>>>>>>>>>>                |
>>>>>>>>>>>>  |   Columns:      85
>>>>>>>>>>>>                |
>>>>>>>>>>>>  |   Depths:       307
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Total Cells:  2348550
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Total size:           3180868 Bytes
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Number of tiles:      336
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Mean tile size:       9466 Bytes
>>>>>>>>>>>>                |
>>>>>>>>>>>>  |   Tile size in memory:  30420 Bytes
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Number of tiles in x, y and  z:   6, 7, 8
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Dimension of a tile in x, y, z:   15, 13, 39
>>>>>>>>>>>>                |
>>>>>>>>>>>>  |
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |        Projection: UTM (zone 17)
>>>>>>>>>>>>                |
>>>>>>>>>>>>  |            N:    4870795    S:    4825715   Res:
>>>>>>>>>>>> 500.88888889              |
>>>>>>>>>>>>  |            E:     651669    W:     609204   Res:
>>>>>>>>>>>> 499.58823529              |
>>>>>>>>>>>>  |            T:        352    B:         45   Res:     1
>>>>>>>>>>>>                |
>>>>>>>>>>>>  |   Range of data:   min = -8.99099445 max = 33.89262009
>>>>>>>>>>>>                |
>>>>>>>>>>>>  |
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Data Source:
>>>>>>>>>>>>                |
>>>>>>>>>>>>  |
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Data Description:
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |    generated by v.vol.rst
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |   Comments:
>>>>>>>>>>>>               |
>>>>>>>>>>>>  |    v.vol.rst input="BHMat1All_May12 at Toronto" wcolumn="cat"
>>>>>>>>>>>> tension=40. \   |
>>>>>>>>>>>>  |    smooth=0.1 segmax=50 npmin=200 npmax=700 dmin=2
>>>>>>>>>>>> wscale=1.0 zscale=1.\   |
>>>>>>>>>>>>  |    0 elevation="BHALL_May24v2"
>>>>>>>>>>>>                |
>>>>>>>>>>>>  |
>>>>>>>>>>>>               |
>>>>>>>>>>>>
>>>>>>>>>>>>  +----------------------------------------------------------------------------+
>>>>>>>>>>>> (Wed May 30 21:11:59 2018) Command finished (0 sec)
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> ----------------------------------------------------------------------------------------
>>>>>>>>>>>>
>>>>>>>>>>>> (Wed May 30 21:10:02 2018)
>>>>>>>>>>>>
>>>>>>>>>>>> r3.colors map=BHALL_May24v2 at Toronto color=aspectcolr
>>>>>>>>>>>>
>>>>>>>>>>>> Color table for raster3d map <BHALL_May24v2 at Toronto> set to
>>>>>>>>>>>> 'aspectcolr'
>>>>>>>>>>>> (Wed May 30 21:10:03 2018) Command finished (0 sec)
>>>>>>>>>>>>
>>>>>>>>>>>> _______________________________________________
>>>>>>>>>>>> grass-user mailing list
>>>>>>>>>>>> grass-user at lists.osgeo.org
>>>>>>>>>>>> https://lists.osgeo.org/mailman/listinfo/grass-user
>>>>>>>>>>>
>>>>>>>>>>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20180602/768aab09/attachment-0001.html>


More information about the grass-user mailing list