[GRASS-user] Problems to show 3d boreholes
Alec Ventura
alecventura at hotmail.com
Tue Jan 19 13:23:53 PST 2016
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 at MapAt 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 at MapAnd at this point i change to 3d and it gives me the follow error: "Loading vector map <test1 at Map> (points) failed No points from vector map <test1 at Map> fall within currentregion"
... 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 grassVery tnks,Alec Ventura.From: carlos.grohmann at gmail.com
Date: Mon, 18 Jan 2016 23:06:03 -0200
Subject: Re: [GRASS-user] Problems to show 3d boreholes
To: alecventura at hotmail.com
CC: grass-user at 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 at 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 at PERMANENT2 - run: v.info map=Pocos_UTM at PERMANENT
ps: i just have imported the .shp file
From: kratochanna at gmail.com
Date: Tue, 12 Jan 2016 21:30:25 -0500
Subject: Re: [GRASS-user] Problems to show 3d boreholes
To: alecventura at hotmail.com
CC: grass-user at lists.osgeo.org
On Tue, Jan 12, 2016 at 7:15 PM, Alec Ventura <alecventura at 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 at Map1>No points from vector map <bore2 at Map1> fall within currentregion
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 at 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 at lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
_______________________________________________
grass-user mailing list
grass-user at 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.comhttp://orcid.org/0000-0001-5073-5572
________________
Can’t stop the signal.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20160119/1782a0a5/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2016-01-18 at 22.49.40.png
Type: image/png
Size: 103490 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20160119/1782a0a5/attachment-0001.png>
More information about the grass-user
mailing list