Re: [GRASS-user] Reclass: Atributte a string
Franz Schiller wrote: > I'm doing a reclass of rasters based on a rule file. In this file I'm doing > a correspondence between values (1,2... 100) to other integer values. > Is it possible to do this but instead to done this to Numbers do this to > Strings just like > 1 = "Agriculture" > 2,3,4= "Forest" You can associate category labels with each category using r.categories. You can't "reclass" to a string. -- Glynn Clements ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] RE: Need help with volumes in NVIZ
Craig, Your instructions gave me success. I think my problem was not putting in a value that was recognized in the data. I now can see isosurfaces. Thanks! Suzanne On 11/11/2010 9:29 AM, Craig Wilson wrote: Suzanne, If you've got to the stage of importing the G3D (volume) file into NVIZ, and if it has the variable (w) attribute running through the volume, isosurfaces can be created by: * Selecting 'Isosurfaces' from the 'Visualization type'; * Clicking 'Add' from the right-hand menu; * Where it prompts you for a digit, enter a value of the variable which is contained within the volume. For example, when I wanted to highlight areas within Loch Lomond which had 0.15mM of dissolved carbon, I entered 15 at this stage (I had multipled all values by 100 at an earlier stage). As long as this number was contained within the original data converted into the G3D file (e.g. if a text file contains 4 columns of x, y, z and the variable attribute, this attribute must be selected when the G3D file is created - from what you said about getting the slices to work, it sounds like this has already been done), the interpolated boundary lines for 0.15mM throughout the lake were then visualized with isosurfaces. In the screen grab below, I performed this for 0.15mM and 0.20mM. _http://img222.imageshack.us/img222/8780/isosurfacesinnviz.jpg_ It's been a few months since I used GRASS, so if those instructions don't work I can go back and figure out precisely how it was done. Craig Wilson ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Using v.colors in WinGRASS
Glynn wrote: > How do people feel about replacing shell > scripts with Python for 6.5? extremely negative. I'm totally against that. please lets just let the stable branch become stable. I don't want to see that change in 6.4.x, so no point to make 6.5 deviate from the stable branch any more than it already has. IMHO 6.5 is only for testing backports and fixes before actually commiting the backport to 6.4. I do not believe it's a fundamental flaw which can not be worked around in the current framework- calling a shell script from another shell script on WinGRASS does work- e.g. r.in.wms calls r.tileset (also a bunch of other private helper shell scripts). [that r.in.wms has other issues is not very relevant to this point] and so I reject the idea to fix a fixable bug by replacing a huge work of well tested code by much less tested replacements. I would like to declare 6.x closed for new major changes ASAP, if that is not already the case. New devel should happen in 7. thanks, Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] raster data output
Take a look at r.what manual page. There you can supply a list of coordinates and the rasters you want to query to get the desired table. On Wed, Nov 10, 2010 at 3:03 PM, Achim Kisseler wrote: > Did you see: > > http://grass.fbk.eu/statsgrass/grass_geostats.html > ? > > Am 10.11.2010 16:23, schrieb Steve Wangen: >> >> Hello list - >> There may be a simple solution to my problem, but I'm a bit new to >> grass and haven't been able to track it down yet, so I was wondering >> if anyone had suggestions. I'm trying to output the cell contents of >> multiple rasters into a single table that I can save to a file with >> the desired result looking like this: >> >> x-coord y-coord raster1 raster2 raster3 ... >> 0 0 1 2 3 ... >> 0 1 0 1 2 ... >> 0 2 0 0 1 ... >> etc... >> >> with the idea of doing some time series analyses, where the subsequent >> rasters are different time steps. Ideally this would come out as one >> table that I could then import into R for analysis. I've had a bit of >> luck outputing the result of r.stats -x to a file using some bash >> scripting and then combining all those files manually, but it seems >> there would be an easier solution. >> >> Any suggestions would be appreciated. >> >> Thanks! >> >> -Steve >> ___ >> grass-user mailing list >> grass-user@lists.osgeo.org >> http://lists.osgeo.org/mailman/listinfo/grass-user > > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.dissolve for lines?
> On 11/05/2010 09:59 PM, Bryan Keith wrote: >>> On 10/19/2010 11:17 PM, Bryan Keith wrote: >>> Hello, I have a vector line file that has lots of two-vertex lines. Many of these two-vertex lines share endpoints and really just make parts of one longer line. I'd like to dissolve these lines into one line when they have the same end point. It looks like v.dissolve only works for areas. Is there somewhere else I should be looking? Thank you. >>> Did you have a look at v.build.polylines ? >>> Additional note: line vectors have a direction. Two line segments with >>> an equivalent node (end point) will *not* connect if they are in >>> opposite directions. In this case you can try the v.edit tool=flip >>> option (on those segments only) >>> >> Thank you for the response. I finally got around to looking at this >> today. I made a simple example and tested with v.build.polylines, and >> it >> works great. Thank you. >> >> However, it turns out that I have many duplicate lines in my dataset. I >> created a simple example where there are duplicate lines, and >> v.build.polylines doesn't work in that case. What do I need to look at >> to >> get rid of duplicate lines? >> >> > > v.clean tool=rmdupl > Might do it if the lines exactly overlap. Yes, that worked in my case. Thanks so much, Micha! Bryan ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Grass freeze in v.digit after some "display attributes"
Hi, I'm experiencieng a strange bug here: with tcl/tk editor or with qgis, after using the "display attributes" tool about 10 to 20 times, GRASS (or qgis) completely freeze.I've tried this:-Uninstal and install everything again.-Many different vectos.-Different databases (dbf, sqlite, pgsql). I have the OSGeo4w version running on Win7. with debug=1 I have these messages for every click with "display attributes". D1/1: Vect_get_field(): field = 1D1/1: Vect_get_dblink(): link = 0 Any clue? Regards. Pablo Torres Carreira ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] RE: Need help with volumes in NVIZ
Suzanne, If you've got to the stage of importing the G3D (volume) file into NVIZ, and if it has the variable (w) attribute running through the volume, isosurfaces can be created by: Selecting 'Isosurfaces' from the 'Visualization type'; Clicking 'Add' from the right-hand menu; Where it prompts you for a digit, enter a value of the variable which is contained within the volume. For example, when I wanted to highlight areas within Loch Lomond which had 0.15mM of dissolved carbon, I entered 15 at this stage (I had multipled all values by 100 at an earlier stage). As long as this number was contained within the original data converted into the G3D file (e.g. if a text file contains 4 columns of x, y, z and the variable attribute, this attribute must be selected when the G3D file is created - from what you said about getting the slices to work, it sounds like this has already been done), the interpolated boundary lines for 0.15mM throughout the lake were then visualized with isosurfaces. In the screen grab below, I performed this for 0.15mM and 0.20mM. http://img222.imageshack.us/img222/8780/isosurfacesinnviz.jpg It's been a few months since I used GRASS, so if those instructions don't work I can go back and figure out precisely how it was done. Craig Wilson ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Difficulties in displaying/metadata reclass maps
Hello all, I have produced a fw of reclass maps and, when I need to display them or even do r.info to them, It takes a lot of time. My question is: - Why is that (since it takes more time than the original file that was reclassified) - And do you advise to create a copy of this map in order to display/get metadata? Thanks Kat ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Reclass: Atributte a string
Greetings I'm doing a reclass of rasters based on a rule file. In this file I'm doing a correspondence between values (1,2... 100) to other integer values. Is it possible to do this but instead to done this to Numbers do this to Strings just like 1 = "Agriculture" 2,3,4= "Forest" ? Thanks Franz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Using v.colors in WinGRASS
2010/11/10 Hamish : >> try r44284, hopefully fixed. > > perhaps buggy: > > - CMDLINE="`basename $0`" > + CMDLINE="`basename "$0"`" > > > becomes "`basename " + $0 + "`" > > afair the " do not nest. > > maybe outer quotes are not needed? > (but regardless I'm fairly sure this is not the cause of the v.colors bug, > just a red herring) it fixes white spaces in the path, e.g. basename: extra operand `(x86)/GRASS-65-SVN/scripts/v.db.addcol' Try `basename --help' for more information. C:/Program Files (x86)/GRASS-65-SVN/scripts/v.db.addcol: line 116: unexpected EOF while looking for matching ``' C:/Program Files (x86)/GRASS-65-SVN/scripts/v.db.addcol: line 125: syntax error: unexpected end of file Martin -- Martin Landa * http://geo.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Using v.colors in WinGRASS
Hi, 2010/11/11 Glynn Clements : >> I have seen something similar with v.rast.stats in wingrass. It seems >> that in wingrass, scripts can not call other scripts, only C modules. >> This might have something to do with the way scripts are called >> through a .bat file. > > How do people feel about replacing shell scripts with Python for 6.5? personally I could imagine this replacement also in GRASS 6.4.1. It's quite big pain to maintain scripts in two languages or better to say there is not enough manpower to do that. Moreover we will have always problems with bash scripts on Windows. Martin -- Martin Landa * http://geo.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Exporting Vector data to SHP (v.out.ogr)- takes a lot of time
Hi, Am 11.11.2010 16:08, schrieb Luis Lisboa: | Type of Map: vector (level: 2) | | | | Number of points: 0 Number of areas: 213.680 | | Number of lines:0 Number of islands: 73.649 | | Number of boundaries: 348.770 Number of faces: 0 | | Number of centroids:213.666 Number of kernels:0 | ...that´ll be a HUGE vector file! So it´s normal, that it takes some minutes. Achim ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Exporting Vector data to SHP (v.out.ogr)- takes a lot of time
Greetings I have produced a vector map (with r.to.vector) with the following metadata: v.info test01 ++ | Layer: teste01 | | Mapset: Regional | | Location:Country| | Database:E:\v3 | | Title: | | Map scale: 1:1 | | Map format: native | | Name of creator: orbit | | Organization: | | Source date: Tue Nov 09 11:11:40 2010 | || | Type of Map: vector (level: 2) | | | | Number of points: 0 Number of areas: 213680 | | Number of lines:0 Number of islands:73649 | | Number of boundaries: 348770 Number of faces: 0 | | Number of centroids:213666 Number of kernels:0 | | | | Map is 3D: No | | Number of dblinks: 1 | | | | Projection: Universal Transverse Mercator (zone 0) | | N: 8832707.83027649S: 7028073.62772369 | | E: 1347458.42001915W:200792.5777576 | | | | Digitization threshold: 0 | | Comments: | | | ++ Now I'm trying to export it to SHAPEFILE but, after 10 minutes the processing bar is not even in 40%. What might be happening? I mean, is it some error with the vector or when I convert to vector (from raster)? Because, taking more than 10 minutes to export to vector it's not normal (I guess) and taking nearly 7 min to convert to vector also... Thanks for the reply Luis L ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] How to drape an aerial image on an DEM map so that it retains its color....
Hello friends I have converted the cotour data to DEM and now looking forward to drape an aerial image of the same area on the DEM that I have got. In order to do this I loaded the DEM map and then draped the imported aerial map (imported using r.in.gdal, but that single jpeg image(satimage.jpeg) got imported into 3 different maps namely satimage.red, satimage.green, satimage.blue) on the DEM. The output view that I get seems okay but its too dark (black and white), though the original aerial image was colourful which I had imported. thanks Nikhil Morajkar Junior Software Engineer, Department of Computer Science and Engineering, Geo-Spatial Information Science & Engineering (GISE), Advance Research Lab, Indian Institute of Technology, Bombay (IITB) Powai, Mumbai -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/How-to-drape-an-aerial-image-on-an-DEM-map-so-that-it-retains-its-color-tp5728128p5728128.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] Using v.colors in WinGRASS
Hi as I said in a previous email, it's still not working. Luisa 2010/11/10 Hamish > > try r44284, hopefully fixed. > > perhaps buggy: > > -CMDLINE="`basename $0`" > +CMDLINE="`basename "$0"`" > > > becomes "`basename " + $0 + "`" > > afair the " do not nest. > > maybe outer quotes are not needed? > (but regardless I'm fairly sure this is not the cause of the v.colors bug, > just a red herring) > > > ?, > Hamish > > > > > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Using v.colors in WinGRASS
Markus Metz wrote: > I have seen something similar with v.rast.stats in wingrass. It seems > that in wingrass, scripts can not call other scripts, only C modules. > This might have something to do with the way scripts are called > through a .bat file. How do people feel about replacing shell scripts with Python for 6.5? -- Glynn Clements ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user