Re: [GRASS-user] problem with reprojection of Hydro1k data
Thank you for your answer Markus. I didn't read the note in the r.proj page. I will try again and then we will see. Stéphanie Markus Neteler wrote: Stéphanie, On Fri, Feb 13, 2009 at 2:18 PM, FAROUX STEPHANIE stephanie.far...@meteo.fr wrote: Hello, First i'd like to know if anybody already used the Hydro1k data from USGS? I try to make global maps of parameters from continental ones. The initial projection is lambert azimutal equal area, with different latitude and longitude of origins for different continents. I defined one location by continent; then i reprojected all maps in a latlon global location, with WGS84 datum. There are gaps between continents, for example between Asia and Europe and between Europe and Africa. I wonder if the problem comes: - from the data which don't perfectly cover the earth ... could be possible, I didn't try with Hydro1k data from USGS. - from my reprojection in grass How did you do that? Say, did you use the v.in.region trick mentioned in the r.proj manual page which helps a lot? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] data base element file
On 16/02/09 00:35, Keleni Raqisia wrote: Im new to Grass. Need to do an image segmentation and which requires me to group and produce a signiture file of image to run the i.smap program. I've created a DBase, Location and a mapset, permanent and also another mapset. I need instructions on how to import this raster files (images) into the mapset , which would allow me to find the name file under each current mapset when running the i.group program or the i.gensignset program. I did a manual copy and paste of this file to the current mapset but this is not use r.in.gdal (File-Import raster maps-Multiple formats using GDAL) Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem with reprojection of Hydro1k data
I did it; results are not exactly the same but it's not enough because gaps between continents are still more important. FAROUX STEPHANIE wrote: Thank you for your answer Markus. I didn't read the note in the r.proj page. I will try again and then we will see. Stéphanie Markus Neteler wrote: Stéphanie, On Fri, Feb 13, 2009 at 2:18 PM, FAROUX STEPHANIE stephanie.far...@meteo.fr wrote: Hello, First i'd like to know if anybody already used the Hydro1k data from USGS? I try to make global maps of parameters from continental ones. The initial projection is lambert azimutal equal area, with different latitude and longitude of origins for different continents. I defined one location by continent; then i reprojected all maps in a latlon global location, with WGS84 datum. There are gaps between continents, for example between Asia and Europe and between Europe and Africa. I wonder if the problem comes: - from the data which don't perfectly cover the earth ... could be possible, I didn't try with Hydro1k data from USGS. - from my reprojection in grass How did you do that? Say, did you use the v.in.region trick mentioned in the r.proj manual page which helps a lot? Markus ___ 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] mossy grass seeds
On 16/02/09 08:52, Markus Neteler wrote: On Mon, Feb 16, 2009 at 8:38 AM, Hamish hamis...@yahoo.com wrote: ... just curious: was GRASS just cheeky naming by a competing US Gov't GIS team or was there tangled roots in the early days? Here are some pointers: * GRASS History II: GRASS Roots by Jim Westervelt (In Proc. Free/Libre and Open Source Software for Geoinformatics: GIS-GRASS Users Conference 2004, Sept. 12-14, Bangkok, Thailand, 2004) And here it is: With a swelling of new customers, the need to maintain the software on multiple computers running at multiple sites, we packaged our 20 programs and called the turn-key solution GRASS the Geographic Resource Analysis Support System. The name continued a series of GIS names based on plants: SAGE, a DOS-based GIS, and MOSS (Mapping Overlay Statistical System), developed by The Bureau of Land Management (BLM). So, yes same as elm/pine etc. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Python scripting
On Fri, 13 Feb 2009 14:14:47 +, Glynn Clements wrote: It's part of the GRASS source code: lib/python/grass.py. It should get installed into $GISBASE/etc/python during the build process, and Init.sh should add that directory to PYTHONPATH. Also, it's only present in 6.4, 6.5 and 7.0, not in earlier versions. thanks. Since GRASS 6.3 integrates better with QGIS an R (at least on my PC) I wonder if I can just copy the file and use it in 6.3!? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.patch, images vs coloured elevation data, aerial photo mosaic
Hi, I have two questions, first a specific one, secondly a more general one. 1) I have some aerial photos that I want to mosaic. I am able to do this with r.patch. However, the photos only cover about half of the region. In the remainder of the region, I would like to show a shaded relief map as a background to give context to the photos. When I include the shaded relief map in the list to patch, it comes out as varying shades of red. My question is therefore, what do I have to do to the shaded relief map to turn it into an 'rgb image' that I can include in r.patch. 2) More generally, what is the best way of converting a good-looking map into a georeferenced image e.g. by converting the elevation data into rgb pixel colour data. Currently I would use d.out.file then gdal_translate, but that is dependent on the screen resolution. Cheers John -- Dr John Stevenson Postdoctoral Research Associate School of Earth, Atmospheric and Environmental Sciences Williamson Building (Room 2.42) University of Manchester Manchester M13 9PL, UK tel. +44(0)161 306 6585; fax. +44(0)161 306 9361; john.steven...@manchester.ac.uk ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem with reprojection of Hydro1k data
I try the v.in.region + v.proj + g.region + r.proj process (explained in the r.proj page) with a region whose reprojection doesn't allow to get a closed boarder: the projection of the raster doesn't work. Is it normal? Thank you FAROUX STEPHANIE wrote: I did it; results are not exactly the same but it's not enough because gaps between continents are still more important. FAROUX STEPHANIE wrote: Thank you for your answer Markus. I didn't read the note in the r.proj page. I will try again and then we will see. Stéphanie Markus Neteler wrote: Stéphanie, On Fri, Feb 13, 2009 at 2:18 PM, FAROUX STEPHANIE stephanie.far...@meteo.fr wrote: Hello, First i'd like to know if anybody already used the Hydro1k data from USGS? I try to make global maps of parameters from continental ones. The initial projection is lambert azimutal equal area, with different latitude and longitude of origins for different continents. I defined one location by continent; then i reprojected all maps in a latlon global location, with WGS84 datum. There are gaps between continents, for example between Asia and Europe and between Europe and Africa. I wonder if the problem comes: - from the data which don't perfectly cover the earth ... could be possible, I didn't try with Hydro1k data from USGS. - from my reprojection in grass How did you do that? Say, did you use the v.in.region trick mentioned in the r.proj manual page which helps a lot? Markus ___ 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 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problem with vdigit in 64 release
hi, I had some trouble with my system. Glibc wouldn't compile and I thought a messed up system might be the cause of my problems with grass. I finally fixed it, but my grass problems haven't changed. will come up with more debug info asap. thanks, kries -- View this message in context: http://n2.nabble.com/Problem-with-vdigit-in-64-release-tp2321670p2334565.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
Re: [GRASS-user] Re: Python scripting
On Mon, Feb 16, 2009 at 11:57 AM, Georg Kaspar ge...@geofs.de wrote: ... thanks. Since GRASS 6.3 integrates better with QGIS an R (at least on my PC) why is that so? The GRASS-R interface was recently updated and should be way faster than previous versions. If you could indicate the problems... Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Python scripting
Georg Kaspar wrote: It's part of the GRASS source code: lib/python/grass.py. It should get installed into $GISBASE/etc/python during the build process, and Init.sh should add that directory to PYTHONPATH. Also, it's only present in 6.4, 6.5 and 7.0, not in earlier versions. thanks. Since GRASS 6.3 integrates better with QGIS an R (at least on my PC) I wonder if I can just copy the file and use it in 6.3!? grass.py should only be needed for the Python scripts which are included in 7.0, and the version of the wxPython GUI in 6.4 and 7.0 (6.3.0 pre-dates grass.py, so the wxPython GUI wasn't using it at that point). Some of the functions in grass.py depend upon features which have been added since 6.3.0. If you have a Python script which needs grass.py, it may also depend upon such features. Also, 6.3.0's Init.sh doesn't set PYTHONPATH, so you will need to do that yourself. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] create raster from line vector
Hello everyone. A simple question, but I'm getting mad with it! I simply have to convert a vector line to raster. I've done it many times in the past, but this time I can't make it work: v.to.rast input=border output=border use=val type=line layer=1 value=1 The output is simply blank (MAX and MIN values are NULL) What stupid error am I doing? Giovanni ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.mapcalc: columns()/modulo - rounding issue, possible bugfeature ?
Hi, I stumbled over this funny behavior of the col() function of r.mapcalc (GRASS6.2 and GRASS6.4, both on Linux). According to The Book, col() returns current col of moving window. I assume that for each row of a given raster layer col() should provide the values 1,2,3,4,5, etc. as it travels through the row. To produce a map of repeating value sequences (for example the numbers 1 2 3 4 ) the module operator can be used on the col()-function: col() operator result 1 mod 4 = 1 2 mod 4 = 2 3 mod 4 = 3 4 mod 4 = 0 5 mod 4 = 1 6 mod 4 = 2 etc this approach was put into mapcalc. The 1 + makes sure that the sequence runs from 1-4 not 0 - 3: r.mapcalc modulo=(1 + ( col() % 4 ) ) A map was produced in a Spearfish location (col() should operate independently from the projection). Here it comes: For some (strange ?)reason, there are ripples in the 1 - 4 sequence. Instead of going 1234 1234 1234 ad infinitum once in a while a number is omitted: 1234 1234 234 1234 134 etc. Has anybody come across this behavior before ? Is this a feature/good-thing ? The phenomenon can be reproduced by using the above mentioned r.mapcalc command, displaying the result as a cell values layer. 1234, 1234, ..., Peter -- Dr. Peter Löwe peter.lo...@gmx.de Jetzt 1 Monat kostenlos! GMX FreeDSL - Telefonanschluss + DSL für nur 17,95 Euro/mtl.!* http://dsl.gmx.de/?ac=OM.AD.PD003K11308T4569a ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] create raster from line vector
Thnaks Nikos. I don't have neither MASKs nor raster values near 0... As you can see in my command I've set the value to 1. I start to think there's something wrong with my GRASS build... 2009/2/16 Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de: On Mon, 2009-02-16 at 15:17 +0100, G. Allegri wrote: Hello everyone. A simple question, but I'm getting mad with it! I simply have to convert a vector line to raster. I've done it many times in the past, but this time I can't make it work: v.to.rast input=border output=border use=val type=line layer=1 value=1 The output is simply blank (MAX and MIN values are NULL) What stupid error am I doing? Giovanni The first thing that fall into my mind: MASK? Perhaps it's a irrelevant thought, but I remember some strange behaviour with floating point values (in your case the pixels that compose your raster line). Could it be they are rounded down to 0? Regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] create raster from line vector
I share my vector, if someone wants to reproduce this problem: http://www.geospatial.it/allegri/bordi.tar.gz EPSG:32632 2009/2/16 G. Allegri gioha...@gmail.com: Thnaks Nikos. I don't have neither MASKs nor raster values near 0... As you can see in my command I've set the value to 1. I start to think there's something wrong with my GRASS build... 2009/2/16 Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de: On Mon, 2009-02-16 at 15:17 +0100, G. Allegri wrote: Hello everyone. A simple question, but I'm getting mad with it! I simply have to convert a vector line to raster. I've done it many times in the past, but this time I can't make it work: v.to.rast input=border output=border use=val type=line layer=1 value=1 The output is simply blank (MAX and MIN values are NULL) What stupid error am I doing? Giovanni The first thing that fall into my mind: MASK? Perhaps it's a irrelevant thought, but I remember some strange behaviour with floating point values (in your case the pixels that compose your raster line). Could it be they are rounded down to 0? Regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] create raster from line vector
On Mon, 2009-02-16 at 15:40 +0100, G. Allegri wrote: I share my vector, if someone wants to reproduce this problem: http://www.geospatial.it/allegri/bordi.tar.gz EPSG:32632 Giovanni, the shapefile's attribute table looks messed up. I just opened it with QGIS and the column names are senseless characters. In addition, all rows are filled with NULL. I don't know it's just a codepage sort of problem or the reason for your empty raster. I will try with GRASS now. If I'll find something useful I'll post back. Kind regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] create raster from line vector
On Mon, 2009-02-16 at 15:52 +0100, Nikos Alexandris wrote: On Mon, 2009-02-16 at 15:40 +0100, G. Allegri wrote: I share my vector, if someone wants to reproduce this problem: http://www.geospatial.it/allegri/bordi.tar.gz EPSG:32632 Giovanni, the shapefile's attribute table looks messed up. I just opened it with QGIS and the column names are senseless characters. In addition, all rows are filled with NULL. I don't know it's just a codepage sort of problem or the reason for your empty raster. I will try with GRASS now. If I'll find something useful I'll post back. Sorry, I can't even read it with ogrinfo. ogrinfo /home/nik/Desktop/in/bordi.shp -al | head -100 ERROR 3: fread(65004) failed on DBF file. ERROR 3: fread(65004) failed on DBF file. [...] INFO: Open of `/home/nik/Desktop/in/bordi.shp' using driver `ESRI Shapefile' successful. Layer name: bordi Geometry: Line String Feature Count: 9 Extent: (1009114.847644, 4446909.136474) - (1312106.717266, 4668313.355348) Layer SRS WKT: PROJCS[WGS_1984_UTM_Zone_32N, GEOGCS[GCS_WGS_1984, DATUM[WGS_1984, SPHEROID[WGS_1984,6378137,298.257223563]], PRIMEM[Greenwich,0], UNIT[Degree,0.017453292519943295]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,9], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[Meter,1]] ��n...@q: String (17.0) L$D=�u: String (207.0) �=?���: String (191.0) ��GMI�FYYY: String (78.0) '6:�+JX�fW`: String (235.0) g��ʚ*;: String (243.0) _*gD88Z��: String (109.0) 2u]「�...@qu�: String (119.0) �'@TN�]: String (223.0) �魗n�676t: String (117.0) +JWv���: String (136.0) \��: String (212.0) o7ʤ��: String (247.0) �_\��: String (126.0) øßøÿÛù¯L: S├⎼␋┼± (207.0) ;4OÓÓÐà4´Þ: S├⎼␋┼± (32.0) B8¬┘┴ì¶%þ¿: S├⎼␋┼± (33.0) Ǹøæ%¬?│: S├⎼␋┼± (200.0) [...] # it continues like that :-( ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.patch, images vs coloured elevation data, aerial photo mosaic
Nikos Alexandris wrote: On Mon, 2009-02-16 at 11:35 +, John Stevenson wrote: Hi, I have two questions, first a specific one, secondly a more general one. 1) I have some aerial photos that I want to mosaic. I am able to do this with r.patch. However, the photos only cover about half of the region. In the remainder of the region, I would like to show a shaded relief map as a background to give context to the photos. When I include the shaded relief map in the list to patch, it comes out as varying shades of red. My question is therefore, what do I have to do to the shaded relief map to turn it into an 'rgb image' that I can include in r.patch. Do you really want to patch the shaed relief map along with the orthophotos in one map? What about r.blend? If you insist on patching, then, I think, you would need to create 3 new versions of your shaded map which will correspond to R(ed), G(reen) and B(lue), rescale them to 0,255 (r.rescale). Then you patch them with the R, G and B mosaic's of the orthophotos and compose (r.composite) an RGB map. I would expect the shaded part to look grey-scaled since you will have the same pixel value in all of the R, G and B layers. But I am 100% sure that it will work. Thanks Nikos, In the end, I made a composite from the original shaded relief: r.composite red=nesja_shade green=nesja_shade blue=nesja_shade levels=32 output=elev_shade_comp then patched it with the photos: r.patch input=p43_trimmed,p44_trimmed,p33_trimmed,elev_shade_comp output=aerial This gave me coloured aerial photos where the data exist, and a grey shaded relief map where they don't. 2) More generally, what is the best way of converting a good-looking map into a georeferenced image e.g. by converting the elevation data into rgb pixel colour data. Currently I would use d.out.file then gdal_translate, but that is dependent on the screen resolution. What about d.out.file in=YourMap out=YourMap resolution=2 # or resolution=4 ? Good point. Cheers John Kind regards, Nikos -- Dr John Stevenson Postdoctoral Research Associate School of Earth, Atmospheric and Environmental Sciences Williamson Building (Room 2.42) University of Manchester Manchester M13 9PL, UK tel. +44(0)161 306 6585; fax. +44(0)161 306 9361; john.steven...@manchester.ac.uk ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] create raster from line vector
??? even the new one that I've uploaded? 2009/2/16 Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de: On Mon, 2009-02-16 at 15:52 +0100, Nikos Alexandris wrote: On Mon, 2009-02-16 at 15:40 +0100, G. Allegri wrote: I share my vector, if someone wants to reproduce this problem: http://www.geospatial.it/allegri/bordi.tar.gz EPSG:32632 Giovanni, the shapefile's attribute table looks messed up. I just opened it with QGIS and the column names are senseless characters. In addition, all rows are filled with NULL. I don't know it's just a codepage sort of problem or the reason for your empty raster. I will try with GRASS now. If I'll find something useful I'll post back. Sorry, I can't even read it with ogrinfo. ogrinfo /home/nik/Desktop/in/bordi.shp -al | head -100 ERROR 3: fread(65004) failed on DBF file. ERROR 3: fread(65004) failed on DBF file. [...] INFO: Open of `/home/nik/Desktop/in/bordi.shp' using driver `ESRI Shapefile' successful. Layer name: bordi Geometry: Line String Feature Count: 9 Extent: (1009114.847644, 4446909.136474) - (1312106.717266, 4668313.355348) Layer SRS WKT: PROJCS[WGS_1984_UTM_Zone_32N, GEOGCS[GCS_WGS_1984, DATUM[WGS_1984, SPHEROID[WGS_1984,6378137,298.257223563]], PRIMEM[Greenwich,0], UNIT[Degree,0.017453292519943295]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,9], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[Meter,1]] ��n...@q: String (17.0) L$D=�u: String (207.0) �=? �� �: String (191.0) ��GMI�FYYY: String (78.0) '6:�+JX�fW`: String (235.0) g ��ʚ *; : String (243.0) _*gD88Z��: String (109.0) 2u]「�...@qu�: String (119.0) �'@TN�]: String (223.0) �魗n�676t: String (117.0) +J W v���: String (136.0) \�� : String (212.0) o7ʤ��: String (247.0) �_\ ��: String (126.0) ˆøßøŒÿÛù¯L: S├⎼␋┼± (207.0) ;4OÓÓÐà 4´Þ: S├⎼␋┼± (32.0) B8¬┘┴šì¶%þ¿: S├⎼␋┼± (33.0) dž¸ ˜øæ%¬?│: S├⎼␋┼± (200.0) [...] # it continues like that :-( ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.patch, images vs coloured elevation data, aerial photo mosaic
On Mon, 2009-02-16 at 14:55 +, John Stevenson wrote: But I am 100% sure that it will work. Oops: But I am *NOT* 100% that it will work is what I wanted to write :-). But I am glad that it worked out. Kind regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] create raster from line vector
Damn! So it was simply a problem with the attribute table. In the meantime I've solved doing a moving average on the rasterized polygon, and extracting only the cat values between the inside and the outside... Not a straightfarword solution :) Thanks Nikos 2009/2/16 Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de: On Mon, 2009-02-16 at 16:00 +0100, Nikos Alexandris wrote: I share my vector, if someone wants to reproduce this problem: http://www.geospatial.it/allegri/bordi.tar.gz EPSG:32632 It works now. No problem. Cheers, Nikos --- history 1 v.in.ogr dsn=/home/nik/Desktop/in/bordi.shp out=testbordi 2 d.mon x0 3 g.region vect=testbordi 4 g.region vect=testbordi -pa 5 d.vect testbordi 6 v.info testbordi 7 v.to.rast testbordi out=testbordirast use=val value=1 8 d.rast -o testbordirast 9 r.info testbodirast 10 r.info testbordirast 11 history # r.info r.info testbordirast -hgr min=1 max=1 north=4668313.355348 south=4446909.136474 east=1312106.717266 west=1009114.847644 Data Source: Vector Map: testbordi in mapset PERMANENT Original scale from vector map: 1:1 Data Description: generated by v.to.rast Comments: v.to.rast input=testbordi output=testbordirast use=val type=p\ oint,line,area layer=1 value=1 rows=4096 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.mapcalc: columns()/modulo - rounding issue, possible bugfeature ?
On 16/02/09 15:31, peter.lo...@gmx.de wrote: Hi, I stumbled over this funny behavior of the col() function of r.mapcalc (GRASS6.2 and GRASS6.4, both on Linux). According to The Book, col() returns current col of moving window. I assume that for each row of a given raster layer col() should provide the values 1,2,3,4,5, etc. as it travels through the row. To produce a map of repeating value sequences (for example the numbers 1 2 3 4 ) the module operator can be used on the col()-function: col() operator result 1 mod 4 = 1 2 mod 4 = 2 3 mod 4 = 3 4 mod 4 = 0 5 mod 4 = 1 6 mod 4 = 2 etc this approach was put into mapcalc. The 1 + makes sure that the sequence runs from 1-4 not 0 - 3: r.mapcalc modulo=(1 + ( col() % 4 ) ) A map was produced in a Spearfish location (col() should operate independently from the projection). Here it comes: For some (strange ?)reason, there are ripples in the 1 - 4 sequence. Instead of going 1234 1234 1234 ad infinitum once in a while a number is omitted: 1234 1234 234 1234 134 etc. Has anybody come across this behavior before ? Is this a feature/good-thing ? After a little testing I come to the (preliminary) conclusion that this depends on whether the number of cols is divisable by 4, i.e. cols % 4 = 0. If this is not the case, then you don't get the complete sequence. So, it depends on your region setting. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] create raster from line vector
On Feb 16, 2009, at 7:40 AM, grass-user-requ...@lists.osgeo.org wrote: Date: Mon, 16 Feb 2009 15:17:27 +0100 From: G. Allegri gioha...@gmail.com Subject: [GRASS-user] create raster from line vector To: grassuser grass-user@lists.osgeo.org Message-ID: e12429640902160617g37ee0a88o85314bcfbf91c...@mail.gmail.com Content-Type: text/plain; charset=ISO-8859-1 Hello everyone. A simple question, but I'm getting mad with it! I simply have to convert a vector line to raster. I've done it many times in the past, but this time I can't make it work: v.to.rast input=border output=border use=val type=line layer=1 value=1 The output is simply blank (MAX and MIN values are NULL) What stupid error am I doing? Giovanni Did you set your computational region to match the vector using g.region? Regardless of what you see in the display, your computational region needs to be set so that GIS operations work on the map or parts of map you want. Michael ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problemas con winGRASS 6.4.0RC3 del instalador OSGeo4W
g.region says: EROR:Invalid ellipsoid '�' in file PROJ_INFO in PERMANENT I removed PROJ_INFO and now grass can be used in gui (gis.m) mode thanks for tip Jarek Markus Neteler pisze: On Sun, Feb 15, 2009 at 8:41 PM, Jarek Jasiewicz jar...@amu.edu.pl wrote: I confirm that error on old latlong location (created over 1 years ago, I didn't use it since that time) due to this error gis.m cannot start, but grass still working in the text mode Addational masage in text console are: Error in startup script: can't read monitor_zooms(1,1,n): no such variable while executing This is typically a GDAL problem. What does g.region -p say? Please test this since you can use GRASS in text mode. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problemas con winGRASS 6.4.0RC3 del instalador OSGeo4W
On Mon, Feb 16, 2009 at 6:05 PM, Jarek Jasiewicz jar...@amu.edu.pl wrote: g.region says: EROR:Invalid ellipsoid '� ' in file PROJ_INFO in PERMANENT How was the PROJ_INFO file generated? Which projection? I removed PROJ_INFO and now grass can be used in gui (gis.m) mode thanks for tip It's no good to remove the file since you lose the spatial reference. Better to actually correct it... Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problem with vdigit in 64 release
Hi, there seems to be a problem with cairo rendering: GUI D3/3: Command(): cmd='g.gisenv', wait=True, returncode=0, alive=False GUI D3/3: Command(): cmd='g.proj -p', wait=True, returncode=0, alive=False GUI D3/3: CreateTool(): tool=498, label=displaymap bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8961768 GUI D3/3: CreateTool(): tool=499, label=rendermap bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8961a28 GUI D3/3: CreateTool(): tool=500, label=erase bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8960c80 GUI D3/3: CreateTool(): tool=501, label=pointer bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x895f980 GUI D3/3: CreateTool(): tool=502, label=query bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x895fc58 GUI D3/3: CreateTool(): tool=503, label=pan bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8960e50 GUI D3/3: CreateTool(): tool=504, label=zoom_in bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8960e40 GUI D3/3: CreateTool(): tool=505, label=zoom_out bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8963038 GUI D3/3: CreateTool(): tool=506, label=zoom_back bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8963008 GUI D3/3: CreateTool(): tool=507, label=zoommenu bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8963090 GUI D3/3: CreateTool(): tool=508, label=analyze bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8960e60 GUI D3/3: CreateTool(): tool=509, label=overlay bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8960e70 GUI D3/3: CreateTool(): tool=510, label=savefile bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8961b60 GUI D3/3: CreateTool(): tool=511, label=printmap bitmap=wx._gdi.Bitmap; proxy of Swig Object of type 'wxBitmap *' at 0x8962dc0 GUI D2/3: MapFrame.__InitDisplay(): GUI D2/3: Map.ChangeMapSize(): width=760, height=545 GUI D3/3: Command(): cmd='g.region -u -g -p -c', wait=True, returncode=0, alive=False GUI D3/3: Map.GetRegion(): {'rows': 12073.0, 'e': 655681.695, 'cells': 101811609.0, 'cols': 8433.0, 'center_northing': 615211.348, 'n': 1218877.2, 'center_easting': 234044.549, 's': 11545.5, 'w': -187592.601, 'ewres': 99.99695244995, 'nsres': 100.00262569} GUI D3/3: BufferedWindow.OnSize(): GUI D2/3: Map.ChangeMapSize(): width=758, height=503 GUI D3/3: BufferedWindow.OnSize(): GUI D2/3: Map.ChangeMapSize(): width=758, height=503 GUI D3/3: BufferedWindow.OnSize(): GUI D2/3: Map.ChangeMapSize(): width=758, height=503 GUI D3/3: BufferedWindow.OnSize(): GUI D2/3: Map.ChangeMapSize(): width=758, height=503 GUI D2/3: Map.ChangeMapSize(): width=758, height=503 GUI D3/3: Map.AdjustRegion(): {'rows': 503.0, 'e': 655681.695, 'cells': 381274.0, 'cols': 758.0, 'center_northing': 615211.348, 'n': 1218877.2, 'center_easting': 234044.549, 's': 11545.5, 'w': -187592.601, 'ewres': 1112.4990765171503, 'nsres': 2400.2618290258447} GUI D3/3: Map.SetRegion(): n-s resol: 2400.261829; n-s resol3: 18864.5578125; rows: 503; north: 1218877.2; t-b resol: 1; zone: 0; bottom: -0.5; rows3: 64; west: -187592.6; top: 0.5; cols: 758; cols3: 128; depths: 1; e-w resol: 1112.499077; proj: 99; e-w resol3: 6588.08046875; east: 655681.7; south: 11545.5; GUI D3/3: Command(): cmd='g.pnmcomp in= mask= opacity= background=255:255:255 width=758 height=503
Re: [GRASS-user] Problemas con winGRASS 6.4.0RC3 del instalador OSGeo4W
Markus Neteler pisze: On Mon, Feb 16, 2009 at 6:05 PM, Jarek Jasiewicz jar...@amu.edu.pl wrote: g.region says: EROR:Invalid ellipsoid '� ' in file PROJ_INFO in PERMANENT How was the PROJ_INFO file generated? Which projection? Sorry I don' rememer it was over one year ago. projection lato long based on WGS84 elipsoid I removed PROJ_INFO and now grass can be used in gui (gis.m) mode thanks for tip It's no good to remove the file since you lose the spatial reference. Better to actually correct it... Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Export specific shapes from a vector file
I need to export specific shapes from a grass vector file to another grass vector file. Also, I need to convert specific polygons from a vector grass file to a raster grass file. How can I do that? Thanks, Leonidas -- View this message in context: http://n2.nabble.com/Export-specific-shapes-from-a-vector-file-tp2337165p2337165.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
Re: [GRASS-user] Export specific shapes from a vector file
On Mon, 2009-02-16 at 12:55 -0800, leonidas wrote: I need to export specific shapes from a grass vector file to another grass vector file. v.extract in=VectorMap out=ExtractedFeatures ... # you can use list= to extract specific categories (the cat column) # or an SQL where clause, e.g. something like where='SomeColumn=SomeString' Also, I need to convert specific polygons from a vector grass file to a raster grass file. I assume by polygons you mean areas. v.to.rast in=VectorMap out=RasterMap type=area ... # use=??? You need to decide which values you will assign to your pixels. How can I do that? Thanks, Leonidas Regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Export specific shapes from a vector file
Thank you Niko, well I understand the first one (using v.extract). -Is it possible to do something like that using v.out.ogr and export to shapefile? -As far as concerns the second question using v.to.rast, how can I export only specific areas using sql where clause or list for specific categories? Thanks again, Leonidas Nikos Alexandris wrote: On Mon, 2009-02-16 at 12:55 -0800, leonidas wrote: I need to export specific shapes from a grass vector file to another grass vector file. v.extract in=VectorMap out=ExtractedFeatures ... # you can use list= to extract specific categories (the cat column) # or an SQL where clause, e.g. something like where='SomeColumn=SomeString' Also, I need to convert specific polygons from a vector grass file to a raster grass file. I assume by polygons you mean areas. v.to.rast in=VectorMap out=RasterMap type=area ... # use=??? You need to decide which values you will assign to your pixels. How can I do that? Thanks, Leonidas Regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://n2.nabble.com/Export-specific-shapes-from-a-vector-file-tp2337165p2337293.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
Re: [GRASS-user] Export specific shapes from a vector file
On Mon, 2009-02-16 at 13:14 -0800, leonidas wrote: Thank you Niko, well I understand the first one (using v.extract). -Is it possible to do something like that using v.out.ogr and export to shapefile? AFAIK, no. The way to do it is: v.extract + v.out.ogr -As far as concerns the second question using v.to.rast, how can I export only specific areas using sql where clause or list for specific categories? Apologies, I overlooked the word specific before. Here (again) I think you would need v.extract + v.to.rast. An example: Let's say you have the CORINE map in whose attribute table you have a column named level3 of type integer (with a lentgh of 3) which contains the CORINE land cover class codes (111,112,113, etc.). And now you need to extract specific classes, that is areas with specific level3 codes (i.e. 111, 112, 122, 123, 141, 142). Then you could do: v.extract in=corine_ellas_v2 out=corine_ellas_v2_subset where='level3=111 or level3=112 or level3=112 or level3=122 or level3=123 or level3=141 or level3=142' --o then v.to.rast in=corine_ellas_v2_subset out=corine_ellas_v2_subset use=attr column=level3 Hope that helps Cheers, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Export specific shapes from a vector file
Yes that was what exactly I need.Thanks again Nikos Alexandris wrote: On Mon, 2009-02-16 at 13:14 -0800, leonidas wrote: Thank you Niko, well I understand the first one (using v.extract). -Is it possible to do something like that using v.out.ogr and export to shapefile? AFAIK, no. The way to do it is: v.extract + v.out.ogr -As far as concerns the second question using v.to.rast, how can I export only specific areas using sql where clause or list for specific categories? Apologies, I overlooked the word specific before. Here (again) I think you would need v.extract + v.to.rast. An example: Let's say you have the CORINE map in whose attribute table you have a column named level3 of type integer (with a lentgh of 3) which contains the CORINE land cover class codes (111,112,113, etc.). And now you need to extract specific classes, that is areas with specific level3 codes (i.e. 111, 112, 122, 123, 141, 142). Then you could do: v.extract in=corine_ellas_v2 out=corine_ellas_v2_subset where='level3=111 or level3=112 or level3=112 or level3=122 or level3=123 or level3=141 or level3=142' --o then v.to.rast in=corine_ellas_v2_subset out=corine_ellas_v2_subset use=attr column=level3 Hope that helps Cheers, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://n2.nabble.com/Export-specific-shapes-from-a-vector-file-tp2337165p2337658.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
Re: [GRASS-user] r.mapcalc: columns()/modulo - rounding issue, possible bugfeature ?
peter.lo...@gmx.de wrote: I stumbled over this funny behavior of the col() function of r.mapcalc (GRASS6.2 and GRASS6.4, both on Linux). According to The Book, col() returns current col of moving window. I assume that for each row of a given raster layer col() should provide the values 1,2,3,4,5, etc. as it travels through the row. To produce a map of repeating value sequences (for example the numbers 1 2 3 4 ) the module operator can be used on the col()-function: col() operator result 1 mod 4 = 1 2 mod 4 = 2 3 mod 4 = 3 4 mod 4 = 0 5 mod 4 = 1 6 mod 4 = 2 etc this approach was put into mapcalc. The 1 + makes sure that the sequence runs from 1-4 not 0 - 3: r.mapcalc modulo=(1 + ( col() % 4 ) ) A map was produced in a Spearfish location (col() should operate independently from the projection). Here it comes: For some (strange ?)reason, there are ripples in the 1 - 4 sequence. Instead of going 1234 1234 1234 ad infinitum once in a while a number is omitted: 1234 1234 234 1234 134 etc. Has anybody come across this behavior before ? Is this a feature/good-thing ? The phenomenon can be reproduced by using the above mentioned r.mapcalc command, displaying the result as a cell values layer. Please be more precise regarding the displaying part. If I run: r.mapcalc modulo = (1 + ( col() % 4 ) ) r.out.ascii modulo | head I see what I expect, not what you describe. OTOH, if you replace r.out.ascii with e.g. d.rast, the data will then be resampled again according to the screen dimensions. If you run d.rast from the command line, the default screen dimensions are 640 x 480. With the default region for Spearfish of 634 x 477, the resampling will be ever so slightly greater than 1:1, resulting in occasional duplicates. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Display output question
Hi all, What's the command-line equivalent to using the JPG export in the display manager? (I would like to script outputs from a number of workspace - GRC - files). Cheers, Richard ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: archaeologist GRASS users - was Thiessen Polygons
Wow, I didn't expect what I thought to be a simple question on how to do something to morph into two separate threads. :) I do have a question for archaeologist out there. Can they suggest some good articles on the use of Thiessen Polygons in archaeology, especially if they deal with prehistoric, monumental architecture? Thanks ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] vector area needs water fill
On Mon, 2009-02-16 at 20:01 -0500, Kurt Springs wrote: Hi folks, I do have another problem I would like some thoughts on. I have vector areas representing water areas (lakes and oceans). I purchased these areas from the Irish and Northern Irish Ordinance Survey. They came as .dxf packages. I extracted them but they would not fill. I opened them in v.digit and found that the vector lines were broken, even missing in places. I did my best to close the vectors. The ocean was just a coastline. I had to modify it by attaching a the bounding box and taking away everything that didn't look like the ocean. Getting these areas to fill even though I think I managed to join most of the nodes. Some I think I really did join, some I know would not join. Is there a simple(ish) way to generate areas that will be I can turn blue and say here are the lakes and oceans? I don't need them for analysis, but I would like to put maps generated into my dissertation and have them be pleasing to the eye. Any help on this will be appreciated. Kurt FWIW, once I received a .dxf file that contained burned areas digitised by GPS-supported devices. After getting the file in GRASS, it looked really bad. As you describe, open boundaries, missing parts of lines. I had to restore it manually. I don't think there is another option than hard-hand-working here. Except if you are lucky and you could use v.clean to close boundaries, add centroids with v.centroid and play around with colors. But I doubt it. Kind regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] vector area needs water fill
On Feb 16, 2009, at 8:37 PM, Nikos Alexandris wrote: On Mon, 2009-02-16 at 20:01 -0500, Kurt Springs wrote: Hi folks, I do have another problem I would like some thoughts on. I have vector areas representing water areas (lakes and oceans). I purchased these areas from the Irish and Northern Irish Ordinance Survey. They came as .dxf packages. I extracted them but they would not fill. I opened them in v.digit and found that the vector lines were broken, even missing in places. I did my best to close the vectors. The ocean was just a coastline. I had to modify it by attaching a the bounding box and taking away everything that didn't look like the ocean. Getting these areas to fill even though I think I managed to join most of the nodes. Some I think I really did join, some I know would not join. Is there a simple(ish) way to generate areas that will be I can turn blue and say here are the lakes and oceans? I don't need them for analysis, but I would like to put maps generated into my dissertation and have them be pleasing to the eye. Any help on this will be appreciated. Kurt FWIW, once I received a .dxf file that contained burned areas digitised by GPS-supported devices. After getting the file in GRASS, it looked really bad. As you describe, open boundaries, missing parts of lines. I had to restore it manually. I don't think there is another option than hard-hand-working here. Except if you are lucky and you could use v.clean to close boundaries, add centroids with v.centroid and play around with colors. But I doubt it. Kind regards, Nikos I am almost sure that I can trace the areas in qgis I made a test with its drawing tools. Is there a way to import that tracing back to GRASS? Kurt ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] how to get instantaneous heading of a line at a pt?
Hi, I have used v.to.points to create a series of points that fall along a line. Now I want to upload to a DB column in the points map the instantaneous heading of the original line at that sample point. any ideas? thanks, Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.generalize for area boundaries?
Hi, I am following the v.generalize tutorial at http://users.ox.ac.uk/~orie1848/tutorial.html (we should move that to the wiki before it disappears) but it doesn't say much about working with areas beyond removing small ones. I have a vector area which has a very steppy boundary, like from r.to.vect with a sawtooth pattern at the cell edges. I want to run a smoothing filter over it to get rid of the jaggy bits. No matter what method I try my output map is always the same as the input map, no vertices are created or destroyed. any ideas how to do this? I know about 'v.clean tool=prune' and Markus Metz's topology-preserving v.simplify (psst- add it to wiki addons) but I'd like to learn more about how to use v.generalize. thanks, Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Thiessen Polygons
Jan Hartmann wrote: With that in mind, if the algorithm you propose would be indeed an approximation to weighted Voronoi polygons, *and* it wouldn't be all to hard to implement (I have no idea about that), would it make sense to propose this as a new RFC for GRASS? Glynn wrote: I mean only that it cannot be done using the approach which r.grow.distance uses, using memory proportional to the number of columns (it uses a number of row buffers, i.e. one-dimensional arrays, with one element per column). [In any case, weighted distances won't produce polygons; at least, not for Euclidean distances. The boundary will only be a straight line if the weights are equal.] However, that doesn't meant that other algorithms wouldn't be feasible, particularly if you're only interested in typical behaviour rather than worst-case behaviour. Also, it may be possible to use the r.grow.distance approach with something other than scaling. An offset would work, optionally combined with a monotonic function of the distance (provided that it's the same for every point). Hi, just some brainstorming ideas for a weighted Voronoi module: input: vector points with weight column output: raster map with weighted Voronoi polygons (or r.to.vect built in) the module would create 2 raster maps: - one with each vector point expanded to a 2.5D hemisphere of influence bubbles centered over it, map value being the bubble height, similar to r.cost or when you lose in the old Missile Command arcade game. - the second map being a categorical raster containing the vector cat which has contributed the maximum bubble at each cell. e.g. canvas_map = all zeros; loop over vector_pts { bubble_height_at_cell = some_calc(); if ( bubble_height_at_cell canvas_map(row,col) ) { canvas_map(row,col) = bubble_height_at_cell; id_map(row,col) = current_vect_cat; } } when done the 2.5D map can be discarded. [low weight points have no area] id_map contains the vector point of note for that area. other ideas: - use r.param.scale to create a feature map and extract all saddle-point boundaries between the bubbles as the voronoi boundaries, (or r.slope.aspect and find areas where slope1 deg then r.thin, r.to.vect) - use some mountain peak prominence algorithm* on the 2.5D map starting at each input vector pt. [*] http://article.gmane.org/gmane.comp.gis.grass.user/19467/ - see v.surf.icw script in addons for other radial basis function ideas for the bubbles beyond the usual IDW 1/distance^2. well, ideas are somewhat abstract/vague, but perhaps something in it. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Seeking simple example of Grass+python through webserver.
Hi folks, I'm looking for a simple example of using grass from within Python through a webserver. PyWPS does what I'm looking for by setting up all the proper grass environment variables, the location, the mapset, a fake grass environment if you will, but it's not as easy to follow this code base as for example the great Grass+PHP examples at http://grass.osgeo.org/wiki/GRASS_and_PHP. Anyone know of some basic python examples like this? Thanks, Tim ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user