[GRASS-user] How we run georectify in command line
Dear all, Any one can tell me the the problem like- georectification in command line Thanks Regards, Swapan Ghosh ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How we run georectify in command line
On 29/08/2012 10:25, SWAPAN GHOSH wrote: Dear all, Any one can tell me the the problem like- "georectification in command line" You can try gdalwarp if you need to reproject, or gdaltransform if you need to just apply a shift with ground control points. Thanks Regards, Swapan Ghosh This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver 052-3665918 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How we run georectify in command line
Actually, I would like to open the Setup for georectification window, how I can do it by writing command in console. Thanks Regards, Swapan On Wed, Aug 29, 2012 at 1:10 PM, Micha Silver mi...@arava.co.il wrote: On 29/08/2012 10:25, SWAPAN GHOSH wrote: Dear all, Any one can tell me the the problem like- georectification in command line You can try gdalwarp if you need to reproject, or gdaltransform if you need to just apply a shift with ground control points. Thanks Regards, Swapan Ghosh This mail was received via Mail-SeCure System. ___ grass-user mailing listgrass-user@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver 052-3665918 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How we run georectify in command line
On Wed, Aug 29, 2012 at 9:25 AM, SWAPAN GHOSH swap.g...@gmail.com wrote: Dear all, Any one can tell me the the problem like- georectification in command line Please see here for the commands: http://grass.osgeo.org/wiki/Image_processing#Geometric_preprocessing.2FGeorectification Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Adjust opacity in command line
Andy wrote: I have been trying to find a way to adjust the opacity (or transparency, to throw another key word in there) of raster and vector layers using the command line. for two rasters you can use the d.shadedmap module for a similar effect. It's just a wrapper around a trick with the d.his module. examples: http://grass.osgeo.org/grass63/screenshots/images/fourmaps.jpg http://grass.osgeo.org/grass63/screenshots/images/spearfish_catchements.png (see raster screenshots on the homepage for scripts that made those) It's a bit limited in what it can be used for, but for the case of a categorical area map (maybe 'v.to.rast type=area') over a shaded relief DEM it works quite nicely. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] datum vs. sphere (grib files)
On Fri, Aug 24, 2012 at 2:19 PM, Margherita Di Leo direg...@gmail.com wrote: Hi, I am currently using some grib files for meteorological data. When I set a location referring to a georeferenced grib file, a default WGS84 is assigned to it, because datum is unknown datum: ** unknown (default: WGS84) ** Actually, the grib files are referred to a sphere, ... could you post the gdalinfo output? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Vector import combo: db.in.ogr/v.in.db data type casting problem
On Mon, Aug 20, 2012 at 3:56 PM, Peter Löwe peter.lo...@gmx.de wrote: hi, here is an issue (probably float/integer typcasting related) issue when importing a large CSV dataset of vector points into GRASS6.4.2. The data set contains hundreds of columns, Note that the default DBF driver does not support an unlimited amount of columns (the specs suggest 128 columns max). so v.in.ascii is not an option (the columns option-string would be ungodly long). Instead, the data is transformed into a DBF Please use SQLite instead, also as DB backend for such operations (ideally, due to recent fixes, the current 6.4.svn). Best, Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Bathymetric data
On Fri, Aug 24, 2012 at 11:12 AM, Peter Löwe peter.lo...@gmx.de wrote: ... The GEBCO data is probably what you are looking for. Google uses it as well: http://www.gebco.net/data_and_products/gridded_bathymetry_data/ Related hints are in the Wiki: http://grass.osgeo.org/wiki/Global_datasets#GEBCO_Bathymetric_Chart Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Vector import combo: db.in.ogr/v.in.db data type casting problem
Hi Markus, Note that the default DBF driver does not support an unlimited amount of columns (the specs suggest 128 columns max). Please use SQLite instead, also as DB backend for such operations (ideally, due to recent fixes, the current 6.4.svn). That's correct. Thanks for following up on this. Peter -- Dr. Peter Löwe peter.lo...@gmx.de ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.fuzzy.system output definition
Hi GRASS users! Maybe anyone has already tried the new GRASS GIS fuzzy inference system (r.fuzzy.system). I try to do my first fuzzy computation and follow the paper Jasiewicz, J. 2011: A new GRASS GIS fuzzy inference system for massive data analysis, Computers Geosciences, 37: 9, 1525-1531. Unfortunately I get an error saying it can't find the output raster map which is indicated in the .map file: r.fuzzy.system maps=path/to/my/test_set.map rules=path/to/my/test.rul family=Zadeh defuz=bisector imp=minimum res=100 output=output ERROR: Raster map output not found my test_set.map file contains beside the other sets also %output $ none {both; 0,20,20,40; linear; 0;1} $ low {both; 20,40,40,60; linear; 0;1} $ medium {both; 40,60,60,80; linear; 0;1} $ high {both; 60,80,80,100; linear; 0;1} The problem is related to the %output in the .map file not the setting for output=output because there is also an error for following case: r.fuzzy.system maps=path/to/my/test_set.map rules=path/to/my/test.rul family=Zadeh defuz=bisector imp=minimum res=100 output=output_ map ERROR: Raster map output not found I am using the GUI interface for the module in GRASS GIS 6.5 on Ubuntu 12.04. Any suggestions what I am doing wrong? Thank you best regards, Johannes ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.stream.order export to shapefile
Hi Markus, Thanks for responding. I still can't seem to get this to work. Here are my steps: r.stream.extract elev=dtm thresh=50 stream_rast=stream_500 stream_vect=stream_500 dir=fdir_500 r.stream.order stream=stream_500 dir=fdir_500 table=stream_order v.db.connect map=stream_500 driv=sqlite table=stream_order key=cat layer=2 Here I get: The table stream_order is now part of vector map stream_500 and may be deleted or overwritten by GRASS modules DBMI-SQLite driver error: Cannot create index: create unique index stream_order_cat on stream_order ( cat ) index stream_order_cat already exists WARNING: Cannot create index Select privileges were granted on the table Next: v.category stream_500 opt=report layer=2 Layer/table: 1/stream_500 type countminmax point248 1237 line 237 1237 boundary 0 0 0 centroid 0 0 0 area 0 0 0 all 485 1237 Layer/table: 2/stream_order type countminmax point248 0 2 line 237 0 1 boundary 0 0 0 centroid 0 0 0 area 0 0 0 all 485 0 2 Why in layer 2 are there only 2 cat values? Shouldn't I get all the cats as in layer 1 when I use key=cat in the v.db.connect?? Do I need to delete and recreate the cats in layer 2 ? And finally, the export: GRASS 6.4.2 (ITM):~/GIS/DEM/LIDAR_EinYahav v.out.ogr -c -e stream_500 type=line dsn=stream_500.shp layer=2 WARNING: 248 point(s) found, but not requested to be exported. Verify 'type' parameter. Exporting 485 geometries... WARNING: 124 features without attributes were written v.out.ogr complete. 237 features written to stream_500 (ESRI_Shapefile). Why am I getting 485 geometries? and why are 124 with no attributes? Many thanks, Micha On 08/28/2012 11:15 PM, Markus Metz wrote: On Tue, Aug 28, 2012 at 9:05 PM, Micha Silvermi...@arava.co.il wrote: A few more details regarding v.out.ogr when layer=2 I see this ticket, possibly similar: http://trac.osgeo.org/grass/ticket/991 I tried to export to PostGIS, and again, all attribute fields are created but all values are NULL. I also tried with the layer 2 database as a dbf file, instead of sqlite. Same result. This was all with GRASS 6.4.2 on scientific linux 6 What does v.category op=report say? Are there any categories in layer 2? If not, there is nothing to export. Markus M Thanks, Micha On 08/28/2012 02:57 PM, Micha Silver wrote: I'm having two problems when exporting the output of r.stream.order to a shapefile. I have attached the created table (saved in my setup in sqlite) to the 'streams' map thru layer 2. Then v.db.select streams layer=2 shows all the details of strahler order, prev_str, etc. First problem: When I try to do v.out.ogr, the column header named next_stream fails because it's 10 characters long, too long for a shapefile dbf. As a result all attributes are missing. Here's a patch to io.c in r.stream.order I put in place to correct this: [micha@SL6 r.stream.order]$ diff -u io.c.orig io.c --- io.c.orig2012-08-28 12:21:46.020275045 +0300 +++ io.c2012-08-28 12:22:17.951382882 +0300 @@ -292,7 +292,7 @@ /* table definition */ char *tab_cat_col_name = cat integer; char *tab_stream = stream integer; -char *tab_next_stream = next_stream integer; +char *tab_next_stream = next_str integer; char *tab_prev_streams; char *tab_strahler = strahler integer; char *tab_horton = horton integer; @@ -300,7 +300,7 @@ char *tab_hack = hack integer; char *tab_length = length double precision; char *tab_cumlength = cum_length double precision; -char *tab_stright = stright double precision; +char *tab_stright = straight double precision; char *tab_fractal = fractal double precision; char *tab_distance = out_dist double precision; char *tab_topo_dim = topo_dim integer; (Also corrects a small typo in the column name stright) After recompiling, then rerunning the r.stream.order addon, the sqlite table now has all column headers= 10 characters suitable to shapefile. Problem 2, where I'm stuck. I export the vector map using : v.out.ogr streams layer=2 dsn=streams.shp The process finishes OK, and the shape is created with all attrib columns, but *all* values are NULL (including the cat). Any ideas what I've done wrong? Thanks, Micha This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver GIS Consultant, Arava Development Co. http://www.surfaces.co.il ___ grass-user mailing
[GRASS-user] Fwd: Re: R + grass7
fyi a question regarding an translation (flags,parameter) table between grass6x and grass7.Helmut Original-Nachricht Datum: Wed, 29 Aug 2012 14:35:37 +0200 (CEST)Von: Roger Bivand An: Helmut Kudrnovsky Betreff: Re: R + grass7 On Wed, 29 Aug 2012, Helmut Kudrnovsky wrote: Dear Roger, Im trying to improve the GRASS GIS-R-coupling on windows (http://trac.osgeo.org/grass/ticket/1149). May I ask you if spgrass6 (http://cran.at.r-project.org/web/packages/spgrass6/index.html) is GRASS GIS 7-aware. Yes, the description on: http://cran.r-project.org/web/packages/spgrass6/index.html uses 6+ to signal that subsequent versions are intended to work. Since the package only uses system(), it will work as long as --interface-description returns an XML file. Plenty of examples and scripts stop working because GRASS 6 and 7 modules have different flags and/or parameter names, but in principle only for that reason. Its hard for me to check - it would be useful to have a translation table for flag and parameter names by module where changes have occurred. Best wishes, Roger -- Roger Bivand Department of Economics, NHH Norwegian School of Economics, Helleveien 30, N-5045 Bergen, Norway. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] [GRASS-dev] Fwd: Re: R + grass7
On Wed, Aug 29, 2012 at 2:50 PM, Helmut Kudrnovsky hel...@web.de wrote: fyi a question regarding an translation (flags,parameter) table between grass6x and grass7. We are maintaining a list of changes here: http://trac.osgeo.org/grass/wiki/Grass7/NewFeatures#Renamedoptions The idea to be least invasive while gaining maximal consistency. I hope that only a few scripts/examples are failing after migrating from GRASS 6 to 7. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Fwd: Re: R + grass7
Hi, You can find a list of changes in option settings and names here: http://trac.osgeo.org/grass/wiki/Grass7/NewFeatures#Renamedoptions Cheers, Paulo On 08/29/2012 02:50 PM, Helmut Kudrnovsky wrote: fyi a question regarding an translation (flags,parameter) table between grass6x and grass7. Helmut Original-Nachricht Datum: Wed, 29 Aug 2012 14:35:37 +0200 (CEST) Von: Roger Bivand An: Helmut Kudrnovsky Betreff: Re: R + grass7 On Wed, 29 Aug 2012, Helmut Kudrnovsky wrote: Dear Roger, I'm trying to improve the GRASS GIS-R-coupling on windows (http://trac.osgeo.org/grass/ticket/1149). https://service.gmx.net/de/cgi/derefer?TYPE=3DEST=http%3A%2F%2Ftrac.osgeo.org%2Fgrass%2Fticket%2F1149%29. May I ask you if spgrass6 (http://cran.at.r-project.org/web/packages/spgrass6/index.html) https://service.gmx.net/de/cgi/derefer?TYPE=3DEST=http%3A%2F%2Fcran.at.r-project.org%2Fweb%2Fpackages%2Fspgrass6%2Findex.html%29 is GRASS GIS 7-aware. Yes, the description on: http://cran.r-project.org/web/packages/spgrass6/index.html https://service.gmx.net/de/cgi/derefer?TYPE=3DEST=http%3A%2F%2Fcran.r-project.org%2Fweb%2Fpackages%2Fspgrass6%2Findex.html uses 6+ to signal that subsequent versions are intended to work. Since the package only uses system(), it will work as long as --interface-description returns an XML file. Plenty of examples and scripts stop working because GRASS 6 and 7 modules have different flags and/or parameter names, but in principle only for that reason. It's hard for me to check - it would be useful to have a translation table for flag and parameter names by module where changes have occurred. Best wishes, Roger -- Roger Bivand Department of Economics, NHH Norwegian School of Economics, Helleveien 30, N-5045 Bergen, Norway. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Script in Grass - Mapcalc and MODIS
Hi, I'm quit new using GRASS and I'm discovering the great possibilities of the program. I'm currently working with rainfall data and I what to add the daily data of a year using r.mapcalc. I'm sure there a quick solution making a script but I'm really beginner in script. Does someone can give me some tutorial or indication? What would be the easier? Bash or python? Does also someone has some experience with MODIS images in hdf format? Is it possible to work with it with GRASS? Thank for your support Lucien - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Lucien Blandenier, PhD student Center for hydrogeology and geothermy University of Neuchâtel Rue Emile-Argand 11 CH-2000 Neuchâtel Switzerland Tel : +41 32 718 26 38 lucien.blanden...@unine.ch [Description : Description :]Please consider the environment before printing this email. inline: image002.png___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Script in Grass - Mapcalc and MODIS
Hi Lucien. Does also someone has some experience with MODIS images in hdf format? Is it possible to work with it with GRASS? A few links: http://grass.osgeo.org/wiki/MODIS http://www.safe.iis.u-tokyo.ac.jp/pdf/MODIS_Image_Processing_using_GRASS_ver1.3.pdf ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user