Re: [GRASS-user] v.generalize failing on some geoms
Hi Markus, Il 12/02/2018 21:46, Markus Metz ha scritto: > The original data in EPSG:3003 and EPSG:3857 are not identical. The > original data in EPSG:3003 are topologically correct, while the original > data in EPSG:3857 are not. interesting. so reprojecting creates errors (tiny gaps, in fact). just to document, the process was importing a shapefile into PostGIS, at the same time reprojecting from 3003 to 3857 > Further on, the original data in EPSG:3857 > contain features that are not present in the original data in EPSG:3003 > which is quite strange. forget about this, polygons added after import > Reprojecting the original data in EPSG:3003 to EPSG:3857 within GRASS > works fine, also with subsequent v.generalize. > > That means that the original data in EPSG:3857 are some reprojected > version of the original original data (which are these?). The > reprojection step, apparently performed on polygons, not GRASS areas > (how did you reproject?), introduced topological errors. Please use > native GRASS v.in.ogr + v.proj to reproject polygons. right; an easy alternative when working with QGIS Processing is to add a snap param, which cures the small gaps. All the best, and thanks again. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html https://www.google.com/trends/explore?date=all=IT=qgis,arcgis ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize failing on some geoms
Il 11/02/2018 18:34, Markus Metz ha scritto: > This implies that the data have been reprojected at some stage somewhere > from one SRS to the other. The most likely explanation for the different > results, particularly the topological errors reported by v.in.ogr, is > that non-topological polygons were reprojected. Can you make these geoms > in EPSG:3857 and EPSG:3003 available for testing? Hi Markus, thanks for your thoughts. Here the data, original and simplified, in two CRS (3003 is simplified correctly, 3857 not): http://www.faunalia.eu/~paolo/test_generalize.tar.xz All the best. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html https://www.google.com/trends/explore?date=all=IT=qgis,arcgis ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize failing on some geoms
Thanks Markus for the hint. Replies below. Il 10/02/2018 22:53, Markus Metz ha scritto: > This seems to be a reprojection problem. > > If you reprojected the vector data in GRASS with v.in.ogr + v.proj, this > is a problem of GRASS, granted that v.in.ogr did not complain about > vector topology. > > If you reprojected non-topological polygons, this is a common problem > when reprojecting non-topological polygons. Please import these polygons > first into GRASS, then reproject within GRASS and run v.generalize on > the reprojected GRASS vector. I'm not explicitly reprojecting it. I tried the command from within QGIS>Processing. Here the relevant commands: g.proj -c proj4="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs" v.in.ogr min_area=0.0001 snap=-1 input="/tmp/processing63c3c94ab569415fba39f41c12e3e48f" layer=1518196447.565 output=tmp1518196448046 --overwrite -o g.region n=5538851.39594 s=5150937.66442 e=1377191.37538 w=1004373.60434 res=100 v.generalize input="tmp1518196448046" method=douglas threshold="1000" look_ahead="7" reduction="50" -l output="output081005b4fafd473bbfa438a964de80f1" --overwrite All the best. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html https://www.google.com/trends/explore?date=all=IT=qgis,arcgis ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize failing on some geoms
Il 09/02/2018 19:25, Markus Neteler ha scritto: > On Fri, Feb 9, 2018 at 6:19 PM, Paolo Cavallini <cavall...@faunalia.it> wrote: >> Hi all, >> v.generalize is failing to simplify some of my geometries when in >> EPSG:3857. Same geoms in EPSG:3003 were simplified correctly. >> Incorrect borders were reported, unclear to me why. > > Please let us know which GRASS GIS version you use. 7.4.0-1 thanks -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html https://www.google.com/trends/explore?date=all=IT=qgis,arcgis ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.generalize failing on some geoms
Hi all, v.generalize is failing to simplify some of my geometries when in EPSG:3857. Same geoms in EPSG:3003 were simplified correctly. Incorrect borders were reported, unclear to me why. All the best, and thanks for any hint. === Avvio dell'algoritmo v.generalize.simplify - Vector based simplification generalization g.proj -c proj4="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs" v.in.ogr min_area=0.0001 snap=-1 input="/tmp/processing63c3c94ab569415fba39f41c12e3e48f" layer=1518196447.565 output=tmp1518196448046 --overwrite -o g.region n=5538851.39594 s=5150937.66442 e=1377191.37538 w=1004373.60434 res=100 v.generalize input="tmp1518196448046" method=douglas threshold="1000" look_ahead="7" reduction="50" -l output="output081005b4fafd473bbfa438a964de80f1" --overwrite v.out.ogr -s -e input=output081005b4fafd473bbfa438a964de80f1 type=auto output="/tmp/processing63c3c94ab569415fba39f41c12e3e48f/3d37e56a7f5f4e17a052af23fcace818" format=ESRI_Shapefile output_layer=output --overwrite Sto avviando GRASS GIS... Executing ... ATTENZIONE: Updating spatial reference with embedded proj4 definition ATTENZIONE: Datum non riconosciuto da GRASS e nessun parametro trovato Default region was updated to the new projection, but if you have multiple mapsets `g.region -d` should be run in each to update the region from the default Informazioni sulla proiezione aggiornate ATTENZIONE: Updating spatial reference with embedded proj4 definition ATTENZIONE: Datum non riconosciuto da GRASS e nessun parametro trovato La verifica della proiezione verrà ignorata Check if OGR layer contains polygons... 100 Creating attribute table for layer ... Importing 17 features (OGR layer )... 0..5..11..17..23..29..35..41..47..52..58..64..70..76..82..88..94..100 - Registrando primitive... 36 primitives registered 58960 vertices registered Numero di nodi: 36 Numero di primitive: 36 Numero di punti: 0 Numero di linee: 0 Numero di confini: 36 Numero di centroidi: 0 Numero di aree: - Numero di isole: - - Cleaning polygons - Si sta rompendo i poligoni... Breaking polygons (pass 1: select break points)... 2..5..8..11..13..16..19..22..25..27..30..33..36..38..41..44..47..50..52..55..58..61..63..66..69..72..75..77..80..83..86..88..91..94..97..100 Breaking polygons (pass 2: break at selected points)... 2..5..8..11..13..16..19..22..25..27..30..33..36..38..41..44..47..50..52..55..58..61..63..66..69..72..75..77..80..83..86..88..91..94..97..100 - Rimozione dei duplicati in corso... 1..3..5..7..9..11..13..15..17..19..21..23..25..27..29..31..33..35..37..39..41..43..45..47..50..52..54..56..58..60..62..64..66..68..70..72..74..76..78..80..82..84..86..88..90..92..94..96..98..100 - Breaking boundaries... 0..2..4..6..8..10..12..14..16..18..20..22..24..26..28..30..32..34..36..38..40..42..44..46..48..50..52..54..56..58..60..62..64..66..68..70..72..74..76..78..80..82..84..86..88..90..92..94..96..98..100 - Rimozione dei duplicati in corso... 1..3..5..7..9..11..13..15..17..19..21..23..25..27..29..31..33..35..37..39..41..43..45..47..49..51..53..55..57..59..61..63..65..67..69..71..73..75..77..79..81..83..85..87..89..91..93..95..97..99..100 - Cleaning boundaries at nodes... 1..3..5..7..9..11..13..15..17..19..21..23..25..27..29..31..33..35..37..39..41..43..45..47..49..51..53..55..57..59..61..63..65..67..69..71..73..75..77..79..81..83..85..87..89..91..93..95..97..99..100 - Breaking boundaries... 0..2..4..6..8..10..12..14..16..18..20..22..24..26..28..30..32..34..36..38..40..42..44..46..48..50..52..54..56..58..60..62..64..66..68..70..72..74..76..78..80..82..84..86..88..90..92..94..96..98..100 - Rimozione dei duplicati in corso... 1..3..5..7..9..11..13..15..17..19..21..23..25..27..29..31..33..35..37..39..41..43..45..47..49..51..53..55..57..59..61..63..65..67..69..71..73..75..77..79..81..83..85..87..89..91..93..95..97..99..100 - Cleaning boundaries at nodes... 1..3..5..7..9..11..13..15..17..19..21..23..25..27..29..31..33..35..37..39..41..43..45..47..49..51..53..55..57..59..61..63..65..67..69..71..73..75..77..79..81..83..85..87..89..91..93..95..97..99..100 - Merging boundaries... 2..5..8..11..14..17..20..23..26..29..32..35..38..41..44..47..50..53..56..59..62..65..68..71..74..77..80..83..86..89..92..95..98..100 - Removing dangles...
Re: [GRASS-user] Broken link in manpage?
Il 09/04/2016 18:32, Martin Landa ha scritto: > 2016-04-09 18:21 GMT+02:00 Paolo Cavallini <cavall...@faunalia.it>: >> The problem is even more serious for Windows users, that cannot >> generally compile themselves. > > I already noted that there daily builds of 7.0.4svn [1]. Martin > > [1] https://wingrass.fsv.cvut.cz/grass70/ yes, but this leaves out the huge majority of users, especially win ones. however, as said I respect your decision, even if I disagree, so I'm not going to insist on this issue. all the best. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Broken link in manpage?
Hi Moritz, Il 09/04/2016 15:38, Moritz Lennert ha scritto: > On 09/04/16 11:57, Paolo Cavallini wrote: > 1) Either release point releases more often, but then we have to be even > stricter than now concerning backports to the release branch. Otherwise > we will undermine the reliability even more. IMHO an important module totally broken, without alternatives, is worth an emergency release. > 2) Integrate selected bug fixes as patches in the Debian package. The problem is even more serious for Windows users, that cannot generally compile themselves. > 3) Provide a private Debian repository with daily packages for stable / > testing / unstable branches of Debian. -1: custom repos require a lot of effort to be properly maintained for all distros, and are often source of problems. > All three represent more workload and so the need for the necessary > humanpower. > > 1) is generally triggered when there is a really critical bug. > Apparently not very many see v.generalize as a function that is critical > enough to trigger a new release... I see the point; however, I respectfully disagree with that. All the best. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Broken link in manpage?
Il 09/04/2016 11:54, Martin Landa ha scritto: > sorry, but we have some release policy. There are daily builds at > least for Windows [1] available. Ma > > [1] https://wingrass.fsv.cvut.cz/grass70/ I respect the choices of the GRASS-PSC. However, this IMHO undermines the Notorious GRASS Reliability(TM). All the best. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Broken link in manpage?
Il 09/04/2016 11:14, Martin Landa ha scritto: > 2016-04-09 11:13 GMT+02:00 Paolo Cavallini <cavall...@faunalia.it>: >>> [1] https://trac.osgeo.org/grass/milestone/7.0.4 >> >> exactly, that's what I'm asking for. > > you don't need to ask for release, it's planned already:-) Ma well, I'm asking because I believe levaing an important command, unique of grass, totally broken for over 2 months is not good for the project. all the best. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Broken link in manpage?
Il 09/04/2016 11:12, Martin Landa ha scritto: > 2016-04-09 11:09 GMT+02:00 Paolo Cavallini <cavall...@faunalia.it>: >> packages for deb+osgeo4w > > before that 7.0.4 must be released :-) see [1]. Ma > > [1] https://trac.osgeo.org/grass/milestone/7.0.4 exactly, that's what I'm asking for. thanks. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Broken link in manpage?
Il 09/04/2016 11:08, Martin Landa ha scritto: > well, but it seems to be already fixed by Markus Metz. I am not sure > what you are asking for. Ma packages for deb+osgeo4w thanks -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Broken link in manpage?
Il 09/04/2016 11:01, Martin Landa ha scritto: > what do you mean by major platform and what bug exactly? Ma to me this means Debian/Ubuntu and osgeo4w - other might add osx. https://trac.osgeo.org/grass/ticket/2929 Thanks. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Broken link in manpage?
Il 09/04/2016 10:47, Martin Landa ha scritto: > thanks, fixed in r68225 (backported also to relbr70). Ma thanks. BTW, any hope to have fixed v.generalize available for major platforms anytime soon? All the best. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Broken link in manpage?
Hi all, in the v.generalize manpage there is apparently a broken link to the tutorial: http://users.ox.ac.uk/~orie1848/tutorial.html Should I open a ticket? Thanks. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize crash
Il 22/02/2016 21:33, Markus Neteler ha scritto: > Sure - here is our release schedule: > https://trac.osgeo.org/grass/milestone/7.0.4 > > In some weeks we'll have the first RC. Some more known bugs to be > eliminated since then! ok, thanks. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize crash
Il 21/02/2016 21:10, Markus Neteler ha scritto: > On Sat, Feb 20, 2016 at 8:50 AM, Paolo Cavallini <cavall...@faunalia.it> > wrote: > ... >>> I have opened a ticket: >>> https://trac.osgeo.org/grass/ticket/2929 > > It got fixed in 7.0.svn - please try (will go into 7.0.4). Hi Markus, that was fast, thanks. When is 7.0.4 scheduled? Given the crash, I believe having fixed packages soon is highly desirable. All the best. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize crash
Thanks Markus. All the best. Il 20 febbraio 2016 08:19:14 CET, Markus Neteler <nete...@osgeo.org> ha scritto: >Hi Paolo, > >On Mon, Feb 15, 2016 at 4:38 PM, Paolo Cavallini ><cavall...@faunalia.it> wrote: >> Hi all, >> I get consistent crashes of v.generalize on Debian sid, official deb, >e.g.: >> >> GRASS 7.0.3 (world):~ > v.generalize input=corine@PERMANENT layer=1 >> type=area type=line,boundary,area method=douglas threshold=5.0 >> look_ahead=7 reduction=50 slide=0.5 angle_thresh=3 degree_thresh=0 >> closeness_thresh=0 betweeness_thresh=0 alpha=1.0 beta=1.0 >iterations=1 >> output=corine_simpl > >Same issue here: > >GRASS 7.0.3svn (nc_spm_08_grass7):~ > v.generalize input=soils_general >layer=1 type=area type=line,boundary,area method=douglas threshold=5.0 >look_ahead=7 reduction=50 slide=0.5 angle_thresh=3 degree_thresh=0 >closeness_thresh=0 betweeness_thresh=0 alpha=1.0 beta=1.0 iterations=1 >output=soils_general_simpl >... >Number of areas: 1428 >Number of isles: 244 >- >Generalization (douglas)... >Using threshold: 5 meters >Segmentation fault (core dumped) > >with "gdb" >Program received signal SIGSEGV, Segmentation fault. >0x77bb2146 in Vect_line_intersection2 () >from >/home/neteler/software/grass70/dist.x86_64-pc-linux-gnu/lib/libgrass_vector.7.0.3svn.so > >I have opened a ticket: >https://trac.osgeo.org/grass/ticket/2929 > >Markus -- Sent from mobile. Sorry for being short http://faunalia.eu___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.generalize crash
Hi all, I get consistent crashes of v.generalize on Debian sid, official deb, e.g.: GRASS 7.0.3 (world):~ > v.generalize input=corine@PERMANENT layer=1 type=area type=line,boundary,area method=douglas threshold=5.0 look_ahead=7 reduction=50 slide=0.5 angle_thresh=3 degree_thresh=0 closeness_thresh=0 betweeness_thresh=0 alpha=1.0 beta=1.0 iterations=1 output=corine_simpl Same with different options and data. Does anyone confirm? Thanks. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.contour failing?
Il 19/10/2015 11:30, Markus Neteler ha scritto: > The problem is to "guess" how much memory resources are available > which change continuously (i.e. RAM is naturally static but its > allocation not). sure - perhaps taking the total memory as an upper bound could mitigate at least the most serious cases. > Just check the "history" file of the mapset. Some previously issued > command will have set the computational region. > A candidate might be a g.region vector=vmap with an additional > (unneeded) res=value call (maybe from "Processing"?). right - interestingly enough, the history does not seem to be visible in the current version of the grass qgis plugin - going to explore this. All the best. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.contour failing?
Hi all. I cannot get r.contour working. It always takes a lot of CPU, then crashes. Anyone has more joy? grass 7.0.1-2 from Debian sid. Thanks. -- Paolo Cavallini - www.faunalia.eu QGIS & PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] [Qgis-developer] QGIS and GRASS 7.0.0
Il 09/03/2015 00:38, Pedro Venâncio ha scritto: It looks really really great! Awesome, chapeau! Looking forward to see the crowdfunding. All the best. -- Paolo Cavallini - www.faunalia.eu QGIS PostGIS courses: http://www.faunalia.eu/training.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] failing to import table
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 02/09/2013 09:07, Moritz Lennert ha scritto: On 31/08/13 13:56, Markus Neteler wrote: Hi, does this persist with the current 6.4.3? I tested with post-6.4.3-release 64 release branch and it does persist. Thanks - should I open a ticket then? All the best. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIkOqIACgkQ/NedwLUzIr4UhACcCZy2KD9axwYj7uuGFY8IXs0W g3MAn17xPMu6O3/92x5Wx/o78BftnTxZ =9W6o -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] failing to import table
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 02/09/2013 11:18, Helmut Kudrnovsky ha scritto: Moritz Lennert wrote Have you tried accessing the attribute table ? Working with it through other modules ? What does v.db.connect -p give you ? attribute table manager opens without any problems, some further tests below Thanks for the tests. So it seems that this is sometimes working, sometimes not, right? What could be the reason? All the best. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIkjoMACgkQ/NedwLUzIr41kACdE3P/G58Ytb7qXIFtmnVad+9D EfUAn2ic0pVlYJxtTlVJYjZQtcSy/XxM =Gvxi -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] failing to import table
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 02/09/2013 15:30, Moritz Lennert ha scritto: I've just tried again with a freshly checked out and compiled grass64 release branch. I cannot reproduce the error. I'm not sure what might have caused it in my previous run, but at least in the current release branch (and thus 6.4.3 probably) it seems to be working. ok, thanks; I'll check again as soon as 6.4.3 hits Debian. All the best. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIkk2MACgkQ/NedwLUzIr4YVQCfca+JAkfRdGaUoKCvPxyQkeql W58AniC2aHf3EAA+gToy93VPI2kkRaIZ =Shjh -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] failing to import table
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 31/08/2013 19:07, Markus Neteler ha scritto: Due to time constraints at my end please test with the current stable version first. Unfortunately it is currently not possible for me right now. If there is anybody with a clean 6.4.3, I'll be happy to provide the file. Thanks a lot. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIjATgACgkQ/NedwLUzIr4P2wCfVjtWyX23EWzBVQsUksXz2xuM lLcAnRAvTjJgBOgjoc7q+OkwTH8cV10c =zibT -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Duplicated CAT during import
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 31/08/2013 19:22, Markus Neteler ha scritto: To my knowledge not possible when importing a SHAPE file since each row creates a new category value in the conversion from Simple Features to topology. Yes, that's also my interpretation. You may check the cat column with v.db.select. v.db.select map=test_import0@alessandra layer=1 columns=cat fs=| cat 1 2 Operazione conclusa con successo So, apparently it is just QGIS that displays different polygons with the same CAT as different rows in the table, but grass correctly see two polygons: correct? All the best, and thanks. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIjBHIACgkQ/NedwLUzIr6wEQCgr34d54ZlcYIilJ195Kcz5a/D h5YAoJq/YVAB9b4z26wCwCVv/sLnqUWW =Das3 -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.overlay not working?
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 31/08/2013 19:24, Markus Neteler ha scritto: Vector map liguria_agricola@alessandra is connected by: layer 1 table AGRICOLA in database cat /home/paolo/Desktop/spezia/grass/spezia/alessandra/dbf/ through driver dbf with key liguria_agricola As you see you have the spurious cat there. A proper connection string looks like this: GRASS 6.4.3svn (nc_spm_08):~/grass70 v.db.connect -p roadsmajor Vector map roadsmajor@PERMANENT is connected by: layer 1 table roadsmajor in database /home/neteler/grassdata/nc_spm_08/PERMANENT/dbf/ through driver dbf with key cat You may reconnect the table with v.db.connect. Done, it works. Unclear how the spurious cat crept in, however. Restarting from scratch: v.in.ogr dsn=/home/paolo/liguria agricola-uso suolo.shp output=liguria_agricola1 snap=-1 min_area=0.0001 -o v.db.connect -p liguria_agricola1 GRASS_INFO_MESSAGE(15222,1): Vector map liguria_agricola1@alessandra is connected by: GRASS_INFO_END(15222,1) layer 1 table agricola-uso in database liguria_agricola1 cat /home/paolo/Desktop/spezia/grass/spezia/alessandra/dbf/ through driver dbf with key suolo v.in.ogr dsn=/home/paolo/liguria_agricola_test.shp output=liguria_agricola2 snap=-1 min_area=0.0001 -o v.db.connect -p liguria_agricola2 GRASS_INFO_MESSAGE(15272,1): Vector map liguria_agricola2@alessandra is connected by: GRASS_INFO_END(15272,1) layer 1 table liguria_agricola2 in database /home/paolo/Desktop/spezia/grass/spezia/alessandra/dbf/ through driver dbf with key cat So, apparently spaces in file name cause the issue. All the best, and many thanks. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIjCMkACgkQ/NedwLUzIr4mlgCglmjb9cURYjMFUO5sd7pCSNqg V20An1kbUh7L2mh9DJ6hI26sm6ITOFHA =+TGL -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Duplicated CAT during import
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 31/08/2013 14:04, Markus Neteler ha scritto: I'm encountering problems with subsequent v.overlay, and I suspect duplicated CATs may be the cause. can you please post the relevant column info output of ogrinfo yourmap? which map? The original before import? Here it is: Layer name: test_import Geometry: Polygon Feature Count: 2 Extent: (1570773.587980, 4879203.434159) - (1579628.292941, 4890270.180298) Layer SRS WKT: PROJCS[Monte_Mario_Italy_zone_1, GEOGCS[GCS_Monte Mario, DATUM[Monte_Mario, SPHEROID[International_1924,6378388,297]], 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,150], PARAMETER[false_northing,0], UNIT[Meter,1]] COD_COMPON: String (6.0) COMPONENTE: String (200.0) COD_SUB_CO: String (10.0) SUB_COMPON: String (200.0) NUM_ELEMEN: String (4.0) NUM_SIST_E: Integer (10.0) NUM_TAV: Real (24.5) NORMA: String (200.0) SCHEDA: String (200.0) ID1: Integer (10.0) OGRFeature(test_import):0 COD_COMPON (String) = LAG COMPONENTE (String) = LIGURIA AGRICOLA COD_SUB_CO (String) = AGR SUB_COMPON (String) = AREE AGRICOLE NUM_ELEMEN (String) = (null) NUM_SIST_E (Integer) = (null) NUM_TAV (Real) = 51.0 NORMA (String) = \\rivu\dtuff\img\ptr\norme\LAG.pdf SCHEDA (String) = (null) ID1 (Integer) = 5 MULTIPOLYGON (... Thanks. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIiIqEACgkQ/NedwLUzIr54EgCgpBce4LzxuHWJhkx/TgwSyL2k VdIAniscwHf88szQsgLxIqtNPXoRFvEi =Zj+c -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Duplicated CAT during import
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 31/08/2013 19:09, Markus Neteler ha scritto: I don't see any existing cat column there. What should be duplicated? sorry I was unclear: I'm talking about duplicate records in the (grass generated) cat columns. All the best, and thanks. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIiJGYACgkQ/NedwLUzIr6BPwCfSTSDIyf4wdDQQzIPW55bkfns s2sAn1NH3bk+AH+v5/QFOCPf18YEkjF8 =21tc -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Duplicated CAT during import
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi all. Importing a shapefile with a simple v.in.ogr dsn=/home/paolo/test_import.shp output=test_import_min snap=-1 min_area=0.0001 -o results in a grass layer with duplicated CAT: is this a correct behaviour? I'm encountering problems with subsequent v.overlay, and I suspect duplicated CATs may be the cause. Thanks for any hint. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIgSWkACgkQ/NedwLUzIr5JIQCglDS55zlC+93V94ehDIXUlvIM IPUAn3LRG81o7V7DugbQy6qcNoTYkhVh =DVP1 -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] failing to import table
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi all. Importing a shapefile with space in the name results in an empty table of attributes. Is this a known problem? Should I file a ticket? GRASS 6.4.2-2 from Debian unstable. Thanks. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIgWTAACgkQ/NedwLUzIr5WjwCdE2LpHK77nWhC0v5+yC2frlcm ITkAn2/3HgE4xBWj0k3F4O+nRfJMzT3b =XaQr -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.overlay not working?
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi all. v.overlay fails on some specific layers: === v.overlay ainput=liguria_agricola@alessandra atype=area alayer=1 binput=uso_suolo_agricolo@alessandra btype=area blayer=1 output=over_xor operator=xor Collecting input attributes... DBMI-DBF driver error: Cannot create dbf database: cat /home/paolo/Desktop/spezia/grass/spezia/alessandra/dbf/ === I suspect this is due to the precence of duplicated CATs in one of the origin layers. So the questions are: * how to rearrange CATs, removing duplicates (but keeping all the geometries)? * does anybody confirm this behaviour? * I think it should fail more gracefully; should I file a ticket? Thanks a lot. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.14 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iEYEARECAAYFAlIgXfwACgkQ/NedwLUzIr6ZcQCfeod3EV5peZlyG9McGjwhp6R+ hK8AoKyZuH+4VAG3FGnFk4sgGpKc5yxb =tgnD -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Getting Started with GRASS Python Scripts
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 24/05/2013 21:13, Martin Lacayo ha scritto: Hello Paolo, Yes, my initial thought was to just use GRASS via SEXTANTE. I tired to do this, but did not have any luck. I posted previously on 14/5/2013 to qgis-users with the subject How to run SEXTANTE algorithms outside of QGIS Python console? but did not get a solution to the problem I encountered, which I thought was related to creating a QGIS instance so that I could make calls to SEXTANTE. If you have any suggestions that would be much appreciated. Hi Martin, better ask on qgis-dev or on IRC, most developers read it. All the best. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.12 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/ iEYEARECAAYFAlGjMj8ACgkQ/NedwLUzIr4HIwCfSk0xqZOnfaFjc24GaHvWF9gP FvYAn3xKhrBdvuUijTwLdgQIqXw3Ly7H =3q/V -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Getting Started with GRASS Python Scripts
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 23/05/2013 22:49, Martin Lacayo ha scritto: My ultimate goal is to decrease our development time for our open source ecosystem service tools: http://invest-natcap.googlecode.com, but more specifically at the moment I am trying to use r.viewshed so I do not have to write my own. Why not using GRASS through SEXTANTE for this? In this case you have nothing to do, just call whatever module you need. All the best. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.12 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/ iEYEARECAAYFAlGfGNUACgkQ/NedwLUzIr5QQQCcCyEDZmgSjTO2gnRC+Jk5jGRf 1esAmwSA0nx2cNbrRNNBkPtcotrdZvHt =NNfh -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Making images colours uniform
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 27/03/2013 11:32, Moritz Lennert ha scritto: I'm not sure if I completely understand what you are looking for, but does r.colors raster= do what you want ? Sorry I was unclear. I have a series of ortophotos, each of them with a different colour blend; when mosaicked, the result is ugly. My objective is to make them look more similar. Thanks. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.12 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/ iEYEARECAAYFAlFSy1MACgkQ/NedwLUzIr4ZIwCgj8R3/V0cejYeXu3lYxh+THnB XLUAnRMR80cB3Bj5LSSh52BDiBph3dll =2KOs -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Making images colours uniform
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi all. Is there a command in GRASS to merge adjacent images so that their colour tabel is more uniform? I remember something, but cannot find it back. All the best, and thanks. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.12 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/ iEYEARECAAYFAlFR0+4ACgkQ/NedwLUzIr7d7QCfR+an86eIIH2Ggm51TkgC2CRf 7GoAoIWviiyksi0lFYZqjZqspvf+js/4 =XWWh -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] [OSGeo-Discuss] GRASS community sprint in Genova, Feb 2nd to Feb 9th, 2013
Il 24/01/2013 21:09, Nikos Alexandris ha scritto: From my very user point of view, given the current participation plan, I 'll have the opportunity to: - meet freeriders Markus Neteler, Massimiliano Canata, Paolo Cavalini honoured by that - unfortunately I'll not be at the CS, arriving only on the 6th. hope to see you all then. enjoy your coding! -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] to move or not move ? grass64 - ? - 7.0
Il 10/01/2013 19:56, Paulo van Breugel ha scritto: Hi Milton I am using GRASS 7.0 most of the time, only switching to 6.4 when I want to use it from QGIS as I don't manage to get QGIS to work with GRASS 7.0 (that will need some work on the gdal-grass plugin I think). I do not think it is much work to update the plugin to use grass7; it would be a nice addition: any volunteer or supporter? All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] UTM zone overlay of world
Il 24/07/2012 10:00, RichardC ha scritto: For reference, I've created a Global UTM zones grid in shapefile format*. It's downloadable at http://www.enviroprojects.org/geospatial-services/gis-resources Thanks Richard for this. Would you mind releasing under a more free licence, for possible inclusion in sample data packages? The current licence: Free for non-commercial use only. Users should acknowledge the EnviroProjects Digital Initiative (EPDI: http://www.enviroprojects.org) as the source. would not allow this. All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Fwd: Re: [Qgis-developer] Sextante: v.buffer angle does not work
FYI. Thanks. Messaggio originale Oggetto: Re: [Qgis-developer] Sextante: v.buffer angle does not work Data: Wed, 30 May 2012 18:33:19 +0200 Mittente: Victor Olaya vola...@gmail.com A: cavall...@faunalia.it CC: qgis-developer qgis-develo...@lists.osgeo.org , I had splitted it, and in fact this algorithm just calls the angle option, but it still complains This is your command line v.buffer input=tmp1338388935286 angle=90 tolerance=0.01 -s -c output=output --overwrite and it says ERRORE: Select a buffer distance/minordistance/angle or column, but not both. But it's only passing an angle, not any of the other options!!!. Or you are not supposed to use tolerance with the angle option? Any GRASS expert can help? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague
Il 28/05/2012 21:17, Markus Metz ha scritto: not to java SEXTANTE. While on it (the blaming), I think it is time for Paolo to apologize for his rude [1] comments [2][3][4][5][6][7] earlier in the thread. I'm more than happy to apologize if I have hurt the sensibility of anybody. However, I think there is some misunderstanding, let's try to clarify: -SEXTANTE has two versions now (which are like 2 different projects) -The QGIS version is still unstable, but is a much better software (cleaner design, different things, etc) -It would be great to coordinate both versions, but it has to be done manually. currently, Victor Olaya is putting time on the QGIS version and not in the gvSIG one, so they might get uncoordinated, which Victor do not see as a terrible problem, anyway. -SEXTANTE for gvSIG has more algorithms and is more prepared for that, while SEXTANTE for QGIS has more backends, and is better suited for it. [2] http://osgeo-org.1560.n6.nabble.com/Report-from-ongoing-GRASS-GIS-Community-Sprint-in-Prague-tp4977136p4977385.html the link is wrong, the original source code is at http://code.google.com/p/sextante/source/browse/trunk/soft/bindings/qgis-plugin/src From the link I provided, you can easily find the source code at: http://plugins.qgis.org/plugins/sextante/ so it should be, according to Victor choice: http://sextante.googlecode.com/svn/trunk/ [3] http://osgeo-org.1560.n6.nabble.com/Report-from-ongoing-GRASS-GIS-Community-Sprint-in-Prague-tp4977136p4977411.html Paolo is appropriating python SEXTANTE as the QGIS Python plugin called sextante. See [2] I really do not understand this (see above). What is Python SEXTANTE? Anyway, sending a link and appropriating code seems different acts. [4] http://osgeo-org.1560.n6.nabble.com/Report-from-ongoing-GRASS-GIS-Community-Sprint-in-Prague-tp4977136p4977430.html Paolo is again appropriating python SEXTANTE (Please do nont send bug reports concerning it to other bugtrackers.) as the QGIS Python plugin called sextante. See [2] Bugs concerning the SEXTANTE QGIS Python plugin should be sent to: http://hub.qgis.org/projects/sextante/issues I do not see how this is appropriating anything. [5] http://osgeo-org.1560.n6.nabble.com/Report-from-ongoing-GRASS-GIS-Community-Sprint-in-Prague-tp4977136p4977473.html Paolo: I can correct you :) 1. Paolo was not addressed. 2. The issue has been resolved, no need for any further correction 3. I, not Paolo, gave an explanation and examples for the issue discussed (there is a complete re-write of SEXTANTE in Python). As I explained in a previous mail, Victor Olaya is the person who can best correct anybody on this. I took the liberty of talking on his behalf just because we are in close contact, and he has not subscribed to this list. [6] http://osgeo-org.1560.n6.nabble.com/Report-from-ongoing-GRASS-GIS-Community-Sprint-in-Prague-tp4977136p4977489.html Paolo: ... and join us. Paolo is not a developer of SEXTANTE, but very active for QGIS. Therefore us does obviously not mean SEXTANTE but QGIS. Paolo is again appropriating python SEXTANTE (Please do nont send bug reports concerning it to other bugtrackers.) as the QGIS Python plugin called sextante. See [2] Yes, us means QGIS. We (QGIS) already decided to integrate SEXTANTE into QGIS master as soon as code freeze is over (we are now releasing QGIS 1.8), so the distinction is blurred. [7] http://osgeo-org.1560.n6.nabble.com/Report-from-ongoing-GRASS-GIS-Community-Sprint-in-Prague-tp4977136p4977500.html At this stage of the python version of SEXTANTE, the main objective should be to get a stable version of python SEXTANTE. That is, a generic python SEXTANTE, not just something that works for QGIS and nothing else (with the SEXTANTE/GRASS interface not even working properly, see main comment). Paolo is again appropriating python SEXTANTE (Please do nont send bug reports concerning it to other bugtrackers.) as the QGIS Python plugin called sextante. See [2] Where do you find a generic python SEXTANTE? Having said that, if this still hurts, I apologize once more: http://www.youtube.com/watch?v=m7mIy97_rlo All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague
Il 26/05/2012 20:59, Benjamin Ducke ha scritto: Please report bugs relating to the SEXTANTE GRASS interface to the SEXTANTE bug tracker: http://bugs.gvsigce.org After logging in, switch the selection under Project: (upper right page area) to SEXTANTE: http://gvsigce.sourceforge.net/wiki/index.php/How_to_report_bugs The QGIS client code is just one of several GIS bindings that SEXTANTE supports, but any problems related to calling GRASS modules from SEXTANTE must be fixed in the shared Java code base. Wrong (no java here): http://hub.qgis.org/projects/sextante/issues That is the bugtracker for the gvSIG edition. I know it's confusing: http://hub.qgis.org/issues/5353 Thanks a lot. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague
Il 28/05/2012 14:33, Benjamin Ducke ha scritto: SEXTANTE and gvSIG CE share a bug tracker at http://bugs.gvsigce.org, a decision that was made together with Victor Olaya a while Hi all. Sorry for the misunderstanding. I am not sure Victor (main SEXTANTE dev, AFAIK) is listening here - if so, he can explain much better. Anyway, I (and I suspect also Marckus) was referring to the QGIS Python plugin called sextante. No shared code with the Java version, AFAIK. All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague
Il 28/05/2012 15:08, Markus Metz ha scritto: There is a dedicated bug tracker for QGIS SEXTANTE, and from a user perspective this is the logical place to report issues for QGIS SEXTANTE, in particular when keeping in mind that the GRASS command interfaces in QGIS SEXTANTE are hand-crafted, see v.buffer.angle.txt, ... How should a user take these differences into account? Report to http://bugs.gvsigce.org and explain that this does not apply to gvSIG SEXTANTE, only to QGIS SEXTANTE? The bugtracker for Sextante as a QGIS plugin is: http://hub.qgis.org/projects/sextante/issues Please do nont send bug reports concerning it to other bugtrackers. Thanks. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague
Il 28/05/2012 16:38, Benjamin Ducke ha scritto: does not contain just another binding for SEXTANTE to QGIS (even though the SVN folder structure suggests this) -- Victor, please correct me if I am wrong here. I can correct you :) Rather, it looks like it contains a complete re-write of SEXTANTE in Python (excluding the original SEXTANTE processing tools, which seem to be only available in their Java versions). Some of the tools are being ported. I am entirely unsure how to further maintain a code base in two different programming languages and how to keep my own work on the GRASS/SAGA/R interfaces in Java synchronized with this. Any ideas welcome (but should probably be discussed on the SEXTANTE user mailing list). As far as the GRASS code sprint is concerned, please consider that many people use SEXTANTE under gvSIG, OpenJUMP, GeoTools and other Java-based GIS at this point, and that they are just as interested in seeing the GRASS interface improved, as QGIS users are. For now, I will monitor the QGIS bug tracker and synchronize any relevant tickets manually with bugs.gvsigce.org. I do not see an alternative to manual sync (quite painful and inefficient, I know). All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague
Il 28/05/2012 17:18, Benjamin Ducke ha scritto: Sorry if I have been blunt in my reaction to this. I have put a lot of thought effort into the SEXTANTE-GRASS interface to get it to work as smoothly as it does now. I would hate to see this work die a slow death for lack of shared resources. I understand your point very well. I may be biased, but I think this new route is much more efficient, so I would suggest you to follow Victor and join us. All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague
Il 28/05/2012 17:56, José Antonio Canalejo Alonso ha scritto: we have put also a lot effort in SEXTANTE and in gvSIG/gvSIG CE. I like a lot these projects and also QGIS and I think it's positiv to see both (SEXTANTE-QGIS and SEXTANTE-gvSIG CE) growing and sharing Knowledge. The situation now is a bit confusing, but much better than 6 months ago: http://sextantegis.blogspot.de/2011/11/kinda-goodbye.html Ola José. I do not think the situation is very confusing for SEXTANTE: the gvSIG plugin is stable, the QGIS one is under active development. All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.to.averline -- script to find averaging line(s) of input vector map
Il 11/05/2012 13:03, Alexander Muriy ha scritto: I wrote shell script to find average line(s) of input vector map. Example picture in includes. quite nice, thanks. this solves a long-standing problem, for which I never found a satisfactory solution. very good for averaging multiple GPS tracks of the same road (I think OSMers will appreciate this). All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.to.averline -- script to find averaging line(s) of input vector map
Il 11/05/2012 13:05, Paolo Cavallini ha scritto: Il 11/05/2012 13:03, Alexander Muriy ha scritto: I wrote shell script to find average line(s) of input vector map. BTW: why not adding it to GRASS trunk? -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GeoPDF into GRASS?
Il 03/05/2012 10:37, Markus Neteler ha scritto: While I have no experience with that format, I see that GDAL supports it: http://www.gdal.org/frmt_pdf.html Stephen, it would be nice if you could provide us some sample files to test if the driver works well. Thanks. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] DTM with contour lines and points
Il 20/04/2012 00:51, Pedro Venâncio ha scritto: I thought a way to get around, which is to convert the contours to equally spaced points and join these points whith the elevation points layer. Is there a more correct and straightforward way to do? your real data are vertices anyway, so I think extracting them and patching them with points is the straightforward way. All the best. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: spgrass6
Il 06/03/2012 22:32, Dylan Beaudette ha scritto: Check out the plot() function from the raster package, or the related rasterViz (?) package. They have some neat stuff in there. Roger already mentioned this, but I'll add a little more. Now that we have Thanks for your comments. Please note however that I was talking about a very small (~1Mb) raster, for which the loading time was really surprising. All the best. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] rasters loaded with r.external miss a colour table?
Il 24/02/2012 09:24, Hamish ha scritto: yes, it is doing that, see main.c line ~ 340. if gdal supplies a color table it uses it, otherwise r.external checks if the map is type GDT_Byte and if so applies the grey255 rules. Or, if some other data type applies G_make_grey_scale_colors() based on the range. I would suggest that this be changed to not create any color table if none is provided. (and so act like other grass modules do in that case and fall back to rainbow rules) ok, thanks a lot for clarifying. all the best. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.clump misses clumps?
Hi all. Just ran r.clump, no options, and it obviously misses some areas (those at the bottom). It only finds 932 areas in my case, 1038 in other machines. Could someone check if it's a general phenomenon, before opening a ticket? Thanks. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.clump misses clumps?
Hi all. Just ran r.clump, no options, and it obviously misses some areas (those at the bottom). It only finds 932 areas in my case, 1038 in other machines. Could someone check if it's a general phenomenon, before opening a ticket? Thanks. Explained: the areas are there, it's just the color table that is missing for higher values. It seems to me a minor bug, right? All the best. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.clump misses clumps?
It seems to me a minor bug, right? All the best. Ticket created (#1598). Thanks. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to create polygon from two selected contour line
Il 25/02/2012 08:45, SWAPAN GHOSH ha scritto: Hello Grassuser, Anybody can please help me To create polygon from two selected contour line using grass comand.. I'm afraid I can't - I tried recently, no joy. Good luck, if you succeed please let us know. All the best. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] rasters loaded with r.external miss a colour table?
Hi all. Importing a dtm raster assign a correct colour table to it. If I load it with r.external, instead, only a greyscale is applied. This happens through the qgis plugin, do not know if it's the same with other interfaces). Is this a known issue? All the best. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] rasters loaded with r.external miss a colour table?
Il 23/02/2012 14:04, Markus Neteler ha scritto: On Thu, Feb 23, 2012 at 12:03 PM, Paolo Cavallinicavall...@faunalia.it wrote: Hi all. Importing a dtm raster assign a correct colour table to it. Which format do you have? For example, GeoTIFF has only limited color table support. adf (the same raster when imported has its nice colormap). Which operating system and which GRASS version? Debian testing, official package: 6.4.1-2 Thanks. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] rasters loaded with r.external miss a colour table?
Il 23/02/2012 14:57, Markus Neteler ha scritto: On Thu, Feb 23, 2012 at 12:03 PM, Paolo Cavallinicavall...@faunalia.it wrote: Hi all. Importing a dtm raster assign a correct colour table to it. If I load it with r.external, instead, only a greyscale is applied. Do you use r.colors? This happens through the qgis plugin, do not know if it's the same with other interfaces). Is this a known issue? I tried in GRASS: # quick and dirty test: GRASS 6.4.3svn (latlong_wgs84):~/data/maps/fao_forest_fra2000/frst_gez Sorry I was unclear: importing the map results in an automatically nicely coloured map, whereas r.external results in a greyscale. Apparently r.in.gdal applies automatically an appropriated colour map, r.external does not. All the best. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] spgrass6
Hi all. I'm working with spgrass6 R module, and I found that spplot() is horribly slow (several *minutes* for a 5k cell raster). Is this a common, known problem, or a local one? Thanks. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] spgrass6
Il 22/02/2012 10:38, Markus Neteler ha scritto: I am afraid that it is a generic R problem. Other example: While GRASS made a PCA in 1-2 min on a set of three radar maps I had to kill R after 20min because I got impatient... OK, I see, thanks. Are R devs aware of this? Should we open tickets? All the best. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] spgrass6
Il 22/02/2012 12:21, Benjamin Ducke ha scritto: As opposed to GRASS, R has not been designed with computational and/or memory efficiency as a priority. Oh! I was not fully aware of this. Maybe your analysis would allow you to run your computations on a representative sample instead of the whole dataset? In our case, the analysis is vert simple: just ssplot() of a small (2k by 2k cells) raster. Thanks. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: spgrass6
Il 22/02/2012 15:54, Roger Bivand ha scritto: Using the improved raster graphics handling for square cells with image() rather than spplot() and useRaster=TRUE - equivalent to image.SpatialGridDataFrame() and useRasterImage=TRUE with the same matrix takes 1.2 seconds on x11/cairo. You didn't say which version of R you are using - the raster graphics facilities have been improved recently. 2.14.1-1 on Debian testing. Did you try using image() instead of spplot() if your cells are square, and if rasterImage() is available in your version of R? Yes, image() is resonably fast, but misses automatic legend, etc. Thanks a lot for your reply. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] two questions: snap and close lines
Sample location here: http://int.faunalia.it/~paolo/test.tar.gz Confirmed by other users. Thanks. Il 27/01/2012 18:53, Markus Metz ha scritto: The only two options I can think of is giving v.edit a try, otherwise hand-editing. Urgh! Too bad, thanks anyway. Why does breaking not work? Maybe you need remove duplicates and clean small angles at nodes. You could try v.clean input=patched@paolo type=line tool=snap,break,rmdupl,rmsa thresh=1,0,0,0 output=patched_1 I'm not quite sure break do not work: new nodes are there, but in subsequent steps they are not treated as closed polygons. Thanks. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] error in example
Hi all. I'm following the examples in: http://casoilresource.lawr.ucdavis.edu/drupal/node/438 but I get an error in: new_data - SpatialGridDataFrame(grd, data=data.frame(k=rep(1,G$cols*G$rows)), proj4string=CRS(x@proj4string@projargs)) Error in validityMethod(object) : unequal number of objects in full grid and data slot Any hint? Thanks a lot. BTW: the link to sample data should be changed to http://mirrors.ibiblio.org/grass/sampledata/spearfish_grass60data-0.3.tar.gz -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.sunmask (position) and r.shaded.relief
Hi all. What is in short, the difference between the two commands? Thanks. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.sunmask (position) and r.shaded.relief
Il 07/02/2012 14:43, Markus Neteler ha scritto: Hope this helps, Markus Thanks a lot! -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize: strange behaviour on win?
Il 30/01/2012 16:35, Markus Metz ha scritto: My previous description was obviously incomplete. Output topology is always correct, independent of the version. In 6.4.1, topologically Thanks Markus, this clarifies the issue a lot. All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize: strange behaviour on win?
Il 28/01/2012 11:36, Martin Landa ha scritto: type=area ... type=boundary ...no good! AFAIK the parser reads this combination as `type=area,boundary` in this case. the results are the same when using ara+boundary or boundary alone. thanks. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] two questions: snap and close lines
Il 26/01/2012 20:41, Markus Metz ha scritto: On Mon, Jan 23, 2012 at 8:00 PM, Paolo Cavallini cavall...@faunalia.it wrote: Hi all. A couple of issues I ran into today: - I have a layer with some good and some bad polygons; how do I import it into GRASS letting the bad ones snapping onto the good ones, and not vice versa? You could try v.edit, otherwise it will be hand-editing with a digitizer. With topological vectors, snapping one set of polygons to another set of polygons does not make sense, because there are no polygons in topological vectors, only shared boundaries. That is probably the reason why this functionality is not available in grass. Hi Markus M, thanks for your reply. Let me explain better my use case: - I have a shapefile, with (almost) adjacent polygons - if I import them, I can snap the almost adjacent lines, and get common boundaries - the problem is that I cannot control which of the two adjacent lines stay fixed, and which one moves to it - I know that of the two polygons, one is good (certified), the other is not; so I want the first to remain fixed, and have the other snap to it. Did I explain myself better? For contour lines, you would need to patch the vector with contour lines with a vector created with v.in.region type=line, then break lines at intersections, convert lines to boundaries, add centroids. This assumes that not closed contour lines have endpoints exactly at the extents of the vector, e.g. if the western edge of the contour vector is 30 and there is an open contour line ending at say 29., this will not work. If the contour lines were produced with r.contour, this should work. In practice: #set the region around contour v.in.region output=region_line type=line v.patch input=contour,region_line output=patched v.clean input=patched@paolo type=line tool=snap thresh=1 output=patched_1 v.clean input=patched@paolo type=line tool=break output=patched_break v.type input=patched_break@paolo output=patched_break_boundaries type=line,boundary v.centroids input=patched_break_boundaries output=polygons It mostly works, but apparently breaking does not (only closed rings become different polygons). In the meantime, I spotted a few minor issues in the qgis-grass interface, I'm going to fix them. Thanks a lot. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize: strange behaviour on win?
Il 26/01/2012 20:53, Markus Metz ha scritto: You would need to provide the input, not the output as an example. https://int.faunalia.it/~paolo/province.zip Considering the detail of the output and the threshold used, it is no surprise that not all boundaries were generalized because that would break topology. If you want to generalize with a threshold of 10, I would recommend to first remove all areas with a size up to 10. This is not the point. The very same input and command work in Linux, leave ungeneralized boundaries in Windows. If you consider that also r.to.vect for lines also behaves differently between the two OSs, I suspect an OS-specific bug here. Any confirmation? Thanks a lot. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize: strange behaviour on win?
Il 27/01/2012 19:02, Markus Metz ha scritto: same output. That is, if Province_general3 has been generated on windows, I do get the same result on Linux, no difference between the OS's. Are you using an older grass version on Linux? aptitude show grass Versione: 6.4.1-1 Thanks. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] two questions: snap and close lines
Il 27/01/2012 18:53, Markus Metz ha scritto: The only two options I can think of is giving v.edit a try, otherwise hand-editing. Urgh! Too bad, thanks anyway. Why does breaking not work? Maybe you need remove duplicates and clean small angles at nodes. You could try v.clean input=patched@paolo type=line tool=snap,break,rmdupl,rmsa thresh=1,0,0,0 output=patched_1 I'm not quite sure break do not work: new nodes are there, but in subsequent steps they are not treated as closed polygons. Thanks. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.generalize: strange behaviour on win?
Hi all. Commands like: v.generalize input=province@Mapset_stefano type=area layer=1 -c type=boundary method=douglas threshold=10 look_ahead=7 reduction=50 slide=0.5 angle_thresh=3 degree_thresh=0 closeness_thresh=0 betweeness_thresh=0 alpha=1.0 beta=1.0 iterations=1 output=Province_general3 on Linux go well, whereas on windows some of the features are not simplified: see https://int.faunalia.it/~paolo/generalize.zip for an example. Anyone confirms? Should I open a bug? All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.water.outlet
Il 04/01/2012 20:55, Micha Silver ha scritto: I guess you mean r.water.outlet? Yes, sorry. 1- The input must be the drainage *direction* map from r.watershed. ok 2- And it's crucial that the easting and northing parameters for the outlet point lay exactly on a stream. Zoom in very close to the streams raster at the outlet point and check for the exact coordinates. Checked, it seems ok. The output maps is all 0s, with only one cell (close to the point I chose, but not on it) with value 1. Thanks. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.water.outlet
Hi all. I tried this command, but I always get a map full of 0. I start from a drainage map from r.watershed, as described. Anyone confirms, or am I doing something wrong? All the best. -- Paolo Cavallini See: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Can't edit GRASS in QGis anymore
Date: Tue, 1 Nov 2011 16:38:16 -0700 From: Roger Andr? ran...@gmail.com Subject: [GRASS-user] Can't edit GRASS in QGis anymore To: GRASS user list grass-user@lists.osgeo.org Message-ID: CAOL+ROGBJeD5AaJk=mwE=qpeoykcidcwnm9rbm1_v8rspb9...@mail.gmail.com Content-Type: text/plain; charset=utf-8 I realize this might be more of a QGis question than a GRASS one, but perhaps the same problem has affected someone here. In QGis when I click on the Edit GRASS Vector Layer button, I no longer can see the GRASS layer and do any editing operations to it. The Symbology pane of the GRASS Edit window shows all of the data types as being un-checked, and I cannot turn them back on. What makes this doubly frustrating is that at one point this worked fine, I fiddled with the settings in the Symbology window, and afterwards I had the problem I'm describing now. I've tried removing and reinstalling QGis and GRASS and I cannot seem to affect this behavior. Any clues greatly appreciated. Reinstalling does not help, removing your preferences (under ~/.qgis) should do. All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.drain
Hi all. I'm testing r.drain, but I cannot get reasonable results over a dtm: the line produced is always much too short, stopping where it seems it should proceed. I exaggerated elevation (*10), as suggested by the manpage, but results were only slightly better. Any suggestion? Thanks. -- http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.drain
Il giorno gio, 24/03/2011 alle 12.47 +0100, Soeren Gebbert ha scritto: Hi Paolo, try to use r.fill.dir on the elevation map with ~3 - 4 runs before you use r.drain. Thanks Soeren, that improved the results considerably. It would be good to have an option to set the level of filling (by recursion, if nothing better is available). All the best. -- http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GRASS plugin
Il giorno mar, 22/03/2011 alle 08.42 +0100, Soeren Gebbert ha scritto: Indeed it is currently not user friendly to install and configure a WPS GRASS GIS environment and GRASS does not have a WPS server included. But GRASS GIS 7 supports WPS server by providing the capability of creating WPS XML process descriptions for every module. And there is the wps-grass-bridge which implements direct GRASS GIS 7 back-end support for PyWPS and ZOO WPS. The 52N WPS server uses the wps-grass-bridge to implement GRASS GIS support too. My hope is that i.e. 52N will provide a graphical installer for there WPS GRASS GIS 7 and Sextante solution to make the deployment and use of a WPS server quite simple. Hi all. Thanks for the suggestion. However, IMHO requiring the installation of a WPS server, plus the configuration of all modules, is a rather heavy requirement, and it will make the installation of the whole beast much more cumbersome. Furthermore, 52N is in java, which will make things only worse. In addition, this does not seem a short term solution, and I think we need something working smoothly *now*. All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: [Qgis-user] GRASS plugin
Il giorno lun, 21/03/2011 alle 14.49 +0100, jr.morre...@enoreth.net ha scritto: I'm willing to help maintain this plugin but as I'm not a coder it will be limited to keeping the modules' UI up to date. Thanks Jean-Roc. Anyone willing to help for core issues? All the best. -- http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] current state of r.li
Hi all. Sorry for crossposting. As some of you know, the r.li suite of GRASS commands allows landscape analyses[0]. Its interface is rather complex, and is still in TclTk, not ported to either wxpython or qgis. As such, it is now more difficult to use than it should be, and it will become unusable when TclTk support will be dropped. The possible solution (thanks Radim) is to rewrite the interface as a qgis python plugin. It should not be a huge work (we provisionally estimate 2-3 weeks). The question is: is there anybody willing to invest either his/her time, or some money, to write such a plugin? We (Faunalia) would be happy to help if necessary. All the best. -- http://www.faunalia.it/pc [0]http://grass.fbk.eu/gdp/html_grass64/r.li.setup.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: QGIS GRASS plugin shadedrelief issue
Date: Thu, 23 Sep 2010 11:06:21 +0300 From: Maris Nartiss maris@gmail.com Subject: Re: [GRASS-user] Re: QGIS GRASS plugin shadedrelief issue To: nunosousa84 nuno_sousa_1...@hotmail.com Sorry for You. You do made backup of important data before messing around with QGIS, right? You can try to launch GRASS and rebuild vector topology with v.build. If it doesn't help - You'r SOL. And don't forget to infrorm QGIS developers how well their GRASS support works. Good luck, Maris. Hi all. Sorry I was unclear: it is a known bug for QGIS to crash in case of an unclean vector. There is no indication that QGIS is more inclined to spoil your vector than any other GRASS GUI. If you stop GRASS while it is writing a vector, problems are likely to occur, with any UI (even from the shell). You can clean your vector also from within QGIS, by simply running v.build.all. And yes, backupping data is always a good idea, regardless of what sw you are using. BTW: Maris, why are you trolling around? We do our best, like everybody, to improve the tools we are using, for the benefit of everybody, and it's no fun to read your polemic notes. All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] sunmask
Thanks. Fixed in https://trac.osgeo.org/qgis/changeset/13427 All the best. On Wed, 5 May 2010 23:52:43 +0200, Markus Neteler nete...@osgeo.org wrote: Yes: the second parameter is missing. Markus -- http://faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] sunmask
Hi all. I'm trying r.sunmask, but I get errors: === r.sunmask elev=t_dtm...@permanent year=2010 month=05 day=05 hour=10 minute=07 timezone=+1 -z output=shadow Using map center coordinates: 1663547.265000 4801864.162000 Neither sun position nor east/north, date/time/timezone definition are complete. === Is the QGIS module incomplete, or? All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize crash
Markus Metz ha scritto: Try another smoothing algorithm, e.g. chaiken. Thanks. Yes, Chaiken and Hermite work. However, IMHO it should not crash, just stop. All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.generalize crash
Markus Metz ha scritto: Right. I'll fix it, it should abort with fatal error: Out of memory. Thanks a lot for this. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.generalize crash
The snake algorithm crashes: ~$ v.generalize input=cont...@paolo type=line layer=1 -c type=line method=snakes alpha=1.0 beta=1.0 output=cont_snake GRASS_INFO_MESSAGE(22465,1): Generalization (snakes)... GRASS_INFO_END(22465,1) GRASS_INFO_PERCENT: 1 GRASS_INFO_PERCENT: 3 GRASS_INFO_PERCENT: 5 GRASS_INFO_PERCENT: 7 Killed Any explanation? All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [Fwd: [GRASS-user] r.cost vs. costallocation (arcgis)]
António M. Rodrigues ha scritto: Allocation surface / grid: a raster whose value of each cell is the code of the nearest target. Example: If the possible targets are hospital 1 and 2, each cell in the allocation surface will contain values 1 or 2, depending on which is the nearest hospital. I was surprised r.cost didn't provide this information (allocation surface), since this is quite standard output when using this methodology with arcinfo (using a function called costdistance); but yes, I'm quite sure there must be a way of doing this in GRASS. If I understand correctly, you just need to calculate several cost maps (you can do that in a loop), and use mapcalculator for this (function MIN). Would this do? All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [Fwd: [GRASS-user] r.cost vs. costallocation (arcgis)]
From: António M. Rodrigues amrodrig...@fcsh.unl.pt To: grass-user@lists.osgeo.org Subject: [GRASS-user] r.cost vs. costallocation (arcgis) Date: Sun, 14 Mar 2010 13:21:14 + Cost distance and costallocation (arcinfo and arcgis) output two surfaces: an accessibility and an allocation surface. For verious reasons, I need both. In GRASS, I had a look at r.cost, r.walk r.watershed and r.drain, trying to find a toll which could create an allocation grid based on a cost surface. From what I could see, none of these tools will do. I find it hard to believe that: 1. No one have had this problem before using GRASS, 2. there is no solution using GRASS. If someone, anyone, could point me in any direction, it would be greatly appreciated. Even something like: dude, just give up, and go back to arcgis... If the solution involves programming a new tool, probably based on r.cost, is there some tutorials online dealing with programming in GRASS? Hi António. Could you please explain us more clearly what do you (rather, ESRI) mean by accessibility and allocation? I would be very surprised if this could not be done with GRASS. All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.watershed fails
Hamish ha scritto: the error message indicates that qgis is trying to load a map named 1000. OK, so it must be something wrong with the View output code in the qgis-grass-plugin. Ticket opened: https://trac.osgeo.org/qgis/ticket/2335 Thanks for help understanding. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] meaning of g.region zoom?
Hamish ha scritto: fyi 'g.region zoom=' starts at the current region settings*, so maybe try: g.region rast=foo g.region zoom=foo * this lets you zoom in on individual clusters within the map There must be something wrong here: g.region zoom=los_175_...@paolo now results in: 10647 paolo 40 0 221m 27m 8096 R 98 2.7 0:50.37 g.region forever - more and more puzzled. This both from qgis and from grass. Anyone can confirm this behaviour? Thanks. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.los crashing
Hi all. While running r.los at increasing values for obs_elev and max_dist, the module starts crshing, e.g.: r.los input=...@permanent coordinate=1676844,4830417 obs_elev=50 max_dist=5000 output=los4 Using maximum distance from the viewing point (meters): 5000.00 Module crashed or killed This both from QGIS plugin and from the console. g.region -p projection: 99 (Transverse Mercator) zone: 0 datum: rome40 ellipsoid: international north: 4834913.78536118 south: 4803364.45415656 west: 1650370.47056664 east: 1689398.62729587 nsres: 24.99947005 ewres: 25.00202225 rows: 1262 cols: 1561 cells: 1969982 Lowering resolution solves the problem, so I guess it's a problem with memory. GRASS 6.4.0~rc5+40109-1 from Debian unstable. Seems the same as http://trac.osgeo.org/grass/ticket/111 , right? Any fix within sight? All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.watershed fails
Hi all. While running: r.watershed elevation=...@permanent threshold=100 basin=watershed2 and loading the output through QGIS, I get repeated errors: /home/paolo/Desktop/Dati/grass/Toscana_corso/paolo/cellhd/1000 is not a valid or recognized raster data source Is this a local problem, or a general phenomenon? Any possible explanation? All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] meaning of g.region zoom?
Hi all. If I zoom to a small map full of nulls, with only a non-nulla area on one side, and zoom back to the full extent, I do not see the difference: === g.region zoom=los_175_...@paolo g.region -p projection: 99 (Transverse Mercator) zone: 0 datum: rome40 ellipsoid: international north: 4834913.78536118 south: 4821192.33020234 west: 1668983.89916058 east: 1685595.88381969 nsres: 100.156607 ewres: 100.07219674 rows: 137 cols: 166 cells: 22742 g.region zoom=...@permanent g.region -p projection: 99 (Transverse Mercator) zone: 0 datum: rome40 ellipsoid: international north: 4834913.78536118 south: 4821192.33020234 west: 1668983.89916058 east: 1685595.88381969 nsres: 100.156607 ewres: 100.07219674 rows: 137 cols: 166 cells: 22742 === Where am I wrong? All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] borken links in web site
Hi all. In the web page: http://grass.osgeo.org/gdp/tutorials.php there are several broken links, e.g. http://home.wlu.edu/~dharbor/gisrs/index.html http://www.geo.unipr.it/~gis/ All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] does the grass support windows 7
Markus Neteler ha scritto: On Tue, Nov 24, 2009 at 8:10 AM, Paolo Cavallini cavall...@faunalia.it wrote: maven apache ha scritto: does the grass support windows 7 system? We are finding that 7 is more or less the same as vista (marketing... ;) ): GRASS used from QGIS plugin crashes at about the same commands as in vista. Did you compile it on Windows 7? are you joking ;) ? No windows here - it's just the results of several courses we are doing, always using the OSGeo4W installer. maybe compiling for vista or 7 solves the problems, no idea - however, the fact that the crashes are often the same in 7 and vista makes me think they are in fact the same system. All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] does the grass support windows 7
Markus Neteler ha scritto: We have installed OSGeo4W-QGIS-GRASS here on XP today. It delivered a version from May 2009 which is about 2000 SVN revisions ago. Is there hope to get a fresh recompile from the SVN 6.4 release branch? Or did we miss anything? AFAIK, you're right - any volunteer to help with OSGeo4w compilation? All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] does the grass support windows 7
maven apache ha scritto: does the grass support windows 7 system? We are finding that 7 is more or less the same as vista (marketing... ;) ): GRASS used from QGIS plugin crashes at about the same commands as in vista. All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user