[GRASS-user] Vector to 2D Raster

Francois_Chartier · May 8, 2018, 12:26am

Hi,

I am trying to build a Raster out of a vector file (points) which consists of soil properties at a given elevation and coordinates XY.

I have attached an example file that I have extracted from the main table which contains all XY Z and properties. I am proceeding with building a raster for each elevation as I think that grass has issues with using 3D vectors with stacked data (vertically) from Borehole logs (that’s a different issue). Then I will combine the 2D raster into one 3D raster.

I am reasonably new to Grass, and so far I have been able to combine multiple DEM together, add attributes to table from different Rasters (DEM, Bedrock surface), but I am unable to create a raster from a simple table, or at least see the raster. In Grass it shows that the raster is in the Tree, but I am unable to show the colors etc. I am not sure what I am doing incorrectly. I have tried several methods such as v.in.ascii, etc.

As I have been trying this for many days, I am sending this email to the mailing list.

My ultimate goal is to build a 3D raster of soil properties over a large area (40 km x 40 km x 200 m depth).

Thanks for your help!

F.

Sent from Mail for Windows 10

	Virus-free. www.avast.com

(attachments)

Elev145Mat1_May7.xlsx (17.7 KB)

AB1 · May 8, 2018, 3:47pm

Hi Francois,

I think you can use directly v.vol.rst directly on a 3D vector (v.in.ascii with -z flag) to interpolate the points and obtain the 3D raster.

Cheers,

···

Andrea Balotti

Il 08/05/2018 02:26, Francois Chartier ha scritto:

Hi,

I am trying to build a Raster out of a vector file (points) which consists of soil properties at a given elevation and coordinates XY.

I have attached an example file that I have extracted from the main table which contains all XY Z and properties. I am proceeding with building a raster for each elevation as I think that grass has issues with using 3D vectors with stacked data (vertically) from Borehole logs (that’s a different issue). Then I will combine the 2D raster into one 3D raster.

I am reasonably new to Grass, and so far I have been able to combine multiple DEM together, add attributes to table from different Rasters (DEM, Bedrock surface), but I am unable to create a raster from a simple table, or at least see the raster. In Grass it shows that the raster is in the Tree, but I am unable to show the colors etc. I am not sure what I am doing incorrectly. I have tried several methods such as v.in.ascii, etc.

As I have been trying this for many days, I am sending this email to the mailing list.

My ultimate goal is to build a 3D raster of soil properties over a large area (40 km x 40 km x 200 m depth).

Thanks for your help!

F.

Sent from Mail for Windows 10

