RE: [GRASS-user] r3.in.ascii - Difficulties understanfing format
Hi all, I still have a question concerning the visualization of 3D raster in Nviz. The region settings are the following: projection: 0 (x,y) zone: 0 north: 7 south: 0 west: 0 east: 7 top:11. bottom: 0. nsres: 1 nsres3: 1 ewres: 1 ewres3: 1 tbres: 1 rows: 7 rows3: 7 cols: 7 cols3: 7 depths: 11 cells: 49 3dcells:539 Type of Map: 3d cell Number of Categories: 0 Data Type:double Rows: 7 Columns: 7 Depths: 11 Total Cells: 539 Projection: x,y (zone 0) N: 7S: 0 Res: 1 E: 7W: 0 Res: 1 T: 11B: 0 Res: 1 Range of data: min = 0 max = 3 But I can't see anything in NVIZ. I don't see the wireframe that is sometimes visible when moving the viewpoint, and when I ry to add a volume/isosurface, it shows nothing. What did I miss here? Thanks, Tom -Original Message- From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Tom van der Putte Sent: woensdag 14 januari 2009 10:23 To: 'Nikos Alexandris' Cc: 'grass-user' Subject: RE: [GRASS-user] r3.in.ascii - Difficulties understanfing format Nikos, Thanks very much for the 101 on rasters in Grass, helped me out big time. No worries on the length of the answer, the more the better ;) Tom -Original Message- From: Nikos Alexandris [mailto:nikos.alexand...@felis.uni-freiburg.de] Sent: dinsdag 13 januari 2009 16:50 To: t...@vdputte.nl Cc: grass-user Subject: Re: [GRASS-user] r3.in.ascii - Difficulties understanfing format Tom, my apologies for writing something long, and repeat within the text some times things already mentioned or known. I did so on purpose -- just to practice and see if I understand the concept myself. Ah, please keep in your posts Cc to the grass-user mailing list. Others will definitely benefit from our mistakes or the correct examples we discuss. On Tue, 2009-01-13 at 15:13 +0100, t...@vdputte.nl wrote: Hi Nikos, Thanks for the quick response! So what I gather from your answer is that the region-definitions are set seperately (I'm guessing in the mapset?) and that you can only display data that falls within the boundaries of this region? Correct. Just to make the big picture clear: Location = One coordinate system Many mapsets with separate region settings for each. In addition, an important and very useful characteristic of GRASS' raster modules (r.*) is that they operate *only* within the defined region. The my problem transforms a bit,in the sense that I want to display a 3D array/raster. Ultimately, this needs to be georeferenced, but for testing purposes I just want to display this raster in Nviz. So how can I create a general Mapset in which I can display a 3D raster with a resolution of 1 (e.g North-South, E-W and T-B all 3 - 0 with a raster of 3x3x3) If you have non-georeferenced spatial data you can create a generic XY location and set the region to match the extent of your imported (in the xy location) raster no-matter how big/small it is. [see [1][2] and search on the mailing list [3] for how to create an xy location -- it's a standard question when one starts to explore data with grass] # set region g.region rast=YourRaster -p # or: g.region rast=YourRaster -pa ## read g.region help for more It should, I think, obtain automatically the required rows columns parameters which should correspond to the columns rows of your raster. If not just grab the resolution of your raster somehow (e.g. gdalinfo YourRaster) and set them manually. I think that south north should match columns rows. east west are irrelevant and, I think, you can set the latter to 0. Hmmm... talking about 3D I think the correct term is voxels, and hence you can set the vertical resolution to 1 as well. But you can experiment with whatever you think you need... The important, in order to *see* your data is here the (x=)columns, (y=)rows and (z=)levels. Now, the *how* they will appear will be affected by the defined x,y,z resolution. Since the units in an XY location are pixels it's correct to set the resolution to 1. A few more notes concerning georeferenced maps/locations While the resolution of the georeferenced map itself is as is (let's say you have the CORINE land cover of some European country, the raster version of 250m pixel size), and refers to the units used by the defined coordinate system, the resolution you define with g.region will affect the output of most raster processing modules. # corine_250m # match region-extent to the corine_250m map g.region rast=corine_250 res=500 -pa # produce a 500m corine map r.mapcalc corine_500m=corine_250m # so simple ;-) I hope you understand my problem. Thanks, Tom [...] I hope I
RE: [GRASS-user] r3.in.ascii - Difficulties understanfing format
I was a little bit to hasty, forgot to mention that these are the region settings, projection: 0 (x,y) zone: 0 north: 7 south: 0 west: 0 east: 7 top:11. bottom: 0. nsres: 1 nsres3: 1 ewres: 1 ewres3: 1 tbres: 1 rows: 7 rows3: 7 cols: 7 cols3: 7 depths: 11 cells: 49 3dcells:539 --- And this the info on the 3D raster I wanted to visualize 3D raster info Type of Map: 3d cell Number of Categories: 0 Data Type:double Rows: 7 Columns: 7 Depths: 11 Total Cells: 539 Projection: x,y (zone 0) N: 7S: 0 Res: 1 E: 7W: 0 Res: 1 T: 11B: 0 Res: 1 Range of data: min = 0 max = 3 Cheers, Tom ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] DEM export to GeoTiff
run r.out.gdal -c the -c flag supresses exporting of these long colortables, they cause problems. -c is not a valid flag ... __ Yahoo! Canada Toolbar: Search from anywhere on the web, and bookmark your favourite sites. Download it now at http://ca.toolbar.yahoo.com.___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Display thematic vector map - WinGRASS 6.3 (Vista)
On 15/01/09 05:22, Michael Barton wrote: I'm pretty sure that this is a Vista issue, though one that has been solved in 6.4. And as the OP is using 6.3, he doesn't have access to that solution, yet. So, we'll have to wait for the version of WinGRASS to decide whether there might be additional issues in Vista. Moritz Michael On Jan 14, 2009, at 7:33 PM, Matt wrote: Ken, You are not alone. I tried the thematic vector example on page 171 of Neteler Mitasova and also encountered the You must open a display monitor message. I'm running 6.3.0-4 on Vista. Matt On Tue, Jan 6, 2009 at 12:51 PM, Moritz Lennert mlenn...@club.worldonline.be wrote: On 06/01/09 17:51, Michael Barton wrote: Thanks for the correction. WinGRASS supports PNG output, not a monitor. d.vect.thematic was rewritten to use direct rendering through PNG, so it should work on Windows. This was done in 6.4, so it is not available in winGRASS, yet, which is based on 6.3. Just replacing the script doesn't work, either as it use the fs= parameter of v.db.connect which wasn't introduced until 6.4. So, another strident call for the a 6.4 winGRASS. Moritz ___ 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
Re: [GRASS-user] DEM export to GeoTiff
Stefan Mietke wrote: run r.out.gdal -c the -c flag supresses exporting of these long colortables, they cause problems.. -c is not a valid flag ... The -c flag for r.out.gdal is new in version 6.4. I don't know what version you use, probably 6.3, there -c is indeed not a valid flag for r.out.gdal. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] DEM CGIAR
Hi all, I downloaded srtm-cgiar data (Lat long datum WGS84). I need to use these geotiff in UTM33 N (datum EUR 50). By using r.proj,the resulting raster shows a relevant offset from correct position (2000 m in N-S direction). It's due to different DATUM? Thanks to All regards Edmondo ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] DEM export to GeoTiff
Stefan Mietke wrote: run r.out.gdal -c the -c flag supresses exporting of these long colortables, they cause problems.. -c is not a valid flag ... The -c flag for r.out.gdal is new in version 6.4. I don't know what version you use, probably 6.3, there -c is indeed not a valid flag for r.out.gdal. Of course I have version 6.3! 6.4 is not stable yet, isn't it? __ Looking for the perfect gift? Give the gift of Flickr! http://www.flickr.com/gift/___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] DEM export to GeoTiff
Stefan Mietke wrote: [...] Of course I have version 6.3! 6.4 is not stable yet, isn't it? On the download page [1] it says that GRASS 6.3 is a testing version and GRASS 6.4 is the next stable version. GRASS 6.4 is available as release candidate 2, with binaries for Mac and Linux. If possible, upgrade to 6.4, it really has a lot of improvements over 6.3. [1] http://grass.osgeo.org/download/index.php ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] DEM CGIAR
On Thu, Jan 15, 2009 at 12:48 PM, Edmondo Elisei edmondo.eli...@gmail.com wrote: Hi all, I downloaded srtm-cgiar data (Lat long datum WGS84). I need to use these geotiff in UTM33 N (datum EUR 50). By using r.proj,the resulting raster shows a relevant offset from correct position (2000 m in N-S direction). It's due to different DATUM? Possibly also the ellipsoid you selected is wrong? please post the output of g.proj -w of your UTM location. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] DEM CGIAR
the output of g.proj for my UTM33 N location: PROJCS[UTM Zone 33, Northern Hemisphere, GEOGCS[international, DATUM[European_Datum_1950, SPHEROID[International_1924,6378388,297], TOWGS84[-87.000,-98.000,-121.000]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,15], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1]] ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] DEM export to GeoTiff
Okay, I can produce now some readable tiffs, but the elevation data is lost finally. I have now three RGB values instead of one elevation datum. Or is it not possible to export DEMs with GRASS at the moment? Kind regards, Stefan __ Looking for the perfect gift? Give the gift of Flickr! http://www.flickr.com/gift/___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Display thematic vector map - WinGRASS 6.3 (Vista)
Probably the d.vector.thematic script can be copied from 6.4 and run in 6.3 with little or no modification. Michael __ C. Michael Barton, Professor of Anthropology Director of Graduate Studies School of Human Evolution Social Change Center for Social Dynamics Complexity Arizona State University Tempe, AZ 85287-2402 USA voice: 480-965-6262; fax: 480-965-7671 www: http://www.public.asu.edu/~cmbarton On Jan 15, 2009, at 3:27 AM, Moritz Lennert wrote: On 15/01/09 05:22, Michael Barton wrote: I'm pretty sure that this is a Vista issue, though one that has been solved in 6.4. And as the OP is using 6.3, he doesn't have access to that solution, yet. So, we'll have to wait for the version of WinGRASS to decide whether there might be additional issues in Vista. Moritz Michael On Jan 14, 2009, at 7:33 PM, Matt wrote: Ken, You are not alone. I tried the thematic vector example on page 171 of Neteler Mitasova and also encountered the You must open a display monitor message. I'm running 6.3.0-4 on Vista. Matt On Tue, Jan 6, 2009 at 12:51 PM, Moritz Lennert mlenn...@club.worldonline.be wrote: On 06/01/09 17:51, Michael Barton wrote: Thanks for the correction. WinGRASS supports PNG output, not a monitor. d.vect.thematic was rewritten to use direct rendering through PNG, so it should work on Windows. This was done in 6.4, so it is not available in winGRASS, yet, which is based on 6.3. Just replacing the script doesn't work, either as it use the fs= parameter of v.db.connect which wasn't introduced until 6.4. So, another strident call for the a 6.4 winGRASS. Moritz ___ 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
Re: [GRASS-user] DEM export to GeoTiff
Finally I solved the problem with the –c flag … r.out.gdal –c mydem mydem.tif This will keep all elevation data in the file. An upgrade to GRASS 6.4 (RC2) is necessary. __ Connect with friends from any web browser - no download required. Try the new Yahoo! Canada Messenger for the Web BETA at http://ca.messenger..yahoo.com/webmessengerpromo.php___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] DEM CGIAR
On Thu, Jan 15, 2009 at 2:41 PM, Edmondo Elisei edmondo.eli...@gmail.com wrote: the output of g.proj for my UTM33 N location: PROJCS[UTM Zone 33, Northern Hemisphere, GEOGCS[international, DATUM[European_Datum_1950, SPHEROID[International_1924,6378388,297], TOWGS84[-87.000,-98.000,-121.000]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,15], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1]] Looks ok to me. If also the LatLong is ok, then probably your other data you compare to are shifted? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] testing
Just a test! ~ Eric. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Any success converting grc files?
Hi gang, Has anyone tried converting tcl/tk grc workspace files to wxgui workspace format recently? I have a grc I've been using with gis.m for some time, and converting it in the wxgui seems to hang and never complete (5 minutes, no success) - or maybe I'm too impatient. ~ Eric. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] r3.in.ascii - Difficulties understanfing format
Tom van der Putte wrote: I still have a question concerning the visualization of 3D raster in Nviz. But I can't see anything in NVIZ. I don't see the wireframe that is sometimes visible when moving the viewpoint, and when I try to add a volume/isosurface, it shows nothing. What did I miss here? when you get to the point where you can see the wireframe box when you move the viewport puck, but then nothing is rendered, what you have to do is add some new iso-surfaces. Try the Custom volume creation example here: http://grass.osgeo.org/wiki/Help_with_3D#Examples Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] starting wxnviz
Hi list, I'm starting to use grass 64 and the GUI has been improved a lot; the new icons are really nice too. Unfortunately, I've been trying to use wxnviz, but haven't been able to do so. I've compiled grass with wxwidgets, and nviz (trad) works fine. I'm attaching the problem at the end. Best, Jaime Traceback (most recent call last): File /usr/local/grass-6.4.0RC2/etc/wxpython/gui_modules/toolbars.py, line 253, in OnSelectTool self.mapdisplay.AddToolbar(nviz) File /usr/local/grass-6.4.0RC2/etc/wxpython/gui_modules/mapdisp.py, line 2769, in AddToolbar import nviz File /usr/local/grass-6.4.0RC2/etc/wxpython/gui_modules/nviz.py, line 49, in module import nviz_mapdisp File /usr/local/grass-6.4.0RC2/etc/wxpython/gui_modules/nviz_mapdisp.py, line 40, in module import grass6_wxnviz as wxnviz ImportError : No module named grass6_wxnviz ¡Todo sobre la Liga Mexicana de fútbol! Estadisticas, resultados, calendario, fotos y más:lt; http://espanol.sports.yahoo.com/___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
