[GRASS-user] i.atcorr whit Aster Data - null value - WARNING: Unknown atmospheric model!....
Hi, I have a problem with i.atcorr with data aster. I tried it with grass 64 and 7, on ubuntu 15.04 64-bit. Maybe I'm wrong to set the parameters of the txt? (Is there anyone can tell me if they are correct? I have attached the file .met)...Otherwise I do not know. /10 - geometrical conditions=ASTER 01 16 15.28 -74.9 -14.8 - month day hh.ddd longitude latitude (hh.ddd is in GMT decimal hours) 2- atmospheric mode=midlatitude summer 1- aerosols model=continental 15 - visibility [km] (aerosol model concentration) -0.121 - mean target elevation above sea level [km] -1000- sensor on board a satellite 72 - aster band 1 (0.480-0.645) / The result I get is always with null values. / | Total Cells: 27293388 | |Projection: UTM (zone 18) | |N: -1605247.26436317S: -1678867.26436317 Res:15 | |E: 552980.14168175W: 469565.14168175 Res:15 | | Range of data:min = -nan max = -nan / In the output window: / wavelength less than 0.25 micron: let's take s(l)=s(0.25) WARNING: Unknown atmospheric model! WARNING: Unknown aerosol model!/ Only in some cases there is a warning wavelength. Here there is the file .met AST_L1B_00301162001152851_20150713051848_23493.met http://osgeo-org.1560.x6.nabble.com/file/n5215412/AST_L1B_00301162001152851_20150713051848_23493.met I've seen other threads such as Yann but I did not understand the problem http://osgeo-org.1560.x6.nabble.com/i-atcorr-with-Aster-td3997467.html Thank you -- View this message in context: http://osgeo-org.1560.x6.nabble.com/i-atcorr-whit-Aster-Data-null-value-WARNING-Unknown-atmospheric-model-tp5215412.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Grib file: strange error
Hi everyone. Using the software zygrib (Also seen here http://grass.osgeo.org/wiki/GRIB) I download the file .grb of the area of interest. When importing the file in GRASS I get many raster ... ok but each raster is stretched along the y What is the problem? http://osgeo-org.1560.n6.nabble.com/file/n4985647/41.png Thanks Gabriel PS: I attach the file http://osgeo-org.1560.n6.nabble.com/file/n4985647/20120703_120237_.grb 20120703_120237_.grb -- View this message in context: http://osgeo-org.1560.n6.nabble.com/Grib-file-strange-error-tp4985647.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] Grib file: strange error
peifer wrote On 03/07/2012 14:22, Gabriele N. wrote: Hi everyone. Using the software zygrib (Also seen here http://grass.osgeo.org/wiki/GRIB) I download the file .grb of the area of interest. When importing the file in GRASS I get many raster ... ok but each raster is stretched along the y It looks like x and y size of the pixels are simply different. See below. Hermann gdalinfo reports: Pixel Size = (0.500,-33.2680001) grib_dump says: jDirectionIncrementInDegrees = 33.268; iDirectionIncrementInDegrees = 0.5; Sure, but if you look the image of the same file of zyGrib all seems ok. When I import the data in GIS there is the problem along the y axis. http://osgeo-org.1560.n6.nabble.com/file/n4985721/20120703_120237_zyGrib.jpg http://osgeo-org.1560.n6.nabble.com/file/n4985721/20120703_120237_zyGrib.jpg 20120703_120237_zyGrib.jpg Gabriele -- View this message in context: http://osgeo-org.1560.n6.nabble.com/Grib-file-strange-error-tp4985647p4985721.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: GrassAddOn - i.pr problem
Ciao. I.pr modules seem to work but I found a problem: i.pr.features I saw this post because I have the same problem i.pr.features -s training=TL_ipr_training_TL_ASTER_2006_02_01_RGB.52 features=aasasssas ERROR: read_training- Can't open file TL_ipr_training_TL_ASTER_2006_02_01_RGB.52 for reading Have you found a solution? Thank you very much Gabriele -- View this message in context: http://osgeo-org.1560.n6.nabble.com/GrassAddOn-i-pr-problem-tp3921153p4359359.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: i.pr -- g.extension and compiling error on ubuntu 11.10 64 bit
hamish-2 wrote Gabriele wrote: I installed but I have the same error. With g.extension I can install other add-ons (eg i.landsat.acca) but not i.pr but are they actually on the disk somewhere, slightly in the wrong place? i.e. if you look in $GRASS_ADDON_PATH/ is there a bin/ dir in there that has them? if so, just move them into the $GRASS_ADDON_PATH dir manually. (the automatic move fails as the module names do not match the directory name, same trouble with gipe, mcda, ...) Ciao Hamish, I do not understand what you said:) but I did the following: I copied the folder i.pr in /usr/local/grass-6.5.svn/bin/ (where there is also i.landsat.acca that I installed with g.extension). Then cd /usr/local/grass-6.5.svn/bin/i.pr/ and make MODULE_TOPDIR=/usr/local/grass-6.5.svn Now it seems to work ... I start to try it. thank you very much Gabriele -- View this message in context: http://osgeo-org.1560.n6.nabble.com/i-pr-g-extension-and-compiling-error-on-ubuntu-11-10-64-bit-tp4352067p4356276.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Re: i.pr -- g.extension and compiling error on ubuntu 11.10 64 bit
margherita wrote I can confirm that installation of i.pr by g.extension seems not to work properly. So, I suggest you to install it by hand following these steps: cd your path/grass6_devel/imagery svn co https://svn.osgeo.org/grass/grass-addons/grass6/imagery/i.pr/ cd i.pr/ make cd .. cd .. make libs make install This worked for me. HTH, madi Ciao Madi, also seems to function well. Now I am beginning to test it. Thanks:) Gabriele -- View this message in context: http://osgeo-org.1560.n6.nabble.com/i-pr-g-extension-and-compiling-error-on-ubuntu-11-10-64-bit-tp4352067p4356496.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: i.pr -- g.extension and compiling error on ubuntu 11.10 64 bit
Thanks hamish. I installed tcl8.5-dev and I followed this guide http://grass.osgeo.org/wiki/Compile_and_Install_Ubuntu and, although I got error ((Errors in: /usr/local/src/grass6_devel/lib/ogsf /usr/local/src/grass6_devel/lib/nviz /usr/local/src/grass6_devel/misc/m.nviz.image /usr/local/src/grass6_devel/visualization/nviz)) seems to work GRASS65. I need i.pr. .. but still wrong. I have done: cd /add-on/i.pr/ make MODULE_TOPDIR=/usr/local/src/grass6_devel/ PRLIB make[1]: ingresso nella directory /home/gab/add-on/i.pr/PRLIB make[1]: uscita dalla directory /home/gab/add-on/i.pr/PRLIB i.pr_blob make[1]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_blob make htmlcmd make[2]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_blob make /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/docs/html/i.pr.blob.html HTMLSRC=/usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/bin/i.pr.blob make[3]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_blob make[3]: /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/docs/html/i.pr.blob.html è aggiornato. make[3]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_blob make[2]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_blob make mancmd make[2]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_blob make /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/man/man1/i.pr.blob.1 MANSRC=/usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/docs/html/i.pr.blob.html make[3]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_blob make[3]: /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/man/man1/i.pr.blob.1 è aggiornato. make[3]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_blob make[2]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_blob make[1]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_blob i.pr_classify make[1]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_classify make htmlcmd make[2]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_classify make /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/docs/html/i.pr.classify.html HTMLSRC=/usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/bin/i.pr.classify make[3]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_classify make[3]: /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/docs/html/i.pr.classify.html è aggiornato. make[3]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_classify make[2]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_classify make mancmd make[2]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_classify make /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/man/man1/i.pr.classify.1 MANSRC=/usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/docs/html/i.pr.classify.html make[3]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_classify make[3]: /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/man/man1/i.pr.classify.1 è aggiornato. make[3]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_classify make[2]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_classify make[1]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_classify i.pr_features make[1]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_features make htmlcmd make[2]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_features make /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/docs/html/i.pr.features.html HTMLSRC=/usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/bin/i.pr.features make[3]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_features make[3]: /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/docs/html/i.pr.features.html è aggiornato. make[3]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_features make[2]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_features make mancmd make[2]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_features make /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/man/man1/i.pr.features.1 MANSRC=/usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/docs/html/i.pr.features.html make[3]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_features make[3]: /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/man/man1/i.pr.features.1 è aggiornato. make[3]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_features make[2]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_features make[1]: uscita dalla directory /home/gab/add-on/i.pr/i.pr_features i.pr_features_additional make[1]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_features_additional make htmlcmd make[2]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_features_additional make /usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/docs/html/i.pr.features_additional.html HTMLSRC=/usr/local/src/grass6_devel/dist.x86_64-unknown-linux-gnu/bin/i.pr.features_additional make[3]: ingresso nella directory /home/gab/add-on/i.pr/i.pr_features_additional make[3]:
[GRASS-user] i.pr -- g.extension and compiling error on ubuntu 11.10 64 bit
Ciao. I have installed GRASS (6.4.1) from the repository and I can not install i.pr (add-on). With g.extension (with the script downloaded from http://trac.osgeo.org/grass/browser/grass/branches/releasebranch_6_4/scripts/g.extension/g.extension) and also g.extension present in my installation of GRASS I have error): (Mon Jan 30 22:57:31 2012) g.extension -s extension=i.pr svnurl=http://svn.osgeo.org/grass/grass-addons/grass6 /GRASS_ADDON_PATH is not defined, installing to ~/.grass6/addons/ Fetching i.pr from GRASS-Addons SVN (be patient)... Ai.pr/i.pr_subsets Ai.pr/i.pr_subsets/main.c Ai.pr/i.pr_subsets/Makefile Ai.pr/i.pr_subsets/old ... . Ai.pr/i.pr_features_selection/Makefile Ai.pr/Makefile Estratta revisione 50573. Compiling i.pr... /usr/lib/grass64/include/Make/Grass.make:418: attenzione: sovrascrittura dei comandi per l'obiettivo /home/gab/stage_ master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin /usr/lib/grass64/include/Make/Grass.make:409: attenzione: comandi obsoleti per l'obiettivo /home/gab/stage_master/GRA SSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin ignorati /usr/lib/grass64/include/Make/Grass.make:418: attenzione: sovrascrittura dei comandi per l'obiettivo /home/gab/stage_ master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin /usr/lib/grass64/include/Make/Grass.make:409: attenzione: comandi obsoleti per l'obiettivo /home/gab/stage_master/GRA SSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin ignorati mkdir -p /home/gab/stage_master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab-laptop/30962.0c entropy.c gcc '-I/usr/lib/grass64/include' '-I/home/gab/.grass6/addons/include' '-I/home/gab/stage_master ... . /usr/lib/grass64/include/Make/Grass.make:418: attenzione: sovrascrittura dei comandi per l'obiettivo /home/gab/stage_ master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin /usr/lib/grass64/include/Make/Grass.make:409: attenzione: comandi obsoleti per l'obiettivo /home/gab/stage_master/GRA SSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin ignorati /bin/sh: /usr/lib/grass64/tools/mkhtml.sh: not found make htmlcmd make[2]: ingresso nella directory /home/gab/stage_master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab-laptop/30962.0/i.pr/i.pr_blob .. /usr/lib/grass64/include/Make/Grass.make:418: attenzione: sovrascrittura dei comandi per l'obiettivo /home/gab/stage_ master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin /usr/lib/grass64/include/Make/Grass.make:409: attenzione: comandi obsoleti per l'obiettivo /home/gab/stage_master/GRA SSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin ignorati /usr/lib/grass64/include/Make/Grass.make:418: attenzione: sovrascrittura dei comandi per l'obiettivo /home/gab/stage_ master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin /usr/lib/grass64/include/Make/Grass.make:409: attenzione: comandi obsoleti per l'obiettivo /home/gab/stage_master/GRA SSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin ignorati sh: Can't open /build/buildd/grass-6.4.1/tools/g.html2man/g.html2man make[5]: *** [/home/gab/stage_master/GRASSDATA/utm_wgs84/dot torato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/man/man1/i.pr.features.1] Errore 127 make[4]: *** [mancmd] Errore 2 make[3]: *** [cmd] Errore 2 /usr/lib/grass64/include/Make/Grass.make:418: attenzione: sovrascrittura dei comandi per l'obiettivo /home/gab/stage_ master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin /usr/lib/grass64/include/Make/Grass.make:409: attenzione: comandi obsoleti per l'obiettivo /home/gab/stage_master/GRA SSDATA/utm_wgs84/dottorato_bari/.tmp/gab- laptop/30962.0/dist.x86_64/bin ignorati make[5]: ingresso nella directory /home/gab/stage_master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab-laptop/30962.0/i.pr/i.pr_blob mkdir -p /home/gab/stage_master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab-laptop/30962.0/dist.x86_64/man/man1 GRASS_PERL=/usr/bin/perl VERSION_NUMBER=6.4.1 sh /build/buildd/grass-6.4.1/tools/g.html2man/g.html2man /home/gab/stage_master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab-laptop/30962.0/dist.x86_64/docs/html/i.pr.blob.html /home/gab/stage_master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab-laptop/30962.0/dist.x86_64/man/man1/i.pr.blob.1 1 make[5]: uscita dalla directory /home/gab/stage_master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab-laptop/30962.0/i.pr/i.pr_blob make[4]: uscita dalla directory /home/gab/stage_master/GRASSDATA/utm_wgs84/dottorato_bari/.tmp/gab-laptop/30962.0/i.pr/i.pr_blob make[3]:
[GRASS-user] Re: i.pr -- g.extension and compiling error on ubuntu 11.10 64 bit
Markus Neteler wrote On Mon, Jan 30, 2012 at 11:16 PM, Gabriele N. lt;gis.gn@gt; wrote: Ciao. I have installed GRASS (6.4.1) from the repository and I can not install i.pr (add-on). With g.extension ... unfortunately 6.4.1 is too old in a sense that on Ubuntu problems are known. They have been resolved (hopefully) since then in SVN and will be official in 6.4.2 (or use already the current 6.4.svn). Markus ___ grass-user mailing list grass-user@.osgeo http://lists.osgeo.org/mailman/listinfo/grass-user Ok, grazie Markus. At this point, I have to compile GRASS, and then install i.pr. but I have the other problem: At this point I tried to compile GRASS 6.5 and then I wanted to install i.pr but ... #configure sudo CFLAGS=-g -Wall ./configure --prefix=/usr/local/src/grass6_devel --with-cxx --with-python --with-widgets --with-wxwidgets --with-freetype=yes --with-ffmpeg=yes --with-postgres=yes --with-sqlite --with-motif=no --with-glw=no --enable-largefile=yes --with-freetype-includes=/usr/include/freetype2 --with-postgres-includes=/usr/include/postgresql --with-opengl-libs=/usr/include/GL --with-proj-share=/usr/share/proj --without-fftw --without-ffmpeg --with-readline=no --with-tcltk-includes=/usr/include/tcl8.5/ --with-gdal --enable-64bit ...and checking whether to use Tcl/Tk... yes checking for location of Tcl/Tk includes... /usr/include/tcl8.5/ configure: error: *** Tcl/Tk includes directory /usr/include/tcl8.5/ does not exist. tcl8.5 is in /usr/share/tcltk/tcl8.5/ I have also tried: sudo ln -s /usr/share/tcltk/tcl8.5/ /usr/include/ but..nothing Thanks for any help -- View this message in context: http://osgeo-org.1560.n6.nabble.com/i-pr-g-extension-and-compiling-error-on-ubuntu-11-10-64-bit-tp4352067p4352165.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Subset several bands of data at one go
Ciao Chetan. I do not know if I understand your problem. I advise you to create a mask (which defines the area in respect of which cut the raster) and call the raster with similar names (eg landsat_1, landsat_2,). At this point you should do something like this: for i in ` g.mlist type=rast pattern=land* ` do r.mapcalc $i = ( $i ) Ciao Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Subset-several-bands-of-data-at-one-go-tp6118506p6134798.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] third raster - as a function of two input raster (r.series?)
Ciao grass users. I have two raster. 1) raster 1: n: 2728570 minimum: -0.791045 maximum: 0.774648 range: 1.56569 mean: 0.138333 mean of absolute values: 0.185515 standard deviation: 0.158265 variance: 0.025048 variation coefficient: 114.409 % sum: 377452.4525757029 2) raster 2: n: 2645080 minimum: -141.448 maximum: 88.4339 range: 229.882 mean: -0.00535481 mean of absolute values: 20.1545 standard deviation: 23.7113 variance: 562.226 variation coefficient: -442804 % sum: -14163.8968614182 I would like in output a third raster that contains only the areas that have high values of the first raster and with high values of the second raster (and with the areas that have low values of the first raster and low values of the second raster). Obviously the scale of values are different. To do this I thought of using r.series or r.statistics but perhaps is not the right solution. Help / advice? Thanks Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/third-raster-as-a-function-of-two-input-raster-r-series-tp6067501p6067501.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: third raster - as a function of two input raster (r.series?)
Ciao Marcello. I thought there was another solution but I adopted the solution proposed by you. Thank you very much Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/third-raster-as-a-function-of-two-input-raster-r-series-tp6067501p6068331.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: import more kml from directory
Ciao. Thanks. Now, although crude, it works. However, I put exit (or exit 0 exit 1) at the end but the script does not finish properly. What do you say? Thanks http://osgeo-org.1803224.n2.nabble.com/file/n6057767/v.import_kml.sh v.import_kml.sh -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/import-more-kml-from-directory-tp6050317p6057767.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: import more kml from directory
Ciao Hamish and thanks. It works ok :) But I have also tried with my small script (attached) that does not work and I did not understand why. Maybe, if you want, you can watch it. Thank you very much. http://osgeo-org.1803224.n2.nabble.com/file/n6052924/v.import_kml.sh v.import_kml.sh -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/import-more-kml-from-directory-tp6050317p6052924.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] import more kml from directory
Ciao Grass users. Probably my question is trivial, but I can not fix. 1) I want to import many files. kml (all in the same folder) and then I would like a merge(with v.patch). I tried to import v.in.ogr all together maybe I'm wrong. 2) Maybe I can do a merge of all files with ogr2ogr? Thanks for any advice -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/import-more-kml-from-directory-tp6050317p6050317.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.mapcalc -- 3 maps
Hy Glynn, thanks for the explanation. Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-mapcalc-3-maps-tp5929440p5952651.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.reclass.area...aborted
For information: I tried r.reclass.area under windows (before I did it from ubuntu) always with grass 6.4. Under widows will not give me any warning but seems that the result is the same. Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-reclass-area-aborted-tp5949232p5952657.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.reclass.area...aborted
Ciao grass users. I have a warning when I launch the command r.reclass.areaaborted. The map, however, is generated but I'm not sure that everything is correct. What do you say? Thanks r.reclass.area input=fc_map_1001_int_fl1_NO_6_bin_fl2fl1ec_recode_meno_10_11___15 output=fc_map_1001_int__min_10pix lesser=62.5 Generating a clumped raster map... Pass 1... Pass 2... r.clump completo. 186967 clumps. Generating a reclass rules file with area size less than or equal to 62.5 hectares... Generating output raster map fc_map_1001_int__min_10pix... Aborted (Fri Jan 21 21:02:34 2011) Comando terminato (197 sec) -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-reclass-area-aborted-tp5949232p5949232.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.mapcalc -- 3 maps
Excellent Maciek and Glynn :) Although the logic is different (in the solution of Maciek I do a check before the NULL value, while in the solution of Glynn I do a check before non-NULL values), I believe that the result is the same. One note: with the solution of Glynn (at least using the shell from QGIS) it works with single quotes ' and not with double quote . gab@gab-laptop:~$ r.mapcalc output_2 = if(!isnull(a),a,if(!isnull(b),b,c)) bash: !isnull: event not found gab@gab-laptop:~$ r.mapcalc 'output_2 = if(!isnull(a),a,if(!isnull(b),b,c))' # OK Thank you very much Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-mapcalc-3-maps-tp5929440p5932658.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.mapcalc -- 3 maps
Ciao grass users. I have three maps A, B, C. I would use r.mapcalc to combine them according to this logic: Put the values of A and if A is null you put the values of B and if B is null you put the values of C. I have tried in various ways without success. .. now I'm trying with a double IF... r.mapcalc output = if((isnull(a)),b,if(isnull(b)),c) .but in this case I do not put the values of A Help / advice? Thanks -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-mapcalc-3-maps-tp5929440p5929440.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.mapcalc -- 3 maps
Hi Carlos and thanks. A and B do not overlap, and I should take the values of A and B also (always). Then.in areas that are empty, I would then put the values of C. I tried withr.mapcalc output = if((isnull(a)),b,if(isnull(b),a),c) but I lose the values of C -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-mapcalc-3-maps-tp5929440p5929564.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.mapcalc -- 3 maps
Hi Sylvain and thanks. I had not thought of r.patch and I think that works (but I have to carefully examine the result) I did so r.patch input = a, b, c = output output_patch Now I want to do it with r.mapcalc:) -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-mapcalc-3-maps-tp5929440p5929594.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Arcgrid - strange values
Hi Daniel, Hi Frank. I tried to export in different ways (as GeoTIFF, as float etc. ..) and import it into arc in different ways. At this point, I think it is a limitation of arc someone knows something? Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Arcgrid-strange-values-tp5802385p5827736.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Importing of raster data doesn't seem to work through GUI
Hy Chetan. Look here, http://grass.osgeo.org/wiki/LANDSAT maybe it can be useful Ciao Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Importing-of-raster-data-doesn-t-seem-to-work-through-GUI-tp5805160p5805396.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Best operating system for GRASS
Hi. My advice is to use linux. In particular, I use Ubuntu (Linux human:)). However GRASS runs well on wind but I think that using Linux has many more advantages. For example, through the repository every time there is an updated version of GRASS you can get it automatically. Morover (if I'm not mistaken) on the system wind is not exactly how it works on linux. Ciao Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Best-operating-system-for-GRASS-tp5798774p5802364.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Arcgrid - strange values
Hi. I am working with colleagues that working on ArcGIS. It 's okay but now I have a problem with esri grid. I am using to export from GRASS r.out.arc, but the values of the cells in arcgis are different from what I read in GRASS. The raster has values with 13 digits. Here's how I export: r.out.arc input=map output=map13.asc dp=13 I tried to import into ArcGIS as integer or as a float (ArcMap - Toolbox - Conversion Tools - To Raster - ASCII to Raster ). Help / advice? Thanks Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Arcgrid-strange-values-tp5802385p5802385.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] reprojection from Projection Mollweide
Hy everyone. I have a problem with a reprojection. I have a global raster (a grid esri) that I imported (by choosing the file w001001.adf) in one location. I created the location with a reference system that has these parameters (in a file prj.adf): Projection Mollweide Units METERS Zunits NO Xshift 0.0 Yshift 0.0 Parameters 6370997.0 6370997.0 0 0 0.0 / * Longitude of projection center So far so ok. Then I tried to reproject the raster in a location with ETRS_1989_LAEA but I get (at least on Italy) a strong localization error. I tried also with WGS84 (EPSG 4326) but nothingand even from the command line making various adjustments to the parameters .. for example: gdalwarp-s_srs proj = moll ellps sphere lon_0 = 0 = x_0 y_0 = 0 = 0 a=6370997 b=6370997 no_defs units=m -t_srs EPSG : 4326 map.tif map_wgs84.tif Does anyone have experience on? Thanks Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/reprojection-from-Projection-Mollweide-tp5797013p5797013.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: reprojection from Projection Mollweide
I found some other posts I'm getting close to the solution (I get good results): gdalwarp -s_srs '+proj=goode +ellps=sphere +lon_0=30E +towgs84=0,0,0,0,0,0,0 +x_0=3335846.22854 +y_0=2' -t_srs '+proj=latlong +ellps=wgs84 towgs84=0,0,0,0,0,0,0' -te 7 33 34 59 map.tif map_wgs84.tif I do not know why, but that's okay with +proj=good and not +proj = moll In addition, by removing the parameter y_0 = 2 I get a more accurate overlay. So by the time I reprojected raster in the WGS84 as follows: gdalwarp -s_srs ' proj=goode ellps=sphere lon_0=30E towgs84=0,0,0,0,0,0,0 x_0=3335846.22854' -t_srs ' proj=latlong ellps=wgs84 towgs84=0,0,0,0,0,0,0' -te 7 33 34 59 map.tif map_wgs84.tif Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/reprojection-from-Projection-Mollweide-tp5797013p5797451.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] file .img and Interrupted Goode Homolosine
Hy everyone. Two problems: 1) I can not import a file erdas why is not supported. gdalinfo image.img ERROR 4: `image.img' not recognised as a supported file format. gdalinfo failed - unable to open 'image.img'. 2) Projection: Interrupted Goode Homolosine I could not find the EPSG code, and then I did not create the location...however, I found the following parameters (but I do not know if they are correct -s_srs '+proj=goode +ellps=sphere +towgs84=0,0,0 +lon_0=100w +x_0=-9487.43' Grass 6.4.0 Thanks for any advice Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/file-img-and-Interrupted-Goode-Homolosine-tp5778028p5778028.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: file .img and Interrupted Goode Homolosine
Hi Martin. Yes, strange ... I tried also with QGIS, and obviously does not work ... : / This file is strange:) .. also because I can not understand how to create the location with the projection Interrupted Goode Homolosine...boooh Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/file-img-and-Interrupted-Goode-Homolosine-tp5778028p5778080.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster -
Glynn Clements wrote: The same result can be achieved with one command: r.mapcalc output = A * 1.0 + B The multiply will cause A to be converted to double to match the other argument (floating-point constants are double unless an explicit f suffix is given), the multiply will yield a double result, and the addition will convert B to double to match this. perfect, I tried it and it actually works:) This speeds up operations. Now I have to combine new maps and then to reclassify the maps, I will use r.recode ... we'll see what happens:) Thank you Glynn Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5772145.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster -
Hy Daniel. I do not think that r.cross can solve the problem. With r.cross I could get a concatenation of label: For example, with values A = 12 B = 3323 C = 1 the result with the label category 12, category 3323, category 1 However, if you know a chance I can prove otherwise. As a result, the value of cell in the raster output should be 12332301 Ciao Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5766350.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster -
yes, maybe with grep and r.cross ... but then you think you can convert the result of r.cross as cell value? Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5766761.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster -
I have a new problem, this time I would like to concatenate two maps A and B. A is with 8 digits, B is with 4 digits. So I do as usual: r.mapcalc output = (( A * 1) ( B )) But I get strange values with negative values. I did several tests but without resolve. Perhaps r.mapcalc may have problems with numbers greater than 10 digits Thanks Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5767794.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster -
Hy Glynn. Thanks for the interesting explanation. I did then so: r.mapcalc A_new = (double(A)) and r.mapcalc B_new = (double(B)) and then r.mapcalc output = ((A_new * 1)+(B_new )) I need to check the result, but it seems that the result is right. Thank you very much Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5768432.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster -
Daniel Victoria wrote: It's not r.mapcalc fault. I imagine it's probably a limitation due to the number of bits in your image. For instance, 8 bit images can have up to 256 values. 16 bit images can have 65,535. So, if you have too many values, you will have problems. Not sure how to avoid this in your case but, maybe you will need to work with categories after all. Are you planing on doing any math operations in the concatenated values or is it more for displaying pourposes? Hi Daniel... Maybe it would be useful only concatenation but I would also have the possibility to apply tools such as r.statistics and other similar that work on the value of cells. Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5768618.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster - to link together
Sorry, I have not explained well. The problem was related to the value of C map. That is, I take the example above. If the value of C map is equal to 1, I get a result - A =12 B=3323 C=1 =12332301 (8 digit number) and instead could be useful 1233231 (number 7-digit). However I checked and I need a value to 8 digits . so ok:) Thank you very much Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5744524.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Values from multiple raster - to link together
Hello everyone. I have to solve a problem already discussed here: http://osgeo-org.1803224.n2.nabble.com/r-mapcalc-concatenates-values-td1884970.html#a1884973 In summary, I have for example 3 maps A, B, C. I need to combine the values of the overlapping cells. For example if I have a cell in the map A with value 11, in the corresponding cell of the map B I have and the corresponding cell of the map C is . In output I would like to have a map that has a value in that cell 11. I tried with r.mapcalc, but without Successfully. Help? Thanks Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5741032.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster - to link together
Excellent Alex:). The maps are as follows: A - value 0, 1, 2, ... 19 B - value 0, 1001, 1002, 1003 3351 C - value from 0, 1, 2, 3, 4 .. 15 So, I did as you suggested: r.mapcalc output = (( A *100)+( B *100)+ C ) right? And if I had decimal numbers or strings instead of numbers? Thanks Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5741794.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: g.extension (e.g. to install i.landsat.toar)
Hi Markus. I have installed grass-dev but now there is another problem (there is a problem of permissions for the folder creation and does not solve if I create the folder before launching g.extension): g.extension i.landsat.toar Fetching i.landsat.toar from GRASS-Addons SVN (be patient)... Ai.landsat.toar/landsat_set.c Ai.landsat.toar/local_proto.h Ai.landsat.toar/main.c Ai.landsat.toar/description.html Ai.landsat.toar/landsat.c Ai.landsat.toar/earth_sun.c Ai.landsat.toar/landsat.h Ai.landsat.toar/landsat_met.c Ai.landsat.toar/Makefile Ai.landsat.toar/earth_sun.h U i.landsat.toar Estratta revisione 44172. Compiling i.landsat.toar... mkdir -p /build/buildd/grass-6.4.0/dist.i686-pc-linux-gnu/include/grass mkdir: impossibile creare la directory `/build/buildd/grass-6.4.0/dist.i686-pc-linux-gnu': Permesso negato make: *** [/build/buildd/grass-6.4.0/dist.i686-pc-linux-gnu/include/grass] Errore 1 ERRORE: Compilation failed, sorry. Please check above error messages. Thanks -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/g-extension-e-g-to-install-i-landsat-toar-tp5662708p5698083.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: g.extension (e.g. to install i.landsat.toar)
Markus, still does not work. I saw some differences compared to the old script. To understand ... which is the part that has been changed to fix this? GRASS 6.4.0 (utm_wgs84):~ g.extension i.landsat.toar Fetching i.landsat.toar from GRASS-Addons SVN (be patient)... Ai.landsat.toar/landsat_set.c Ai.landsat.toar/local_proto.h Ai.landsat.toar/main.c Ai.landsat.toar/description.html Ai.landsat.toar/landsat.c Ai.landsat.toar/earth_sun.c Ai.landsat.toar/landsat.h Ai.landsat.toar/landsat_met.c Ai.landsat.toar/Makefile Ai.landsat.toar/earth_sun.h U i.landsat.toar Estratta revisione 44172. Compiling i.landsat.toar... mkdir -p /build/buildd/grass-6.4.0/bin.i686-pc-linux-gnu mkdir: impossibile creare la directory `/build': Permesso negato make: *** [/build/buildd/grass-6.4.0/bin.i686-pc-linux-gnu] Errore 1 ERRORE: Compilation failed, sorry. Please check above error messages. GRASS 6.4.0 (utm_wgs84):~ g.extension i.landsat.acca Fetching i.landsat.acca from GRASS-Addons SVN (be patient)... Ai.landsat.acca/tools.c Ai.landsat.acca/local_proto.h Ai.landsat.acca/main.c Ai.landsat.acca/description.html Ai.landsat.acca/algorithm.c Ai.landsat.acca/Makefile U i.landsat.acca Estratta revisione 44172. Compiling i.landsat.acca... mkdir -p /build/buildd/grass-6.4.0/bin.i686-pc-linux-gnu mkdir: impossibile creare la directory `/build': Permesso negato make: *** [/build/buildd/grass-6.4.0/bin.i686-pc-linux-gnu] Errore 1 ERRORE: Compilation failed, sorry. Please check above error messages. -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/g-extension-e-g-to-install-i-landsat-toar-tp5662708p5699112.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: g.extension (e.g. to install i.landsat.toar)
I have the same problem (ubuntu 10.10, with repo ubuntugis) GRASS 6.4.0 (utm_wgs84):~ g.extension i.landsat.toar Fetching i.landsat.toar from GRASS-Addons SVN (be patient)... Ai.landsat.toar/landsat_set.c Ai.landsat.toar/local_proto.h Ai.landsat.toar/main.c Ai.landsat.toar/description.html Ai.landsat.toar/landsat.c Ai.landsat.toar/earth_sun.c Ai.landsat.toar/landsat.h Ai.landsat.toar/landsat_met.c Ai.landsat.toar/Makefile Ai.landsat.toar/earth_sun.h U i.landsat.toar Estratta revisione 44145. Compiling i.landsat.toar... Makefile:8: /usr/lib/grass64/include/Make/Module.make: Nessun file o directory make: *** Nessuna regola per generare l'obiettivo «/usr/lib/grass64/include/Make/Module.make». Stop. ERRORE: Compilation failed, sorry. Please check above error messages. /usr/lib/grass64$ ls AUTHORS.gz bwidget docsetcGPL.TXT locale bin COPYING driver fonts lib scripts Thanks -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/g-extension-e-g-to-install-i-landsat-toar-tp5662708p5693486.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: R, Moran index, getis e ord.....
Thanks Robbie. I try to take a look. Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/R-Moran-index-getis-e-ord-tp5588509p5599570.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] R, Moran index, getis e ord.....
Hello everyone. I have a raster which is the result of calculating an index of remote sensing (with decimals). I would like to calculate the Moran's index and the index of Getis and Ord To do this I installed R (following http://grass.fbk.eu/statsgrass/grass6_r_install.html). Then I installed the packages install.packages(spdep,dependencies=TRUE) install.packages(maptools,dependencies=TRUE) Then from the shell: R Then load these libraries library(spgrass6) library(spdep) library(maptools) library(boot) library(Matrix) library(lattice) I read the raster index_water_2009-readRAST6(index_water_2009) At this point I stopped. I saw on the mailing list and some tutorials. I can not understand the syntax to use. It only works on vector? Someone can give me some help? Thanks -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/R-Moran-index-getis-e-ord-tp5588509p5588509.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] link image or application in DB sqlite
Hi to everyone. I have a vector of points and I want to connect each point to a photo and a link (to open a photo or a web page with the browser firefox). I have the attributes in SQLite. With sqlite manager I added a column 'image' of type blob. I uploaded an image but does not work. When I use the identifier ..: GRASS 6.5.svn (32633):/home/gab/GRASSDATA WARNING: SQLite driver: unable to parse decltype: BLOB WARNING: SQLite driver: unable to parse decltype: BLOB WARNING: SQLite driver: column 'image', SQLite type 4 is not supported WARNING: SQLite driver: unable to parse decltype: BLOB Alternatively, I tried to do a field by putting the path of the images but I do not know if the identifier can open the image. I should get it to work with QGIS. In addition, it also works on windows XP?Otherwise I can use a virtual machine. Help / suggestions? Thanks Gabriele PS Grass 65 ubuntu 8.10 -- View this message in context: http://n2.nabble.com/link-image-or-application-in-DB-sqlite-tp2594524p2594524.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] link one to many-sqlite
Hi. I'm working with grass64 on ubuntu 8.10. I have created a sqlite db. 1) I have a shape of municipalities that have problems of topology. In addition, some municipalities have more non-adjacent areas. To not have the municipalities 'split' I tried to import without topological corrections but I do not think is the right thing. v.dissolve does not work if the polygons are not adjacent. Right? 2) In addition to the linked table, I want to connect another table but it has many records for each municipalities and then to do a link one to many (see picture example in attachment). link_table.gif Should I put this table on layer 2? or create a copy of the shape and try to connect this table? I have attached the table with v.db.connect on layer 2 and it seems ok, but if I query the poligon (identify in display mode) opens the table only on layer 1. Tips / help? Thanks Gabriele -- View this message in context: http://n2.nabble.com/link-one-to-many-sqlite-tp2541600p2541600.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] GRASS 7 and SQLITE o POSTGRES
Hi Markus, Benjamin Markus Neteler wrote: Last Friday I met the SpatiaLite developer at the Italian GRASS/GFOSS meeting (http://gfoss2009.crs4.it/) which was a great conference btw. He is interested in linking GRASS and SpatiaLite, so maybe later this year new features are added to put GRASS vector and also raster data into SpatiaLite. Markus this is really a good news :-) Benjamin Ducke wrote: That would be simply great! With SpatiaLite support evolving in GRASS, QGIS and (perhaps) gvSIG, there would finally be a sane way for direct exchange of vector data between open source GIS without crippling it by going through abominations like Shapefiles. Ben could be an asset important to work with open source. Thanks Gabriele - Original Message - From: Markus Neteler nete...@osgeo.org To: Moritz Lennert mlenn...@club.worldonline.be Cc: GRASS user list grass-user@lists.osgeo.org, Gabriele N. gis...@libero.it Sent: Sunday, March 1, 2009 9:10:14 PM GMT +00:00 GMT Britain, Ireland, Portugal Subject: Re: [GRASS-user] GRASS 7 and SQLITE o POSTGRES On Sat, Feb 21, 2009 at 1:42 PM, Moritz Lennert mlenn...@club.worldonline.be wrote: ... The SQLite support in GRASS concerns the attribute management, not geometries. To access geometries in a SpatiaLite, you have to go through v.external. Last Friday I met the SpatiaLite developer at the Italian GRASS/GFOSS meeting (http://gfoss2009.crs4.it/) which was a great conference btw. He is interested in linking GRASS and SpatiaLite, so maybe later this year new features are added to put GRASS vector and also raster data into SpatiaLite. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Files attached to this email may be in ISO 26300 format (OASIS Open Document Format). If you have difficulty opening them, please visit http://iso26300.info for more information. ___ 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/GRASS-7-and-SQLITE-o-POSTGRES-tp2360247p2408210.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] GRASS 7 and SQLITE o POSTGRES
Given the limitations that exist with using GRASS and Postgres I do not know how. There is someone who manages a GIS (with others), based mainly on GRASS? 1) If yes, how do? how organized these data? how to manage them? 2) If I put everything on postgres, To do spatial analysis with GRASS I am forced to import data from Postgres. right? Thanks for any advice. Moritz Lennert wrote: On 21/02/09 14:25, Gabriele N. wrote: Hi Moritz. I have read here http://grass.osgeo.org/wiki/GRASS_7_ideas_collection Database establish SQLite as default DBMI driver (DBF is too limited). done May 2008. That's because I thought that SQLite support in GRASS 7 concerns the attribute management and the geometries. No, only attribute management. Moritz ___ 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/GRASS-7-and-SQLITE-o-POSTGRES-tp2360247p2405193.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] GRASS 7 and SQLITE o POSTGRES
Hi Benjamin, thanks for the reply. I read something about spatial lite. I work mainly with GRASS and thought to use it with SQLite because GRASS 7 should be able to work directly on the SQLite for the geometries and the attributes (in short, without v.exeternal). I think, but I do not know if it is correct, that implementing a DB in SQLite might be possible to access data with GRASS 7, QGIS, gvSIG, etc... (my colleagues still do not use GRASS) but maybe it's only my hope :-) With postgres/postgis I should always connect with v.external and I could not do spatial analysis using the tools of GRASS ... right? Thank you very much Gabriele Benjamin Ducke wrote: Hi Gabriele, A spatial database extension allows you to store geographic features as part of the database itself. There are many advantages to this, especially if you work with huge vector datasets. SQLite has a spatial extension equivalent to PostGIS. It is called SpatiaLite. It has the same functionality as PostGIS. The latest version of OGR already has some basic support for SQLite, so hopefully in the near future all open source GIS will support SQLite/SpatiaLite as a datasource. However, if you are planning to to set up a spatial data infrastructure for collaborative work, you should probably choose PostgreSQL/PostGIS as your data backend, because it supports user access control and is currently better supported by open source GIS. Ben - Original Message - From: Gabriele N. gis...@libero.it To: grass-user@lists.osgeo.org Sent: Friday, February 20, 2009 5:38:48 PM GMT +00:00 GMT Britain, Ireland, Portugal Subject: [GRASS-user] GRASS 7 and SQLITE o POSTGRES Hello. I work with colleagues sharing the files (shapes .. etc) that are on a NAS server. I use ubuntu 8.10 (with GRASS locally) and my colleagues use windows. I would do this: - A GIS with all these data - Building a database (sqlite or postgres/postgis/ on NAS) - An excellent organization of the data - To enable access to data with other software (such as QGIS and gvSIG) GRASS 7 will default SQLite and therefore should be different to GRASS 6. In the sense that now GRASS 6 connects to the DB (postgres and / or sqlite) but works differently than the dbf ... or wrong? SQLite does not have the spatial component as postgis with postgres? What does this? Can you tell me the steps to follow? What do you recommend? Thanks Gabriele -- View this message in context: http://n2.nabble.com/GRASS-7-and-SQLITE-o-POSTGRES-tp2360247p2360247.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 -- Files attached to this email may be in ISO 26300 format (OASIS Open Document Format). If you have difficulty opening them, please visit http://iso26300.info for more information. ___ 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/GRASS-7-and-SQLITE-o-POSTGRES-tp2360247p2363407.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] GRASS 7 and SQLITE o POSTGRES
Hello. I work with colleagues sharing the files (shapes .. etc) that are on a NAS server. I use ubuntu 8.10 (with GRASS locally) and my colleagues use windows. I would do this: - A GIS with all these data - Building a database (sqlite or postgres/postgis/ on NAS) - An excellent organization of the data - To enable access to data with other software (such as QGIS and gvSIG) GRASS 7 will default SQLite and therefore should be different to GRASS 6. In the sense that now GRASS 6 connects to the DB (postgres and / or sqlite) but works differently than the dbf ... or wrong? SQLite does not have the spatial component as postgis with postgres? What does this? Can you tell me the steps to follow? What do you recommend? Thanks Gabriele -- View this message in context: http://n2.nabble.com/GRASS-7-and-SQLITE-o-POSTGRES-tp2360247p2360247.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] script by bash to Python
Hi. All my scripts in bash work with tcltk. Now, I take as example my 2 scripts: v.to.dem.sh and v.impact.sh. With WxPython all the scripts were not working with grass 63. With Grass64, instead, some scripts works (es v.impact.sh) and some other scripts does not works (es.v.to.dem.sh). But I do not know what may be the characteristic that makes v.impact.sh work and that does not work v.to.dem.sh. I hope to have been clearer. Thanks :-) Gabriele PS The error is always the same. Here's an example: GRASS 6.4.svn (gb2):~ v.to.dem.sh Traceback (most recent call last): File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/menuform.py, line 1687, in module GrassGUIApp( grassTask( sys.argv[1] ) ).MainLoop() File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/menuform.py, line 193, in __init__ xml.sax.parseString( getInterfaceDescription( grassModule ) , processTask( self ) ) File /usr/lib/python2.5/xml/sax/__init__.py, line 49, in parseString parser.parse(inpsrc) File /usr/lib/python2.5/xml/sax/expatreader.py, line 107, in parse xmlreader.IncrementalParser.parse(self, source) File /usr/lib/python2.5/xml/sax/xmlreader.py, line 123, in parse self.feed(buffer) File /usr/lib/python2.5/xml/sax/expatreader.py, line 211, in feed self._err_handler.fatalError(exc) File /usr/lib/python2.5/xml/sax/handler.py, line 38, in fatalError raise exception xml.sax._exceptions.SAXParseException: unknown:21:73: not well-formed (invalid token) [Raster MASK present] GRASS 6.4.svn (gb2):~ hamish_b wrote: Gabriele wrote: I speak again of this post. With grass64 svn the situation has improved. Some scripts in bash run well and others do not. which ones work, which ones don't? what does others do not run well mean? The same Python GUI error, or..? But I do not understand what are the differences between a script that works and a script that does not work. I do not see the differences. need more details please. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://www.nabble.com/script-by-bash-to-Python-tp18292175p19645026.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] script by bash to Python
A step forward :-) I think that the letters with the accent in the code create problems. I replaced the letters è with e ' and now the script v.to.dem.sh work. Still does not go well: give me an error message when I launch (picture attached) and does not work the STOP button but gives a result in output. Unfortunately, doing the same with other scripts I did not have the same effect. I try to understand what other anomalies can be. Thank you. I communicate further progress. Gabriele http://www.nabble.com/file/p19657956/Screenshot-Error%2Bin%2Bcommand%2Bexecution.png Screenshot-Error+in+command+execution.png Glynn Clements wrote: Gabriele N. wrote: All my scripts in bash work with tcltk. Now, I take as example my 2 scripts: v.to.dem.sh and v.impact.sh. With WxPython all the scripts were not working with grass 63. With Grass64, instead, some scripts works (es v.impact.sh) and some other scripts does not works (es.v.to.dem.sh). But I do not know what may be the characteristic that makes v.impact.sh work and that does not work v.to.dem.sh. I hope to have been clearer. Thanks :-) Gabriele PS The error is always the same. Here's an example: GRASS 6.4.svn (gb2):~ v.to.dem.sh Traceback (most recent call last): File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/menuform.py, line 1687, in module GrassGUIApp( grassTask( sys.argv[1] ) ).MainLoop() File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/menuform.py, line 193, in __init__ xml.sax.parseString( getInterfaceDescription( grassModule ) , processTask( self ) ) File /usr/lib/python2.5/xml/sax/__init__.py, line 49, in parseString parser.parse(inpsrc) File /usr/lib/python2.5/xml/sax/expatreader.py, line 107, in parse xmlreader.IncrementalParser.parse(self, source) File /usr/lib/python2.5/xml/sax/xmlreader.py, line 123, in parse self.feed(buffer) File /usr/lib/python2.5/xml/sax/expatreader.py, line 211, in feed self._err_handler.fatalError(exc) File /usr/lib/python2.5/xml/sax/handler.py, line 38, in fatalError raise exception xml.sax._exceptions.SAXParseException: unknown:21:73: not well-formed (invalid token) This indicates the the output from v.to.dem.sh --interface-description isn't valid XML. This is likely to be caused by the contents of the option descriptions at the top of the script (the #% lines). Does one of the entries have unknown as its value? If so, which field (gisprompt, label, etc) is it in? I suspect that the code which generates the XML data isn't always converting , and to lt;, gt; and amp; when it should. You can work around the problem by using e.g. (unknown), but we should really fix this in the code which generates the XML. -- Glynn Clements [EMAIL PROTECTED] ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://www.nabble.com/script-by-bash-to-Python-tp18292175p19657956.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] script by bash to Python
Martin Landa-2 wrote: Hi, 2008/7/5 Gabriele N. [EMAIL PROTECTED]: I use grass63 with wxpython and I do not know how to use some scripts that I have in bash. The graphical interface opens but when I try to launch the script I get the following error: GRASS 6.3.0 (gb2):~ v.topo_5 $ v.topo_5 percorso=/home/gab/prova.shp nome_vett=foglio num=3 campo=foglio Traceback (most recent call last): File /usr/lib/grass/etc/wxpython/gui_modules/menuform.py, line 763, in OnRun self.goutput.RunCmd(cmd) File /usr/lib/grass/etc/wxpython/gui_modules/goutput.py, line 209, in RunCmd generalCmd = subprocess.Popen(cmdlist, NameError: global name 'subprocess' is not defined If I run the script with grass -tcltk works well. Upgrading to grass64 (devbr6) should fix it. Martin I speak again of this post. With grass64 svn the situation has improved. Some scripts in bash run well and others do not. But I do not understand what are the differences between a script that works and a script that does not work. I do not see the differences. Thanks -- View this message in context: http://www.nabble.com/script-by-bash-to-Python-tp18292175p19630269.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] script by bash to Python
Hy Glynn. Ah ... I work on ubuntu 8.04. But what could be the commands in the script that GRASS not 'recognize'? Compared to version 63 with the 64 now run more scripts (always bash). All the scripts still work well in tcltk. ? Thanks Gabriele Glynn Clements wrote: Gabriele N. wrote: I use grass63 with wxpython and I do not know how to use some scripts that I have in bash. The graphical interface opens but when I try to launch the script I get the following error: GRASS 6.3.0 (gb2):~ v.topo_5 $ v.topo_5 percorso=/home/gab/prova.shp nome_vett=foglio num=3 campo=foglio Traceback (most recent call last): File /usr/lib/grass/etc/wxpython/gui_modules/menuform.py, line 763, in OnRun self.goutput.RunCmd(cmd) File /usr/lib/grass/etc/wxpython/gui_modules/goutput.py, line 209, in RunCmd generalCmd = subprocess.Popen(cmdlist, NameError: global name 'subprocess' is not defined If I run the script with grass -tcltk works well. Upgrading to grass64 (devbr6) should fix it. I speak again of this post. With grass64 svn the situation has improved. Some scripts in bash run well and others do not. But I do not understand what are the differences between a script that works and a script that does not work. I do not see the differences. The above error is due to a bug in 6.3.0. The bug is only triggered when running a command which isn't recognised as a GRASS command. Recognised commands are those corresponding to files in the bin, scripts and etc/gui/scripts directories, with any .exe or .bat extension removed. It's also possible that individual scripts have bugs which are specific to Windows. GRASS probably doesn't get a fraction of the testing on Windows that it does on Linux and MacOSX. -- Glynn Clements [EMAIL PROTECTED] ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://www.nabble.com/script-by-bash-to-Python-tp18292175p19637676.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.mapcalc concatenates values?
Hi everyone. I have two questions 1) I have some raster (with the same format and are converted from vector). The cells have the value of CAT. So the value = CAT With r.mapcalc I would define a map of output that I give, as the value for each cell, indications of the different CAT. For example, if I have the raster A that has value = 1 and raster B that has value = 2, I will have raster C output value = 1,2 If I have 3 maps in input, value = 1,2,3 etc. ... I made several attempts but I do not know if r.mapcalc concatenates the values. 2) if r.mapcalc is not possible, I tried to make the transition to vector. The passages that I made are: - r.to.vect maps raster input (es rasterA and rasterB) - V.overlay and I get in dbf, in a column the value of rasterA and in another column the value of rasterB. I've created a third column to concatenate the 2 columns but it seems to me that the operator '| |' does not work with dbf. How do I concatenate the values in 2 columns? Thank you very much Gabriele -- View this message in context: http://www.nabble.com/r.mapcalc-concatenates-values--tp19392795p19392795.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] grass64 and qgis-plugin-grass
Hello. Following this recent post: http://www.nabble.com/Re%3A-grass6.4-develbrunch-revision-32726%3A-undefined-reference-to-%60G_no_gisinit%27-td18947949.html, I installed so grass64 (ubuntu 8.04): I've just built gdal 1.5.2 (without grass) ./configure --with-ecw=yes --without-grass and grass (--with-gdal) CFLAGS=-O3 -mcpu=i686 -Wall LDFLAGS=-s ./configure --with-cxx --with-python --with-widgets --with-wxwidgets --with-freetype=yes --with-ffmpeg=yes --with-postgres=yes --with-gdal --with-sqlite=no --with-motif=no --with-glw=no --enable-largefile=yes --with-tcltk-includes=/usr/include/tcl8.4 --with-freetype-includes=/usr/include/freetype2 --with-postgres-libs=/usr/include/postgresql/libpq --with-postgres-includes=/usr/include/postgresql --with-opengl-libs=/usr/include/GL --with-proj-share=/usr/share/proj --without-fftw --without-ffmpeg --with-readline=no instead of --with-gdal had to make --with-gdal=/usr/local/bin/gdal-config ? make make install and then, gdal-grass plugin 1.4.3 (gdal-grass-1.4.3) ../configure --with-gdal=/usr/local/bin/gdal-config --with-grass=/usr/local/grass-6.4.svn make make install -- --- Now I have a problem with QGIS. If I open qgis from grass63 (installed by the binaries) all is well. If I open QGIS (installed by the binaries) from GRASS64 : Loaded / usr / lib / qgis / libgrassplugin.so qgis: symbol lookup error: / usr/lib/libqgisgrass.so.0.11: undefined symbol: G_no_gisinit If I remove qgis-plugin-grass QGIS is fine. How can I fix? Thank you Gabriele -- View this message in context: http://www.nabble.com/grass64-and-qgis-plugin-grass-tp19030535p19030535.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] v.out.ogr multiple
Thanks Moritz. I have made some improvements. I give you an example and the new script. Now I export all vectors with name that starts with foglio_1 But I do not know how to select a folder where to export the shapes. Now, only works if I enter the path by hand or select from a graphical user interface a file in the folder (prova). But it begins to work :-) Still suggestions? Thanks Gabriele http://www.nabble.com/file/p18359390/Schermata-v.out.ogr_multiplo.sh.png Schermata-v.out.ogr_multiplo.sh.png http://www.nabble.com/file/p18359390/v.out.ogr_multiplo.sh v.out.ogr_multiplo.sh Moritz Lennert-2 wrote: On 09/07/08 10:42, Gabriele N. wrote: Hi Moritz. Moritz Lennert-2 wrote: On 08/07/08 22:33, Gabriele N. wrote: Hello. I want to create a script bash (GUI) to export many vector in shapefile having the opportunity to decide which export layer (0, 1, 2) and type=line,area ecc.. Does this mean interactive decision layer by layer, or that the module should allow for entering layer and type as an option ? I would also like to know if you can use the characters as' * 'to be able to export all instance vector starting with' nome_vector '. There is already some scripts? If not, I started to write something in bash but I am in difficulty. It will have to be something like: for map in `g.mlist type=vect pattern=$GIS_OPT_PATTERN do v.out.ogr $map type=$GIS_OPT_TYPE layer=$GIS_OPT_LAYER done where GIS_OPT_PATTERN would be something like 'nome_vector*' If you need to automatically decide which layers and types to export, you might be able to do some grep/awk magic with v.category -g $map option=report. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user I would like to insert layer and type as an option. Attached is what I have done (GUI and script). However, there are errors. http://www.nabble.com/file/p18356967/Schermata-v.out.ogr_multiplo.sh.png Schermata-v.out.ogr_multiplo.sh.png http://www.nabble.com/file/p18356967/v.out.ogr_multiplo.sh v.out.ogr_multiplo.sh Where I wrong? It would be easier to help you if you gave us more information about the errors... But from looking at your script, I already see this: for map in g.mlist type=vect pattern=$GIS_OPT_PATTERN should be: for map in `g.mlist type=vect pattern=$GIS_OPT_PATTERN` note the backticks (`) and then v.out.ogr input=${map} ... should be v.out.ogr input=$map Moritz Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://www.nabble.com/v.out.ogr-multiple-tp18348459p18359390.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.out.ogr multiple
Hello. I want to create a script bash (GUI) to export many vector in shapefile having the opportunity to decide which export layer (0, 1, 2) and type=line,area ecc.. I would also like to know if you can use the characters as' * 'to be able to export all instance vector starting with' nome_vector '. There is already some scripts? If not, I started to write something in bash but I am in difficulty. Thank you Gabriele -- View this message in context: http://www.nabble.com/v.out.ogr-multiple-tp18348459p18348459.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] script by bash to Python
Hello 1) I use grass63 with wxpython and I do not know how to use some scripts that I have in bash. The graphical interface opens but when I try to launch the script I get the following error: GRASS 6.3.0 (gb2):~ v.topo_5 $ v.topo_5 percorso=/home/gab/prova.shp nome_vett=foglio num=3 campo=foglio Traceback (most recent call last): File /usr/lib/grass/etc/wxpython/gui_modules/menuform.py, line 763, in OnRun self.goutput.RunCmd(cmd) File /usr/lib/grass/etc/wxpython/gui_modules/goutput.py, line 209, in RunCmd generalCmd = subprocess.Popen(cmdlist, NameError: global name 'subprocess' is not defined If I run the script with grass -tcltk works well. 2) In this regard, as I do reopen grass with wxpython? The first few times I feel I get The LOCATION gab does not exist. Please create first. Thank you very much -- View this message in context: http://www.nabble.com/script-by-bash-to-Python-tp18292175p18292175.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] script by bash to Python
Hi Martin. I like that. Thank you very much Gabriele Martin Landa-2 wrote: Hi, Upgrading to grass64 (devbr6) should fix it. generally speaking, if you want to test wxGUI and report its bugs, upgrade to grass64.svn is necessary (wxpython code in grass63 is quite out-dated, just a few serious bugs fixes were backported for 6.3.1). Martin -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa * ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://www.nabble.com/script-by-bash-to-Python-tp18292175p18295088.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] Enclose the polygons within 50 metres
Hi Christian Indeed this was doing. I made the buffer (but on many polygons gives a bit of problems) and now , as suggested by you, I can go to raster with v.to.rast. Or should I try to dissolve the buffer by working again with the vector. Thanks Gabriele Christian Schwartze-2 wrote: Hi, don't know if I unterstood you right, but what about buffering your polygons (v.buffer with buffer=50m), i.e. the buildings. So you get a map of overlapping areas where features not more far away than 50 meters. Afterwards you can use v.to.rast for rasterizing the buffer zones and describing small towns for you case. Regards, Christian. Am Sonntag, den 15.06.2008, 08:43 -0700 schrieb Gabriele N.: Hi all I have 2 vector polygonal; one representing the buildings and another that represents the infrastructure. I think to put them together with v.patch. I have to create another vector polygonal, which is to enclose all the polygons that are close to each other within 50 metres. I need to identify small towns. How can I do? 1) v.patch 2) v.distance? 3)? Thank you very much Gabriele This mail was sent through http://webmail.uni-jena.de ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://www.nabble.com/Enclose-the-polygons-within-50-metres-tp17851388p17919567.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Enclose the polygons within 50 metres
Hi all I have 2 vector polygonal; one representing the buildings and another that represents the infrastructure. I think to put them together with v.patch. I have to create another vector polygonal, which is to enclose all the polygons that are close to each other within 50 metres. I need to identify small towns. How can I do? 1) v.patch 2) v.distance? 3)? Thank you very much Gabriele -- View this message in context: http://www.nabble.com/Enclose-the-polygons-within-50-metres-tp17851388p17851388.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] dissolve and multipolygon(1) convert polygon to points(centroid) (2)
Hello all. I have two problems. 1) v.dissolve? it is possible to dissolve the polygons (nonadjacent) who have the same attribute obtaining a multipolygon? 2) v.type I noticed that to get points (representing the centroid) from polygons must change the order by hand ... type = point, centroid ..with.. type = centroid, point. From GUI as you can get. v.type [EMAIL PROTECTED] [EMAIL PROTECTED] type=point,centroid That exchange v.type [EMAIL PROTECTED] [EMAIL PROTECTED] type=centroid,point Thanks Gabriele ubuntu 7.10 -grass6.3.ORC6 -- View this message in context: http://www.nabble.com/dissolve-and-multipolygon%281%29--convert-polygon-to-points%28centroid%29--%282%29-tp16286855p16286855.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] center line with v.generalize ?
hamish_b wrote: Gabriele wrote: I have a theme of a polygonal rivers. I have to locate the center line as a line that identifies the river. Obviously the edges of the river are not parallel and can not use v.parallel. .. I would locate the midline of the river. So I would like to transform the polygon with the middle line. Maris: Just a quick idea - if vector map contains only rivers, You could try out v.to.rast and then r.thin till You get something similar to river and then convert back to vectors with r.to.vect. Gabriele: I have already tried with r.thin etc. .. but sometimes I need a little too pixels (because the rivers is very narrow and contorted) to avoid errors of approximation. The problem is that choosing a region with high resolution happens ' out of memory error '. I have tried to solve this problem in a way similar to Maris's suggestion. I was looking at fjords not rivers so my width was 500-1000m and I could use the method at 10m resolution without region size problems. For a long-thin river you might be able to do the processing piece-by- piece in a moving region window then patch all the center line parts together. I am more interested in river-width than river mile, I had hoped to create a center line then v.to.rast that line and for each cell make a line normal to the center-line and measure its distance (ie shore to shore distance normal to center-line), and output a profile cross section (2 * r.transect + bathy DEM). River width is interesting for things like atmosphere-water surface coupling + heat/gas exchange, photo-reactive dissolved organic chemistry exposure to sunlight, etc. -- a thinner section of the river/fjord will expose less top-1m volume and (depending on depth changes) may speed up the flow further reducing available reaction time*spatial exposure. Anyway I experimented with a few things, what I ended with AFAICR was creating a land raster MASK with v.to.rast and the coastline, then running r.cost to find distance to nearest shore. I then ran a combination of r.slope.aspect + r.mapcalc slope5 and 'r.param.scale param=feature' to look for ridges in the cost map. Then r.thin + r.to.vect. The result was a nice start, but there were some problems I still don't know the answers to. For one thing the presence of islands in the channel split the distance in two. I guess you could use v.extract and/or v.dissolve to remove all islands before r.cost, then subtract island width from the r.transect width later on? A second problem was what to do when you came to a place where the channel forked. here is an interesting link: Dynamic Segmentation and Thiessen Polygons: A Solution to the River Mile Problem [using Thiessen polygons] by William W. Hargrove, Richard F. Winterfield, Daniel A. Levine: http://research.esd.ornl.gov/CRERP/DOCS/RIVERMI/P114.HTM here are some old grass mailing list links, I though there was something more recent, but don't see it now. the gmane threading seems a bit broken, may nabble does better? http://thread.gmane.org/gmane.comp.gis.grass.user/12308 http://article.gmane.org/gmane.comp.gis.grass.user/12332 http://article.gmane.org/gmane.comp.gis.grass.user/12346 http://thread.gmane.org/gmane.comp.gis.grass.user/12352 http://img93.imageshack.us/my.php?image=schermorast73iv.png be careful as river-mile can be a fractal problem. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user Your information is very interesting. I will try to take inspiration from all of the advice that I read. Many thanks Gabriele -- View this message in context: http://www.nabble.com/center-line-with-v.generalize---tp15063073p15084050.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] center line with v.generalize ?.
Hi list -- View this message in context: http://www.nabble.com/center-line-with-v.generalize--.-tp15062970p15062970.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] center line with v.generalize ?
Hi list. I have a theme of a polygonal rivers. I have to locate the center line as a line that identifies the river. Obviously the edges of the river are not parallel and can not use v.parallel. I am making attempts with v.generalize but unable to solve. Suggestions? Thanks Gabriele -- View this message in context: http://www.nabble.com/center-line-with-v.generalize---tp15063073p15063073.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] center line with v.generalize ?
I would locate the midline of the river. So I would like to transform the polygon with the middle line. Thanks -- View this message in context: http://www.nabble.com/center-line-with-v.generalize---tp15063073p15063622.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] v.coordinate.sh - improve this script?
Moritz Lennert-2 wrote: On 24/01/08 14:12, Gabriele N. wrote: Hi list. I made a script (attached) for the definition of interactive text file (containing the coordinates input and output) to be loaded later in v.transform. Could you explain how this is different from the existing File-Georectify tool in the gis.m GUI ? Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user I do not know well the File- Georectify tool in the gis.m GUI and are not sure if I can even assign coordinates output with the click. However, I now better study (recently I have switched to GRASS), and if I was wrong... I apologize. Thanks -- View this message in context: http://www.nabble.com/v.coordinate.sh---improve-this-script--tp15063748p15065908.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