First i want to tnks all the folks who spare some time to help the others.
I do some tests after your e-mail Carlos and i find out in which step is the problem:
After i import my .shp file i execute the commands: (# on left line is a python command and @ is a command console command)
import sys, os, numpy, argparse
import grass.script as g
import grass.script.setup as gsetup
g.run_command(‘v.buffer’, input=‘Pocos_UTM’, output=‘mapaComBuffer’, type=‘point’, distance=‘1000’, flags = ‘t’)
g.run_command(‘v.extract’,input=mapWithBuffer, cats=1, output=‘tmp’,type=‘area’, quiet=True, overwrite=True)
@ d.vect map=tmp@Map
At this point i change to 3d and everthing is ok so i continue:
g.run_command(‘v.extrude’, input=‘tmp’, output=“test1”, zshift=300,height=1000,overwrite=True)
@ d.vect map=test1@Map
And at this point i change to 3d and it gives me the follow error:
“Loading vector map test1@Map (points) failed
No points from vector map test1@Map fall within current
region”
… And while i was writing this e-mail i notice that the borehole appeared o.O, so it seens it have worked!!
By the way, how do you put the surface? you said that you have made a crude interpolation, i don’t understand. Sorry i still have too little knowledge with grass
Very tnks,
Alec Ventura.
From: carlos.grohmann@gmail.com
Date: Mon, 18 Jan 2016 23:06:03 -0200
Subject: Re: [GRASS-user] Problems to show 3d boreholes
To: alecventura@hotmail.com
CC: grass-user@lists.osgeo.org
Hello Alec
I did some tests with your file and it worked fine.
Importing the vector into a UTM22south region indeed results in ‘zone invalid’ from v.info, but it works. (even with the ‘override dataset projection’ option)
(note: the data is on UTM zone 22J - Santa Catarina State)
In the 3D view, initially, the symbol size was really small (100 meters - map units) and the extent of the points is quite large 212x367 kilometers) so it was hard to see the vector points. But they are there. In the ‘Appearance’ tab, you can increase this size.
I attached a screenshot of the point (size=1000) and a crude interpolation of the PROF_BASE column.
Given the spatial distribution of the points, I suggest you create a mask of the convex hull (v.hull + v.to.rast + r.mask) of your points and restrict the interpolation to this area.
You might also want to double check the encoding on the DBF file. I didn’t got any of the Portuguese special characters correctly (I tend to run away from using characters like ‘ç’ ‘ã’, etc).
best
Carlos
On Mon, Jan 18, 2016 at 10:09 PM, Alec Ventura <alecventura@hotmail.com> wrote:
So, i as able to find some time to work here again…
I have created a new location using the shape file as source for it, and when i choose the file, it successfully recognize the location as UTM 22s.
But right after GRASS had opened and i just imported the shape to my viewport and see the metadata of shape, it again says : Projection: UTM (zone invalid)
My steps:
created a location based on shape file attached,
1 - run: d.vect map=Pocos_UTM@PERMANENT
2 - run: v.info map=Pocos_UTM@PERMANENT
ps: i just have imported the .shp file
From: kratochanna@gmail.com
Date: Tue, 12 Jan 2016 21:30:25 -0500
Subject: Re: [GRASS-user] Problems to show 3d boreholes
To: alecventura@hotmail.com
CC: grass-user@lists.osgeo.org
On Tue, Jan 12, 2016 at 7:15 PM, Alec Ventura <alecventura@hotmail.com> wrote:
Im trying to draw some 3d boreholes, i have based my attempts on gene solution:
http://osgeo-org.1560.x6.nabble.com/Automatic-3D-geological-boreholes-representation-automate-v-extrude-from-a-table-my-solution-in-Pythn-td4978801.html
I believe im on the right path, but by some reason when i change to 3d it gave me the fallow message:
Error loading vector map bore1@Map1
No points from vector map bore2@Map1 fall within current
region
And it happens for all boreholes, when i access theyr metadata it gives:
| Map format: native |
| Type of map: vector (level: 2) |
|
| Number of points: 0 Number of centroids: 0 |
| Number of lines: 0 Number of boundaries: 0 |
| Number of areas: 0 Number of islands: 0 |
| Number of faces: 25 Number of kernels: 1 |
| Number of volumes: 0 Number of holes: 0 |
|
| Map is 3D: Yes |
| Number of dblinks: 1 |
|
| Projection: UTM (zone invalid) |
|
| N: 7064064.24981919 S: 7062068.91228081 |
| E: 371187.626007 W: 369196.94006096 |
| B: 148 T: 496 |
|
| Digitization threshold: 0 |
| Comment: |
| Generated by v.extrude from vector map tmp@Map1 |
| ±---------------------------------------------------------------------------+ |
Im using UTM 22s, but on metadata it says UTM (zone invalid).
But on 2d, it appears on viewport.
Any ideas?
Tnks.
Perhaps you could share the data to test? You can use v.pack.
Best,
Anna
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
–
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -
http://carlosgrohmann.com
http://orcid.org/0000-0001-5073-5572
Can’t stop the signal.