Re: [GRASS-user] Reclass: Atributte a string
Hello Glynn Thanks for the reply. About using r.categories, can I apply it to a raster map and then export it as a normal raster? Thanks Franz On Fri, Nov 12, 2010 at 6:18 AM, Glynn Clements gl...@gclements.plus.comwrote: 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 gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Exporting Vector data to SHP (v.out.ogr)- takes a lot of time
The problem is that, I left it this night processing and, when I returned it has freezed. So, v.out.ogr is not being able to export vector files produced by r.to.vect... As I said, maybe it's related with some option in r.to.vect that is messing this up. Or something related with tables. Can anyone give me some tips on this? Thanks Luis ...that´ll be a HUGE vector file! So it´s normal, that it takes some minutes. On Thu, Nov 11, 2010 at 3:08 PM, Luis Lisboa luislisboa1...@gmail.comwrote: 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
Re: [GRASS-user] Re: Need some help in converting a contour map to DEM map...
of the DEM after exagerating it.Now I am 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.jpg) got imported into 3 different maps namely satimage.red, satimage.green, satimage.blue, therefore is just draped only one of it on the DEM). The output view that I get seems okay but its too dark (black and white-- please see the attachment), though the original aerial image was colourful which I had imported. Why is it happening so?? Can any one help me with this. Hi Nikhil, take a look at http://grass.fbk.eu/grass64/manuals/html64_user/imageryintro.html it might give you a good introduction. especially this part: Visualizing (true) color composites To quickly combine the first three channels to a near natural color image, the GRASS command d.rgb can be used or the graphical GIS manager (gis.m). It assigns each channel to a color which is then mixed while displayed. With a bit more work of tuning the grey scales of the channels, nearly perfect colors can be achieved. Channel histograms can be shown with d.histogram. Regards, Achim ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Exporting Vector data to SHP (v.out.ogr)- takes a lot of time
Hello Achim... Not really. It didn't work Let me start from the beggining: 1- My orioginal raster file has the following: | Type of Map: raster Number of Categories: 11 | | Data Type:CELL | | Rows: 5381 | | Columns: 4481 | | Total Cells: 24112261 | |Projection: UTM (zone -24) | |N: 9659623.5017993S: 8045266.26306668 Res: 300.01063719 | |E: 1077874.87767815W: -266394.50364563 Res: 299.993167 | | Range of data:min = 1 max = 11 | It's a landcover map 2- I convert to vectorusing the following expression r.to.vect -v input=teste01_l...@regional output=teste05 feature=area and it returns: Extracting areas... Building topology for vector map teste05... Registering primiti 907000907256 primitives registered 3527832 vertices registered Building areas... 321179 areas built 85875 isles built Attaching islands... Attaching centroids... Number of nodes: 671952 Number of primitives: 907256 Number of points: 0 Number of lines: 0 Number of boundaries: 586103 Number of centroids: 321153 Number of areas: 321179 Number of isles: 85875 Number of areas without centroid: 26 r.to.vect complete. My vector is: || | Type of Map: vector (level: 2) | | | | Number of points: 0 Number of areas: 321179 | | Number of lines:0 Number of islands:85875 | | Number of boundaries: 586103 Number of faces: 0 | | Number of centroids:321153 Number of kernels:0 | | | | Map is 3D: No | | Number of dblinks: 1 | | | | Projection: Universal Transverse Mercator (zone 0) | | N: 9659623.5017993S: 8045266.26306668 | | E: 911978.65632755W: -266394.50364563 | | | | Digitization threshold: 0 | | Comments: | | | ++ this is a vector composed by areas (as far as I wanted and I believe) 3- I reduced the region in order to export just a piece 4- I exported but, when I open it in QGIS I get an empty vector file 5- Then again: v.out.ogr -c input=test...@regional type=area dsn=c:\teste07.shp Exporting 321179 areas (may take some time)... And after 20 minutes it seems freezed It this too much polygons? On Fri, Nov 12, 2010 at 9:06 AM, Achim Kisseler a...@jupiter.uni-freiburg.de wrote: Am 12.11.2010 09:40, schrieb Luis Lisboa: The problem is that, I left it this night processing and, when I returned it has freezed. So, v.out.ogr is not being able to export vector files produced by r.to.vect... As I said, maybe it's related with some option in r.to.vect that is messing this up. Or something related with tables. Can anyone give me some tips on this? I guess it freezes because it uses so much memory. Please check if v.out.ogr works with a small vector map. E.g. shrink region before r.to.vect or choose only some features from the vector map to export. ? Achim Thanks Luis ...that´ll be a HUGE vector file! So it´s normal, that it takes some minutes. On Thu, Nov 11, 2010 at 3:08 PM, Luis Lisboa luislisboa1...@gmail.com mailto:luislisboa1...@gmail.com wrote: Greetings I have produced a vector map (with r.to.vector) with the following metadata: v.info http://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 | |
[GRASS-user] current state of r.li
Hi all. Sorry for crossposting. As some of you know, the r.li suite of GRASS commands allows landscape analyses[0]. Its interface is rather complex, and is still in TclTk, not ported to either wxpython or qgis. As such, it is now more difficult to use than it should be, and it will become unusable when TclTk support will be dropped. The possible solution (thanks Radim) is to rewrite the interface as a qgis python plugin. It should not be a huge work (we provisionally estimate 2-3 weeks). The question is: is there anybody willing to invest either his/her time, or some money, to write such a plugin? We (Faunalia) would be happy to help if necessary. All the best. -- http://www.faunalia.it/pc [0]http://grass.fbk.eu/gdp/html_grass64/r.li.setup.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Need some help in converting a contour map to DEM map...
Hello sir Thanks for the link. I finally got the coloured aerial image draped on the DEM map and which can be viewed in nviz. Since the aerial image was getting imported in an rgb files (3 different files) I first needed to combine them into one raster map, which could be done using r.composite. which takes in these 3 files (.red, .green, .blue) and one output map name of users choice. Thanks once again for your help sir. 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/Need-some-help-in-converting-a-contour-map-to-DEM-map-tp5720519p5732236.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] How to delete deadends?
Dear group, I have a lot of lines that are linked to another forming various nets. Several of those lines are deadends. Now I search for a option to delete dead-ends within those networks. It would be perfect to delete the lines, or to search for the nodes that link to more than one line. I was searching in v.generalize, v.net, v.to.point, v.to.db but won't find the proper approach. Doas anybody of you have an idea? Thanks for your feedback, Patrick Schirmer ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to delete deadends?
On Fri, Nov 12, 2010 at 2:07 PM, Patrick_schirmer patrick_...@gmx.ch wrote: Dear group, I have a lot of lines that are linked to another forming various nets. Several of those lines are deadends. Now I search for a option to delete dead-ends within those networks. It would be perfect to delete the lines, or to search for the nodes that link to more than one line. I was searching in v.generalize, v.net, v.to.point, v.to.db but won't find the proper approach. Try to remove dangles with v.clean tool=rmdangle If threshold is 0, all dangles will be removed, otherwise only dangles shorter than threshold. Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] i.smap and NULL (Landsat 7)
Dear Grass List, I am trying to do radiometric classifications on Landsat 7 SLC-Off imagery (i.e. it has lots of NULL data). My question is regarding i.smap. The documentation says: The module i.smap does not support MASKed or NULL cells. Therefore it might be necessary to create a copy of the classification results using e.g. r.mapcalc. Does this mean i.smap will not work correctly with Landsat 7 SLC-Off imagery? When I try i.gensigset and i.smap I get results, but I'm not sure if I can trust them due to the disclaimer. Thanks for any help or advice, Nick ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to delete deadends?
Patrick_schirmer wrote: I have a lot of lines that are linked to another forming various nets. Several of those lines are deadends. Now I search for a option to delete dead-ends within those networks. It would be perfect to delete the lines, or to search for the nodes that link to more than one line. I was searching in v.generalize, v.net, v.to.point, v.to.db but won't find the proper approach. not a v.net specific thing, but in general for vector maps you can use v.clean to remove (or isolate) line dangles with the 'rmdangle' (or 'chdangle') tools. maybe that helps, Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Exporting Vector data to SHP (v.out.ogr)- takes a lot of time
Luis wrote: 1- My orioginal raster file has the following: ... Registering primiti 907000 907256 primitives registered 900k is a lot, but not crazy. 3527832 vertices registered Building areas... 321179 areas built 300k areas is not small, but it is not huge either. ... Exporting 321179 areas (may take some time)... And after 20 minutes it seems freezed wait longer... maybe hours or overnight. export to shapefile with v.out.ogr is known to be quite slow. are you sure it is frozen? check in process monitor (top on mac or linux, ctrl-alt-del on Windows) if cpu usage by v.out.ogr is 100%, and how much memory it is using. if it is frozen cpu will be 0%, if it is still working cpu will be 100%. It this too much polygons? perhaps. maybe resample raster to a coarser grid or try running v.build.polylines on the vector map prior to export. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Reclass: Atributte a string
Franz Schiller wrote: Thanks for the reply. About using r.categories, can I apply it to a raster map and then export it as a normal raster? category labels are attached to an existing map; there's nothing to export. It might help if you explain what your ultimate goal is. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user