Re: [GRASS-user] Re: grass-user Digest, Vol 49, Issue 27
Is projecting to NH state plane and option? Along the lines of using something more localized. Mark On May 14, 2010, at 2:06 PM, Kurt Springs ferret_b...@mac.com wrote: Thanks Rich and Dylan I downloaded the pdf of document #1395. At the moment I am leaning toward Lambert Conic Conformal (1SP) since it seems to use Lat/Long of Natural Origin, in case I need to use a GPS. If I am reading you right Latitude and longitude don't even come into the equation, just the projection. I've been looking at the website http://www.dmap.co.uk/ utmworld.htm. I was mistaken it was 18 and 19T that NH falls in. However, it appears to be just the western most sliver. However, if I don't have to figure out the conversion, so much the better. Kurt On May 14, 2010, at 12:00 PM, grass-user-requ...@lists.osgeo.org wrote: Message: 2 Date: Fri, 14 May 2010 08:10:20 -0700 From: Dylan Beaudette dylan.beaude...@gmail.com Subject: Re: [GRASS-user] Latitude/Longitude vs UTM To: Rich Shepard rshep...@appl-ecosys.com Cc: GRASS user list grass-user@lists.osgeo.org Message-ID: aanlktimnbx-ed-ocxmdgtm_pw5ngzm48cbabs98tw...@mail.gmail.com Content-Type: text/plain; charset=ISO-8859-1 On Fri, May 14, 2010 at 6:22 AM, Rich Shepard rshep...@appl-ecosys.com wrote: On Thu, 13 May 2010, Kurt Springs wrote: This was interesting in that it told me that r.topidx could not be run with latitude and longitude and I had to convert to UTM. I was wondering if this is the answer to the problem and I just had to convert to UTM. Kurt, Lat/Long represents geographic coordinates, not a projection of location on a mathematial model of the earth. UTM is the Universal Transverse Mercador projection that we see on most printed (or computer displayed) maps of the earth. There is documentation within the GRASS Web site that provides a good explanation of the differences. GRASS modules work on geographic projections, not just coordinates. There is a USGS technical report from the mid-1980s that's the standard on projections. While it is becoming more rare to locatate, see if you can find a copy. I think that Rich is referring to this USGS document, #1395 http://pubs.er.usgs.gov/usgspubs/pp/pp1395 Definitely worth the price if you want to become an expert in map projections. One other question. New Hampshire appears to fall within two UTM zones (19T and 20T) Is there a way for a maps set to contain two UTM zones? Yes. Don't use UTM. In this case use a regional projection that suits your needs: 1) navigation -- use a conformal projection 2) area statistics -- use an equal-area projection ... etc ... Variations on the Albers or Lambert (conformal) conic projections work quite well for large regions that are wider than tall, but for such as small state should be just fine. We use an Albers equal-area projection to house soil survey data for the lower 48 states. Interesting. NH is a tall, narrow state so one would assume it would be within a single zone. Regardless, yes there is a way to reproject locations in one zone on the other, but it's non trivial and I've not done it. I wouldn't recommend it. The desirable properties of the UTM system (i.e. the fairly good compromise between distortion, preservation of angles, and preservation of area) only occur within a zone's boundaries. The farther you move from the central meridian of the UTM zone, the more distortion you will encounter-- therefore 'projecting' UTM z10 data into UTM z11 is technically possible, but not a great idea. Oregon is primarily in Zone 10, but the eastern edge (I don't recall the distance within the state) is in Zone 11. The available DEM and hydrologic data were reprojected from 11 to 10 by the supplying agency. Hmm... Dylan Rich ___ 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] dems from coordinate lists
If I follow correctly, instead of v.to.rast, you need to interpolate a raster DEM from the points. v.surf.rst produces nice results, but there are other interpolation modules as well in the raster category. Best regards, Mark On May 13, 2010, at 10:26 AM, Hanlie Pretorius hanlie.pretor...@gmail.com wrote: Hi, I've obtained DEMs in text files with columns X, Y and Z at 25m spacing. The first three entries in the text file are: - X,Y,Z 99550,2.9883e+06,1473.47 99550,2.98828e+06,1473.57 99550,2.98825e+06,1473.63 - To me it seems the easiest way to import these is the following: 1. v.in.ascii 2. Set the region to the extents of new vector file and the resolution to 25m. 3. v.to.rast. This 'works' but I get horizontal strips of no data in my raster DEM at 25m spacings. Also, when I look at the raster and the vector layer together, the vector points are not always on the edges of the raster cells. Can someone perhaps help me to fix this? g.region: projection: 99 (Transverse Mercator) zone: 0 datum: ** unknown (default: WGS84) ** ellipsoid: wgs84 north: 2988300 south: 2959850 west: 74300 east: 99550 nsres: 25 ewres: 25 rows: 1138 cols: 1010 cells: 1149380 v.info: +-- --- --- + | Layer: dem_2628cc_...@c83 | Mapset: C83 | Location:sa_tm_19deg_E | Database:C:\Hanlie\grassdata | Title: | Map scale: 1:1 | Map format: native | Name of creator: Administrator | Organization: | Source date: Thu May 13 16:06:07 2010 |-- --- --- | | Type of Map: vector (level: 2) | | Number of points: 1151529 Number of areas: 0 | Number of lines:0 Number of islands:0 | Number of boundaries: 0 Number of faces: 0 | Number of centroids:0 Number of kernels:0 | | Map is 3D: Yes | Number of dblinks: 0 | | Projection: Transverse Mercator | N: 2988300S: 2959850 | E: 99550W: 74300 | B: 1429.79T:1740.2 | | Digitization threshold: 0 | Comments: | +-- --- --- + r.info: +-- --- --- + | Layer:dem_2628cc_...@c83 Date: Thu May 13 16:15:50 2010 | Mapset: C83Login of Creator: Administrator | Location: sa_tm_19deg_E | DataBase: C:\Hanlie\grassdata | Title:Categories ( dem_2628cc_25m ) | Timestamp: none |-- --- --- | | | | Type of Map: raster Number of Categories: 0 | Data Type:DCELL | Rows: 1138 | Columns: 1010 | Total Cells: 1149380 |Projection: Transverse Mercator |N:2988300S:2959850 Res:25 |E: 99550W: 74300 Res:25 | Range of data:min = 1429.85 max = 1739.41 | | Data Source: |Vector Map: dem_2628cc_...@c83 in mapset C83 |Original scale from vector map: 1:1 | | Data Description: |generated by v.to.rast | | Comments: |v.to.rast input=dem_2628cc_...@c83 output=dem_2628cc_25m use=z\ | type=point,line,area layer=1 value=1 rows=4096 | +-- --- --- + Thanks Hanlie ___ 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] r.watershed: pretty flat areas
Great news. It sounds like this will be helpful with the very flat and wetland rich terrain in Florida, US. Mark On May 12, 2010, at 5:51 AM, Markus Metz markus.metz.gisw...@googlemail.com wrote: Hi all, r.watershed in trunk r42236 has a new option to beautify flat areas, activated with a new -b flag, works with both SFD and MFD. The method is adapted from Garbrecht Martz (1997). The justification for the new option is that some regard straight flow lines in flat areas as an issue, e.g. Nardi et al (2008), Hydrogeomorphic properties of simulated drainage patterns using digital elevation models: the flat area issue http://dx.doi.org/10.1623/hysj.53.6.1176 For r.watershed, this is not that much of an issue because it does not create flat areas by sink filling. For naturally flat areas like lakes and broad rivers, it is not possible to determine the correct drainage direction only from a DEM, information about river and lake bottom would be needed to more accurately determine flow direction, but DEMs only show water surface (apart from bathymetry). Therefore I chose to describe the new option as Beautify flat areas and not as correct flat areas. The new option is currently only implemented for the all in memory mode, and it requires more memory and more time. Contrary to other implementations, this new one modifies elevation values only temporarily to adjust drainage direction, flow is accumulated using the original elevation values. I'm not really convinced about this option, but since methods to modify flat areas appear regularly in the literature, I added one. If there is no strong support for the new option, I would remove it again. I guess it's mostly useful to produce pretty maps. Markus M ___ 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] Re: Extract coordinates of vertices/nodes
If a polygon has no other touching it, it should have 2 nodes and at least 2 verticies. When polygons touch other polygons, topology comes into play. So when polygon B is adjacent to, and touching polygon A, that makes a node because it is a shared line segment between polygon A and B. The lines break with nodes when other features connect, to create topologically correct, non-overlapping polygons. Mark On Apr 30, 2010, at 8:15 AM, Sophie Leguedois sophie.legued...@ensaia.inpl-nancy.fr wrote: Thank you very much for your useful answers. The commands v.to.point and v.out.ascii are doing the trick! About vertices and nodes, I noticed that you can have several nodes for the boundary limit of a polygon. Why so? Is it beacuse the nature of the neighbour polygon is different? Regards, So.L. - Sophie Leguédois Researcher - Inra (French National Institute for Agricultural Research) Laboratoire Sols et Environnement (UMR 1120 LSE) Ensaia - 2 avenue de la Forêt de Haye BP 172 - 54505 Vandoeuvre-lès-Nancy - France Tel: +33 (0)3 83 59 57 62 http://www.lse.inpl-nancy.fr/ -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Extract-coordinates-of-vertices-nodes-tp4973423p4985298.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 mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] specify tunnel/bridge in roadmap
If the bridge and tunnel are valid route options, can you remove the node and make them topologically incorrect? Meaning, crossing lines with no node intersection so there is no option to turn. I'm not sure if module needs topologically correct lines/network. Alternatively, if the bridge or tunnel is not a valid route, can you break the network connection by creating a gap in the line/network for either the bridge or tunnel so one route is not an option? Mark On Apr 29, 2010, at 3:51 PM, Robbie Heremans robbie.herem...@gmail.com wrote: Dear, In order to perform (shortest path) network analysis, how can I specify that the crossing of 2 lines in a roadmap is not an intersection but a tunnel/bridge (impossible to turn left or right)? Is v.digit the right way to do this (deleting en changing lines/ vertices)? Robbie ___ 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] r.to.vect Produces Too Many Areas
I haven't used the python GUI much. How about v.digit to graphically delete the blocks? Mark On Feb 22, 2010, at 4:33 PM, Rich Shepard rshep...@appl-ecosys.com wrote: On Mon, 22 Feb 2010, M S wrote: Since now a vector, can you select and delete those extraneous blocks? Mark, I've tried, by selecting both 'delete features' and 'delete nodes', but when I redraw the display they're still there. The wxPython GUI help web page doesn't tell me just how to do this. I exit and confirm the deletions, but they're still there. Probably user error. Rich ___ 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] Raster Map subset
There is probably a more eloquent way, but I have taken a vector watershed boundary, converted to a raster, and set as mask. Then used a r.mapcalc expression like newMap=oldMap+0. Mark On Feb 6, 2010, at 3:23 AM, Micha Silver mi...@arava.co.il wrote: Rich Shepard wrote: On Fri, 5 Feb 2010, stephen sefick wrote: Is there a way to get a subset of a larger raster map. I have zoomed into an area and would like to make a raster map only containing that area. Can I do this? I am going to use this as input into R for some topmodel runs. Stephen, Yes, you can. This subject has come up several times in the past month or so so the archives have a lot of options for you. One to look at closely is v.in.region. After zooming to the area you want, run v.in.region and save that vector map. You can then use it with the g.region command to limit the displayed bounds of your raster map. You can use r.patch to make a new raster map that's the size of the newly created vector map. I'm not sure this is entirely correct. v.in.region saves the current *computational* region to a vector, not the area zoomed to. Setting the desired region first with g.region, I think, is the correct way to go. Then r.mapcalc can create a new raster, limited to the new smaller region. Are you aware of r.topmodel within GRASS? Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Micha Silver Arava Development Co. +972-52-3665918 http://surfaces.co.il ___ 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] Need help with Grass GIS!!
Would you please post the full command usage? Thanks Mark On Jan 14, 2010, at 4:56 PM, sarah moore hopefading2...@yahoo.com wrote: To the Grass GIS community: I am new to Grass GIS and ultimately I want to connect attribute tables. First, however, I need get data from Vector works (on Mac) to a Grass GIS 6.4 (on PC). I have tried these steps: File, Import Vector, DXF, but I get error messages stating that I need to rename the vector map, in which I have tried. Thank you, Sarah ___ 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] help, interpolation of vector results in one value?
Is your region set to the extent of the vector points? Is the right attribute column being spicified to interpolate? Mark On Jan 9, 2010, at 2:34 PM, James Kebinger jkebin...@gmail.com wrote: Hello, I've imported a vector of points into GRASS - it's approximately a grid with ~24k points according to GRASS. I'm trying to interpolate it, but what ever methods I try (v.surf.idw or v.surf.bspline), the resulting raster is just one point, approximately in the middle of the grid with the average value of the z values. I am able to see the points rendered as a vector field, so the data looks like it's been imported correctly, but I can't figure out this interpolation step. Thanks! ___ 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] r.shade.relief: Results Not As Expected
Have you tried a different azimuth? Sometimes certain values can cause features to look inverted for given terrain conditions. Mark On Jan 4, 2010, at 8:32 PM, Rich Shepard rshep...@appl-ecosys.com wrote: On Tue, 5 Jan 2010, Glynn Clements wrote: Run r.slope.aspect on the interpolated map, and check for grid-like artifacts. Glynn, Did this. Yes, there is a grid-like appearance to the resulting map. How good is the underlying data? If it's noisy, it may need to be filtered before interpolating. How do I check for noisy data? I have all the basic information about the source data and derived maps. Also: which version of GRASS are you using? 6.5svn. Thanks, Rich ___ 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] Re: jgrass source code?
Jarek: just a note that the stream moules thus far are *awesome*. I love the fine tuning through input parameters. The stream extraction and sub-basins have worked great. Anxious to try the other modules, and see the new updates coming! Great work and thanks for the contribution! Mark On Jan 2, 2010, at 2:08 PM, Jarek Jasiewicz jar...@amu.edu.pl wrote: Margherita Di Leo pisze: Rich Shepard wrote: On Fri, 1 Jan 2010, Margherita Di Leo wrote: did you try to ask in jgrass mailing list? [1] Margherita, No. Just curiosity.. which module of jgrass in your opinion grass needs? Those not in the GRASS distribution. The suggestion was made by others on the GRASS mail list, and I do hydrologic modeling, so I decided to look into the Hi Note that big part of that work is done or almost done, see r.stream.extract, r.stream.order, r.stream.distance and r.stream.stats modules in add-ons. In next first days I will send major update of r.stream.order and r.stream basins with new functionality. Maybe it could be better you if join us and help to improve or expand that what is, instead of creating it from scratch? Also note that python (and java too) is not good tool for hydrological modeling (especially artificial code created with converters) because python is a scripting language but hydrological modeling needs loots of RAM and loots of processing time Jarek ___ 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.patch: Correct Syntax?
v.select using the polygon to make selection. Mark On Dec 31, 2009, at 10:08 AM, Rich Shepard rshep...@appl-ecosys.com wrote: On Wed, 30 Dec 2009, Hamish wrote: Here is an example using the Spearfish sample dataset: # preform the overlay operation v.overlay ain=soils bin=zoom_box op=and out=soils_zoom_box Hamish, The map types I need to clip include areas, lines, and points. The v.overlay man page tells me that the atype options are line and area. What do I do for points (e.g., dams)? Is there a different module to clip point maps to a smaller region? Thanks, Rich ___ 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.overlay: An Observation About Time
v.select can be a nice first pass to filter back a substantial amount of unnecessary data. Mark On Dec 30, 2009, at 2:33 PM, Rich Shepard rshep...@appl-ecosys.com wrote: I have a map of water bodies for the entire state. Yesterday, when I re-imported it (correctly), it took about 3 hours to complete; most of the time was waiting for it to 'break boundaries'. Today I'm clipping it to a very small region, one drainage basin and parts of the surrounding basins, yet it's taking about the same amount of time. Yes, there are a lot of nodes, vertices, and primitives, but I want only those within the area of the overlaid map. So I'm surprised it's taking as long as it is. This is an extensive original map: Building topology for vector map bl_lakes... Registering primitives... 220819 primitives registered 4589198 vertices registered Number of nodes: 158848 Number of primitives: 220819 Number of points: 0 Number of lines: 0 Number of boundaries: 220819 Number of centroids: 0 Number of areas: - Number of isles: - Breaking lines... Removing duplicates... Cleaning boundaries at nodes... Breaking lines... This is not a complaint, but an observation. Rich ___ 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] New version of r.stream.order
This sounds great. Thanks for the contribution. I have been wanting to use these stream analysis tools, but ran into a packaging problem under kubuntu, and am still trying to sort it out. I seem to recall it having to do with the grass-dev package in ubuntugis repo. Mark On Nov 9, 2009, at 9:47 AM, Jarosław Jasiewicz jar...@amu.edu.pl wrote: Hi all! I updated r.stream order to the new version. In new version Horton and Hack do not require accum map. It use the longest stream length (calculated internally) to determine both Hack and Horton. But accum it is still leave as an option. So now Hack stream order indicate the longest (previously most accumulated) stream in the catchment. Due to the module do not require to read huge accum maps for horton and hack orders it works almost twice faster than previously. I also added r.stream.del (see more in description) but its functionality is now integrated with r.stream.extract so I add it only as an option for r.watershed. It will be probably depreciated in the future. regards Jarek Jasiewicz ___ 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] having problems importing shapefiles
Sometimes I get those import errors, and usually delete a problematic column or not use a table. I figured it was a problematic dbf, but perhaps not? Mark On Oct 28, 2009, at 5:17 PM, Carlos Grohmann carlos.grohm...@gmail.com wrote: Hi all, I'm facing some serious issues importing shapefiles (latest svn). I always end up with a buffer overflow, no matter the size of the file (error msg at the end of this) I did managed to use v.external, but if I try to reproject the linked shp, it comes without any table.. I also noticed that sometimes, using wxgui, I see python eating about 40% of my cpu... this looks bad Just some more notes on wxgui: v.overlay: window title set at v.overla, and in v.external, the fiels for the dns path is set to be too long, it should be smaller. On a side note: I know that programming the new GUI is a lot of work, but every time I have to use r.mapcalculator on wxgui, I just give up the annoying clicking and usually just start the old tcltk GUI so I can run the stuff as I am sued to it, like if (A x, null(), B). Can we have something like that in wx??? please?? :) tks Carlos GRASS 6.4.0RC5 (newLocation):~ *** buffer overflow detected ***: v.in.ogr terminated === Backtrace: = /lib/tls/i686/cmov/libc.so.6(__fortify_fail+0x48)[0x264de8] /lib/tls/i686/cmov/libc.so.6[0x263e20] /lib/tls/i686/cmov/libc.so.6[0x263779] /lib/tls/i686/cmov/libc.so.6(__snprintf_chk+0x34)[0x263664] /usr/lib/libgdal1.5.0.so.1(_ZN10OGRFeature16GetFieldAsStringEi+0x34e) [0x135e3be] /usr/lib/libgdal1.5.0.so.1(OGR_F_GetFieldAsString+0x28)[0x135f1c8] v.in.ogr(main+0x21e8)[0x804d9d8] /lib/tls/i686/cmov/libc.so.6(__libc_start_main+0xe6)[0x19ab56] v.in.ogr[0x804aab1] === Memory map: 0011-0011d000 r-xp 08:01 321277 /usr/lib/grass64/lib/libgrass_dbmibase.6. 4.0RC5.so 0011d000-0011e000 r--p c000 08:01 321277 /usr/lib/grass64/lib/libgrass_dbmibase.6. 4.0RC5.so 0011e000-0011f000 rw-p d000 08:01 321277 /usr/lib/grass64/lib/libgrass_dbmibase.6. 4.0RC5.so 0011f000-00151000 r-xp 08:01 280413 /usr/lib/libproj.so. 0.5.5 00151000-00152000 r--p 00032000 08:01 280413 /usr/lib/libproj.so. 0.5.5 00152000-00154000 rw-p 00033000 08:01 280413 /usr/lib/libproj.so. 0.5.5 00154000-00169000 r-xp 08:01 279041 /usr/lib/ libgeos_c.so.1.5.0 00169000-0016a000 r--p 00014000 08:01 279041 /usr/lib/ libgeos_c.so.1.5.0 0016a000-0016b000 rw-p 00015000 08:01 279041 /usr/lib/ libgeos_c.so.1.5.0 0016b000-0016d000 r-xp 08:01 768/lib/ libkeyutils-1.2.so 0016d000-0016e000 r--p 1000 08:01 768/lib/ libkeyutils-1.2.so 0016e000-0016f000 rw-p 2000 08:01 768/lib/ libkeyutils-1.2.so 00171000-00182000 r-xp 08:01 321293 /usr/lib/grass64/lib/libgrass_dig2.6.4.0R C5.so 00182000-00183000 r--p 0001 08:01 321293 /usr/lib/grass64/lib/libgrass_dig2.6.4.0R C5.so 00183000-00184000 rw-p 00011000 08:01 321293 /usr/lib/grass64/lib/libgrass_dig2.6.4.0R C5.so 00184000-002c2000 r-xp 08:01 2505 /lib/tls/i686/cmov/libc-2.10.1.so 002c2000-002c4000 r--p 0013e000 08:01 2505 /lib/tls/i686/cmov/libc-2.10.1.so 002c4000-002c5000 rw-p 0014 08:01 2505 /lib/tls/i686/cmov/libc-2.10.1.so 002c5000-002c8000 rw-p 00:00 0 002c8000-0032c000 r-xp 08:01 3913 /usr/lib/libodbc.so. 1.0.0 0032c000-0032d000 r--p 00063000 08:01 3913 /usr/lib/libodbc.so. 1.0.0 0032d000-00331000 rw-p 00064000 08:01 3913 /usr/lib/libodbc.so. 1.0.0 00331000-00334000 r-xp 08:01 758/lib/libgpg- error.so.0.4.0 00334000-00335000 r--p 2000 08:01 758/lib/libgpg- error.so.0.4.0 00335000-00336000 rw-p 3000 08:01 758/lib/libgpg- error.so.0.4.0 00337000-00338000 r-xp 00:00 0 [vdso] 00338000-00381000 r-xp 08:01 5305 /usr/lib/ libjasper.so.1.0.0 00381000-00382000 r--p 00048000 08:01 5305 /usr/lib/ libjasper.so.1.0.0 00382000-00385000 rw-p 00049000 08:01 5305 /usr/lib/ libjasper.so.1.0.0 00385000-0038b000 rw-p 00:00 0 0038b000-003ad000 r-xp 08:01 279373 /usr/lib/ libmfhdf.so.4.1r4 003ad000-003ae000 r--p 00021000 08:01 279373 /usr/lib/ libmfhdf.so.4.1r4 003ae000-003af000 rw-p 00022000 08:01 279373 /usr/lib/ libmfhdf.so.4.1r4 003af000-003b rw-p 00:00 0 003b-0042c000 r-xp 08:01 279372 /usr/lib/libdf.so. 4.1r4 0042c000-0042d000 r--p 0007c000 08:01 279372 /usr/lib/libdf.so. 4.1r4 0042d000-0042f000 rw-p 0007d000 08:01 279372 /usr/lib/libdf.so. 4.1r4 0042f000-00457000 rw-p 00:00 0 00457000-00473000 r-xp 08:01 289756 /usr/lib/libogdi.so. 3.2 00473000-00474000 r--p 0001b000 08:01 289756 /usr/lib/libogdi.so. 3.2 00474000-00475000 rw-p 0001c000 08:01 289756 /usr/lib/libogdi.so. 3.2 00475000-00497000 r-xp 08:01 5307 /usr/lib/libjpeg.so. 62.0.0 00497000-00498000 r--p 00021000 08:01 5307 /usr/lib/libjpeg.so.
Re: [GRASS-user] Help creating dissolved multipolygon
Is there output for the dissolve command? Mark On Oct 20, 2009, at 4:37 PM, Roger André ran...@gmail.com wrote: Hi all, I'm having trouble creating a single, dissolved feature and exporting it as a multipolygon. Here is what I've done: - layer I started with: v.info -t global_mask nodes=5448 points=0 lines=0 boundaries=2727 centroids=2723 areas=2723 islands=2721 faces=0 kernels=0 primitives=5450 map3d=0 - added dissolve field: v.db.addcol global_mask col=diss int v.db.update global_mask col=diss value=1 - check that diss field is there: v.info -c global_mask Displaying column types/names for database connection of layer 1: INTEGER|cat DOUBLE PRECISION|cat_ INTEGER|diss - dissolve: v.dissolve global_mask out=new_mask col=diss --o - inspect new layer: v.info -t new_mask nodes=5448 points=0 lines=0 boundaries=2727 centroids=2723 areas=2723 islands=2721 faces=0 kernels=0 primitives=5450 map3d=0 - export layer as shapefile: v.out.ogr -p input=new_mask type=area dsn=/tmp olayer=new_mask - check new shapefile: ogrinfo -summary /tmp new_mask INFO: Open of `/tmp' using driver `ESRI Shapefile' successful. Layer name: new_mask Geometry: Polygon Feature Count: 2723 Extent: (-180.00, -90.00) - (180.00, 83.623596) Layer SRS WKT: GEOGCS[GCS_WGS_1984, DATUM[WGS_1984, SPHEROID[WGS_1984,6378137,298.257223563]], PRIMEM[Greenwich,0], UNIT[Degree,0.017453292519943295]] cat: Real (11.0) As you can see, I still have 2723 features, when I expect to have 1. What am I doing wrong? Roger ___ 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] r.basins.fill usage
Thanks to all for the help and guidence. Mark On Oct 16, 2009, at 6:20 AM, Jarosław Jasiewicz jar...@amu.edu.pl wrote: Hi try r.stream.basins from grass add-ons it has options (-c and -l flags) to generate subbasin without any manual work. Look also into tutorial on grass wiki about it (wiki - dokuments -r.stream.* Jarek Achim Kisseler pisze: Based on this excerpt from the manual, it seems the main external/perimeter watershed boundary, as well as the internal ridges need inputted? The raster map layer depicting ridges should include the ridge which defines the perimeter of the watershed. The coded stream network can be generated as part of the r.watershed program, but the map layer of ridges will need to be created by hand My question is really about the need to manually digitize ridges inside the watershed boundary delineated by r.watershed. As opposed to running r.watershed with a smaller threshold to get these smaller basin delineations within a given study area. If one has to manually digitize the internal ridges in a watershed, I wondered what the point of the module was. It obviously has application, but it seems I am missing something obvious. Hope that is more clear. Now its much clearer to me. Thanks. In my opinion the description ao r.basins.fill is now very detailed. I don't really understand it. So maybe someone else can help and maybe we could find out and enhance the description: I think this one: r.basins.fill generates a raster map layer depicting subbasins, based on input raster map layers for the coded stream network (where each channel segment has been coded with a unique category value).. is clear. Input is: raster map X from r.watershed stream=X Not clear to me is: ...and for the ridges within a given watershed. which is explained as: The raster map layer depicting ridges should include the ridge which defines the perimeter of the watershed. What should the parameter mean? How should that be stored in a raster-map? So try out: - take the stream-network and run r.mapcalc: r.mapcalc perimeter=if(X,5,null) Resulting map is a rastermap with value 5, where stream-network is not null. The five is raster-cells or map units or what else? - run r.basins.fill with that and see what happens. I guess: from each network-segment the subbasin is created which has a minimum distance of 5 (cells?) to the stream. So its something like a buffer. ...in contrast to: ...but the map layer of ridges will need to be created by hand (for example, through digitizing done in v.digit) = the ridges-map here is a vector-file! Achim PS: --%-- DESCRIPTION r.basins.fill generates a raster map layer depicting subbasins, based on input raster map layers for the coded stream network (where each channel segment has been coded with a unique category value) and for the ridges within a given watershed. The raster map layer depicting ridges should include the ridge which defines the perimeter of the watershed. The coded stream network can be generated as part of the r.watershed program, but the map layer of ridges will need to be created by hand (for example, through digitizing done in v.digit). The resulting output raster map layer will code the subbasins with category values matching those of the channel segments passing through them. A user-supplied number of passes through the data is made in an attempt to fill in these subbasins. If the resulting map layer from this program appears to have holes within a subbasin, the program should be rerun with a higher number of passes. --%-- ___ 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] r.basins.fill usage
The man page seems to imply that internal divided need digitized? Much thanks, Mark On Oct 15, 2009, at 9:27 AM, Achim Kisseler a...@jupiter.uni- freiburg.de wrote: I'm not sure if I understood you right: the border of a watershed is unique, if you see the basin starting from outlet of a river into the ocean. But there are unlimited subbasins (up to the resolution of the raster)! So you have to define a characteristic for the watersheds, eg. a basin size from the threshold of r.watershed. Greets, Achim M S schrieb: Upon first glance, r.basins.fill seems like a great tool to generate subbasins within a watershed. Although I must not be understanding something obvious and essential. If one has to digitize in the internal ridges in a given watershed, why run the ridge input through the module if one has essentially just delineated the subbasins ridges? Is there an advantage to using a smaller threshold on r.watershed to get smaller basin delineations (ceasing when thin strip basins get introduced)? Thanks for any feedback, Mark ___ 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] [GRASS-dev] error on gui startup - solved
Hi all: just a summary and new twist with GRASS GUI error on Kubuntu, which was solved thanks to Markus. On my kubuntu 9.04, the rgb.txt file was in another location than the previous solution, but locate found it in /usr/share/vim/vim72/rgb.txt I put that in the /etc/X11 directory and rebooted and all is good. On upgrading a system to 9.04, the file was still in /etc/X11, on a new install it was not there. Mark On Wed, Aug 5, 2009 at 10:10 AM, Markus Neteler nete...@osgeo.org wrote: On Wed, Aug 5, 2009 at 2:47 PM, M Smsei...@gmail.com wrote: I did a fresh install of Kubuntu 9.04. Then installed GRASS and QGIS (and associated dependencies with apt-get). When I launch GRASS I get the following error(s): bnow...@mapping:~$ grass -tcltk Cleaning up temporary files ... Starting GRASS ... Error in startup script: unknown color name black Please check: [GRASS-user] Tcl/Tk gui and nviz errors http://lists.osgeo.org/pipermail/grass-user/2009-July/thread.html#51777 for a solution. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Flood simulation
Under map type conversions, you can convert raster to vector. Mark On Aug 1, 2009, at 4:36 PM, Pavel Iacovlev iacovlev.pa...@gmail.com wrote: Good day all, I want to simulate a simple flood situation. My steps are 1. Download the elevation data 2. Use r.lake to generate the raster image of the flood area No problem with this, works like a charm. Now I would like to get the vector variant of the flood area so I can run intersects with my other vector data to see what is flooded and whats not. Any tips with what tool or combination of commands I can achieve this ? -- http://iap.md, The future is open ___ 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] Building topologically correct stream network
I'm not sure why there would be disconnected segments in the accumulation output. Mark On Jul 20, 2009, at 3:35 PM, stephen sefick ssef...@gmail.com wrote: I have a low relief area in the costal plain of georgia that I am trying to produce a topologically correct stream network. There are many breaks in the flow accumulation grid and I would like to delineate the streams and correct the topology after r.thin and r.to.vect have been preformed so that I can then use v.strahler to order the streams. Thanks for any help and I will provide and data that would be helpful. thanks, -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ 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] Site Suitability like Arc Spatial Analyst
If you have polygon data, convert each one to raster with weighted value. Then use the map calculator to add together for a weighted map. Mark On Jun 4, 2009, at 8:25 PM, Frank Aragona fr...@agroinnovations.com wrote: Grass list, Anybody ever done a site suitability analysis exactly as ArcGIS Spatial Analyst would do? Ranking different variables per priority, inputing different layers, etc.? I'm quite sure this is possible using GRASS, I just can't really imagine how it would be done. Please point me in the right direction. Thanks in advance, -- Frank Aragona ___ 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] edges in basin map from r.watershed
IMO the esri term of a filled DEM being hydrologically correct is a misnomer. Natural terrain has real depressions that impact surface water flow. Big depressional wetlands can retain water and release via groundwater or evapotranspiration. I like how r.watershed acomodates known depressions and handles the flow as interception. Also, then one can use other tools to find problematic areas in a raw DEM. I modeled an internally drained basin using known depressions in GRASS, and it worked fantastic. One problematic example for esri is modeling an internally drained or sink-watershed. These have no surface water outlet. If one filled it to get flow out of a 9 sq. mile watershed, the esri analyses are then meaningless. Which is one reason why filling a dem just to get an esri module to work and calling the DEM hydrologically correct is a misnomer an a limitation. Mark On Feb 13, 2009, at 11:09 AM, Markus Metz markus.metz.gisw...@googlemail.com wrote: Christian Schwartze wrote: Dear GRASS users, with r.watershed I get strange basin boundaries for some areas und I'm not able to give account of it. Attached you can find that part of the basin map which looks curiously. I means the sharp-edged regions... Whats the reason? This is most probably a flat area (no slope). Flow direction, flow accumulation, stream segments and basins can not reasonably be calculated for flat areas, these are regarded as missing information and some assumption has been made by the algorithm. What could help is to use a raster DEM as input that is *not* filled, some would say not hydrologically correct, but r.watershed works better with the raw, not filled DEM. What could also help, if this does not work or it really is a flat area, is r.watershed of grass7 with multiple flow direction. Note that the result may look nicer, but it still holds true that drainage direction (and therefore all other output) has to be estimated for flat areas. The A * Search of r.watershed is doing a pretty good job, and multiple flow accumulation can improve it a bit more, within limits. Basis is an Arc Info .adf raster file for DEM data. I think Arc Info wants a depressionless, filled, hydrologically correct DEM. r.watershed does explicitely not want such a DEM, it wants a raw DEM with depressions, not filled, and not hydrologically correct. I hope that helps, Markus M ___ 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] Re: understanding r.watershed
R.terraflow outputs catchments around sinks in the dem. R.watershed should be making huge delineations in comparison. Mark On Feb 2, 2009, at 11:40 AM, Markus Metz markus.metz.gisw...@googlemail.com wrote: Georg Kaspar wrote: If these lakes have an outflow, i.e. water is leaving these lakes, the results will be more realistic when you omit the depression input to r.watershed and only use the (not filled) DEM. If you are talking about the basins output having NULL values around these depression, this is because your basin threshold value was too high, set it to a lower value. Check the current flow accumulation output for a reasonable threshold. so, my accumulation map contains values from -144714 to 58920 with a majority of cells between 0-100. what would be an appropriate threshold value? You can try standard deviation of flow accumulation (r.univar -g). I already tried 5000 and it looked similar to the output I received from r.terraflow, The basins of r.watershed look similar to the basins of r.terraflow? but when running r.watershed with depression input I still receive those null()-areas... ... in the basins output I assume. If you really want to treat lakes as real depressions, i.e. water is not supposed to leave these lakes, then there will always be NULL areas around the lakes unless you set the threshold to something much smaller than 100, but then the stream segments and basins become meaningless because there are too many... by the way, this is what my region looks like: GRASS 6.2.3 Rather use 6.4.0RC3 instead of 6.2.3, you should get a couple of nice surprises... [...] cells: 466400 Not that many cells, should take just a few seconds to get the basins (with 6.4.0RC3). ___ 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] Tips for setting up an new FOSS-GEO-linux-box
That is really cool. Thanks for putting this together! Mark On Nov 2, 2008, at 9:13 PM, Nikos Alexandris [EMAIL PROTECTED] wrote: I've started the attempt to collect in a wiki-page tips for setting up a GIS workstation. http://wiki.osgeo.org/wiki/GIS_workstation_setup_tips It certainly needs modifications. Kind regards, Nikos ___ 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] How to install Grass into Ubuntu
sudo apt-get install grass Mark On Oct 19, 2008, at 11:28 AM, Van PD Tri [EMAIL PROTECTED] wrote: Hi to all, I am a new user of GRASS as well as Ubuntu. I faced troubles on understanding the guidelines from GRASS intrusion. Could any one help me? preferly step by step. Many thanks for your help. Best regards, Tri Van This message was sent using IMP, the Internet Messaging Program. ___ 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] (no subject)
Can your CAD system read shapefiles? Mark ~ iPhone'd On Sep 29, 2008, at 3:33 AM, Patrick Schirmer [EMAIL PROTECTED] wrote: Dear GRASS-users, does someone of you know how to export a GRASS-polygone layer to CAD polygones? I used GRASS6.3.0 to create closed polygones from a 'unclean' dwg file and want to export this GIS-file to go on drawing within a CAD- programm (we're an architecure/urban planner office and my collegues are not used to use GIS). I tried v.out.dxf , but this does explode the polygones to single lines, so it ain't working. Dwg might be a solution but there's no v.out.dwg in my GRASS, converters from dxf to dwg are not free of license and the registration in opendwg also is not free of charge. Thanks for your help, Patrick Schirmer -- GMX startet ShortView.de. Hier findest Du Leute mit Deinen Interessen! Jetzt dabei sein: http://www.shortview.de/[EMAIL PROTECTED] ___ 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] Re: grass
I have used this module a lot. Any more details? Maybe input parameters and steps leading up to this error please? Thanks. Mark ~ iPhone'd On Aug 31, 2008, at 5:22 AM, Markus Neteler [EMAIL PROTECTED] wrote: On Fri, Aug 29, 2008 at 12:48 PM, Nikos Alexandris [EMAIL PROTECTED] wrote: Dear GRASS friends, A friend has some problem with r.terraflow for which I have no experience. On Tue, 2008-08-26 at 16:09 +0200, Christian Herold wrote: Assertion failed: nrows * ncols == str-stream_len(), file grass2str.h, line 144 Any suggestions? Yes: tell us which GRASS version and which operating system... Markus ___ 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