Virus-free. www.avast.com
_______________________________________________
grass-user mailing list
[grass-user@lists.osgeo.org](mailto:grass-user@lists.osgeo.org)
[https://lists.osgeo.org/mailman/listinfo/grass-user](https://lists.osgeo.org/mailman/listinfo/grass-user)

neteler · May 8, 2018, 9:42pm

On Tue, May 8, 2018 at 5:47 PM, <balotti.and@gmail.com> wrote:

Hi Francois,

I think you can use directly v.vol.rst directly on a 3D vector (v.in.ascii with -z flag) to interpolate the points and obtain the 3D raster.

Good hint - I have added above to this section:

https://grass.osgeo.org/grass74/manuals/raster3dintro.html#3d-raster-import

Best,
Markus

Francois_Chartier · May 8, 2018, 11:28pm

Thanks – I will try and let you know the results.

Sent from Mail for Windows 10

From: Markus Neteler
Sent: May 8, 2018 5:42 PM
To: Andrea Balotti
Cc: GRASS user list
Subject: Re: [GRASS-user] Vector to 2D Raster

On Tue, May 8, 2018 at 5:47 PM, balotti.and@gmail.com wrote:

Hi Francois,

I think you can use directly v.vol.rst directly on a 3D vector (v.in.ascii with -z flag) to interpolate the points and obtain the 3D raster.

Good hint - I have added above to this section:

https://grass.osgeo.org/grass74/manuals/raster3dintro.html#3d-raster-import

Best,

Markus

grass-user mailing list

grass-user@lists.osgeo.org

https://lists.osgeo.org/mailman/listinfo/grass-user

	Virus-free. www.avast.com

Francois_Chartier · May 12, 2018, 6:50pm

Here is the error message that i get:

v.vol.rst input=BHMat1All@Toronto dmin=0.5 elevation=may11

ERROR: Unable to get layer info for vector map

(Sat May 12 14:46:07 2018) Command finished (0 sec)

I am not sure what ‘Layer’ info is this referring to. I have 4 columns, X, Y, Z, Property (numerical). 70,000 points all together.

Best

Francois

Sent from Mail for Windows 10

From: Francois Chartier
Sent: May 8, 2018 7:28 PM
To: Markus Neteler; Andrea Balotti
Cc: GRASS user list
Subject: RE: [GRASS-user] Vector to 2D Raster

Thanks – I will try and let you know the results.

Sent from Mail for Windows 10

From: Markus Neteler
Sent: May 8, 2018 5:42 PM
To: Andrea Balotti
Cc: GRASS user list
Subject: Re: [GRASS-user] Vector to 2D Raster

On Tue, May 8, 2018 at 5:47 PM, balotti.and@gmail.com wrote:

Hi Francois,

I think you can use directly v.vol.rst directly on a 3D vector (v.in.ascii with -z flag) to interpolate the points and obtain the 3D raster.

Good hint - I have added above to this section:

https://grass.osgeo.org/grass74/manuals/raster3dintro.html#3d-raster-import

Best,

Markus

grass-user mailing list

grass-user@lists.osgeo.org

https://lists.osgeo.org/mailman/listinfo/grass-user

	Virus-free. www.avast.com

hellik · May 12, 2018, 7:26pm

v.vol.rst input=BHMat1All@Toronto dmin=0.5 elevation=may11
ERROR: Unable to get layer info for vector map
(Sat May 12 14:46:07 2018) Command finished (0 sec)

could you post here the result of

v.info map=BHMat1All@Toronto

and

v.info -c map=BHMat1All@Toronto

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

Francois_Chartier · May 12, 2018, 7:33pm

v.info message

(Sat May 12 15:29:10 2018)

v.info map=BHMat1All@Toronto

±---------------------------------------------------------------------------+

| Name: BHMat1All |

| Mapset: Toronto |

| Location: Toronto |

| Database: C:\Users\ \Documents\grassdata |

| Title: |

| Map scale: 1:1 |

| Name of creator: |

| Organization: |

| Source date: Sun May 06 14:53:16 2018 |

| Timestamp (first layer): none |

|----------------------------------------------------------------------------|

| Map format: native |

|----------------------------------------------------------------------------|

| Type of map: vector (level: 2) |

| |

| Number of points: 76887 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: 0 |

| |

| Projection: UTM (zone 17) |

| |

| N: 4870795 S: 4825715 |

| E: 651669 W: 609204 |

| B: 45 T: 352 |

| |

| Digitization threshold: 0 |

| Comment: |

| |

±---------------------------------------------------------------------------+

(Sat May 12 15:29:10 2018) Command finished (0 sec)

And

v.info -c map=BHMat1All@Toronto

ERROR: Database connection for map BHMat1All@Toronto is not defined in DB

file

C:\Users\ >

It would seem that it needs to be connected to DB?

Regards,

Sent from Mail for Windows 10

From: Helmut Kudrnovsky
Sent: May 12, 2018 3:26 PM
To: grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] Vector to 2D Raster

v.vol.rst input=BHMat1All@Toronto dmin=0.5 elevation=may11

ERROR: Unable to get layer info for vector map

(Sat May 12 14:46:07 2018) Command finished (0 sec)

could you post here the result of

v.info map=BHMat1All@Toronto

and

v.info -c map=BHMat1All@Toronto

best regards

Helmut

–

Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

grass-user mailing list

grass-user@lists.osgeo.org

https://lists.osgeo.org/mailman/listinfo/grass-user

	Virus-free. www.avast.com

hellik · May 12, 2018, 7:42pm

v.info -c map=BHMat1All@Toronto

ERROR: Database connection for map <BHMat1All@Toronto> is not defined in DB

