[GRASS-user] mapcalc combining raster and vector

Francois Chartier fra.chartier at gmail.com
Thu May 9 19:25:57 PDT 2019


Location:        Toronto                                                   |
 | Database:        C:\Users\hello\Documents\grassdata
  |
 | Title:
   |
 | Map scale:       1:1
   |
 | Name of creator: hello
   |
 | Organization:
  |
 | Source date:     Thu May 09 19:17:57 2019
  |
 | Timestamp (first layer): none
  |
 |----------------------------------------------------------------------------|
 | Map format:      native
  |
 |----------------------------------------------------------------------------|
 |   Type of map: vector (level: 2)
   |
 |
  |
 |   Number of points:       316866          Number of centroids:  0
  |
 |   Number of lines:        0               Number of boundaries: 0
  |
 |   Number of areas:        0               Number of islands:    0
  |
 |   Number of faces:        0               Number of kernels:    0
  |
 |   Number of volumes:      0               Number of holes:      0
  |
 |
  |
 |   Map is 3D:              Yes
  |
 |   Number of dblinks:      1
  |
 |
  |
 |   Projection: UTM (zone 17)
  |
 |
  |
 |               N:           4891804    S:           4789063
   |
 |               E:            701767    W:            578337
   |
 |               B:         -272.2963    T:         472.80127
   |
 |
  |
 |   Digitization threshold: 0
  |
 |   Comment:
   |
 |
  |
 +----------------------------------------------------------------------------+
(Thu May 09 19:22:41 2019) Command finished (0 sec)


V.DB.SELECT

cat,int_1,int_2,dbl_1,dbl_2,dbl_3,str_1,dbl_4,dbl_5,dbl_6,str_2,dbl_7
1,11105041,17,580165.3,4873100,279.89767,SAND,277.7976701,279.59767,1.21,WWR,5e-006
2,11105041,17,580165.3,4873100,279.89767,SAND,275.5976698,277.7976701,1.21,WWR,5e-006
3,11105041,17,580165.3,4873100,279.89767,SAND,269.1976702,275.5976698,1.21,WWR,5e-006
4,11105041,17,580165.3,4873100,279.89767,SAND,277.7976701,278,1.21,WWR,5e-006
5,11105041,17,580165.3,4873100,279.89767,SAND,277.7976701,279,1.21,WWR,5e-006
6,11105041,17,580165.3,4873100,279.89767,SAND,275.5976698,276,1.21,WWR,5e-006
7,11105041,17,580165.3,4873100,279.89767,SAND,275.5976698,277,1.21,WWR,5e-006

Le jeu. 9 mai 2019 à 04:59, Francois Chartier <fra.chartier at gmail.com> a
écrit :

> sure i will send that later.  thks
>
> On Thu., May 9, 2019, 04:00 Micha Silver, <tsvibar at gmail.com> wrote:
>
>>
>> On 08/05/2019 16:10, Moritz Lennert wrote:
>>
>> On 8/05/19 15:03, Francois Chartier wrote:
>>
>> Hi Micha,
>>
>> I have soil samples from boreholes where i have data points at every
>> meter interval and where soil type changes.  i would like to create a
>> weighted average of soil property for lets say the first 1 m below the DEM,
>> so that i can evaluate potential for infiltration across a watershed.
>>
>>
>> You can use v.what.rast to read the DEM value into the attribute table of
>> your points. Then you can do whatever calculations you need based on that.
>>
>> Moritz
>>
>>
>> After Moritz's answer, It seems that your question is a data aggregation
>> problem. If we can see the actual structure of your data (i.e. v.info -c
>> or v.db.select then more specific responses might be possible. Best
>> would be a small subset of the actual data.
>>
>>
>>
>>
>> On Wed., May 8, 2019, 08:17 Moritz Lennert, <mlennert at club.worldonline.be
>> <mailto:mlennert at club.worldonline.be> <mlennert at club.worldonline.be>>
>> wrote:
>>
>>     On 8/05/19 14:00, Micha Silver wrote:
>>      > You will most likely have to convert the points to raster. But
>>     before
>>      > that, can you explain what you mean by "located for 10 m from
>>     DEM", "10
>>      > to 20".
>>
>>     Or, if you mean zones of equal height, you can transform them to
>> vector
>>     polygons and use v.vect.stats on these polygons and the points.
>>
>>     Moritz
>>
>>      >
>>      >
>>      > On 08/05/2019 13:31, Francois Chartier wrote:
>>      >> Hi
>>      >>
>>      >> I would like to use mapcalc with types of input a DEM (raster)
>>     and xyz
>>      >> property (vector).  basically, I would like to do a weighted
>>     average
>>      >> of data points located for 10 m from DEM, then 10 to 20 etc.
>>      >>
>>      >> _______________________________________________
>>      >> grass-user mailing list
>>      >> grass-user at lists.osgeo.org <mailto:grass-user at lists.osgeo.org>
>> <grass-user at lists.osgeo.org>
>>      >> https://lists.osgeo.org/mailman/listinfo/grass-user
>>      >
>>      > --
>>      > Micha Silver
>>      > Ben Gurion Univ.
>>      > Sde Boker, Remote Sensing Lab
>>      > +972-523-665918
>>      >
>>      >
>>      > _______________________________________________
>>      > grass-user mailing list
>>      > grass-user at lists.osgeo.org <mailto:grass-user at lists.osgeo.org>
>> <grass-user at lists.osgeo.org>
>>      > https://lists.osgeo.org/mailman/listinfo/grass-user
>>      >
>>
>>
>>
>>
>> _______________________________________________
>> grass-user mailing list
>> grass-user at lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>>
>> --
>> Micha Silver
>> Ben Gurion Univ.
>> Sde Boker, Remote Sensing Lab
>> cell: +972-523-665918
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190509/f18d14f9/attachment.html>


More information about the grass-user mailing list