Re: [GRASS-user] Move vertex in v.digit releasebranch
Hi Martin, Martin Landa wrote: Hi, 2008/12/21 Casagrande Luca thedokto...@gmail.com: Hello everybody! I am testing the releasebranch with the wxpython interface and I found this little problem: - open roads from spearfish and start edit it - select the move vertex tool - after you select a vertex with the left click, if the button is pressed again by mistake, the display freeze. should be fixed in r35035 (relbr64). Thanks, it's fine now. I'll keep testing.. Luca -- View this message in context: http://n2.nabble.com/Move-vertex-in-v.digit-releasebranch-tp1886429p2102368.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Help with reprojection
Hi, I'm new to GRASS, and am having trouble with reprojection. I've spent two days trying to figure this out, and I feel that I'm close, but reprojection is still failing. Could anyone clue me in what I'm missing? The input data set is Natural Earth II map of the world http://www.nacis.org/data/NE2/NE2_WDB1.jpg 16200x8100 JPG in plate carrée projection. I need to reproject it into Lambert azimuthal equal-area projection centered on N0 E15. I'm using WinGRASS 6.3.0. First, I've created GRASS location geographic with following region and projection: g.region -p projection: 99 (Equidistant Cylindrical (Plate Caree)) zone: 0 datum: wgs84 ellipsoid: wgs84 north: 8100 south: 0 west: 0 east: 16200 nsres: 1 ewres: 1 rows: 8100 cols: 16200 cells: 13122 g.proj -p datumtrans=0 PROJ_INFO- name : Equidistant Cylindrical (Plate Caree) datum : wgs84 towgs84: 0.000,0.000,0.000 proj : eqc ellps : wgs84 a : 6378137.00 es : 0.0066943800 f : 298.2572235630 lat_0 : 0.00 lat_ts : 0.00 lon_0 : 0.00 -PROJ_UNITS unit : meter units : meters meters : 1.0 I then imported the JPG map, which got split into color channels and composed the channels together into single map. After, I created a new location lambert: g.region -p projection: 99 (Lambert Azimuthal Equal Area) zone: 0 datum: wgs84 ellipsoid: wgs84 north: 8100 south: 0 west: 0 east: 16200 nsres: 1 ewres: 1 rows: 8100 cols: 16200 cells: 13122 g.proj -p datumtrans=0 -PROJ_INFO- name : Lambert Azimuthal Equal Area datum : wgs84 towgs84: 0.000,0.000,0.000 proj : laea ellps : wgs84 a : 6378137.00 es : 0.0066943800 f : 298.2572235630 lat_0 : 0.00 lon_0 : 15.00 x_0: 0.00 y_0: 0.00 -PROJ_UNITS unit : meter units : meters meters : 1.0 Finally, when I try reprojection I get Input raster map is outside current region error: r.proj input=NE2 location=geographic mapset=PERMANENT method=cubic Input Projection Parameters: +proj=eqc +lat_0=0.00 +lat_ts=0.00 +lon_0=0.00 +a=6378137 +rf=298.257223563 +no_defs +towgs84=0.000,0.000,0.000 Input Unit Factor: 1 Output Projection Parameters: +proj=laea +lat_0=0.00 +lon_0=15.00 +x_0=0.00 +y_0=0.00 +a=6378137 +rf=298.257223563 +no_defs +towgs84=0.000,0.000,0.000 Output Unit Factor: 1 Input raster map is outside current region How can I make this work? Thanks, Alex ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.quantile: what it exaclty do?
On Fri, Jan 2, 2009 at 4:55 AM, Jarek Jasiewicz jar...@amu.edu.pl wrote: Hi I tirend used r.quantile on SRTM data for Poland (min value about to 0, max value about to 2500 m a.s.l) r.quantile input=pol...@permanent quantiles=4 percentiles=0.001,0.01,0,1,0.25,0.50,0.75,0.90,0.99,0.999 bins=100 and I recived: 0:0.00:0.00 1:0.001000:0.00 2:0.01:0.00 3:0.25:0.00 4:0.50:1.00 5:0.75:2.00 6:0.90:3.00 7:0.99:3.00 8:0.999000:3.00 9:1.00:3.00 It looks like the results are divided by 1000 and rounded to the nearest integer. Did I something wrong? Jarek I thought that the user supplies one of [quantiles=] | [percentiles=] ... could that be related to the odd output. I have verified (with R) that r.quantile can compute correct quantiles... however, on my system I need to set bins=100 or so, as the higher i set 'bins' odd things would happen. Glynn should know for sure. Dylan ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] RE: NVIZ problem
Many thanks for the suggestions. I tried for more hours than I care to admit to solve this without success. I eventually removed the HD4870 and exchanged it for my old nvidia 8800GTS. GRASS now works in intrepid64 (although it is 6.2!). I could not get 6.4 to compile but I probably gave up too easily (it was late last night) and I will try again shortly. Interesting GRASS also works on my Vista64 bit machine with the HD4870 catalyst drivers installed (although installation of the catalyst drivers was not a simple matter!). Why does working in 64bit seem to cause so many driver issues! Both of my machines have 4Gb ram and I simply want/need to use it. I always knew Vista was a pain but this has also shaken my trust in linux. Anyhow I now have both systems working so I am happy - for now :). Thanks again for all the help. Darrel Prof Darrel Maddy Chair of Quaternary Sciences Newcastle University ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.quantile: what it exaclty do?
Dylan Beaudette pisze: On Fri, Jan 2, 2009 at 4:55 AM, Jarek Jasiewicz jar...@amu.edu.pl wrote: Hi I tirend used r.quantile on SRTM data for Poland (min value about to 0, max value about to 2500 m a.s.l) r.quantile input=pol...@permanent quantiles=4 percentiles=0.001,0.01,0,1,0.25,0.50,0.75,0.90,0.99,0.999 bins=100 and I recived: 0:0.00:0.00 1:0.001000:0.00 2:0.01:0.00 3:0.25:0.00 4:0.50:1.00 5:0.75:2.00 6:0.90:3.00 7:0.99:3.00 8:0.999000:3.00 9:1.00:3.00 It looks like the results are divided by 1000 and rounded to the nearest integer. Did I something wrong? Jarek I thought that the user supplies one of [quantiles=] | [percentiles=] ... could that be related to the odd output. I have verified (with R) that r.quantile can compute correct quantiles... however, on my system I need to set bins=100 or so, as the higher i set 'bins' odd things would happen. Glynn should know for sure. Dylan Thanks it really seems that the problem was in default number of bins ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Display thematic vector map - WinGRASS 6.3 (Vista)
Hello I'm a GRASS GIS novice and have recently downloaded and successfully installed GRASS 6.3.0 directly onto MS-Windows Vista. I have been playing with the GIS manager GUI and have run into a problem trying to create thematic Vector maps and display them in the map display. I am using the North Carolina data set (example data set provided) and a variety of vector maps to try and create graduated color, point and line maps to get a feel of how it works. I have been adding the thematic layer and specifying all the parameters correctly in the GUI. The actual binning of the numeric data is occurring properly because a new window will pop-up with the correct legend for the parameters specified, however nothing will display in the map display window. I do get a message in the command line shell saying No such monitor exists 'PNG' . I've tried to run d.vect.thematic through the command line shell and get an error saying must open display monitor. I'm not sure what the problem is. Can anyone help? Thanks K. Brookes Toronto, Canada _ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Volumes of intersections
Is it possible to use grass to calculate the volume of a part of a 3-D vector object? What I mean is if I have a mountain, for example, that is defines in one vector map. Then, in another map, I define an area (a plane for example). This represents a slice through the mountain. Is it possible to calculate the volume of the mountain that lies above the second plane? (how much rock is above the slice through it?) It seems like it should be doable, but I am not really sure how to get going. Any suggestions? Thanks, --Adam ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: GRASS-user] Help with reprojection
On Jan 2, 2009, at 4:38 PM, grass-user-requ...@lists.osgeo.org wrote: Message: 2 Date: Fri, 2 Jan 2009 21:04:44 +0200 From: Alex Bernstein pofi...@gmail.com Subject: [GRASS-user] Help with reprojection To: grass-user@lists.osgeo.org Message-ID: bf69017a0901021104n29f77efal2206f53128b08...@mail.gmail.com Content-Type: text/plain; charset=ISO-8859-1 Hi, I'm new to GRASS, and am having trouble with reprojection. I've spent two days trying to figure this out, and I feel that I'm close, but reprojection is still failing. Could anyone clue me in what I'm missing? The input data set is Natural Earth II map of the world http://www.nacis.org/data/NE2/NE2_WDB1.jpg 16200x8100 JPG in plate carrée projection. I need to reproject it into Lambert azimuthal equal-area projection centered on N0 E15. I'm using WinGRASS 6.3.0. First, I've created GRASS location geographic with following region and projection: g.region -p projection: 99 (Equidistant Cylindrical (Plate Caree)) zone: 0 datum: wgs84 ellipsoid: wgs84 north: 8100 south: 0 west: 0 east: 16200 nsres: 1 ewres: 1 rows: 8100 cols: 16200 cells: 13122 g.proj -p datumtrans=0 PROJ_INFO- name : Equidistant Cylindrical (Plate Caree) datum : wgs84 towgs84: 0.000,0.000,0.000 proj : eqc ellps : wgs84 a : 6378137.00 es : 0.0066943800 f : 298.2572235630 lat_0 : 0.00 lat_ts : 0.00 lon_0 : 0.00 -PROJ_UNITS unit : meter units : meters meters : 1.0 I then imported the JPG map, which got split into color channels and composed the channels together into single map. After, I created a new location lambert: g.region -p projection: 99 (Lambert Azimuthal Equal Area) zone: 0 datum: wgs84 ellipsoid: wgs84 north: 8100 south: 0 west: 0 east: 16200 nsres: 1 ewres: 1 rows: 8100 cols: 16200 cells: 13122 g.proj -p datumtrans=0 -PROJ_INFO- name : Lambert Azimuthal Equal Area datum : wgs84 towgs84: 0.000,0.000,0.000 proj : laea ellps : wgs84 a : 6378137.00 es : 0.0066943800 f : 298.2572235630 lat_0 : 0.00 lon_0 : 15.00 x_0: 0.00 y_0: 0.00 -PROJ_UNITS unit : meter units : meters meters : 1.0 Finally, when I try reprojection I get Input raster map is outside current region error: r.proj input=NE2 location=geographic mapset=PERMANENT method=cubic Input Projection Parameters: +proj=eqc +lat_0=0.00 +lat_ts=0.00 +lon_0=0.00 +a=6378137 +rf=298.257223563 +no_defs +towgs84=0.000,0.000,0.000 Input Unit Factor: 1 Output Projection Parameters: +proj=laea +lat_0=0.00 +lon_0=15.00 +x_0=0.00 +y_0=0.00 +a=6378137 +rf=298.257223563 +no_defs +towgs84=0.000,0.000,0.000 Output Unit Factor: 1 Input raster map is outside current region How can I make this work? Your problem is that this map is an image, not a GIS map. It is in what is called an xy projection in GRASS. That is, it is in no projection at all. The coordinates of each cell are simply the xy coordinates from the corner, not their position in the world. In a 'real' GIS map, each cell would have some kind of real-world geographic coordinates (latlon, UTM, etc). To get from a graphic image like this to geospatial data, you need to georectify the map. It would be a good idea for you to get an intro GIS book so that you can understand what is going on with this if it is unfamiliar to you. The steps you need to follow in GRASS are... 1) create an xy location 2) import your jpg map into the xy location 3) open GRASS in your Lambert location 4) start the georectifier and georectify the map from the xy to Lambert location. Looking at the jpg, it has no latlon or other geographic markers. This will make it difficult to georectify accurately. If you are looking for nice geospatial data of world topography, it might be better for you to download some data rather than trying to georectify a jpg image. GTOPO30 or similar scale (1km or larger) topography can be downloaded from a variety of places. Check the GRASS site to links to geospatial data. Also check out the GLCF http://glcf.umiacs.umd.edu/ . For global-scale imagery, try MODIS satellite data. The good thing about these data sets is that they are already georectified and can be imported directly into GRASS as GIS data. You can then easily reproject to something else. Michael ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Volumes of intersections
I don't know the latest devlopments on the 3d vector system, but I would convert the vector to rasters (2d or 3d) and then use map algebra... giovanni 2009/1/3 Adam Dershowitz adershow...@exponent.com: Is it possible to use grass to calculate the volume of a part of a 3-D vector object? What I mean is if I have a mountain, for example, that is defines in one vector map. Then, in another map, I define an area (a plane for example). This represents a slice through the mountain. Is it possible to calculate the volume of the mountain that lies above the second plane? (how much rock is above the slice through it?) It seems like it should be doable, but I am not really sure how to get going. Any suggestions? Thanks, --Adam ___ 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] r.quantile: what it exaclty do?
On Fri, Jan 2, 2009 at 8:40 PM, Dylan Beaudette dylan.beaude...@gmail.com wrote: On Fri, Jan 2, 2009 at 4:55 AM, Jarek Jasiewicz jar...@amu.edu.pl wrote: Hi I tirend used r.quantile on SRTM data for Poland (min value about to 0, max value about to 2500 m a.s.l) r.quantile input=pol...@permanent quantiles=4 percentiles=0.001,0.01,0,1,0.25,0.50,0.75,0.90,0.99,0.999 bins=100 and I recived: 0:0.00:0.00 1:0.001000:0.00 2:0.01:0.00 3:0.25:0.00 4:0.50:1.00 5:0.75:2.00 6:0.90:3.00 7:0.99:3.00 8:0.999000:3.00 9:1.00:3.00 It looks like the results are divided by 1000 and rounded to the nearest integer. Did I something wrong? Jarek I thought that the user supplies one of [quantiles=] | [percentiles=] ... could that be related to the odd output. I have verified (with R) that r.quantile can compute correct quantiles... however, on my system I need to set bins=100 or so, as the higher i set 'bins' odd things would happen. Glynn should know for sure. Can you please file a bug report on this? Preferable with a Spearfish/NC reproducible example? Thanks Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] RE: NVIZ problem
On Fri, Jan 2, 2009 at 9:17 PM, Darrel Maddy darrel.ma...@newcastle.ac.uk wrote: Many thanks for the suggestions. I tried for more hours than I care to admit to solve this without success. I eventually removed the HD4870 and exchanged it for my old nvidia 8800GTS. GRASS now works in intrepid64 (although it is 6.2!). I could not get 6.4 to compile but I probably gave up too easily (it was late last night) and I will try again shortly. Please let us know what failed. Maybe we can quickly help. Interesting GRASS also works on my Vista64 bit machine with the HD4870 catalyst drivers installed (although installation of the catalyst drivers was not a simple matter!). Why does working in 64bit seem to cause so many driver issues! Both of my machines have 4Gb ram and I simply want/need to use it. I always knew Vista was a pain but this has also shaken my trust in linux. Perhaps an Ubuntu issue here? I am using 64bit Linux without troubles for years meanwhile. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user