file

just before you mentioned:

"I have 4 columns, X, Y, Z, Property (numerical)."

have you tried the example of the v.vol.rst-manual?

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

Francois_Chartier · May 12, 2018, 7:57pm

I have reimported the data set and have linked it to a new db as i assume this is the issue.

I rerun v.vol.rst and i have to define ‘cat’ as the name of col that contains the w-values attribute to interpolate. These match the attribute i want to interpolate from the db.

Now the issue is the data is not in the region…zero points in the region.

Sent from Mail for Windows 10

From: Helmut Kudrnovsky
Sent: May 12, 2018 3:43 PM
To: grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] Vector to 2D Raster

v.info -c map=BHMat1All@Toronto

ERROR: Database connection for map BHMat1All@Toronto is not defined in DB

file

just before you mentioned:

“I have 4 columns, X, Y, Z, Property (numerical).”

have you tried the example of the v.vol.rst-manual?

best regards

Helmut

–

Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

grass-user mailing list

grass-user@lists.osgeo.org

https://lists.osgeo.org/mailman/listinfo/grass-user

	Virus-free. www.avast.com

hellik · May 12, 2018, 8:05pm

Now the issue is the data is not in the region....zero points in the region.

you have to set the computational region:

https://grass.osgeo.org/grass75/manuals/g.region.html

have a look into the v.vol.rst-manual how to set also the 3D region

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

Francois_Chartier · May 12, 2018, 8:23pm

Bare with me. I guessed the issue was for the region was not defined with the 3d region. So I ran g.region and define the region in NS, EW, and Top and bottom. Now V.vol.rst is running…crushing numbers.

Sent from Mail for Windows 10

From: Francois Chartier
Sent: May 12, 2018 3:57 PM
To: Helmut Kudrnovsky; grass-user@lists.osgeo.org
Subject: RE: [GRASS-user] Vector to 2D Raster

I have reimported the data set and have linked it to a new db as i assume this is the issue.

I rerun v.vol.rst and i have to define ‘cat’ as the name of col that contains the w-values attribute to interpolate. These match the attribute i want to interpolate from the db.

Now the issue is the data is not in the region…zero points in the region.

Sent from Mail for Windows 10

From: Helmut Kudrnovsky
Sent: May 12, 2018 3:43 PM
To: grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] Vector to 2D Raster

v.info -c map=BHMat1All@Toronto

ERROR: Database connection for map BHMat1All@Toronto is not defined in DB

file

just before you mentioned:

“I have 4 columns, X, Y, Z, Property (numerical).”

have you tried the example of the v.vol.rst-manual?

best regards

Helmut

–

Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

grass-user mailing list

grass-user@lists.osgeo.org

https://lists.osgeo.org/mailman/listinfo/grass-user

	Virus-free. www.avast.com

Francois_Chartier · May 12, 2018, 9:22pm

Not enough disk space for temp file… too much data i guess for my computer. Is there a possibility to write temp files to an external hard drive?

Sent from Mail for Windows 10

From: Francois Chartier
Sent: May 12, 2018 4:23 PM
To: Helmut Kudrnovsky; grass-user@lists.osgeo.org
Subject: RE: [GRASS-user] Vector to 2D Raster

Bare with me. I guessed the issue was for the region was not defined with the 3d region. So I ran g.region and define the region in NS, EW, and Top and bottom. Now V.vol.rst is running…crushing numbers.

Sent from Mail for Windows 10

From: Francois Chartier
Sent: May 12, 2018 3:57 PM
To: Helmut Kudrnovsky; grass-user@lists.osgeo.org
Subject: RE: [GRASS-user] Vector to 2D Raster

I have reimported the data set and have linked it to a new db as i assume this is the issue.

I rerun v.vol.rst and i have to define ‘cat’ as the name of col that contains the w-values attribute to interpolate. These match the attribute i want to interpolate from the db.

Now the issue is the data is not in the region…zero points in the region.

Sent from Mail for Windows 10

From: Helmut Kudrnovsky
Sent: May 12, 2018 3:43 PM
To: grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] Vector to 2D Raster