[GRASS-user] working with shapefiles
Hi all, my problem is about topology and attribute storing : (I work on Grass 6.3-cvs, connected to a postgresql database) * I start from a line shapefile representing polygons boundaries, * import it with v.in.ogr, and the option 'boundaries' ; In order to build a topologically correct map, grass generates nodes at each boundary intersection, and it is actually what I want... The problem is elsewhere : the source shapefile, non-topological, contains certain lines that are not splitted at intersections. So, let's say the original shape contains 340 lines, the imported map then will contain 385 boundaries. But the attribute table contains 340 records (340 categories). For some reasons (exporting the table in an ascii format) I would like to have 385 records in the attribute table, that is in a way to duplicate certain records. What would you suggest me to do ? (my temporary and dirty solution was to export the boundaries in a line shapefile, then re-import it... not very clean) Thank you ! Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] rectifying digital images
Le samedi 05 janvier 2008 à 13:44 +0100, Markus Neteler a écrit : On Jan 5, 2008 1:35 PM, Markus Metz [EMAIL PROTECTED] wrote: ... Anyway, digital cameras are apparently not supported in i.ortho.photo. Why this? I did so, using a handheld digital camera: M. Neteler, D. Grasso, I. Michelazzi, L. Miori, S. Merler, and C. Furlanello, 2005. An integrated toolbox for image registration, fusion and classification. International Journal of Geoinformatics. Special Issue on FOSS/GRASS 2004 GIS-IDEAS 2004, 1(1), pp. 51-61, March 2005. PDF: http://ojs.info.gscc.osaka-cu.ac.jp/IJG/index.php/IJG/article/view/103/42 I heard of some tools that can add fiducial marks to digital photographs You can use the photo corners and calculate things yourself using jhead or other tools. You mean pointing corners of the image as four fiducials, then assigning them the dimensions of the sensor, in a coordinate system centered on the estimated principal point ? If you do that with i.ortho.photo what value would you give to the focal length, the one provided by the exif ? Vincent 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
Re: [GRASS-user] export e00
Hi Roland, unfortunately v.out.e00 is no longer maintained (formerly available at version 5). Some time ago I had initiated a thread about this subject on the french grass list (I just notice the whole archive of this list disappeared... does anybody know more about it ?) and Michel Wurtz answered me he did not adapt his module in the sixth version. sorry, Vincent Maybe the simplest way would be to try scripting an ungenerate Le lundi 07 janvier 2008 à 15:39 +0100, Roland Kaiser a écrit : Dear all! Please educate me: How can I export to ARC/INFO ungenerate format or e00 in GRASS 6.3 cvs. Thanks Roland ___ 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] export e00
Roland, You can still find the message in Google cache : http://209.85.129.104/search?q=cache:EjfdqnQ9BZYJ:www.linux-nantes.fr.eu.org/pipermail/grass-fr/2006-September/001006.html+michel+wurtz+bain+e00hl=frct=clnkcd=1gl=fr (in french). VB Le lundi 07 janvier 2008 à 15:39 +0100, Roland Kaiser a écrit : Dear all! Please educate me: How can I export to ARC/INFO ungenerate format or e00 in GRASS 6.3 cvs. Thanks Roland ___ 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] Vector drawings in PDF - converting them
Within Illustrator maybe you can build a script that adds intermediate points to curves (don't really remember, but it sounds like 'densify paths', or 'add anchor points'), attended that vertices will be recognized when imported Vincent Le jeudi 17 janvier 2008 à 09:03 +0100, Andreas Neumann a écrit : ok, the problems are with the curves. The ogr2ogr tool has a converter from Interlis to other ogrformats. Interlis supports arcs (defined by three points). They also calculate new vertices along the arcs to be compatible to other formats that don't support curves. One can define how many Every SVG viewer has code on how to display curves (arcs, quadratic and cubic spline curves). But I don't know if you can dig up and work with that code directly. Many of them probably forward their drawing commands to underlying graphics libraries like Cairo, AGG, OpenVG, Java2D, etc. As you said, you might also want to ask the Batik developers on the Batik list. These people are usually friendly and often help where they can. Andreas Andreas, Thanks for you comments. I agree completely. The curves are definitely where the problem lies. There would need to be a clever algorithm that converted them to polylines with some sort of configurable resolution. However, my guess is the code used convert SVG to raster formats like TIF and PNG has to do something similar. I began looking into using the facilities in Batik to accomplish this. The PNGTranscoder looks to be a good place to start. Haven't gotten very far in this past weekend when I began to look into it, no blatant roadblocks yet. --Kurt Andreas Neumann wrote: Hi, I don't have an immediate solution to your problem, however, I would like to discuss the use of the SVG format. Also, are you using curves in your original data? SVG would be a great format for transforming non-GIS vector maps into a GIS format. But it is not so surprising that most GIS only export and don't import SVG. SVG is usually a presentation format, not a transfer format for GIS data. Also, SVG has a lot of features that can't be easily transfered to GIS - think about elliptical arcs, cubic and quadratic spline curves. Those would have to be transfered into the OGC geometry models, where support for curves is more or less in its infancy. A number of GIS software allows the export of SVG, not always in a very form, though: Postgis (ok), Grass (did not try that), Mapserver (haven't tried that), ESRI (crappy, they seem to have some sort of resolution for their vector export). But, I think it would be great to have svg support incorporated into the ogr tools. It would make sense, also since FME and other OS and commercial GIS support SVG. Andreas I have a TON of vector drawings available to me in PDF format that I'm georeferencing. So far, I'm doing it by converting them to PNGs and using r.to.vect to massage them. This is really tedious and the data isn't as clean as a straight vector conversion would be. In Illustrator, I can see that all the vector math is there. That is, the PDF isn't simply a wrapper around raster data. I've tried opening the PDFs in Illustrator and exporting them as DXFs, but Grass ignores a lot of the data when I do this...they look very different in Grass once imported. I can export them as SVG from Illustrator, but surprisingly, there doesn't appear to be a single open-source tool available out there that converts SVG into some mainstream GIS file format. Does anyone have experience doing this type of conversion? Pointers welcome. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Orthorectification of Historic Aerial Photographs
May it be a stupid suggestion, but did you try to contact anyone at the USGS ? perhaps they maintain data from these old times... (here in France IGN is always helpful, and can provide answers to such historical requests) Vincent. Le jeudi 10 avril 2008 à 18:33 +0200, Markus Neteler a écrit : On Thu, Apr 10, 2008 at 5:29 PM, Shaun Melissa Busler [EMAIL PROTECTED] wrote: Unfortunately, the focal length of the camera is not on the photo border. The photo must be too old. Do you have any recommendations on what focal length I should try? It was taken in 1939 for the USDA Agricultural Adjustment Administration. You could try 150mm or 305mm. In fact, i.ortho.photo is quite robust, just try and hope that it converges. (and/or: Maybe someone on the PROJ4 list knows.) 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
[GRASS-user] Re: Post-processing ps.map output with inkscape
Yes, with ps2pdf the combination Grass/Inkscape works great... By the way (perhaps off subject, but it may be interesting for those who can't give up with Illustrator) : Inkscape 0.46 now handles transparency between layers. In the past it was not implemented which was a big limitation for clean graphical map processing. VB ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] how to calculate areas
to my mind, one has to be cautious with this module, and be aware of how surfaces are estimated. It appears that results show a big variability given the current region resolution. (just try r.surf.area on a sloppy plane [generated with r.plane], knowing its inclination and thus its actual surface in advance) Bye, Vincent Le mardi 16 septembre 2008 à 02:28 -0700, Alex Mandel a écrit : Alex Mandel wrote: Maciej Tomaszczyk wrote: hello I have DEM of some part of mountains chain and i would like to calculate a real area of this (area in 3d not plane area). I wonder if Grass can do it. best maciej Maciej Tomaszczyk Polish Geological Institute Geological Mapping Department http://www.pgi.gov.pl http://www.opengeology.pl Just to clarify, you mean surface area as opposed to footprint of the mountain range. So the sum of all the surfaces of the 3d model? Alex The solution to that problem would be r.surf.area which should give you both the 2D and 3D area. You'll have to actually try this as I just read about it in Markus and Helena's book. Alex ___ 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] Display window export
Has anybody had this problem : I try to export a raster image of the Display window (tcltk interface), with the button on the top bar : 'Export display to graphics file JPG* very high resolution (300% your current resolution)' The file generated through this command is actually 3 x the display resolution, but features (e.g. vector lines, raster maps, etc.) are drawn at the screen low resolution, i.e. 3 times coarser than the .jpg image resolution ! in other words 9 contiguous pixels have the same value. What does this suggest you ? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] interesting article on TIN generation
Dylan, As I spend a couple of days computing a huge lidar dataset (first split the original file in 145 pieces ! then performing a loop with v.surf.rst, and rearrange data in a single raster), it is of great interest for us. Perhaps could this code be turned into grass module ? well, it's beyond my skills, but... Le lundi 22 septembre 2008 à 15:53 -0700, Dylan Beaudette a écrit : Thought this would be of interest: http://www.cs.unc.edu/~isenburg/papers/ilsst-tin2dem-06.pdf Cheers, Dylan ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] interesting article on TIN generation
Le mardi 23 septembre 2008 à 01:17 -0700, Hamish a écrit : I still look for a nice method to paste together overlapping splines cleanly. Maybe first generating /sufficiently/ overlapping tiles, then adjust adjacent ones in the middle of the overlap (ok, not very clean) IM(V)HO the fact that TINs are used so much to make raster surfaces is an artifact kludge from a history of other well known GIS which was historically strong with vector data but weak with raster data Yes, it's hard to ged rid of some old /educational/ reflexes. I agree ! but anyway, TINGRID was very useful. If speed is a real limiting factor I'd suggest trying v.surf.idw(2), or perhaps r.surf.nnbathy*. Very useful module if one needs to generate a non-interpolated surface, fitting exactly a given set of points. Could it become a default command in Grass ? Perhaps could this code be turned into grass module ? well, it's beyond my skills, but... I could have sworn there was a TIN module on the wiki addons page, but don't see it now. Maybe it was lurking in the R-interface tutorials or somewhere? (The grass wiki as-setup doesn't search for 3 letter words, which is a bit of a problem) Maybe you confuse with the TIN how-to ? using Paraview and Meshlab http://grass.osgeo.org/wiki/HOWTO_create_3D_TIN 3c, Hamish [*] Is there any thoughts on moving r.surf.nnbathy into the main source? It requires an external dependency to use, but so do many other scripts. To me it's a valuable addition to the available quiver of interpolation methods; a nice compromise between IDW and splines. Before doing that I think to change it to be v.surf.nnbathy (its first step is r.to.vect). ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] interesting article on TIN generation
Sorry I forgot the PS [*] Is there any thoughts on moving r.surf.nnbathy into the main source? It requires an external dependency to use, but so do many other scripts. To me it's a valuable addition to the available quiver of interpolation methods; a nice compromise between IDW and splines. Before doing that I think to change it to be v.surf.nnbathy (its first step is r.to.vect). Yes indeed it would be faster when starting from a set of vector points. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] nnbathy
Given your system, don't know if it can help, but the Natural Neighbours interpolation library is available here : http://packages.ubuntu.com/hardy/libs/libcsiro0 Le mardi 23 septembre 2008 à 12:18 +0100, John Stevenson a écrit : Hi, I am keen to try r.surf.nnbathy, but I cannot find the nn libraries online. The link given in the instructions: http://www.marine.csiro.au/~sak007 is dead. Does anyone know where I can get them from now? Cheers John ___ 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] nnbathy
Sorry for the wrong way, and thank you Morritz for contacting P. Sakov ! VB Le mardi 23 septembre 2008 à 14:50 +0200, G. Allegri a écrit : You're right Moritz. I thought it contained also the executable... 2008/9/23 Moritz Lennert [EMAIL PROTECTED]: On 23/09/08 14:07, G. Allegri wrote: I could use the one from the Ubuntu repositories. Just out of curiosity: how does this work ? IIUC, the Ubuntu library version does not contain the nnbathy executable which r.surf.nnbathy calls, but only a C-library. How do you use this for the script ? Moritz Thanks 2008/9/23 Moritz Lennert [EMAIL PROTECTED]: On 23/09/08 13:49, Moritz Lennert wrote: On 23/09/08 13:18, John Stevenson wrote: Hi, I am keen to try r.surf.nnbathy, but I cannot find the nn libraries online. The link given in the instructions: http://www.marine.csiro.au/~sak007 is dead. Does anyone know where I can get them from now? http://web.archive.org/web/20080203105504/http://www.marine.csiro.au/~sak007/ I have also put it on my machine, in case the above link doesn't work anymore: http://geog-pc40.ulb.ac.be/grass/nnbathy/ As a follow-up: I have sent a mail to Pavel Sakov to find out where the code is available now. Moritz ___ 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] r.digit. no zoom/pan option while digitizing ?
Nagesh, me too I think that r.digit is not handy ; I suggest you to digitize your mask contour with v.digit. Then you can generate the corresponding raster map with v.to.rast, and apply r.mask to the latter. Bye, VB Le mercredi 24 septembre 2008 à 11:27 +0530, Nagesh Bhatkar a écrit : Hello Everyone, I was looking to create a raster mask by hand using r.digit module. The geocoded map is large and zooming is essential. There's no zoom/pan option while digitizing. I have to zoom into part of the raster map prior to running the module r.digit using d.zoom. I have to exit the r.digit module after each segment of the mask is created. And then reset region to fit the entire map. And then display the mask created and then select the next region to be zoomed in. Create many masks and then patch them using r.patch. I have gone through the archives, It says, One can't force a running r.digit program to change the region it uses internally. So one could change the region stored in the WIND file using g.region, but that won't affect r.digit. Is there any other way to create a raster mask! Regards Nagesh ___ 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] Importing polygon maps with overlapping features
Hello there, yes this can be seen as a limitation of grass, or more precisely of the topological data model. But I think one must remind of the strict border that should lie between geometry and a semantic contents. With a relational database structure you can of course overlap various characteristics of a single area : on one hand you have a set of adjacent -- that is topologically clean -- polygons, on the other hand a table or a set of related tables holding the complexity of attributes ; the cat value of a centroïd has a univoque link with a row of a table, then you can define associations of polygons (ajacent, included, non-adjacent). Grass does not provide tools for directly manipulating let's say 'multi-polygon' objects, but the ability to link a database to simple geometry allows the handling of complex objects. __ (a bit out of context : for those who worked with ArcInfo, it is sometimes frustrating not to have in Grass a feature class like REGION-SUBCLASS, which is to my mind a far better solution than multipolygons (the OGC standard geographic features)) Yours, Vincent. Le mercredi 07 janvier 2009 à 14:47 +0200, Maris Nartiss a écrit : Hello, GRASS vector model is advanced, but sometimes it fails for simple usage. It's one of best features is also it's point of failure. Vector areas support in GRASS is built around bogous assumption, that areas can not overlap. Such assumption holds true for many vector usages (i.e. property boundaries don't overlap), but fails for other vector usage patterns. Let's assume one GRASS user wants to create vector area map with suitable animal habitat areas. Does gay wolf and brown bear habitat areas may overlap in real world? Yes they can. Can they overlap in GRASS? No. User is forced to adapt semantics (habitat area) to data storage limitations (one map per species). Probably this is not the best example, simply I could not make better one fast. Sorry for language and trolling, Maris. 2009/1/7, Benjamin Ducke benjamin.du...@oxfordarch.co.uk: The only trouble this gives me in practice arises when I need to import data created in non-topological GIS (e.g. ArcView Shapefiles) that contains overlapping polygons. Granted, those should not exist in the first place, but bad data quality and faulty topology is a constant reality in my field of work. With overlapping polygons, centroids cannot always be related to exactly one polygon, so the topology building fails for those cases and attribute data does not get attached correctly. Cheers, Ben ___ 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] v.digit under wxGUI
Hello, Compiling grass64 on Ubuntu Hardy for 64bit system, I got this error when running v.digit under wxpython interface : /usr/local/lib/libgdi.so: undefined symbol: Py_InitModule4 Already pointed in this thread http://lists.osgeo.org/pipermail/grass-dev/2008-May/037763.html but I can't catch any solution in the links proposed by Markus. Does anyone have a solution ? or should I give up on my 64bit architecture ? Thank you, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.digit under wxGUI
OK, so in order to fully enjoy the wxpython GUI on an 64 bit architecture (especially the wxdigit module) one should first read carefully the README file located in the your_grass_source_location/gui/wxpython directory. The section named 7 - VECTOR DIGITIZER tells you to create a symbolic link to a wx library : sudo ln -s `locate _gdi_.so` /usr/local/lib/libgdi.so The command locate _gdi_.so may points towards the 32 bit library and ignore (?) the 64 bit libs path (/usr/lib64/...), so it fails. In my configuration I had to type : sudo ln -s /usr/lib64/python2.5/site-packages\ /wx-2.8-gtk2-unicode/wx/_gdi_.so /usr/local/lib/libgdi.so And it worked ! Le dimanche 18 janvier 2009 à 18:59 +0100, Markus Neteler a écrit : Excellent. Please post a short note to the list so that the next one in the trap finds it :) thanks Markus On Sun, Jan 18, 2009 at 6:21 PM, Vincent Bain b...@toraval.fr wrote: Yes ! you guessed my mistake, the link was pointing to the lib32 lib... Now it works fine. Thank you very much, Vincent. Le dimanche 18 janvier 2009 à 17:59 +0100, Markus Neteler a écrit : Vincent, On Sun, Jan 18, 2009 at 10:35 AM, Vincent Bain b...@toraval.fr wrote: Hello, Compiling grass64 on Ubuntu Hardy for 64bit system, I got this error when running v.digit under wxpython interface : /usr/local/lib/libgdi.so: undefined symbol: Py_InitModule4 Already pointed in this thread http://lists.osgeo.org/pipermail/grass-dev/2008-May/037763.html but I can't catch any solution in the links proposed by Markus. Does anyone have a solution ? or should I give up on my 64bit architecture ? No need to give up on 64bit, I am using it like many others here. Please post the result of file /usr/local/lib/libgdi.so I get /usr/local/lib/libgdi.so: symbolic link to `/usr/lib64/python2.5/site-packages/wx-2.8-gtk2-unicode/wx/_gdi_.so' As you see it links to the 64bit version of _gdi_.so. Perhaps you linked the 32bit version? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] interactive multiple selection
Thank you Markus, I tried several methods like d.extract, or d.where | v.what map=zonpi test and extracting cat values from the test file. Now I am wondering if it is worth (let's better say realistic given my skills) to try and write a custom module, inspired from v.digatt. As I already had a little experience in perl and wxwidget (frontend for data capture in postgres database) I would appreciate later on to involve in writing a python add-on, but like everyone, my principal ennemy is the lack of time ! Beside data storing, treatment, analysis, I think it would be great to improve data *production* tools, and it's what I would like to contribute. Vincent. Le mercredi 21 janvier 2009 à 09:32 +0100, Markus Neteler a écrit : On Wed, Jan 21, 2009 at 9:25 AM, Vincent Bain b...@toraval.fr wrote: Hello, considering a vector map, do you know a simple way to perform multiple selections (e.g. highlighting a set of polygons) interactively either on the display (tcltk or wxpython) or through a X monitor ? the aim in the end is to update database records on the basis of a cat list. I tried the addon v.digatt but I can't run it correctly (I think a possible reason is my vector is linked to a pg db driver, and the module was tested only with dbf driver). Coudl this help? d.extract help Description: Selects and extracts vectors with mouse into new vector map. It's not perfect but maybe a starter. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] interactive multiple selection
Very much related to this, please see my suggestions for updates to the GRASS query functions. Most should be relatively minor changes, but would make querying much easier. And it would make it easier to implement GUI wrappers for interactive querying. http://grass.osgeo.org/wiki/GRASS_6.3_Feature_Plan Updates to vector querying OK, good perspective ! Concerning the GUI wrapping of commands such as v.what or v.select, will it imply the ability to query directly within the display window ? Vincent. Thanks much for your interest. Michael ___ 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] interactive multiple selection
Yes. That's the idea. Martin Landa already has tried to implement this to some extent in the new wxPython GUI. But it involves a lot of workarounds for inadequacies in the GRASS vector querying. These updates to querying modules would make this easier and more robust. Michael Nice, Thank you all for that huge amount of work ! Vincent. ___ 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 other projection
Hello Mohammed, could you be a bit more explicit ? what is your problem ? where is located your data ? and so on... No one can really guess what you need. VB Le jeudi 22 janvier 2009 à 07:15 -0800, Mohammed Rashad a écrit : help we with using conical equal area projection for maps in GRASS ___ 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] wxGUI, display
Another suggestion, perhaps quite superficial : it would be interesting to reload the windows sizes and positions (layer manager, display, terminal). Vincent Le jeudi 22 janvier 2009 à 17:16 +0100, Martin Landa a écrit : Hi, 2009/1/22 Moritz Lennert mlenn...@club.worldonline.be: 1) saving/opening the entire workspace, i.e. including the information of which Displays are open, and 2) saving/loading the content of a particular display and so make it possible to load this information into another display ? this would require a bit more effort to implement it. Please create for this enhancement ticket on trac. I am afraid it will be not implemented soon, so it's better to have in trac. Martin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: using other projection
Perhaps a good way to proceed is to read some documentation for a first approach. This document https://subversion.gdf-hannover.de/dl.php?download=gdf_grass60_v1.2_en.pdf can help you, especially chapters 4 to 6. In a few words, you have to create a location according to your projection parameters (type, datum, ellipsoïd, etc.), another location in XY coordinate without geocoding in which you import your rasterized map. Then you can proceed georeferencing from local location towards target location. Good luck, Vincent. Le jeudi 22 janvier 2009 à 07:46 -0800, Mohammed Rashad a écrit : i need to use conical equal area because my maps are done using it. i scanned the maps and want to create a conical equal area (Albers equal area) with two standard parallels. so please send the steps in georeferencing those maps into the above projection ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] wxGUI, display
Yes it is (formerly I believe it was not, sorry). Vincent Le jeudi 22 janvier 2009 à 18:16 +0100, Martin Landa a écrit : Hi, 2009/1/22 Vincent Bain b...@toraval.fr: Another suggestion, perhaps quite superficial : it would be interesting to reload the windows sizes and positions (layer manager, display, terminal). it's already part of workspace, window position and their size should be reloaded. Martin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] georectify tool in wxGUI
Martin, updated at revision 35542 but the error remains. Perhaps am I doing something wrong ? Vincent Le jeudi 22 janvier 2009 à 20:15 +0100, Martin Landa a écrit : Hi, 2009/1/22 Vincent Bain b...@toraval.fr: self.parent.goutput.WriteLog(_('POINTS file %s saved') % self.file['points']) File /usr/local/grass-6.4.0svn/etc/wxpython/gui_modules/g output.py, line 235, in WriteLog line += diff * ' ' UnboundLocalError : local variable 'diff' referenced before assignment should be fixed in r35541. Martin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] georectify tool in wxGUI
Yes, now it works ! Thank you, Vincent. Le jeudi 22 janvier 2009 à 21:23 +0100, Martin Landa a écrit : Hi, 2009/1/22 Vincent Bain b...@toraval.fr: Martin, updated at revision 35542 but the error remains. Perhaps am I doing something wrong ? in relbr64 fixed in r35543. Martin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] d.extract strange behavior
Hi list ! Following upon a previous post concerning multiple selections on vector maps, I explored the capabilites of d.extract. I have a vector map1 linked to a pg database. When running : d.extract in=map1 out=map2 I select a set of boundaries (or centroids, or anything else). Map2 is properly created with the previously pointed boundaries, but when I look into the newly created linked table, it contains ALL the records of the source map, not only the few selected ones. Another troubling thing is the way selections operate : some objects cannot be catched, others are selected though falling out of the box... Hope it can help for future enhancement, Bye, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] grass startup cmd line
Hi list, as I'm quite a lazy boy, and am currently working on a single grass project, I would appreciate to launch grass from a terminal so that it starts with a given GUI, a given gisbase/location/mapset and a given workspace... This command grass64 -wxpython my_gisbase/my_location/my_mapset works, but is it possible to specify a default workspace to load i.e. grass64 -wxpython my_gisbase/my_location/my_mapset my_workspace.gxw Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass startup cmd line
Thank you, and I fear it is not possible to launch this command in the child process from the parent (e.g. grass64 -wxpython ... g.gui workspace=my_workspace.gxw) (OK, I purged the former grass-list address) Vincent Le mardi 10 février 2009 à 10:56 +0100, Martin Landa a écrit : Hi, 2009/2/10 Vincent Bain b...@toraval.fr: as I'm quite a lazy boy, and am currently working on a single grass project, I would appreciate to launch grass from a terminal so that it starts with a given GUI, a given gisbase/location/mapset and a given workspace... This command grass64 -wxpython my_gisbase/my_location/my_mapset works, but is it possible to specify a default workspace to load i.e. grass64 -wxpython my_gisbase/my_location/my_mapset my_workspace.gxw this is currently possible only from GRASS cmd, e.g. g.gui workspace=my_workspace.gxw (assuming that wxpython is default GUI) Martin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] g.findfile
Hi, Using grass6.4.0svn (revision 35838), the result of g.findfile appears not to run properly for raster query : create a raster map test, run g.findfile element=raster file=test. Looks like the file does not exist : name= mapset= fullname= file= (Tested successfully with a vector map) What am I doing wrong ? I remember the command worked in previous versions. Thanks, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] g.findfile
Thank you Martin ! (for my info, is it a recent change or was the previous keyword really 'raster' ?) Vincent Le mardi 10 février 2009 à 14:39 +0100, Martin Landa a écrit : Hi, 2009/2/10 Vincent Bain b...@toraval.fr: Using grass6.4.0svn (revision 35838), the result of g.findfile appears not to run properly for raster query : create a raster map test, run g.findfile element=raster file=test. Looks like the file does not exist : name= mapset= fullname= file= (Tested successfully with a vector map) What am I doing wrong ? I remember the command worked in previous versions. use element=cell to find raster map Martin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass startup cmd line
Thank you Hamish, As you suggest, I wrote a shell script containing the instruction g.gui workspace=... After declaring the GRASS_BATCH_JOB variable, and launching grass in text mode, the wxpython GUI starts and my_workspace loads correctly, but the problem is I can't get the control on the process ; I guess grass is waiting for instructions from /within/ the shell script ! Otherwise I tried to pipe the command like this : echo 'g.gui gui=wxpython workspace=my_workspace.gxw' | grass64 -text my_gisbase/my_location/my_mapset There again gui and workspace load properly but the grass prompt does not appear. Anyway, it's a minor problem... but it could be interesting to improve the capabilities of this concept of workspace by allowing to launch grass only by e.g. double-clicking on a .gxw icon. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass startup cmd line
Hello Moritz, unfortunately, I tried it too, but I got bailed out, with several error statements that don't really spark off anything to me : Traceback (most recent call last): File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 70, in module import gui_modules.preferences as preferences File /usr/local/grass-6.4.0svn/etc/wxpython/gui_modules/preferences.py, line 799, in module globalSettings = Settings() File /usr/local/grass-6.4.0svn/etc/wxpython/gui_modules/preferences.py, line 488, in __init__ self.ReadSettingsFile() File /usr/local/grass-6.4.0svn/etc/wxpython/gui_modules/preferences.py, line 558, in ReadSettingsFile mapset_file = os.path.join(gisdbase, location_name, mapset_name, self.fileName) File /usr/lib/python2.5/posixpath.py, line 60, in join if b.startswith('/'): AttributeError: 'NoneType' object has no attribute 'startswith' I did not go further in this direction, am I wrong ? Thank you, Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] db.select and SELECT COUNT query
Hi, has anybody ever had this kind of message ? Working on grass6.4.0svn, considering a database connection through the pg driver ; in order to count records matching an attribute (attr) value I execute this command on a table (table1) linked to a vector map : echo select count(*) from table1 where attr=201 | db.select The result should be 3. And db.select returns : ATTENTION : column 'count' : type int8 (bigint) is stored as integer (4 bytes) some data may be damaged count 3 In the end, it returns the result, but what does this message mean ? (Of course the same instruction given at the psql prompt raises no special warning). Any idea ? Thank you in advance, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] db.select and SELECT COUNT query
Sorry for this /stupid/ question (might be time to go to sleep !) I realize it must be related to a truncature operated between pg and grass or something like it. As it was annoying in a recursive query script, I turned the command in this : echo select count(*) from table1 where attr=201\ | db.select 2/dev/null Good evening ;-) Vincent Le mercredi 11 février 2009 à 21:36 +0100, Vincent Bain a écrit : Hi, has anybody ever had this kind of message ? Working on grass6.4.0svn, considering a database connection through the pg driver ; in order to count records matching an attribute (attr) value I execute this command on a table (table1) linked to a vector map : echo select count(*) from table1 where attr=201 | db.select The result should be 3. And db.select returns : ATTENTION : column 'count' : type int8 (bigint) is stored as integer (4 bytes) some data may be damaged count 3 In the end, it returns the result, but what does this message mean ? (Of course the same instruction given at the psql prompt raises no special warning). Any idea ? Thank you in advance, Vincent. ___ 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] face to DEM?
Hello Adam, hope I understood what you mean to do. If I had to cope with your problem, I would import the source file as points, then run v.to.rast, and r.surf.nnbathy with the l interpolation method. The risk for this solution is the triangulation performed by r.surf.nnbathy be different from your source TIN... Good luck, Vincent. Le samedi 14 février 2009 à 16:42 -0800, Adam Dershowitz a écrit : I have some 3-d vector faces. They were defined in a text file like this: F 10 1 2 1 2 2 2 ... ... etc. I imported them into grass like this: v.in.ascii -zn input=faces.txt out=faces format=standard and all seems fine. I can see the faces, and if I click on points around the edges I can get a height as well. I can't figure out how to generate a raster DEM from these faces. Essentially I would like to know the approximate elevation at each raster point inside the faces. I tried v.to.rast but then I get Column parameter missing For v.to.rast3 I get Unable to get layer info for vector map I also tried v.extrude but I see all 0s (ie Number of areas: 0 etc.) Did I do something wrong on my import? Did I miss some other way to do the conversion? I want to be able to do some mapcalcs with the heights of an these faces compared to another raster map, but I can't figure out how to convert these faces to something that I can work with. Thanks much, --Adam ___ 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] error saving workspace
Hello, currently working on grass64.svn (v.35911), I cannot save my workspace anymore under wxpython GUI. An error window pops up with the following message : Writing current settings to workspace file failed (global name 'val' is not defined). Am I doing something wrong ? Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] error saving workspace
Indeed,:-) Thank you. Vincent Le mardi 17 février 2009 à 20:39 +0100, Martin Landa a écrit : Hi, 2009/2/17 Vincent Bain b...@toraval.fr: currently working on grass64.svn (v.35911), I cannot save my workspace anymore under wxpython GUI. An error window pops up with the following message : Writing current settings to workspace file failed (global name 'val' is not defined). fixed in r35912. Martin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.digit - new wiki page for migrators
Hi, it may not be the right place to develop proposals, but if I could piggyback on what Jan said, I think his suggestion is very pertinent to those who -- as I said in a recent previous post -- really produce vector data ! one has to be able to work with both hands (one on the keyboard, the other on the mouse or the graphic pad). It quickly becomes much more efficient than clicking around all the time ;-) Another enhancement I would really expect is the ability to have such shortcuts for dynamically panning and zooming, while digitizing. Sorry if I only can give suggestions, but few coding solutions... Thanks, Vincent Le mercredi 18 février 2009 à 14:53 +0100, Jan Hartmann a écrit : A question from a very old GRASS-hand (I have used it since we got Internet in Amsterdam in 1991, having downloaded and compiled it, with the complete X-Window system, from a 9600-baud telephone line. Had even to compile the gcc-compiler, since native AIX cc didn't like GRASS): In the old, pre-6 v.digit, all commands were given by keyboard shortcuts. That took some exercise, but once you had it in your fingers (literally), you could digitize much faster than with the present click-button interface. I would guess it isn't all too difficult to add the old keyboard shortcuts to the system. Any idea from the developers whether this could be done? I still have GRASS 5.7 running for digitizing only. Beats every other package. Jan Dr. J. Hartmann Department of Geography University of Amsterdam Micha Silver wrote: I've posted a new page on the wiki [1] entitled Digitizing Area Features. It's targeted at new GRASS users, migrating from other non-topological software. I'd appreciate any comments from more experienced users as to correctness, completeness, or clarity. Thanks, Micha [1] http://grass.osgeo.org/wiki/Digitizing_Area_Features ___ 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.digit - new wiki page for migrators
Not yet very familiar with grass wiki, but perhaps here : http://grass.osgeo.org/wiki/GRASS_6.4_Feature_Plan#Wishlist Vincent Le mercredi 18 février 2009 à 15:55 +0100, Jan Hartmann a écrit : Yes, exactly, especially the panning and zooming. If this is not the place to place suggestions for enhancements, where can I put them? Vincent Bain wrote: Hi, it may not be the right place to develop proposals, but if I could piggyback on what Jan said, I think his suggestion is very pertinent to those who -- as I said in a recent previous post -- really produce vector data ! one has to be able to work with both hands (one on the keyboard, the other on the mouse or the graphic pad). It quickly becomes much more efficient than clicking around all the time ;-) Another enhancement I would really expect is the ability to have such shortcuts for dynamically panning and zooming, while digitizing. Sorry if I only can give suggestions, but few coding solutions... Thanks, Vincent Le mercredi 18 février 2009 à 14:53 +0100, Jan Hartmann a écrit : A question from a very old GRASS-hand (I have used it since we got Internet in Amsterdam in 1991, having downloaded and compiled it, with the complete X-Window system, from a 9600-baud telephone line. Had even to compile the gcc-compiler, since native AIX cc didn't like GRASS): In the old, pre-6 v.digit, all commands were given by keyboard shortcuts. That took some exercise, but once you had it in your fingers (literally), you could digitize much faster than with the present click-button interface. I would guess it isn't all too difficult to add the old keyboard shortcuts to the system. Any idea from the developers whether this could be done? I still have GRASS 5.7 running for digitizing only. Beats every other package. Jan Dr. J. Hartmann Department of Geography University of Amsterdam Micha Silver wrote: I've posted a new page on the wiki [1] entitled Digitizing Area Features. It's targeted at new GRASS users, migrating from other non-topological software. I'd appreciate any comments from more experienced users as to correctness, completeness, or clarity. Thanks, Micha [1] http://grass.osgeo.org/wiki/Digitizing_Area_Features ___ 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 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] output tables
Which db driver do you use ? which variables do you need to export (elevation, length, etc.)? Have a look at v.to.db. Good luck, Vincent Le samedi 28 février 2009 à 09:26 +, Matthew Mulbrandon a écrit : Hello, I just want to output all my data in file tenn_county_all into a comma delimited file with all variables. I have done v.out.ascii format=standard it works but I only get x,y coord and the ID but I need all variables. -84.28851567|35.99395611|1 -84.25999586|36.01053534|2 -84.25146063|35.95974559|3 -84.24508926|36.01258551|4 -84.20325976|36.00941242|5 Thanks Matthew __ From: mulb...@hotmail.com To: grass-user@lists.osgeo.org Subject: bounding box grid and polygon layer Date: Wed, 28 Jan 2009 10:27:54 + Hello, I want to give census tracts categories values ( layer tenn_county_all) based on there centroids location in a grid layer that I have made (tenn_county_grid). The grid is nubered 1-42 and also as letters i.e. row A colunm G. I was planning to do something like if centroid of tract X in grid box 12 write out category 12 for that tract X. I don't want to do this by hand but I am not sure how or if I can query between two layers. I am not attached to any particular method just need something that will do the job. I looked up info on the web and found bounding box in v.edit so I think that is a way to go but have not seen anything on how to do this. Thanks, Matthew __ Windows Live™: E-mail. Chat. Share. Get more ways to connect. See how it works. __ Windows Live™ Hotmail®:…more than just e-mail. Check it out. ___ 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 update vector layers?
Hi Thybério, it is actually v.patch that you need to run for this purpose (AFAIU your problem). v.patch allows to merge several maps with attribute preservation (-e flag), but it only works for layer 1 features. I suggest you to check if each linked table has the same columns definition. Instead of creating a new merge map, you can choose the append file to existing file option (-a flag), then be carefull with the advice given in the man page : When using the -a flag, the user has to make sure that the features in the different maps added to the output map do not have overlapping category numbers, unless identical category numbers reflect identical attributes, otherwise the attributes of the added maps are lost. After you perform v.patch don't forget to clean your vector (e.g. removing duplicate features with v.clean) Good luck, Vincent Le jeudi 05 mars 2009 à 21:38 -0300, Thybério Luna Freire a écrit : I have a location with sugarcane areas. They are in a vector layer and each area have a line in a table in a pg database. There are new areas being planted in the farm and i have them in a shape file. I imported the new vectors and i want put them into the layer sugarcane to update it. My question is: how can i do it without interfering in the table structure and without to need create a new table? I tried v.overlay, but it adds new columns to the new table and a don't want i new table. Then, i tried v.patch, but this module is not apropriated to do this. It looks like very simple to do, but simply don't know how. Help me!! Thybério Luna - Brazil ___ 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 update vector layers?
If you are sure category values don't overlap, then you can use the -a flag. But in this case I would NOT give both arguments -a and -e; i.e. given that the sugarcane map already has a related table, the features from sugarcane_update will populate the table with new rows (it is the reason why cat numbers, which tag geometrical features, should be distinct from previous rows in the sugarcane table)... Else you can try to create a new map : v.patch -e in=sugarcane,sugarcane_update out=merge, If you suspect duplicate objets (e.g. adjacent areas boundaries), don't forget to run v.clean with the tool rmdupl. Vincent. Le samedi 07 mars 2009 à 16:03 -0300, Thybério Luna Freire a écrit : Hi Vincent, I've tried the following command: v.patch -ae input=sugarcane,sugarcane_update output=sugarcane --o it returns: ERROR: Output vector map sugarcane is used as input I'm using the -e and -a flags and overwriting the old vector layer, but the message is the same. I'm doing anything wrong? Both vector maps have exactly the same columns and there aren't identical attributes: Cat (integer) Cane_id(Integer) thanks, Thyberio ___ 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 update vector layers?
Having a look at the man page of the v.patch command, I was surprised by a semantic ambiguity : saying v.patch Create[s] a new vector map layer by combining other vector map layers sounds quite odd, shouldn't we say :Create a new vector map by combining other vector maps, given that only layer 1 is supported ? well, maybe there is no confusion and I misunderstand something ? In the current french translation of wxpython GUI, there are several confusions between these notions of maps and layers*, which IMHO should be clearly distinguished if we want to show a coherent vector data model. Bye, VB __ * I know I should keep working on translation files, but for the moment I'm really over-loaded. Le samedi 07 mars 2009 à 21:02 +0100, Markus Neteler a écrit : On Sat, Mar 7, 2009 at 8:03 PM, Thybério Luna Freire thybe...@gmail.com wrote: Hi Vincent, I've tried the following command: v.patch -ae input=sugarcane,sugarcane_update output=sugarcane --o it returns: ERROR: Output vector map sugarcane is used as input You have to write the output to a new map name. 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
Re: [GRASS-user] How update vector layers?
Thybério, as Markus said, you cannot proceed this way because your output map is already involved in the patch process. So using the -a option, you must not include the target map. Try : v.patch -a -e --overwrite input=sugarcane_update output=sugarcane To work safely regarding your database structure, I would suggest you to test on a copy of map sugarcane (again : particularly be sure of non-overlaping cat values...) Keep trying, you'll get it ! Vincent Le samedi 07 mars 2009 à 21:46 -0300, Thybério Luna Freire a écrit : The great problem in create a new map is that there is a database with rules of integrity. If the table sugarcane is linked to other table by primary/foreign keys i can't simply delete it and put the new table in your place. This is my problem. I have to mantain the table sugarcane and keep updating it, when necessary. thyberio ___ 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] import DWG files ?
Hi Grégoire, it is still possible to compile grass with the the option --with-opendwg, then specify where is located your ad2.h and other files needed by v.in.dwg. You have to subscribe to Open Design Alliance (ODA) to obtain the opendwg toolkit ; AFAIR you have to pay $250 to register (formerly free...). Yours, Vincent Le mercredi 01 avril 2009 à 17:50 +0200, Grégoire Pigeon a écrit : Hi all, is it still possible to import dwg files in grass6.2.3 I remember that I did that with grass6.0.1 and the OpenAlliance library. best regards, Grégoire ___ 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] concurrent versions of grass
Hi, till recently I used to work on a Hardy-Ubuntu GNU/linux. I was able to run various grass versions installed on the system. Now on Debian I cannot install let's say 640_rc2 if 6.5.svn is already installed. The error returned by the system deals with the replacement of file /usr/local/bin/gem6, which belongs to the other grass version and cannot be overwritten. Has anyone experienced that ? (am I explicit enough ?) BTW, I use checkinstall instead of make install. Bye, Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] concurrent versions of grass
Perhaps, but it could be interesting to check whether the error is reproductible, or it is due to a configuration mess on my system... Thank you, VB Le mercredi 01 avril 2009 à 22:18 +0200, Markus Neteler a écrit : On Wed, Apr 1, 2009 at 10:02 PM, Vincent Bain b...@toraval.fr wrote: Hi, till recently I used to work on a Hardy-Ubuntu GNU/linux. I was able to run various grass versions installed on the system. Now on Debian I cannot install let's say 640_rc2 if 6.5.svn is already installed. The error returned by the system deals with the replacement of file /usr/local/bin/gem6, which belongs to the other grass version and cannot be overwritten. Mhh, perhaps we have to rename gem6 - gem6x with x=0..n? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] concurrent versions of grass
In my case, what if I voluntarily remove gem6 before installing a second version of Grass ? will the previous grass install be damaged ? Le mercredi 01 avril 2009 à 22:28 +0200, Markus Neteler a écrit : On Wed, Apr 1, 2009 at 10:24 PM, Vincent Bain b...@toraval.fr wrote: Perhaps, but it could be interesting to check whether the error is reproductible, Yes: because in all GRASS 6 versions it is called gem6 (unfortunately). So it must conflict... If the others agree I'll start to add the second version number starting with GRASS 6.5. Then behaviour doesn't change for 6.4. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] concurrent versions of grass
Thank you, I'll try it VB Le mercredi 01 avril 2009 à 22:53 +0200, Markus Neteler a écrit : Unless you don't use GEM then you are fine to remove/rename it. Markus On Wed, Apr 1, 2009 at 10:33 PM, Vincent Bain b...@toraval.fr wrote: In my case, what if I voluntarily remove gem6 before installing a second version of Grass ? will the previous grass install be damaged ? Le mercredi 01 avril 2009 à 22:28 +0200, Markus Neteler a écrit : On Wed, Apr 1, 2009 at 10:24 PM, Vincent Bain b...@toraval.fr wrote: Perhaps, but it could be interesting to check whether the error is reproductible, Yes: because in all GRASS 6 versions it is called gem6 (unfortunately). So it must conflict... If the others agree I'll start to add the second version number starting with GRASS 6.5. Then behaviour doesn't change for 6.4. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] working on a mounted external disk
Hello, has anybody ever had problems accessing a gisdbase located on an external disk ? In my case, I can open a grass session in a distant location/mapset : grass65 /mnt/path_to_my_disk/my_location/topo From within this place, everything works well, OK. If I start grass from my home directory and try to reproject data from the distant directory : grass65 /home/vincent/my_gisdbase/my_location/my_mapset r.proj input=dem50 location=/mnt/path_to_my_disk/my_location mapset=topo output=dem50 Then an error occurs : ERREUR :Mapset topo in input location /mnt/path_to_my_disk/my_location/ - not found whereas the disk is properly mounted there and accessible... Any idea on why grass seems not to see the mounted directory ? Thank you Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] working on a mounted external disk
OH, that simple ! shame on me :-( Sorry for being so dazed... A detail concerning r.proj in that context of distinct gisdbases : if the distant location has the same name as the current (i.e. in my example both locations are called 'my_location'), then the command returns a comprehensible error : ERROR: You have to use a different location for input than the current It would be nice if the dbase argument was checked /before/ the location argument, in order to admit same location names, not necessarily meaning that we project data in an identical projection system. Vincent Le mercredi 08 avril 2009 à 13:41 +0200, Moritz Lennert a écrit : this should be: r.proj input=dem50 location=my_location mapset=topo dbase=/mnt/path_to_my_disk/ output=dem50 If you give it location=/mnt/path_to_my_disk/my_location it thinks that the location is called /mnt/path_to_my_disk/my_location as the error message indicates: ERREUR :Mapset topo in input location /mnt/path_to_my_disk/my_location/ - not found Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] working on a mounted external disk
Thank you Moritz, I tried to specify the path to the target database location this way : r.proj input=dem50 location=/mnt/path_to_my_disk/my_location mapset=topo dbase=/mnt/path_to_my_disk/ output=dem50 unfortunately it still returns the same error. Well, maybe I should change the title of this post, because I notice now that the problem is the same with resident paths on the machine (apparently nothing to deal with mounted directory). Vincent Le mercredi 08 avril 2009 à 10:48 +0200, Moritz Lennert a écrit : You need to use the dbase= option of r.proj to direct it to another database. In your case it is looking for a location called /mnt/path_to_my_disk/my_location/ in your current gisdatabase. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] working on a mounted external disk
In raster/r.proj/main.c, line 208 I wrote : if (strcmp(inlocation-answer, G_location()) == 0 !indbase-answer || strcmp(indbase-answer, G_gisdbase())== 0)) G_fatal_error(_(Input and output locations can not be the same)); G_get_window(outcellhd); Then compiled and installed again, but the error remains, and what's more puzzling to me is that the message string is different from what I have in the modified file (You have to use a different location for input than the current)... What have I done wrong ? VB Le mercredi 08 avril 2009 à 14:32 +0200, Moritz Lennert a écrit : On 08/04/09 14:13, Vincent Bain wrote: It would be nice if the dbase argument was checked /before/ the location argument, in order to admit same location names, not necessarily meaning that we project data in an identical projection system. Try this (untested): Index: raster/r.proj/main.c === --- raster/r.proj/main.c (révision 36609) +++ raster/r.proj/main.c (copie de travail) @@ -219,7 +219,7 @@ setname = imapset-answer ? imapset-answer : G_store(G_mapset()); -if (strcmp(inlocation-answer, G_location()) == 0) +if (strcmp(inlocation-answer, G_location()) == 0 (!indbase-answer || strcmp(indbase-answer, G_gisdbase()) == 0)) G_fatal_error(_(You have to use a different location for input than the current)); G_get_window(outcellhd); Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] working on a mounted external disk
Apparently a fix for this issue done by Markus a long time ago in grass6, but this got lost, when (IIUC) r.proj.seg became r.proj in grass7. OK, I've got it : I made changes to raster/r.proj/main.c instead of raster/r.proj.seg/main.c Now it works fine with the right main.c file ! Thank you, Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] error in wx GUI Orthorectification module
Hi, on grass65 (rev 36698), in the wx GUI, calling 'Imagery Ortho photo rectification' results in the freezing of the application. Closing the X monitor allows to escape from the error. Opening a monitor prior to the call results in the same freezing (what looks strange to me is that the command launches a new monitor X1 while X0 is already started). Unfortunately no message appears neither on the console, nor on the messages manager. On the other hand, no problem running i.ortho.photo from the command line. Any idea ? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Multiple to one question
Hello Adam, maybe another solution in this case would be a set of 2 tables : * one linking to the geometry, that is containing nothing but cat values, * another one, containing a cat column (related to the geometric table) and different data columns corresponding to your sampling. Does this help ? VB Le lundi 13 avril 2009 à 14:23 -0700, Adam Dershowitz a écrit : I am trying to set up a new project in Grass, and I have a question about the best approach. I have different vector locations, and at each one there were multiple samples taken. At the moment I have each sample as a row in a data base. My question is how best to put this data into a set of vector points. I believe that I can do it in either of two ways (of not others). 1) I can create a vector point at each location, then I think that I can have multiple cats for that object. So I think I can do cat=1,3,6 for a given location. Will that work OK? 2) I can just create different vector objects, that happen to be at the identical location, and have each one point to a different cat. If the above is not clear, here is a bit more detailed example. At location A there was a sample collected on 1/1 with a value of 2.1, on 2/2 with a value of 2.2 and on 3/3 with a value of 3.3 The above data is already 3 rows in a database. I want to be able to display data about point A (say, average value or things like that). Should I just create a vector point A and then do cat=1,2,3 or should I create 3 different vector points at A, each one having a different cat? Any guidance about the benefits or limitations each approach (or any other approach to consider) would be greatly appreciated. Thanks, --Adam ___ 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] Multiple to one question
Hmmm, not sure you understood what I suggested (or maybe I did not catch your point of view !). Let's consider you add two measures a week for point A (cat=1), and 1 measure monthly for point B (cat=2), then after one year you have : * in table1 (cat integer): 2 recors (cat=1 and cat=2) * in table2 (cat integer, mes float): 116 records (104 records with cat=1 and 12 records with cat=2). Where is there a problem for you ? Hope this helps, VB Le mardi 14 avril 2009 à 09:38 -0700, Adam Dershowitz, Ph.D., P.E. a écrit : On Apr 14, 2009, at 4:06 AM, Moritz Lennert wrote: On 14/04/09 08:37, Vincent Bain wrote: Hello Adam, maybe another solution in this case would be a set of 2 tables : * one linking to the geometry, that is containing nothing but cat values, * another one, containing a cat column (related to the geometric table) and different data columns corresponding to your sampling. I think that if all you want is calculate some means or similar across dates and then display the results, Vincent's solution is the easiest. But you could also use layers [1]: layer 1 = January round of sampling layer 2 = February round of sampling etc. You would have to give each point a category value in each layer (cf v.category) and then either create separate tables for each period linking each to one of the layers or at least create some obvious cat values (i.e. 100s for January, 200s for February, etc) and link on single table to all the layers, but with different cat values in each layer. Moritz [1] See Vector object categories and attribute management on http://grass.osgeo.org/grass64/manuals/html64_user/vectorintro.html for a quick introduction Thanks, But, the problem with both of these approaches, columns, and layers, (Vincent or Moritz version) is that I don't have consistent times for each site. So, at site A I might have 5 samples, once a month and at site B I have 2 samples, one each year, and site C I have a few spread over a few years. So both of those approaches essentially need to have a column, or layer, for each possible time of sampling. But that is not really appropriate for the quasi-random times of the samples. Does this help ? VB Le lundi 13 avril 2009 à 14:23 -0700, Adam Dershowitz a écrit : I am trying to set up a new project in Grass, and I have a question about the best approach. I have different vector locations, and at each one there were multiple samples taken. At the moment I have each sample as a row in a data base. My question is how best to put this data into a set of vector points. I believe that I can do it in either of two ways (of not others). 1) I can create a vector point at each location, then I think that I can have multiple cats for that object. So I think I can do cat=1,3,6 for a given location. Will that work OK? 2) I can just create different vector objects, that happen to be at the identical location, and have each one point to a different cat. If the above is not clear, here is a bit more detailed example. At location A there was a sample collected on 1/1 with a value of 2.1, on 2/2 with a value of 2.2 and on 3/3 with a value of 3.3 The above data is already 3 rows in a database. I want to be able to display data about point A (say, average value or things like that). Should I just create a vector point A and then do cat=1,2,3 or should I create 3 different vector points at A, each one having a different cat? Any guidance about the benefits or limitations each approach (or any other approach to consider) would be greatly appreciated. Thanks, --Adam ___ 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] traveling salesman problem in air
Sorry if I'm beside the point, but I insist, being myself interested in this topic... My suggestion of a delaunay triangle net to join polygons centroids seems to be totally beside the point :-( I would just like to know why. Probably sth I did not catch ? Thank you VB Le mercredi 15 avril 2009 à 18:04 +0200, Moritz Lennert a écrit : On 15/04/09 09:46, Martina Schäfer wrote: Yes, the centroids should be fine. it actually did not work to run v.net.visibility when I tried but I'll give it another try with the centroids. Keep us posted. This is an interesting application of these modules... Moritz ___ 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] Rotate a map display
Jérôme, If you need to warp a map for a strict display purpose, IMHO it's not to be performed from within grass given that this operation makes no sense geographically. Don't know what the context is, but if you just have to rotate an image output of a map maybe you'd better look towards image manipulation tools (e.g. imagemagick and the -rotate option, which can easily be integrated in a script process). Let us know a bit more about your needs ! Vincent Le mercredi 15 avril 2009 à 10:22 +0200, Jérôme - GeoRezo.net a écrit : Hello, I want to make a map which follow a mobile. So I will have a set of raster maps and a polyline which represent a the path. One main feature is that map could follow the mobile's route. The y-axis could be the mobile'direction. My question is how can I do this. Best regards, ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] import dxf 3d polylines
Hi list, I need to import a dxf file containing 3d contourlines. I was convinced v.in.dxf was able to do that, unfortunately, it seems it does not. My contour lines are imported flat, ie with elevation set to 0... and no specific field provided with altitude. Reading the man page I see one can import 3d faces, but nothing about 3d polylines. Has anybody else dealt with this situation ? would there be an alternative such as multiple formats transformations with tools like Meshlab, Paraview, blender...? Thank you, Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] import dxf 3d polylines
Hello Nick, the thing is v.in.dxf does not propose this option. VB Le samedi 18 avril 2009 à 15:08 -0500, Nick Cahill a écrit : Did you import with the create 3d vector checked? I haven't used the v.in.dxf recently, but importing Microstation .dgn files, I'm always forgetting to check that box. Cheers, Nick Cahill On Apr 18, 2009, at 2:08 PM, Vincent Bain wrote: Hi list, I need to import a dxf file containing 3d contourlines. I was convinced v.in.dxf was able to do that, unfortunately, it seems it does not. My contour lines are imported flat, ie with elevation set to 0... and no specific field provided with altitude. Reading the man page I see one can import 3d faces, but nothing about 3d polylines. Has anybody else dealt with this situation ? would there be an alternative such as multiple formats transformations with tools like Meshlab, Paraview, blender...? Thank you, Vincent ___ 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] import dxf 3d polylines
Yes, I use grass65 (rev 36785). Now I remember I formerly used to compile grass with opendwg ; I would possibly perform this kind of import with v.in.dwg on .dwg files. (Anyway, now it's over because last year I gave up paying for my ODA membership). You can try importing this 3d dxf file : http://www.toraval.fr/telec/TEST.DXF (layers containing contourlines are named cl_m and cl_s) Vincent Le dimanche 19 avril 2009 à 10:09 +0200, Markus Neteler a écrit : You didn't say which GRASS version you use. There were recently fixes implemented in v.in.dxf. http://trac.osgeo.org/grass/wiki/Release/6.4.0RC4-News Perhaps you can also point us to a data set to try? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] import dxf 3d polylines
It works ! thank you very much Huidae, VB Le dimanche 19 avril 2009 à 23:35 -0400, Huidae Cho a écrit : Vincent, It's fixed in 36804. Huidae On Sun, Apr 19, 2009 at 10:39:31AM +0200, Markus Neteler wrote: Vincent, I can reproduce the problem. Please file it as bug in trac: http://trac.osgeo.org/grass/ including the relevant info. Markus On Sun, Apr 19, 2009 at 10:35 AM, Vincent Bain b...@toraval.fr wrote: Yes, I use grass65 (rev 36785). Now I remember I formerly used to compile grass with opendwg ; I would possibly perform this kind of import with v.in.dwg on .dwg files. (Anyway, now it's over because last year I gave up paying for my ODA membership). You can try importing this 3d dxf file : http://www.toraval.fr/telec/TEST.DXF (layers containing contourlines are named cl_m and cl_s) Vincent Le dimanche 19 avril 2009 à 10:09 +0200, Markus Neteler a écrit : You didn't say which GRASS version you use. There were recently fixes implemented in v.in.dxf. http://trac.osgeo.org/grass/wiki/Release/6.4.0RC4-News Perhaps you can also point us to a data set to try? 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
[GRASS-user] wx GUI profile analysis tool
Hello all, the profile analysis tool available at wx map display is definitely usefull. As a suggestion, it would be more powerfull if we could set, in the settings dialog box, an optional fixed ratio scale (1:1). In the case raster values are homogenous to distances (elevation values) then topographic transects would represent genuine slope values. Is it feasible ? Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] wx GUI profile analysis tool
Hello Moritz, IIUC what you want, I think you can approach this by setting the min and max values of both axes to the same values in the Profile plot settings and resizing the plot window in order to have a square grid. That's it, I agree it's a matter of comfort, but... I am lazy ! ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] problem with multiple selection on layer manager
Hi, trying to drop a set of maps from my display, I noticed the multiple map selection is not possible within the layer manager window through Ctrl +click or Shift+click. But it is possible to perform from the keyboard : Ctrl+up/down arrows. Oh, I forgot : I work on grass65 svn version. Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] wx python memory allocation
Hello list, for I've been testing grass70 python interface I've noticed a crash frequently related to memory consumption : the python process requires an /increasing/ amount of memory while e.g. wandering on the map display. The problem is obvious when the display copes with raster data. After several zooming/panning on a 10 000 000 cells image, memory allocated to the python process quickly runs beyond 200/300 Mo, and never decreases eventhough I set the display back to nothing... Could it be related to anything like a pyramid building process that does not work properly ? is the problem already well known or is it just normal, regarding the size of the displayed raster map ? Thank you in advance, Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GRASS project movement
Perhaps useful to raise the particular status of db data : if you use dbf driver, the directory will be saved, but if you have an external db connection, beware of data loss... Le vendredi 11 septembre 2009 à 12:42 +0200, Markus Neteler a écrit : On Fri, Sep 11, 2009 at 9:11 AM, christian Brandt chrisbrand...@hotmail.com wrote: Dear list, I have a question in terms of GRASS project movement. I would like to archive GRASS projects so that I can still use it but hence it does not take needed space. So if I package the complete GRASS GIS data directory does it also take the total information including projection, coordinate system... or do I have to start with a new project? (I did not use GRASS for a long time, so...) You may want to package per location. Locations are independent from each other (so just put each location directory in a different package). 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
[GRASS-user] db.execute transactions
Hi List, working on grass65 svn, with the pg dbms driver, I am trying to run several slq statements from a shell script. These statements are generated within a loop and I would like them to be executed in a transaction. I see it might not be possible this way : echo BEGIN TRANSACTION; | db.execute echo myfirststatement; | db.execute echo mysecondstatement | db.execute ... echo COMMIT; | db.execute as it returns WARNING: there is no transaction in progress Should I dump all the commands in a temporary file $tmp then run db.execute input=$tmp ? or is there a better way ? Thank you, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] grass70 and display monitor
Dear grass users, my question could be insane but I am wondering if there is an equivalent command to d.mon in grass70. In several man pages I saw a reference to d.frame but it seems not be implemented yet. Or will grass70 give up with x monitors ? considering one still can launch grass in text mode, I guess there might be a way to open a graphical device. thank you for your commentaries on this question, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass70 and display monitor
Thank you Michael, I found the thread... Yours, Vincent Le mercredi 02 décembre 2009 à 10:37 -0700, Michael Barton a écrit : There is no use of xmonitors in GRASS 7. There is a recent and ongoing discussion of alternatives on the Dev list. Look for the Shell scripts thread. Michael C. Michael Barton Director, Center for Social Dynamics Complexity Professor of Anthropology, School of Human Evolution Social Change Arizona State University Phone: 480-965-6262 Fax: 480-965-7671 www: www.public.asu.edu/~cmbarton, http://csdc.asu.edu On Dec 2, 2009, at 10:34 AM, grass-user-requ...@lists.osgeo.org wrote: Date: Wed, 02 Dec 2009 18:22:15 +0100 From: Vincent Bain b...@toraval.fr Subject: [GRASS-user] grass70 and display monitor To: GRASS user list grass-user@lists.osgeo.org Message-ID: 1259774535.11312.12.ca...@vincent-desktop Content-Type: text/plain Dear grass users, my question could be insane but I am wondering if there is an equivalent command to d.mon in grass70. In several man pages I saw a reference to d.frame but it seems not be implemented yet. Or will grass70 give up with x monitors ? considering one still can launch grass in text mode, I guess there might be a way to open a graphical device. thank you for your commentaries on this question, Vincent. ___ 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] grass70 and display monitor
Thank you for your answers and advices. Glynn, with your method based on displaying grass files on the fly in a viewer, I guess all d.* interactive commands will fail, won't they ? Deep down I want to know how much I will have to adapt some shell scripts I wrote for grass6* when my users migrate to grass70. My feeling on this point with care for durability of the application - and I think grass is mostly used by people in search of customization - is it would be nice that be maintained these fundamentals... To my mind it is one of the most important benefits of open source apps. On the other hand I am not aware enough of the constraints that lead to drop x monitors. Bye, Vincent Le jeudi 03 décembre 2009 à 06:00 +, Glynn Clements a écrit : Vincent Bain wrote: my question could be insane but I am wondering if there is an equivalent command to d.mon in grass70. In several man pages I saw a reference to d.frame but it seems not be implemented yet. Or will grass70 give up with x monitors ? considering one still can launch grass in text mode, I guess there might be a way to open a graphical device. GRASS 7.0 does not support monitors. The various display drivers are libraries against which the display library is linked. There are no persistent monitor processes as in 6.x, meaning that there is no persistence of state between d.* commands, so no d.frame, d.font etc. All rendering parameters are set via environment variables. Some of these are listed in the variables manual page, while others are listed in the manual pages for the various drivers. You can approximate the pre-7.0 workflow using an image viewer which automatically refreshes the display whenever the file changes. For X, you can use the ximgview program, e.g.: export GRASS_PNGFILE=map.bmp d.erase ximgview export GRASS_PNG_MAPPED=TRUE export GRASS_PNG_READ=TRUE # more d.* commands ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass70 and display monitor
Hi John and others, (and for easy reference as well to whom it may concern) I think the discussion that best corresponds to my question is here : http://lists.osgeo.org/pipermail/grass-dev/2009-August/045399.html Vincent. Le mercredi 02 décembre 2009 à 10:47 -0800, John C. Tull a écrit : And for easy reference, here is the thread: http://lists.osgeo.org/pipermail/grass-dev/2009-November/047135.html Cheers, John On Dec 2, 2009, at 9:46 AM, Vincent Bain wrote: Thank you Michael, I found the thread... Yours, Vincent Le mercredi 02 décembre 2009 à 10:37 -0700, Michael Barton a écrit : There is no use of xmonitors in GRASS 7. There is a recent and ongoing discussion of alternatives on the Dev list. Look for the Shell scripts thread. Michael C. Michael Barton Director, Center for Social Dynamics Complexity Professor of Anthropology, School of Human Evolution Social Change Arizona State University Phone: 480-965-6262 Fax: 480-965-7671 www: www.public.asu.edu/~cmbarton, http://csdc.asu.edu On Dec 2, 2009, at 10:34 AM, grass-user-requ...@lists.osgeo.org wrote: Date: Wed, 02 Dec 2009 18:22:15 +0100 From: Vincent Bain b...@toraval.fr Subject: [GRASS-user] grass70 and display monitor To: GRASS user list grass-user@lists.osgeo.org Message-ID: 1259774535.11312.12.ca...@vincent-desktop Content-Type: text/plain Dear grass users, my question could be insane but I am wondering if there is an equivalent command to d.mon in grass70. In several man pages I saw a reference to d.frame but it seems not be implemented yet. Or will grass70 give up with x monitors ? considering one still can launch grass in text mode, I guess there might be a way to open a graphical device. thank you for your commentaries on this question, Vincent. ___ 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] grass70 and display monitor
Thanks for your reply, May I be considered bothersome, but was it really impossible to achieve (like in grass65) the development of a wxGUI (which I really enjoy too) while keeping the complete set of historical command line environment ? Le jeudi 03 décembre 2009 à 14:40 +, Glynn Clements a écrit : Vincent Bain wrote: Glynn, with your method based on displaying grass files on the fly in a viewer, I guess all d.* interactive commands will fail, won't they ? There are no interactive d.* commands in 7.0. The 7.0 display architecture doesn't have any facility to query a mouse. Any modules which depended upon such functionality have either been removed entirely or have had the corresponding functionality removed. Deep down I want to know how much I will have to adapt some shell scripts I wrote for grass6* when my users migrate to grass70. My feeling on this point with care for durability of the application - and I think grass is mostly used by people in search of customization - is it would be nice that be maintained these fundamentals... To my mind it is one of the most important benefits of open source apps. This isn't going to change. If you want interaction, you need to either extend the GUI, or leverage existing functionality (e.g. by using the digitiser to create a vector map or the georectifier to create GCPs, then using this data as input). On the other hand I am not aware enough of the constraints that lead to drop x monitors. The interactive features made it impractical to build a decent GUI on top of the display commands. Additionally, the monitor-based architecture meant that it took an excessive amount of work to implement relatively minor improvements, as well as making debugging difficult. Oh, and none of it worked on Windows. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass70 and display monitor
Le vendredi 04 décembre 2009 à 14:16 -0700, Michael Barton a écrit : Markus, This is helpful. Much more so than simply those that ask 'why can't we do things the way we did'. Michael, as far as I am concerned by your remark, I just wish to distinguish my attitude from that of a potential pure fastidious/demanding consumer : I just keep very cautious from giving advices when I feel like being a too modest contributor... To finish off with my contribution, I initially bewailed the loss of some functionalities that straightly threatened some home-made pieces of code I use daily and intensively. I rencently worked for a customer who needed to upgrade a series a ArcInfo AML routines I wrote for him some years ago for AI7.2 and which did not work properly now on AI9. I praised him that developing his future personalized solutions on GRASS would be a guarantee of long-lasting, reliabilty, and so on... oops ! Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass70 and display monitor
Hamish, you're right, anyone can still freely use former but stable versions of Grass (like Jan does with grass5 for digitizing) for his particular purposes. Looking towards future I just figure out it's high time I began learning python ;-) Anyway, thank you all for this interesting discussion, Vincent. Le samedi 05 décembre 2009 à 02:51 -0800, Hamish a écrit : Vincent wrote: I praised him that developing his future personalized solutions on GRASS would be a guarantee of long-lasting, reliabilty, and so on... oops ! backwards compatibility since grass 6.0.0 in March 2005 is no mistake. it suggests that any code written for grass 7 will work for a similar time frame. and there is nothing stopping anyone from running grass 6.x well into the future, or continue with grass 5.4 for that matter. to limit ourselves from any change would guarantee obsolescence. the great trick is to compress all that good change into a single instant which are far appart, rather than to allow continual change with no stability. regards, Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass70 and display monitor
Thank you Michael and others, I really appreciate the time you allow to clarify these points. I am convinced the discussion will have been benefic to many other users. Yours, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] DEM creation from stereo pairs
Michael, don't know if it matches your needs as it is not an automatic process but I remember a long time ago I used e-foto, formerly only available on windows. It was a very rough photogrammetry station which allowed to digitize contour lines on the basis of stereo couples. Nowadays I see the project is still alive and was ported to Qt : http://www.efoto.eng.uerj.br Hope it can help, Bye, Vincent Le vendredi 11 décembre 2009 à 21:02 -0700, Michael Barton a écrit : Can anyone recommend open source tools for creating a DEM from stereo image pairs? Thanks Michael Barton C. Michael Barton Director, Center for Social Dynamics Complexity Professor of Anthropology, School of Human Evolution Social Change Arizona State University Phone: 480-965-6262 Fax: 480-965-7671 www: www.public.asu.edu/~cmbarton, http://csdc.asu.edu ___ 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] DEM creation from stereo pairs
Compiling from Qt seems quite easy on Linux : after unzipping the source code archive, from within each module directory, you just have to run qmake-qt4 to compile .pro files, then make install in order to create (local) executables. I tried with 'exterior' module, it works fine. Vincent. Le samedi 12 décembre 2009 à 08:31 -0700, Michael Barton a écrit : Thanks. I'm looking into this. To run on Mac and Linux, one needs to compile this from within QT. I haven't done that and don't know what is required. Michael C. Michael Barton Director, Center for Social Dynamics Complexity Professor of Anthropology, School of Human Evolution Social Change Arizona State University Phone: 480-965-6262 Fax: 480-965-7671 www: www.public.asu.edu/~cmbarton, http://csdc.asu.edu On Dec 12, 2009, at 2:17 AM, Vincent Bain wrote: Michael, don't know if it matches your needs as it is not an automatic process but I remember a long time ago I used e-foto, formerly only available on windows. It was a very rough photogrammetry station which allowed to digitize contour lines on the basis of stereo couples. Nowadays I see the project is still alive and was ported to Qt : http://www.efoto.eng.uerj.br Hope it can help, Bye, Vincent Le vendredi 11 décembre 2009 à 21:02 -0700, Michael Barton a écrit : Can anyone recommend open source tools for creating a DEM from stereo image pairs? Thanks Michael Barton C. Michael Barton Director, Center for Social Dynamics Complexity Professor of Anthropology, School of Human Evolution Social Change Arizona State University Phone: 480-965-6262 Fax: 480-965-7671 www: www.public.asu.edu/~cmbarton, http://csdc.asu.edu ___ 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] DEM creation from stereo pairs
and don't know what is required. Oh, on Debian-based distros, qmake-qt4 should be in the libqt4-dev package Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] ps.map Advice
If you really need to have full control on your virtual printer output, have a look at gs arguments : http://www.ghostscript.com/doc/7.07/Use.htm you might calculate what resolution gives you 300dpi by measuring the width of the printed map box in inches and multiplying by the desired dpi. any thing more than that is wasted disk space, Hamish Well, in the context of high quality press output it's a bit more delicate ; depends on what image you produce (bw, grayscale, color), what printing device it is (inkjet printer, laser printer, offset press... each having different lineature values and different frame types), and paper quality. Then, as you say, 300 dpi : * is most often enough for common color map printing * should be increased to 600 dpi for bw images e.g. text+thin black contour lines on a white bg to be sent on a laser printer, andbut should be considered the least acceptable resolution for high quality paper production. I agree that it's useless beyond 600 dpi. Yours, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] ps.map and line style
Hello, in a map composition I need to draw contours with two particular linestyles: - a line bedecked with crosses: xxxx - a double line: = I am wondering if there is a simple way to specify e.g. a postscript linear pattern that applies to vlines, the same way we can with vareas. Any idea ? Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] ps.map and line style
Hi John, thanks for the trick, I'll give it a try, though I need two black lines with /transparency/ in between. Vincent. Le mardi 12 janvier 2010 à 10:10 +, John A Stevenson a écrit : Vincent Bain wrote: Hello, in a map composition I need to draw contours with two particular linestyles: - a line bedecked with crosses: xxxx - a double line: = Hi Vincent, I have plotted a double line like this: vlines my_lines color yellow width 1.5 label Chirp lines lpos 4 style 000111 end vlines my_lines color black width 3 label Chirp lines lpos 4 style 000111 end The top entry in the file plots on top in the map. This gives a black dotted line with a yellow middle. Specifying the same legend position recreates the pattern in the legend. Later John ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] ps.map and line style
Yes Hamish, this is what I was planning to do in the end. If I go this way, I'll probably have to implement an 'angle' attribute in order to align my point symbols on the line direction... Will give it a try too, thank you, Vincent. Le mardi 12 janvier 2010 à 02:57 -0800, Hamish a écrit : Vincent wrote: in a map composition I need to draw contours with two particular linestyles: - a line bedecked with crosses: xxxx I suspect for that you might have to use the v.lrs modules or v.to.points to make points along the line at a regular interval, then display them over the top with a vpoints command. Perhaps there is another way though. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] working on Oracle Spatial Data
Hello list, could anyone help me finding ressources on the way Grass should be built, configured, etc. in order to connect to data stored in an Oracle database ? I wonder if gdal by default compiles with Oracle support, how one must set the ODBC connection, which grass environment variables have to be specifically set, how concurrent access is managed by v.external during layers edition, and so on... Thanks, I keep grabing the web on my side! Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] working on Oracle Spatial Data
Ok, yet it's part of a groundwork for me. Going further I'll try to feed the wiki page... As a basis, I found a doc page on Oracle website, written by Helena Markus : http://www.oracle.com/technology/pub/articles/mitasova-grass.html Vincent. Le mardi 19 janvier 2010 à 01:17 -0800, Hamish a écrit : Vincent wrote: could anyone help me finding ressources on the way Grass should be built, configured, etc. in order to connect to data stored in an Oracle database ? I wonder if gdal by default compiles with Oracle support, ( check for support with `ogrinfo --formats` ) how one must set the ODBC connection, which grass environment variables have to be specifically set, how concurrent access is managed by v.external during layers edition, and so on... there is a place-holder for that wiki page, but so far it is mostly empty. http://grass.osgeo.org/wiki/Data_formats#Oracle http://grass.osgeo.org/wiki/Oracle_Spatial If you can help document whatever you find and add it to that page :) Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] to separate a polygon
AFAIU : * your network was digitized as a single polygon, contouring streets/roads and crossroads in one go ; * you would like to turn it into a linear reference network. Am I right ? If no special splitting was initially performed on boundaries to distinguish between connected street sections I fear you'll have to do things manually... Vincent Le mardi 09 février 2010 à 16:03 +0100, Yvan Paillé a écrit : My road network was digitized into a polygon. Initially, I need to separate into two distinct parts by half. In a second time, I need to cut each polygon based on the streets and of belonging to adjacent parcels. As if I would define a sidewalk. Only the first step is already causing me problems. Thank you Yvan Vincent Bain a écrit : Yvan, could you please state your intention ? and clarify how your network is represented : are roads delimited by polygon boundaries and crossroads by nodes ? or was the network digitized as a polygon (I mean giving roads thickness) ? First, about data type conversion, you can have a look at v.type. Yours, Vincent. Le mardi 09 février 2010 à 15:14 +0100, Yvan Paillé a écrit : Hello, I would like to separate a polygon that represents a road network into several polygons or polylines. I do not know what functions to use to display vertices and nodes. Thank you in advance for answers that guide me. Yvan ___ 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] clean nodes in topology
Yvan, have a look at v.build.polylines Vincent. Le jeudi 18 février 2010 à 10:50 +0100, Yvan Paillé a écrit : Hello, I have a map with more polygons. The problem is that polygons have multiple nodes distributed randomly and no located at the intersecting polygons. Do you know a function to clean the topology without using manual editing? Thank you ___ 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] graphically selecting multiple features
Hello Kurt, may it help you, I wrote a /very draft/ add-on that could help you : http://grass.osgeo.org/wiki/GRASS_AddOns#v.selmany The name of the module should remind you of old AI stuff... it's a shell script, so not working with grass70, and only tested with postgresql database connection (it will probably fail with dbf driver). Grass has to be launched in text mode (type grass64 -text from your shell prompt), and there must be an open monitor (d.mon). Good luck ! Vincent. Le jeudi 18 février 2010 à 16:53 -0700, Kurtis Blaikie-Birkigt a écrit : Hello, I'm new to GRASS, coming from an Arc* background, and often find that I'm approaching things from the wrong direction. I'm using GRASS 6.4 with the wxPython GUI. I would like to graphically select multiple features from a large vector file and edit an attribute of those features. (I have a poorly classified polygon shapefile of hydrology and want to separate the rivers from the lakes). How can I go about this. With v.what.vect I seem to be able to only edit one feature at a time. Any suggestions will be appreciated. Thanks, Kurt ___ 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] Rectify an image
Nikos, as Maciek said, your image seems to need warping, so I fear you have no other solution than georectify it. Vincent Le lundi 22 février 2010 à 11:03 +, Nikos Dumakis a écrit : How can r.region put my image in place? Because it's a geostationary data so the image corners don't have the same coordinates (e.g. Top Left and Top Right don't have the same latitude). But I do have some GCP's for the image. So r.region doesn't help me (as far as I can see). Nikos On Sun, Feb 21, 2010 at 1:35 PM, Maciej Sieczka msiec...@sieczka.org wrote: Nikos Dumakis pisze: Greetings I have an image without any GCP's but, teorethically it's in WGS84. I have all pixel positions outside that file. My question is, what should I do in GRASS to rectify/project this image? I thought in two solutions If the image doesn't require warping, but is just missing a georereference, use r.region to put in place. Maciek -- Maciej Sieczka http://www.sieczka.org ___ 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] translucency and svg output
Pete, sorry if I'm off board : do you absolutely need to generate a /vector/ output ? Vincent. Le mercredi 24 mars 2010 à 20:15 -0700, pete davidson a écrit : On Wed, Mar 24, 2010 at 7:49 PM, pete davidson caiti...@gmail.com wrote: On Wed, Mar 24, 2010 at 7:42 PM, Glynn Clements gl...@gclements.plus.com wrote: pete davidson wrote: I'm producing maps in svg format using the cairo driver and d.out.file. One thing I'd really like to be able to do is set one or more of the vector layers to be partially translucent. After rummaging through this list, I'm beginning to wonder if its possible at all, since translucency in the guis seems to be handled by g.pnmcomp, but perhaps someone with more experience knows a way to do it. It isn't. The graphics API doesn't support an alpha component. Apart from anything else, it can't easily be implemented for PostScript output, as PostScript doesn't support translucency. Faking it with patterns is possible, but it starts getting complicated if you want more than one translucent layer. The only workaround I've managed to find so far is to export each layer as an individual svg, then use something like inkscape to rejoin them, setting the translucency/transparency of each layer as I do it. Which is fine for a map with only a few layers, but gets tedious when dealing with maps with dozens of layers. SVG is XML, so it shouldn't be especially hard to automate the process using generic XML-processing tools. Ohh, I hadn't thought of editing the svg/xml directly. Thanks for the suggestion, I'll see what I can achieve from there. And then we discover that the 'svg' produced by cairo using the above script is actually just a png in an svg wrapper.. At least for my particular install of Grass. That's frustrating. Ahh well, back to the drawing board. Peter ___ 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] grass 7 and ubuntu start up
Hi Matthew, just an idea : did you check if grass70 exists in your path and if so, is it executable (you can type sth like ls -l /usr/bin/grass* in a term) Yours, Vincent. Le jeudi 25 mars 2010 à 12:24 +, Matthew Mulbrandon a écrit : Hello, Just installed grass 7.0 using the latest snapshot binary march 20 (grass7.0.svn-i686-pc-linux-gnu-20_03_2010) and ubuntu 9.10. I type in grass70 in the shell and nothing happens. No error but also no gui. I try some gui commands and it says command not found. Nothing really works but the program seems to activate. I use grass 6.4 often but am at a loss as to what the problem is. Thanks, Matthew grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user __ The New Busy is not the old busy. Search, chat and e-mail from your inbox. Get started. ___ 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] db.select driver issue
Richard, FYI the problem may be located at the driver level: with 'pg' driver this query works fine. Bye, Vincent Le dimanche 28 mars 2010 à 08:54 +1100, Richard Chirgwin a écrit : In a previous thread, I mentioned trouble with Grass-GIS's database driver on complex queries, but couldn't remember an example. I've bumped into one this morning. In MySQL, this query works fine: select sum(population) from australia in Grass-GIS: echo select sum(population) from Australia | db.select -c returns: DBMI-DBF driver error: SQL parser error: syntax error, unexpected '(', expecting FROM processing '(' in statement: select sum(population) from Australia Error in db_open_select_cursor() Richard ___ 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] d.vect option for dotted or hatched line display
Hello list, is there anywhere a way to choose a linestyle (e.g. dots or hatches) for displaying vector maps ? Did not find such option on the d.vect man page. I know it is achievable in the ps module, but is it possible on the grass display ? Thank you, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] maps France and Germany
AFAIK, this data level in France is not free. The only free administrative data I know is available here : http://professionnels.ign.fr/DISPLAY/000/528/175/5281750/GEOFLADept_FR_Corse_AV_L93.zip as a shapefile. Or may you have a look at this address : http://download.geofabrik.de/osm/europe/ though I did not have time to visualize data. Good luck, Vincent Le mercredi 12 mai 2010 à 06:23 +0200, José Miguel Barrios a écrit : Hi list; Can someone indicate me where can I request simple vector layers of Germany and France at municipality level? Thanks! Miguel ___ 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] Really dumb script question
Hello Richard, just a suggestion, though I'm not a bash expert : k=10 while [ $k -le 250 ]; do i=0 while [ $i -lt $k ]; do image=`echo Background_$i.png` your stack of operations on $image; ((i++)); done k=$((k+10)); done hmmm, perhaps it's not clean at interval borders, but it may be a good beginning ? Good luck, Vincent. Le lundi 28 juin 2010 à 21:17 +1000, Richard Chirgwin a écrit : It's like this. I'm good with Grass-GIS and not good with complicated scripts. An r.in.wms process downloaded too many tiles, and broke r.in.gdalwarp. So, I would like to break up the downloaded tiles into groups of ten, run r.in.gdalwarp on the groups, then patch the grouped tiles together at the end. I have: Background_ _ 0 to 251.png ie - Background__0.png to Background__251.png What I want to do is for i in group of ten tiles; do r.in.gdalwarp input=grouped filenames output=subset method=nearest s_srs=EPSG:4326 done But how do I break the set of 251 tiles into groups for this? Richard ___ 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] Really dumb script question
Oops, checking for your reply, I notice an error in my code at line 3 : replace i=0 with i=$((k-10)) Yours, Vincent. Le lundi 28 juin 2010 à 14:06 +0200, Vincent Bain a écrit : Hello Richard, just a suggestion, though I'm not a bash expert : k=10 while [ $k -le 250 ]; do i=0 while [ $i -lt $k ]; do image=`echo Background_$i.png` your stack of operations on $image; ((i++)); done k=$((k+10)); done hmmm, perhaps it's not clean at interval borders, but it may be a good beginning ? Good luck, Vincent. Le lundi 28 juin 2010 à 21:17 +1000, Richard Chirgwin a écrit : It's like this. I'm good with Grass-GIS and not good with complicated scripts. An r.in.wms process downloaded too many tiles, and broke r.in.gdalwarp. So, I would like to break up the downloaded tiles into groups of ten, run r.in.gdalwarp on the groups, then patch the grouped tiles together at the end. I have: Background_ _ 0 to 251.png ie - Background__0.png to Background__251.png What I want to do is for i in group of ten tiles; do r.in.gdalwarp input=grouped filenames output=subset method=nearest s_srs=EPSG:4326 done But how do I break the set of 251 tiles into groups for this? Richard ___ 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] create a mask to exclude plate area
Bonsoir Étienne, Je crois comprendre que ton problème se rapporte à l'identification de zones sur des critères géomorphologiques. Peut-être pourrais-tu partir sur un traitement avec r.mapcalc ou des outils comme r.reclass pour lisser ta topo ou du moins discrétiser ta carte des pentes, afin d'identifier des zones homogènes. Puis transformer ces zones en polygones (r.to.vect) et faire des ratios entre périmètre et superficie... I guess your problem is related to a geomorphological identification process. I suggest you to try and make a kind of slopes reclassification in order to discretize and aggregate homogenous slope areas, then turn these areas into a polygon vector map. After that you can compute ratios between perimeters, areas, X/Y extension, and so on... Amuse-toi bien Enjoy it ! Vincent. Le lundi 19 juillet 2010 à 21:39 +0200, Etienne DELAY (deust 2010) a écrit : How? Which command(s) exactly? Using r.mapcalc? There is also r.mask [1] Also, tres-important, did you define the region of your interest using g.region[2]? Now I have to exclude the plate areas (eg trays Spanish) but I can not see any use for this test. I thought r.mapcalc place an order with a request to both area and on slope. That is right. r.mapcalc (or most of the raster modules) will ignore the MASKed-out cells by the raster map named MASK. Radom example grepped from the grass-user list: --%--- Re: How To use a raster MASK correctly? --17/10/2009 20:09 --%--- Figured it out. Used the following r.mapcalc expression: r.mapcalc new_hillshade=if(europe_dem 100, 180, europe_hillshade) --%--- Cordialement, Nikos --- [1] http://grass.osgeo.org/grass65/manuals/html65_user/r.mask.html [2] http://grass.osgeo.org/grass65/manuals/html65_user/g.region.html Hello Nikos sorry for the French, I need to think in my language :-) So I have no problem using r.mapcalc for 1st and 2nd criteria, and to assemble the two masks. The region is also parameterize, but now, it has become a 3rd test, I must arrive in isolating the flat portions of the mountain in Spain, if not that is reflected in part authorize my mask, but not know structural problems in they vineyards ... I want to use something like: r.mapcalc Maskplate= if((Area10ha slope=10%), 1, null()) mais l'utilisation de aire ne me satisfait pas vraiment parce quelle peut changer en longueur et largeur... but the use of area does not really satisfy me because what can change in length and width ... Is there a solution would be to make visual grade, to define the geomorphological unite? thank you very much ;-) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user