[GRASS-user] 3D spatial query

Dear list,

I am new to both this list and GRASS< so please bear with me.

What I would like to achieve is a 3D "spatial query".

I have a 3D points vector layer (imported from a shapefile) which represents archaeological finds data. In addition to this I have volume layers (3D rasters) for individual deposits.
I would like to select and either extract or update the finds data according to deposit, but cannot find a command that does this.

I have managed to successfully run v.what.rast - but then realised that this only works in 2D.

Any ideas would be greatly appreciated!

Many thanks and all best,

Anna

--
Anna Kathrin Hodgkinson MA MSt (Oxon) AIFA EES
Supervisor Geomatics

Oxford Archaeology: Exploring the Human Journey
Mobile: +44 (0)7906 638308 (personal)
Switchboard: +44 (0)1524 541 000
http://thehumanjourney.net

--
Anna Kathrin Hodgkinson MA MSt (Oxon) AIFA EES
Supervisor Geomatics

Oxford Archaeology: Exploring the Human Journey
Mobile: +44 (0)7906 638308 (personal)
Switchboard: +44 (0)1524 541 000
http://thehumanjourney.net

--
Anna Kathrin Hodgkinson MA MSt (Oxon) AIFA EES
Supervisor Geomatics

Oxford Archaeology: Exploring the Human Journey
Mobile: +44 (0)7906 638308 (personal)
Switchboard: +44 (0)1524 541 000
http://thehumanjourney.net

------
Files attached to this email may be in ISO 26300 format (OASIS Open Document Format). If you have difficulty opening them, please visit http://iso26300.info for more information.

On Thu, Oct 7, 2010 at 2:07 PM, Anna Hodgkinson
<anna.hodgkinson@thehumanjourney.net> wrote:

Dear list,

I am new to both this list and GRASS< so please bear with me.

Welcome here!

What I would like to achieve is a 3D "spatial query".

I have a 3D points vector layer (imported from a shapefile) which represents archaeological finds data. In addition to this I have volume layers (3D rasters) for individual deposits.
I would like to select and either extract or update the finds data according to deposit, but cannot find a command that does this.

I have managed to successfully run v.what.rast - but then realised that this only works in 2D.

You would need a new v.what.rast3d module. Perhaps the code of
r3.stats could be helpful for
that in order to see how to obtain voxel info (I thought we had a
r3.what but that's apparently
not the case).

Any ideas would be greatly appreciated!

Some programming seems to be needed for this.

Markus

Many thanks and all best,

Anna

--
Anna Kathrin Hodgkinson MA MSt (Oxon) AIFA EES
Supervisor Geomatics

...

Dear Markus,

Many thanks for your help!

An r3.what was indeed what I was hoping for. Is there any chance of a command such as this being developed any time soon? It would be very useful!

Thank you very much,

Anna Hodgkinson

----- Original Message -----
From: "Markus Neteler" <neteler@osgeo.org>
To: "Anna Hodgkinson" <anna.hodgkinson@thehumanjourney.net>
Cc: "grass-user" <grass-user@lists.osgeo.org>
Sent: Friday, 8 October, 2010 1:01:45 PM
Subject: Re: [GRASS-user] 3D spatial query

On Thu, Oct 7, 2010 at 2:07 PM, Anna Hodgkinson
<anna.hodgkinson@thehumanjourney.net> wrote:

Dear list,

I am new to both this list and GRASS< so please bear with me.

Welcome here!

What I would like to achieve is a 3D "spatial query".

I have a 3D points vector layer (imported from a shapefile) which
represents archaeological finds data. In addition to this I have
volume layers (3D rasters) for individual deposits.
I would like to select and either extract or update the finds data
according to deposit, but cannot find a command that does this.

I have managed to successfully run v.what.rast - but then realised
that this only works in 2D.

You would need a new v.what.rast3d module. Perhaps the code of
r3.stats could be helpful for
that in order to see how to obtain voxel info (I thought we had a
r3.what but that's apparently
not the case).

Any ideas would be greatly appreciated!

Some programming seems to be needed for this.

Markus

Many thanks and all best,

Anna

-- Anna Kathrin Hodgkinson MA MSt (Oxon) AIFA EES
Supervisor Geomatics

...

--
Anna Kathrin Hodgkinson MA MSt (Oxon) AIFA EES
Supervisor Geomatics

Oxford Archaeology: Exploring the Human Journey
Mobile: +44 (0)7906 638308 (personal)
Switchboard: +44 (0)1524 541 000
http://thehumanjourney.net

------
Files attached to this email may be in ISO 26300 format (OASIS Open Document Format). If you have difficulty opening them, please visit http://iso26300.info for more information.

On Sun, Oct 10, 2010 at 7:03 PM, Anna Hodgkinson
<anna.hodgkinson@thehumanjourney.net> wrote:

Dear Markus,

Many thanks for your help!

An r3.what was indeed what I was hoping for. Is there any chance of a command
such as this being developed any time soon? It would be very useful!

A chance is always there :slight_smile:
Please submit an enhancement ticket here:

http://trac.osgeo.org/grass/

Cheers
Markus