[GRASS-user] doubt about picture link
Is there any solution to give link to picture in databse -- vipin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] opacity through command line
Hi everybody, how can I set the opacity of a raster map through command line (not through GUI interface)? thanks as usual, -- Alberto Pettazzi MeteoGalicia - Departamento de Climatologa y Observacin Consellera de Medio Ambiente e Desenvolvemento Sostible rea Central, Local: L-31-C. Polgono de Fontias 15703 Santiago de Compostela. A Corua Telfono: +34-981-957463 e-mail: [EMAIL PROTECTED] ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] How to write scripts in GRASS
I´d very grateful if anyone can point any tutorial on how to write scripts in Python to perform routine tasks. I remember to have seen anyone at GRASS web site, but I can´t find it. Thanks. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to write scripts in GRASS
http://grass.osgeo.org/wiki/GRASS_and_Python ...maybe your starting point. Regards, Christian. Zitat von Gabriel Messner [EMAIL PROTECTED]: I´d very grateful if anyone can point any tutorial on how to write scripts in Python to perform routine tasks. I remember to have seen anyone at GRASS web site, but I can´t find it. Thanks. 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
Re: [GRASS-user] Different GIS_ERROR_LOG files for different instances of GRASS?
On Wed, Oct 8, 2008 at 12:53 AM, Glynn Clements [EMAIL PROTECTED] wrote: Markus Neteler wrote: I am running simulations which are using GRASS for some calculations, and I would like to use GIS_ERROR_LOG as a log file. How can I have different GIS_ERROR_LOG for different instances of GRASS, so that each instance has it's own log file? Best might be to change the name GIS_ERROR_LOG which is defined in lib/gis/error.c into a variable. ? Done in 7.0 (r33740). Backported to 6.4.svn: r33753 If set, $GIS_ERROR_LOG should be the absolute path to the log file (a relative path will be interpreted relative to the process' cwd, not the cwd at the point you set the variable). If not set, $HOME/GIS_ERROR_LOG is used instead. ok The previous behaviour of writing to both $HOME/GIS_ERROR_LOG and $GISBASE/GIS_ERROR_LOG has been eliminated; only one log file will be written. That's not entirely clear to me. Say, lib/init/variables.html may need a fix for that. As before, the file will only be used if it already exists. Also, the current directory is obtained using getcwd() rather than popen(pwd), so it should work on Windows. Perfect, thanks. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Export Centroids
On 07/10/08 17:48, gridcell wrote: Is there a way to export centroids (v.out.ogr) from a GRASS layer? I have a layer that was originally overlapping polygons in a shapefile. I imported the shapefile into GRASS and build/clean geometry, resulting in 4 polygons. I want to export the centroids of these polygons, and was hoping to have returned the polygon's attributes and a point geometry representing the centroid (4 features). I tried the following but the command never get any results v.out.ogr [EMAIL PROTECTED] type=point,centroid 'dsn=PG: host=localhost user=postgres dbname=test' olayer=blk layer=1 format=PostgreSQL I get the following in my output... 6 Boundary(ies) found, but not requested to be exported. Verify 'type' parameter. Exporting 9 points/lines... ERROR 1: INSERT command for new feature failed. ERROR: geometry requires more points Command: INSERT INTO blk (wkb_geometry , cat, comments) VALUES (GeomFromEWKT('SRID=32768;LINESTRING (1361422.480504325358197 586604.790284495218657)'::TEXT) , 1, 'Polygon 1') ERROR 1: INSERT command for new feature failed. ERROR: geometry requires more points This seems to be an issue with ogr's conversion to PostGIS. You could try exporting to a shapefile and importing that to PostGIS... Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] How do you list/remobe subgroups of a group?
Why can't we list subgroups of a specific group? Even more important, I think, is that there is no option to remove subgroups. Is it ok to remove subgroups manually? (e.g. rm -r ~/grassdb/project/location/mapset/group/somegroup/somesubgroup) Is there a reason for not having these options? Thank you, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Using r.mapcalc to transform '*' values to a float
Glynn, Thank you very much for the help; I've tried the methods you suggested and they work fine. The latter, using r.out.ascii, is what makes the most sense for me. Thanks again, Tom Glynn Clements wrote: Thomas Adams wrote: I have a GRASS location that spans the entire U.S. I plan on re-projecting raster maps into the location from a GRASS location with a much smaller domain. When I do this, the values outside of the projected map region have a value '*' when I query the new raster outside of my area of interest using the GRASS GUI. What I need to do, ultimately, is to use r.out.ascii to write out an ascii grid that covers the full U.S. region regardless of the extent of the actual data. However, before I do this I want to use r.mapcalc to convert the '*' values (are they 'null'?) to 9.999e+20, which is a the missing value for the software that needs to ingest the ascii file output from GRASS. How do I treat the '*' values in r.mapcalc? The '*' values are nulls. You can convert them with r.mapcalc with: r.mapcalc 'newmap = if(isnull(oldmap),9.999e+20,oldmap)' Or you can modify the existing map in place with: r.null themap null=9.999e+20 Or you can leave them as null but replace them in the export process with: r.out.ascii themap null=9.999e+20 -- Thomas E Adams National Weather Service Ohio River Forecast Center 1901 South State Route 134 Wilmington, OH 45177 EMAIL: [EMAIL PROTECTED] VOICE: 937-383-0528 FAX:937-383-0033 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How do you list/remobe subgroups of a group?
On Wed, Oct 8, 2008 at 2:22 PM, Nikos Alexandris [EMAIL PROTECTED] wrote: Why can't we list subgroups of a specific group? Even more important, I think, is that there is no option to remove subgroups. Is it ok to remove subgroups manually? (e.g. rm -r ~/grassdb/project/location/mapset/group/somegroup/somesubgroup) Is there a reason for not having these options? Are you using an oldish GRASS version? http://grass.osgeo.org/grass64/manuals/html64_user/i.group.html i.group [-rlg] group=name [subgroup=string] [input=name[,name,...]] [--verbose] [--quiet] Flags: -r Remove selected files from specified group -l List files from specified (sub)group (fancy) -g List files from specified (sub)group (shell script style) Or does it *fail*? Best Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Export Centroids
I've already tried to shapefile. It builds a shapefile and tries to add linestrings. I'm guessing its a bug in v.out.ogr. I think I have a temporary work around until this is fixed. I use v.type centroid to point and then use v.out.ogr. Seems to work but I'm still testing. gridcell wrote: Is there a way to export centroids (v.out.ogr) from a GRASS layer? I have a layer that was originally overlapping polygons in a shapefile. I imported the shapefile into GRASS and build/clean geometry, resulting in 4 polygons. I want to export the centroids of these polygons, and was hoping to have returned the polygon's attributes and a point geometry representing the centroid (4 features). I tried the following but the command never get any results v.out.ogr [EMAIL PROTECTED] type=point,centroid 'dsn=PG: host=localhost user=postgres dbname=test' olayer=blk layer=1 format=PostgreSQL I get the following in my output... 6 Boundary(ies) found, but not requested to be exported. Verify 'type' parameter. Exporting 9 points/lines... ERROR 1: INSERT command for new feature failed. ERROR: geometry requires more points Command: INSERT INTO blk (wkb_geometry , cat, comments) VALUES (GeomFromEWKT('SRID=32768;LINESTRING (1361422.480504325358197 586604.790284495218657)'::TEXT) , 1, 'Polygon 1') ERROR 1: INSERT command for new feature failed. ERROR: geometry requires more points Command: INSERT INTO blk (wkb_geometry , cat, comments) VALUES (GeomFromEWKT('SRID=32768;LINESTRING (1361632.899105372838676 586307.90450647752732)'::TEXT) , 2, 'Polygon 2') ERROR 1: INSERT command for new feature failed. ERROR: geometry requires more points Command: INSERT INTO blk (wkb_geometry , cat, comments) VALUES (GeomFromEWKT('SRID=32768;LINESTRING (1361538.573653233936056 586446.606141098076478)'::TEXT) , 1, 'Polygon 1') ERROR 1: INSERT command for new feature failed. ERROR: geometry requires more points Command: INSERT INTO blk (wkb_geometry , cat, comments) VALUES (GeomFromEWKT('SRID=32768;LINESTRING (1361538.573653233936056 586446.606141098076478)'::TEXT) , 2, 'Polygon 2') 4 features written http://www.nabble.com/file/p19861164/blocks.png -- View this message in context: http://www.nabble.com/Export-Centroids-tp19861164p19879641.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] Different GIS_ERROR_LOG files for different instances of GRASS?
These are more and more reasons, to use 6.4 instead. Thanks a million for all your suggestions, Rainer On Wed, Oct 8, 2008 at 10:34 AM, Markus Neteler [EMAIL PROTECTED] wrote: On Wed, Oct 8, 2008 at 12:53 AM, Glynn Clements [EMAIL PROTECTED] wrote: Markus Neteler wrote: I am running simulations which are using GRASS for some calculations, and I would like to use GIS_ERROR_LOG as a log file. How can I have different GIS_ERROR_LOG for different instances of GRASS, so that each instance has it's own log file? Best might be to change the name GIS_ERROR_LOG which is defined in lib/gis/error.c into a variable. ? Done in 7.0 (r33740). Backported to 6.4.svn: r33753 If set, $GIS_ERROR_LOG should be the absolute path to the log file (a relative path will be interpreted relative to the process' cwd, not the cwd at the point you set the variable). If not set, $HOME/GIS_ERROR_LOG is used instead. ok The previous behaviour of writing to both $HOME/GIS_ERROR_LOG and $GISBASE/GIS_ERROR_LOG has been eliminated; only one log file will be written. That's not entirely clear to me. Say, lib/init/variables.html may need a fix for that. As before, the file will only be used if it already exists. Also, the current directory is obtained using getcwd() rather than popen(pwd), so it should work on Windows. Perfect, thanks. Markus -- Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT), Dipl. Phys. (Germany) Centre of Excellence for Invasion Biology Faculty of Science Natural Sciences Building Private Bag X1 University of Stellenbosch Matieland 7602 South Africa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Runing grass inside Msys console
Hi, I just installed WinGrass (great job by the way) but in order to run GRASS inside the MSYS console I have to explicitly go to the directory were grass63 is at. IOW, the grass63 startup script is not in the path environment. Is it possible to put the startup script in the path environment, along with other utilities like gdalinfo? Cheers Daniel ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Problem with the signatures from i.class
Hi, i'm using i.class to generate a signature to run on i.maxlik. The module read the signature with no problem, but, the class map result shows only one class and the area that its cover is bigger than the region that i've set in i.class. Thats a large red rectangle. I tried various standard deviations but the result is always the same. what can be happening? Sorry, i'm brazilian and my english is not very good! thanks Thybério Luna Freire ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] flip raster
On Tue, 2008-10-07 at 01:37 +0100, Glynn Clements wrote: José María Michia wrote: I've imported a NetCDF file (ETOPO1 model). The resulted raster appears flipped vertically. How can I fix this? I try with mapcalc, without success: - neighborhood modifier: not accept computations in offset parameter, like map[0,total_rows-row()] - I dont know how to query map in arbitrary coordinates, like map(x(),-1*y()) There isn't any mechanism in r.mapcalc to achieve this. I'm not sure that it's possible to do it entirely within GRASS. Flipping geo-referenced data doesn't really make much sense. Either figure out how to import it with the correct orientation, or export it, flip the exported data, then re-import. Glynn and All, excuse me for hijacking the post. I am looking for a way (just for the fun of the game or for philosophical re-search) to rotate geotiffs at 180 degrees (i.e. flip vertically and horizontally). I've posted about this in the gdal-dev list [1]. Is there a way to accomplish a 180 deg. rotation using the listgeo/geotifcp command line tools? Thank you, Nikos [1] http://lists.osgeo.org/pipermail/gdal-dev/2008-October/018548.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] flip raster
2008/10/8 Nikos Alexandris [EMAIL PROTECTED]: On Tue, 2008-10-07 at 01:37 +0100, Glynn Clements wrote: José María Michia wrote: I've imported a NetCDF file (ETOPO1 model). The resulted raster appears flipped vertically. How can I fix this? I try with mapcalc, without success: - neighborhood modifier: not accept computations in offset parameter, like map[0,total_rows-row()] - I dont know how to query map in arbitrary coordinates, like map(x(),-1*y()) There isn't any mechanism in r.mapcalc to achieve this. I'm not sure that it's possible to do it entirely within GRASS. Flipping geo-referenced data doesn't really make much sense. Either figure out how to import it with the correct orientation, or export it, flip the exported data, then re-import. Glynn and All, excuse me for hijacking the post. I am looking for a way (just for the fun of the game or for philosophical re-search) to rotate geotiffs at 180 degrees (i.e. flip vertically and horizontally). I've posted about this in the gdal-dev list [1]. Is there a way to accomplish a 180 deg. rotation using the listgeo/geotifcp command line tools? I'm not sure, but maybe this can be useful: Alexander Schulze suggested me that I use the library to invest CDO latitudes: http://www.mpimet.mpg.de/fileadmin/software/cdo/ Something like: cdo invertlat filein fileout1 cdo invertlon fileout1 fileout2 If you can use netCDF format, look at this: http://nco.wiki.sourceforge.net/reverse a dimension Saludos José María Thank you, Nikos [1] http://lists.osgeo.org/pipermail/gdal-dev/2008-October/018548.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] doubt about picture link
Hi, well, you can always convert your database from DBF format to PostgreSQL and use their BLOB type. Then you can store the images directly to database. If this is not what you are looking for, just add new column of type STRING (varchar) to the database: v.db.addcol map=your_vector column=image varchar(256) then you can store the link directly. jachym 2008/10/8 vipin poulose [EMAIL PROTECTED]: Is there any solution to give link to picture in databse -- vipin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Jachym Cepicky e-mail: jachym.cepicky gmail com URL: http://les-ejk.cz GPG: http://www.les-ejk.cz/pgp/jachym_cepicky-gpg.pub ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] flip raster
Nikos: Is there a way to accomplish a 180 deg. rotation using the listgeo/geotifcp command line tools? try saving meta data to a file as detailed in the libGeotiff FAQ, then tifftopnm | pnmflip | pnmtotiff and reattaching metadata as shown in the FAQ entry. http://www.remotesensing.org/geotiff/faq.html#preserve_metadata Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Using GRASS with SPOT 5 Images
I've been searching the network for a GRASS tutorial/Example on how to use SPOT 5 images in grass. Basically I have a GeoTIFF images and DTED data. I need to create a scenery for a flight simulator wuith uses a Lat/long projection. Can anyone assists me in this regard. Linda Malinga Software Engineer. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user