[GRASSLIST:10487] Re: Question about v.vol.rst in GRASS 6.1

[Markus Neteler]

Tom,

the settings are not right:

v.info ...
 |            B: 99.000        T: 1155.000                                    |

g.region -p3
top:        1.00000000
bottom:     0.00000000


You see that the 3D region is defined from 0 to 1 meter while
3D points are "flying" from 99 to 1155 meters.

You will have to use
       t   Value for the top edge
       b   Value for the bottom edge
    res3   3D grid resolution (north-south, east-west and top-bottom)
of g.region to get it right.

Markus


On Tue, Feb 21, 2006 at 11:36:51AM -0500, Thomas Adams wrote:
> Markus,
> 
> I did as you suggested ? it looks OK to me, although I'm not sure. 
> Please see the attachment.
> 
> Tom
> 
> 
> Markus Neteler wrote:
> 
> >On Tue, Feb 21, 2006 at 11:13:52AM -0500, Thomas Adams wrote:
> > 
> >
> >>List:
> >>
> >>I am trying to use v.vol.rst in GRASS 6.1cvs on a Linux system. First, 
> >>using the command:
> >>
> >>v.vol.rst -c input=sta_obs_temp_20050731_08 wcolumn=temp 
> >>maskmap=gtopo30.eastern.us segmax=700 tension=10 cvdev=cvdev_map
> >>
> >>I get the following:
> >>
> >>118 records selected from table
> >>WARNING: some points outside of region -- will ignore...
> >>WARNING: strip exists with insufficient data
> >>Processing all selected output files will require 0 bytes of disk space 
> >>for temp files
> >>
> >>WARNING: There are points outside specified 2D/3D region--ignored 118 
> >>points (total points: 118)
> >>ERROR: zero points in the given region!
> >>ERROR: input failed!
> >>
> >>I just don't understand what's going on; all the points appear when I 
> >>display them. I'm able to query the points and perform other analyses, 
> >>such as, v.surf.idw.
> >>
> >>Regards,
> >>   
> >>
> >
> >Tom,
> >
> >please check with 
> >
> >v.info sta_obs_temp_20050731_08
> >
> >top and bottom of the data (since we are using 3D data here).
> >
> >Then check with
> >g.region -p3
> >
> >if the volume region settings are defined appropriately
> >("There are points outside specified 2D/3D region--ignored 118").
> >
> >I assume that you have to define the volume properly.
> >Please suggest a better error message.
> >
> >Regards
> >
> >Markus
> >
> > 
> >
> 
> 
> -- 
> Thomas E Adams
> National Weather Service
> Ohio River Forecast Center
> 1901 South State Route 134
> Wilmington, OH 45177
> 
> EMAIL:	thomas.adams at noaa.gov
> 
> VOICE:	937-383-0528
> FAX:	937-383-0033
> 

> GRASS-GRID > v.info sta_obs_temp_20050731_08
>  +----------------------------------------------------------------------------+
>  | Layer:    sta_obs_temp_20050731_08       Organization:                     |
>  | Mapset:   tea                            Source Date:                      |
>  | Location: mm5                            Name of creator:                  |
>  | Database: /grass/data                                                      |
>  | Title:                                                                     |
>  | Map Scale:  1:1                                                            |
>  | Map format: native                                                         |
>  |----------------------------------------------------------------------------|
>  |   Type of Map:  Vector (level: 2)                                          |
>  |                                                                            |
>  |   Number of points:       118             Number of areas:      0          |
>  |   Number of lines:        0               Number of islands:    0          |
>  |   Number of boundaries:   0               Number of faces:      0          |
>  |   Number of centroids:    0               Number of kernels:    0          |
>  |                                                                            |
>  |   Map is 3D:              1                                                |
>  |   Number of dblinks:      1                                                |
>  |                                                                            |
>  |   Projection: Lambert Conformal Conic (zone 0)                             |
>  |            N: 398185.942    S: -342056.941                                 |
>  |            E: 632334.575    W: -364736.996                                 |
>  |            B: 99.000        T: 1155.000                                    |
>  |                                                                            |
>  |   Digitize threshold: 0.00000                                              |
>  |   Comments:                                                                |
>  |                                                                            |
>  +----------------------------------------------------------------------------+
> 
> Mapset <tea> in Location <mm5>
> GRASS-GRID > g.region -p3
> projection: 99 (Lambert Conformal Conic)
> zone:       0
> datum:      wgs84
> ellipsoid:  wgs84
> north:      486000
> south:      -486000
> west:       -351000
> east:       675000
> top:        1.00000000
> bottom:     0.00000000
> nsres:      9000
> nsres3:     1
> ewres:      9000
> ewres3:     1
> tbres:      1
> rows:       108
> rows3:      972000
> cols:       114
> cols3:      1026000
> depths:     1
> 


-- 
Markus Neteler     <neteler itc it>       http://mpa.itc.it
ITC-irst -  Centro per la Ricerca Scientifica e Tecnologica
MPBA - Predictive Models for Biol. & Environ. Data Analysis
Via Sommarive, 18        -       38050 Povo (Trento), Italy