mapcalc combining raster and vector

11 messages
Open this post in threaded view
|

mapcalc combining raster and vector

 HiI 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 [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user
Open this post in threaded view
|

Re: mapcalc combining raster and vector

 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".  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 [hidden email] 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 [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user
Open this post in threaded view
|

Re: mapcalc combining raster and vector

 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 >> [hidden email] >> 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 > [hidden email] > https://lists.osgeo.org/mailman/listinfo/grass-user> _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user
Open this post in threaded view
|

Re: mapcalc combining raster and vector

 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. On Wed., May 8, 2019, 08:17 Moritz Lennert, <[hidden email]> 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 >> [hidden email] >> 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 > [hidden email] > https://lists.osgeo.org/mailman/listinfo/grass-user > _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user
Open this post in threaded view
|

Re: mapcalc combining raster and vector

 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 > > On Wed., May 8, 2019, 08:17 Moritz Lennert, > <[hidden email] > 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 >      >> [hidden email] >      >> 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 >      > [hidden email] >      > https://lists.osgeo.org/mailman/listinfo/grass-user>      > > > _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user
Open this post in threaded view
|

Re: mapcalc combining raster and vector

 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, <[hidden email] [hidden email]> 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      >> [hidden email] [hidden email]      >> 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      > [hidden email] [hidden email]      > https://lists.osgeo.org/mailman/listinfo/grass-user      > _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user ```-- Micha Silver Ben Gurion Univ. Sde Boker, Remote Sensing Lab cell: +972-523-665918``` _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user
Open this post in threaded view
|

Re: mapcalc combining raster and vector

 sure i will send that later.  thksOn Thu., May 9, 2019, 04:00 Micha Silver, <[hidden email]> 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, <[hidden email] [hidden email]> 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      >> [hidden email] [hidden email]      >> 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      > [hidden email] [hidden email]      > https://lists.osgeo.org/mailman/listinfo/grass-user      > _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user ```-- Micha Silver Ben Gurion Univ. Sde Boker, Remote Sensing Lab cell: +972-523-665918``` _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user
Open this post in threaded view
|

Re: mapcalc combining raster and vector

 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.SELECTcat,int_1,int_2,dbl_1,dbl_2,dbl_3,str_1,dbl_4,dbl_5,dbl_6,str_2,dbl_71,11105041,17,580165.3,4873100,279.89767,SAND,277.7976701,279.59767,1.21,WWR,5e-0062,11105041,17,580165.3,4873100,279.89767,SAND,275.5976698,277.7976701,1.21,WWR,5e-0063,11105041,17,580165.3,4873100,279.89767,SAND,269.1976702,275.5976698,1.21,WWR,5e-0064,11105041,17,580165.3,4873100,279.89767,SAND,277.7976701,278,1.21,WWR,5e-0065,11105041,17,580165.3,4873100,279.89767,SAND,277.7976701,279,1.21,WWR,5e-0066,11105041,17,580165.3,4873100,279.89767,SAND,275.5976698,276,1.21,WWR,5e-0067,11105041,17,580165.3,4873100,279.89767,SAND,275.5976698,277,1.21,WWR,5e-006Le jeu. 9 mai 2019 à 04:59, Francois Chartier <[hidden email]> a écrit :sure i will send that later.  thksOn Thu., May 9, 2019, 04:00 Micha Silver, <[hidden email]> 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, <[hidden email] [hidden email]> 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      >> [hidden email] [hidden email]      >> 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      > [hidden email] [hidden email]      > https://lists.osgeo.org/mailman/listinfo/grass-user      > _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user ```-- Micha Silver Ben Gurion Univ. Sde Boker, Remote Sensing Lab cell: +972-523-665918``` _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user
Open this post in threaded view
|

Re: mapcalc combining raster and vector

 OK, please clarify what each column reprepsents. On 10/05/2019 5:25, Francois Chartier wrote: 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 cat, bore_id??, ??, easting?, northing?, depth?, soil_type, ??, ??,  bulk_density??, WWR ?, hydraulic_conduct?? 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 <[hidden email]> a écrit : sure i will send that later.  thks On Thu., May 9, 2019, 04:00 Micha Silver, <[hidden email]> 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, <[hidden email] [hidden email]> 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      >> [hidden email] [hidden email]      >> 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      > [hidden email] [hidden email]      > https://lists.osgeo.org/mailman/listinfo/grass-user      > _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user ```-- Micha Silver Ben Gurion Univ. Sde Boker, Remote Sensing Lab cell: +972-523-665918``` ```-- Micha Silver Ben Gurion Univ. Sde Boker, Remote Sensing Lab cell: +972-523-665918``` _______________________________________________ grass-user mailing list [hidden email] https://lists.osgeo.org/mailman/listinfo/grass-user
Open this post in threaded view
|

Re: mapcalc combining raster and vector

Sorry i realised later that in my current table i didnt rename the columns. here are the extract from excel (300 K pts).
i am using the TopDepthRasterDEM_MASL as the Z elevation, and the attribute i am interested in are Texture in mm and K in m/sec

 BHID Zone EAST83 NORTH83 GTADEM MAT1 BOTDepthRasterDEM_MASL TOPDepthRasterDEM_MASL TEXTURE_mm BH_Type K 11105041 17 580165.3 4873100 279.8977 SAND 277.7976701 279.59767 1.21 WWR 0.0000050 11105041 17 580165.3 4873100 279.8977 SAND 275.5976698 277.7976701 1.21 WWR 0.0000050 11105041 17 580165.3 4873100 279.8977 SAND 269.1976702 275.5976698 1.21 WWR 0.0000050 11105041 17 580165.3 4873100 279.8977 SAND 277.7976701 278 1.21 WWR 0.0000050 11105041 17 580165.3 4873100 279.8977 SAND 277.7976701 279 1.21 WWR 0.0000050 11105041 17 580165.3 4873100 279.8977 SAND 275.5976698 276 1.21 WWR 0.0000050

Le jeu. 9 mai 2019 à 22:45, Micha Silver <[hidden email]> a écrit :

OK, please clarify what each column reprepsents.

On 10/05/2019 5:25, Francois Chartier wrote:

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
cat, bore_id??, ??, easting?, northing?, depth?, soil_type, ??, ??,  bulk_density??, WWR ?, hydraulic_conduct??
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 <[hidden email]> a écrit :
sure i will send that later.  thks

On Thu., May 9, 2019, 04:00 Micha Silver, <[hidden email]> 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, <[hidden email] [hidden email]> 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
>> [hidden email] [hidden email]
>> 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
> [hidden email] [hidden email]
> https://lists.osgeo.org/mailman/listinfo/grass-user
>

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
```--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918```
```--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918```

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Open this post in threaded view
|

Re: mapcalc combining raster and vector

On 11/05/2019 3:16, Francois Chartier wrote:
Sorry i realised later that in my current table i didnt rename the columns. here are the extract from excel (300 K pts).
i am using the TopDepthRasterDEM_MASL as the Z elevation, and the attribute i am interested in are Texture in mm and K in m/sec

 BHID Zone EAST83 NORTH83 GTADEM MAT1 BOTDepthRasterDEM_MASL TOPDepthRasterDEM_MASL TEXTURE_mm BH_Type K 11105041 17 580165.3 4873100 279.8977 SAND 277.7976701 279.59767 1.21 WWR 0.0000050 11105041 17 580165.3 4873100 279.8977 SAND 275.5976698 277.7976701 1.21 WWR 0.0000050 11105041 17 580165.3 4873100 279.8977 SAND 269.1976702 275.5976698 1.21 WWR 0.0000050 11105041 17 580165.3 4873100 279.8977 SAND 277.7976701 278 1.21 WWR 0.0000050 11105041 17 580165.3 4873100 279.8977 SAND 277.7976701 279 1.21 WWR 0.0000050 11105041 17 580165.3 4873100 279.8977 SAND 275.5976698 276 1.21 WWR 0.0000050

First you can start with the suggestion by Moritz, to get the surface elevation from your DTM uisng v.what.rast.

Then, calculate the depth (from the surface) of each sample with a simple v.db.update expression.

Finally an SQL query in db.execute will get a table of aggregated texture and K for depth ranges that you choose. This table will be piped into v.in.ascii to get a new point vector map for each depth range. (You can run the db.execute command without v.in.ascii to check that you're getting the results you expect)

Here are the commands:

# 1.Surface elevation

v.db.addcolumn map=borehole_samples column="depth DOUBLE, surface_elev DOUBLE"

v.what.rast map=borehole_samples raster=dtm column=surface_elev

# 2.Get Depth of samples

v.db.update map=borehole_samples column=depth query_column="surface_elev - TOPDepthRasterDEM_MASL"

# 3. Now an aggregate query for the texture and K

# depth from 0 to 10 m

db.execute -c sql="SELECT BHID as bhid, EAST83 AS x, NORTH83 AS y, AVG(TEXTURE_mm) as avg_texture, AVG(K) as avg_K WHERE depth>=0 and depth<10 GROUP BY bhid,x,y" | v.in.ascii input=- output=soil_characteristics_0 x=2 y=3 columns="bhid INTEGER, x DOUBLE, y DOUBLE, avg_texture DOUBLE, avg_K DOUBLE"

# depth from 10 to 20

db.execute -c sql="SELECT BHID as bhid, EAST83 AS x, NORTH83 AS y, AVG(TEXTURE_mm) as avg_texture, AVG(K) as avg_K WHERE depth>=10 and depth<20 GROUP BY bhid,x,y" | v.in.ascii input=- output=soil_characteristics_10 x=2 y=3 columns="bhid INTEGER, x DOUBLE, y DOUBLE, avg_texture DOUBLE, avg_K DOUBLE"

# and so on for all the depth ranges you need.

Le jeu. 9 mai 2019 à 22:45, Micha Silver <[hidden email]> a écrit :

OK, please clarify what each column reprepsents.

On 10/05/2019 5:25, Francois Chartier wrote:

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
cat, bore_id??, ??, easting?, northing?, depth?, soil_type, ??, ??,  bulk_density??, WWR ?, hydraulic_conduct??
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 <[hidden email]> a écrit :
sure i will send that later.  thks

On Thu., May 9, 2019, 04:00 Micha Silver, <[hidden email]> 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, <[hidden email] [hidden email]> 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
>> [hidden email] [hidden email]
>> 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
> [hidden email] [hidden email]
> https://lists.osgeo.org/mailman/listinfo/grass-user
>

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
```--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918```
```--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918```
```--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918```

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user