Re: [GRASS-user] GeoPackage and multiple tables
You might need to run VACUUM from gdal tools: ogrinfo -sql "VACUUM" pal.gpkg Regards Saber On Sun, 17 Mar 2019 at 17:16, Ken Mankoff wrote: > > On 2019-03-17 at 10:03 -0700, Ken Mankoff wrote... > > After learning from Markus & Moritz I see that the GPKG files are as > > they should be - vector features are not in triplicate. I think there > > are some strange QGIS import issues where it imports everything 3x but > > this is not a GRASS issue. > > Actually, I am still confused about something with the file sizes. The > data/metadata all appear correct. But if I export a three-layer vector w/ > three tables 1x, the GPGK file is 12 MB and contains 1 table. It is missing > two tables. If I then export each of the databases associated with each of > the layers, they are 0.2 MB. If export the vector 3x, once per layer, with > update, as advised earlier in this thread, the resulting GPKG is ~34 MB, or > ~3x as large as I would expect it to be. > > v.out.ogr -c input=pal output=pal_all.gpkg --o # 12 MB, but only layer1 > table > db.out.ogr input=pal output=l1.csv layer=1 # 7.5k > db.out.ogr input=pal output=l2.csv layer=2 # 124K > db.out.ogr input=pal output=l3.csv layer=3 # 85k > > v.out.ogr -c input=pal output_layer=tbl2_points output=pal.gpkg layer=1 --o > v.out.ogr -c input=pal output_layer=tbl2_areas output=pal.gpkg layer=2 -u > --o > v.out.ogr -c input=pal output_layer=tbl2_lines output=pal.gpkg layer=3 -u > --o > # pal.gpg is 34 MB, with 3 tables. > > Why does adding two additional (small) tables increase the size of the > resulting geopackage by 100s of times the size of the tables? > > -k. > > ___ > grass-user mailing list > grass-user@lists.osgeo.org > https://lists.osgeo.org/mailman/listinfo/grass-user -- Saber Razmjooei www.lutraconsulting.co.uk +44 (0)7568 129733 ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] 3D unstructured Grid
The purpose of MDAL is slightly different from 3D meshes such as PLY. In 3D mesh, you store just elevation (either face or corners) and some type of materials. But in MDAL, the aim is to store quantities on top of the mesh (e.g. a mesh for ground elevation and water depth on top over a period of time). So in essence, it is more like 2.5D (+time). Cheers Saber On Tue, 2 Oct 2018 at 13:40, Carlos Henrique Grohmann de Carvalho < carlos.grohm...@gmail.com> wrote: > That seems very interesting! Any plans to support PLY? > > Carlos > > On Tue, Oct 2, 2018 at 9:32 AM, Saber Razmjooei < > saber.razmjo...@lutraconsulting.co.uk> wrote: > >> Hi, >> >> Sorry to hijack the discussion related to GRASS. But, we are working a >> new library to handle mesh data. MDAL is currently supported in QGIS. But >> the idea is to make it generic so other OS GIS application can handle the >> mesh layer through a generic library similar to GDAL. Here is more info >> about MDAL: >> https://lutraconsulting.github.io/mdal-web/ >> >> You can find some example of the formats supported there too. >> >> Regards >> Saber >> >> >> On Tue, 2 Oct 2018 at 03:01, Vaclav Petras wrote: >> >>> On Sat, Sep 29, 2018 at 4:45 PM Francois Chartier < >>> fra.chart...@gmail.com> wrote: >>> >>>> Can Grass GIS create a 3D unstructured grid based on a 3D point data >>>> set ? >>>> >>> >>> All data and visualizations at the following website were created in >>> GRASS GIS including a 3D regular grid (3D raster in GRASS GIS) interpolated >>> from a set of 3D points (vector which is 3D, i.e. with z, in GRASS GIS). >>> Points were soil properties measurement at different depths, 3D raster is >>> then continuous representation of a selected property. >>> >>> https://geospatial.ncsu.edu/geoforall/multidim-geovis.html >>> >>> Module used for interpolation was v.vol.rst: >>> >>> https://grass.osgeo.org/grass74/manuals/v.vol.rst.html >>> >>> Vaclav >>> ___ >>> grass-user mailing list >>> grass-user@lists.osgeo.org >>> https://lists.osgeo.org/mailman/listinfo/grass-user >> >> >> >> -- >> Saber Razmjooei >> www.lutraconsulting.co.uk >> +44 (0)7568 129733 >> >> _______ >> grass-user mailing list >> grass-user@lists.osgeo.org >> https://lists.osgeo.org/mailman/listinfo/grass-user >> > > > > -- > Prof. Carlos Henrique Grohmann > Institute of Energy and Environment - Univ. of São Paulo, Brazil > - Digital Terrain Analysis | GIS | Remote Sensing - > > http://carlosgrohmann.com > http://orcid.org/-0001-5073-5572 > > Can’t stop the signal. > -- Saber Razmjooei www.lutraconsulting.co.uk +44 (0)7568 129733 ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] 3D unstructured Grid
Hi, Sorry to hijack the discussion related to GRASS. But, we are working a new library to handle mesh data. MDAL is currently supported in QGIS. But the idea is to make it generic so other OS GIS application can handle the mesh layer through a generic library similar to GDAL. Here is more info about MDAL: https://lutraconsulting.github.io/mdal-web/ You can find some example of the formats supported there too. Regards Saber On Tue, 2 Oct 2018 at 03:01, Vaclav Petras wrote: > On Sat, Sep 29, 2018 at 4:45 PM Francois Chartier > wrote: > >> Can Grass GIS create a 3D unstructured grid based on a 3D point data set ? >> > > All data and visualizations at the following website were created in GRASS > GIS including a 3D regular grid (3D raster in GRASS GIS) interpolated from > a set of 3D points (vector which is 3D, i.e. with z, in GRASS GIS). Points > were soil properties measurement at different depths, 3D raster is then > continuous representation of a selected property. > > https://geospatial.ncsu.edu/geoforall/multidim-geovis.html > > Module used for interpolation was v.vol.rst: > > https://grass.osgeo.org/grass74/manuals/v.vol.rst.html > > Vaclav > ___ > grass-user mailing list > grass-user@lists.osgeo.org > https://lists.osgeo.org/mailman/listinfo/grass-user -- Saber Razmjooei www.lutraconsulting.co.uk +44 (0)7568 129733 ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] QGIS 3D crowd-funding
Dear GRASS users, This is related to QGIS, but as there is a large overlap between GRASS and QGIS users, I thought this might be of your interest. QGIS 3 is out and amongst the exciting new features, we have introduced QGIS 3D thanks to the grant from QGIS.org We have decided to run a crowd-funding campaign with a comprehensive proposal to add more features: https://www.lutraconsulting.co.uk/crowdfunding/more-qgis-3d/ Any help from FOSS GIS community is highly appreciated. Regards Saber -- Saber Razmjooei www.lutraconsulting.co.uk ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] file formats from Ordinance Survey maps
Hi Jim, What OS data are you requesting? In general, it is best to have the files in GML as DWG is propietry format and needs extra effort to be read in GRASS. If you are working with OS MasterMap, the GML contains nested attributes and best to convert them to a standard format. We have developed a tool in QGIS to convert GML/GZ files directly to PostGIS (or Shapefile): http://www.lutraconsulting.co.uk/products/ostranslator-ii/ http://www.lutraconsulting.co.uk/products/ostranslator/ You can then use the data in GRASS easily. For other OS products in GML format, you can directly read them in GRASS as the attribute tables are less complex compared with OSMM. Hope it helps, Cheers, Saber From: grass-user [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of James Maas (MED) Sent: 19 January 2016 07:15 To: grass-user@lists.osgeo.org Subject: [GRASS-user] file formats from Ordinance Survey maps I'm attempting to get some maps and use them in GRASS, which gives me a choice of formats to download from. Can anyone tell me if any one of these formats is better and will allow direct import into GRASS 7.0.2. without transformation. For one type of topography map I can get either GML2, DWG or File Geodatabase. For another type of Terrain type map I can get either Shape, GML3 or DWG. I have googled and looked at the choices provided by Dr. Jim Maas Mathematical Modelling in Epidemiology of Infectious Disease Rm 2.13 MED Norwich Medical School University of East Anglia Norwich Research Park, NR4 7TJ U.K. E: j.m...@uea.ac.uk <https://mail.google.com/mail/?view=cm=1=1=j.m...@uea.ac.uk> P: +44 (0)1603 91325 -- This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error please notify the system manager. This message contains confidential information and is intended only for the individual named. If you are not the named addressee you should not disseminate, distribute or copy this e-mail. Please notify the sender immediately by e-mail if you have received this e-mail by mistake and delete this e-mail from your system. If you are not the intended recipient you are notified that disclosing, copying, distributing or taking any action in reliance on the contents of this information is strictly prohibited. Whilst reasonable care has been taken to avoid virus transmission, no responsibility for viruses is taken and it is your responsibility to carry out such checks as you feel appropriate. If this email contains a quote or offer to sell products, carry out work or perform services then our standard terms and conditions (which can be found at http://www.lutraconsulting.co.uk/downloads/Lutra%20Consulting%20Standard%20Terms%20and%20Conditions.pdf shall apply unless explicitly stated otherwise. Saber Razmjooei and Peter Wells trading as Lutra Consulting.___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] 64 bit Python in GRASS 6.4.2
Hi Andy, You can try this: http://download.osgeo.org/osgeo4w/osgeo4w-setup-x86_64.exe Cheers, Saber From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Orr, Andrew B. Sent: 26 September 2014 15:46 To: 'grass-user@lists.osgeo.org' Subject: [GRASS-user] 64 bit Python in GRASS 6.4.2 Hello all, I've been working on a r.cost and r.drain python script that is using very large rasters. I've been running into some memory issues with this script due to the size of the rasters. I read about setting the amount of memory that the map can be stored in with the percent_memory=xx argument for r.cost, however the performance of running the script is becoming a problem (it takes too long to run for our goal). I then started watching the Task Manager and noticed it was peaked out at less ~4 gigs of ram and I realized it's only a 32 bit version of Python installed with GRASS. My question is: is there a compatible 64 bit version of Python for GRASS 6.4.x, and if so, how does one install it? Thanks for your help! -Andy -- This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error please notify the system manager. This message contains confidential information and is intended only for the individual named. If you are not the named addressee you should not disseminate, distribute or copy this e-mail. Please notify the sender immediately by e-mail if you have received this e-mail by mistake and delete this e-mail from your system. If you are not the intended recipient you are notified that disclosing, copying, distributing or taking any action in reliance on the contents of this information is strictly prohibited. Whilst reasonable care has been taken to avoid virus transmission, no responsibility for viruses is taken and it is your responsibility to carry out such checks as you feel appropriate. If this email contains a quote or offer to sell products, carry out work or perform services then our standard terms and conditions (which can be found at http://www.lutraconsulting.co.uk/downloads/Lutra%20Consulting%20Standard%20Terms%20and%20Conditions.pdf shall apply unless explicitly stated otherwise. Saber Razmjooei and Peter Wells trading as Lutra Consulting.___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.lidar.growing.exe stops
Hi Simon, Do you get any error in the command window? Does your computer freeze or slowing down before the message? Cheers, Saber From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Crutchley, Simon Sent: 19 June 2014 12:21 To: 'grass-user@lists.osgeo.org' Subject: [GRASS-user] v.lidar.growing.exe stops I've created a file using v.lidar.edgedetection but when I run v.lidar.growing it runs for a while then I get a Windows message saying it has stopped working and forces closure. I get no error in GRASS. This is my command line - v.lidar.growing input=veg_edge2@PERMANENT output=lidar_edge_1 first=D149660_Petx@PERMANENT. Any suggestions? Simon This e-mail (and any attachments) is confidential and may contain personal views which are not the views of English Heritage unless specifically stated. If you have received it in error, please delete it from your system and notify the sender immediately. Do not use, copy or disclose the information in any way nor act in reliance on it. Any information sent to English Heritage may become publicly available. Portico: your gateway to information on sites in the National Heritage Collection; have a look and tell us what you think. http://www.english-heritage.org.uk/professional/archives-and-collections/por tico/ -- This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error please notify the system manager. This message contains confidential information and is intended only for the individual named. If you are not the named addressee you should not disseminate, distribute or copy this e-mail. Please notify the sender immediately by e-mail if you have received this e-mail by mistake and delete this e-mail from your system. If you are not the intended recipient you are notified that disclosing, copying, distributing or taking any action in reliance on the contents of this information is strictly prohibited. Whilst reasonable care has been taken to avoid virus transmission, no responsibility for viruses is taken and it is your responsibility to carry out such checks as you feel appropriate. If this email contains a quote or offer to sell products, carry out work or perform services then our standard terms and conditions (which can be found at http://www.lutraconsulting.co.uk/downloads/Lutra%20Consulting%20Standard%20Terms%20and%20Conditions.pdf shall apply unless explicitly stated otherwise. Saber Razmjooei and Peter Wells trading as Lutra Consulting.___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problem with Python GUI Grass 6
Diane, Have you tried to use GRASS plugin in QGIS? Cheers, Saber Saber Razmjooei Lutra Consulting www.lutraconsulting.co.uk http://www.lutraconsulting.co.uk/ Tel: +44 (0)1444 848012 From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Diane Palmer Sent: 30 January 2014 10:13 To: grass-user@lists.osgeo.org Subject: [GRASS-user] Problem with Python GUI Grass 6 Dear List Members I am very new to Grass and am having a few problems getting started. I’m running it on Windows 7 Enterprise, 64 bit. 1) I’ve tried the version 6.4.3 installed with QGIS. a) In the Python GUI, when I try to add a vector or raster layer (Spearfish data, Quickstart tutorial), I don’t get a drop-down list. It adds an empty file with the message “Unable to fetch interface description for command ‘d.vect’/'d.rast'”. b) When I try to import a shapefile of my own, it does import it, but I can only see the attribute table, not the map. (I did remember to zoom to it and check the colour visibility).The message is “Command ‘d.vect map=UK_outline@eldp’ failed”. 2) I installed http://grass.osgeo.org/grass64/binary/mswindows/native/WinGRASS-6.4.3-1-Setup.exe WinGRASS-6.4.3-1-Setup.exe and accepted the Microsoft Runtime libraries option. The Python GUI flickers in the toolbar and then crashes. Fortunately the old TclTK GUI works. So I have got a version which works but I’d really appreciate some suggestions on how I can get the attractive Python GUI up and running. Many thanks. Diane Palmer Research Associate CREST (Centre for Renewable Energy Systems Technology) School of Electronic, Electrical and Systems Engineering Loughborough University e-mail: d.pal...@lboro.ac.uk T: +44 (0) 1509 635 604 -- This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error please notify the system manager. This message contains confidential information and is intended only for the individual named. If you are not the named addressee you should not disseminate, distribute or copy this e-mail. Please notify the sender immediately by e-mail if you have received this e-mail by mistake and delete this e-mail from your system. If you are not the intended recipient you are notified that disclosing, copying, distributing or taking any action in reliance on the contents of this information is strictly prohibited. Whilst reasonable care has been taken to avoid virus transmission, no responsibility for viruses is taken and it is your responsibility to carry out such checks as you feel appropriate. Saber Razmjooei and Peter Wells trading as Lut ra Consulting. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] unable to start grass 7
Are you using an Anti-Virus software at all? Cheers Saber From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Nataniel Therrien Sent: 04 June 2013 20:45 To: Markus Neteler Cc: GRASS user list Subject: Re: [GRASS-user] unable to start grass 7 Hi, - winGRASS standalone. - Windows 7 - I can't see the error messages. When I click the software does not start, no error messages. Thank you for your help. On Tue, Jun 4, 2013 at 3:39 PM, Markus Neteler nete...@osgeo.org wrote: (please keep the list in cc) On Tue, Jun 4, 2013 at 9:38 PM, Nataniel Therrien nataniel.therr...@gmail.com wrote: Hi, - winGRASS standalone. - Windows 7 - I can't see the error messages. When I click the software does not start, no error messages. So the first suggestion of this page may apply: http://grasswiki.osgeo.org/wiki/WinGRASS_errors Markus -- This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error please notify the system manager. This message contains confidential information and is intended only for the individual named. If you are not the named addressee you should not disseminate, distribute or copy this e-mail. Please notify the sender immediately by e-mail if you have received this e-mail by mistake and delete this e-mail from your system. If you are not the intended recipient you are notified that disclosing, copying, distributing or taking any action in reliance on the contents of this information is strictly prohibited. Whilst reasonable care has been taken to avoid virus transmission, no responsibility for viruses is taken and it is your responsibility to carry out such checks as you feel appropriate. Saber Razmjooei and Peter Wells trading as Lutra Consulting.___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] unable to start grass 7
Check your AV log to see if it has flagged/blocked/deleted any of your installation files (usually exe files). Cheers Saber -Original Message- From: Nataniel Therrien nataniel.therr...@gmail.com To: Saber Razmjooei saber.razmjo...@lutraconsulting.co.uk Cc: 'Markus Neteler' nete...@osgeo.org, 'GRASS user list' grass-user@lists.osgeo.org Subject: Re: [GRASS-user] unable to start grass 7 Date: Tue, 4 Jun 2013 16:45:51 -0400 Yes, there is an anti-virus installed on the computer Le 2013-06-04 à 15:46, Saber Razmjooei a écrit : Are you using an Anti-Virus software at all? Cheers Saber From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Nataniel Therrien Sent: 04 June 2013 20:45 To: Markus Neteler Cc: GRASS user list Subject: Re: [GRASS-user] unable to start grass 7 Hi, - winGRASS standalone. - Windows 7 - I can't see the error messages. When I click the software does not start, no error messages. Thank you for your help. On Tue, Jun 4, 2013 at 3:39 PM, Markus Neteler nete...@osgeo.org wrote: (please keep the list in cc) On Tue, Jun 4, 2013 at 9:38 PM, Nataniel Therrien nataniel.therr...@gmail.com wrote: Hi, - winGRASS standalone. - Windows 7 - I can't see the error messages. When I click the software does not start, no error messages. So the first suggestion of this page may apply: http://grasswiki.osgeo.org/wiki/WinGRASS_errors Markus __ This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error please notify the system manager. This message contains confidential information and is intended only for the individual named. If you are not the named addressee you should not disseminate, distribute or copy this e-mail. Please notify the sender immediately by e-mail if you have received this e-mail by mistake and delete this e-mail from your system. If you are not the intended recipient you are notified that disclosing, copying, distributing or taking any action in reliance on the contents of this information is strictly prohibited. Whilst reasonable care has been taken to avoid virus transmission, no responsibility for viruses is taken and it is your responsibility to carry out such checks as you feel appropriate. Saber Razmjooei and Peter Wells trading as Lutra Consulting. --- Nataniel Therrien nataniel.therr...@gmail.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Using PVGIS Raster Files in GRASS
You can use r.colors or if you are using graphical interface: Raster Manage colors Color tables (r.colors) to change the color. Cheers Saber -Original Message- From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of James Allan Sent: 03 June 2013 16:21 To: grass-user@lists.osgeo.org Subject: [GRASS-user] Using PVGIS Raster Files in GRASS Hi All, I am very new to GRASS and I need to use the r.sun function to determine the incident solar radiation on a surface. I have downloaded the PVGIS dataset from http://re.jrc.ec.europa.eu/pvgis/download/solar_radiation_cmsaf_download.htm l and have tried importing it into the layer manager as a grid ascii file. I was successfully able to import the file but now all I see is a purple square in the map display. Please see attached. I was expecting to see a map of Europe and Africa with the corresponding radiation data. If you have any experience with solar radiation files or importing these types of files into GRASS I would really appreciate your help. Best Regards James Allan -- This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error please notify the system manager. This message contains confidential information and is intended only for the individual named. If you are not the named addressee you should not disseminate, distribute or copy this e-mail. Please notify the sender immediately by e-mail if you have received this e-mail by mistake and delete this e-mail from your system. If you are not the intended recipient you are notified that disclosing, copying, distributing or taking any action in reliance on the contents of this information is strictly prohibited. Whilst reasonable care has been taken to avoid virus transmission, no responsibility for viruses is taken and it is your responsibility to carry out such checks as you feel appropriate. Saber Razmjooei and Peter Wells trading as Lutra Consulting. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GRASS 6.4 How to extract point-values from multiple rasters
Alternatively you can do it using this script (python-gdal): http://gis-lab.info/qa/extract-values-rasters-eng.html Saber Razmjooei Lutra Consulting http://www.lutraconsulting.co.uk/ www.lutraconsulting.co.uk -Original Message- From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Markus Metz Sent: 04 April 2013 15:51 To: BLANDENIER Lucien Cc: grass [grass-user@lists.osgeo.org] Subject: Re: [GRASS-user] GRASS 6.4 How to extract point-values from multiple rasters On Thu, Apr 4, 2013 at 4:37 PM, BLANDENIER Lucien mailto:lucien.blanden...@unine.ch lucien.blanden...@unine.ch wrote: Dear GRASS community, I have to extract isolated value from several raster. I've found this script but I do not understand how to run it. Would someone have an other script with more comments (notably what are the input I have to give, what are the grass or python function/module) or explain me how to run this script. The GRASS module r.what [0] extracts raster values at a given location. As input you need to give the name of the raster map and the coordinates as east,north [0] http://grass.osgeo.org/grass64/manuals/r.what.html http://grass.osgeo.org/grass64/manuals/r.what.html Markus M I would also be interested by a bash script if someone has such a script. Thank you for your support. Kind regards Lucien ___ grass-user mailing list mailto:grass-user@lists.osgeo.org grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list mailto:grass-user@lists.osgeo.org grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user http://lists.osgeo.org/mailman/listinfo/grass-user -- This email and any files transmitted with it are confidential and intended solely for the use of the individual or entity to whom they are addressed. If you have received this email in error please notify the system manager. This message contains confidential information and is intended only for the individual named. If you are not the named addressee you should not disseminate, distribute or copy this e-mail. Please notify the sender immediately by e-mail if you have received this e-mail by mistake and delete this e-mail from your system. If you are not the intended recipient you are notified that disclosing, copying, distributing or taking any action in reliance on the contents of this information is strictly prohibited. Whilst reasonable care has been taken to avoid virus transmission, no responsibility for viruses is taken and it is your responsibility to carry out such checks as you feel appropriate. Saber Razmjooei and Peter Wells trading as Lutra Consulting.___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] NVIZ Wire Mesh
Have you looked into using Paraview to do that? Cheers Saber -Original Message- From: Anna Kratochvílová kratocha...@gmail.com To: Paul Shapley p.shap...@gmail.com Cc: grass-user@lists.osgeo.org Subject: Re: [GRASS-user] NVIZ Wire Mesh Date: Thu, 26 Jan 2012 15:29:51 +0100 On Thu, Jan 26, 2012 at 2:09 PM, Paul Shapley p.shap...@gmail.com wrote: Hello, I want to have the NVIZ wire mesh opaque but don't know how best to do this. I don't want to be transparent. Does anyone know how this can be done? Sorry but this is not possible. You can create a ticket to support this behaviour but I'm afraid that nobody will work on it now. Anna Thanks, -- Paul J. Shapley ___ 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] GRASS GIS development visualization from 1999 to 2011
Excellent work...both the video and development! Cheers Saber Hi, I have created a video of the GRASS GIS development from 1999 to 2011: http://grass.osgeo.org/wiki/GRASS_development_visualization Enjoy Markus PS: It is on youtube, http://www.youtube.com/watch?v=suyDqmGXoWk The high-res version is also available in some minutes, see Wiki page. ___ 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] Draping of Raster images
This page might be helpful: http://grass.osgeo.org/wiki/Help_with_3D Cheers Saber Dear all, I was wondering whether it is possible to use GRASS to drape raster images over existing 3D surfaces or vector layers. I know ArcGIS is capable of this, is seems a fairly straight-forward procedure there (http://webhelp.esri.com/arcgisdesktop/9.3/tutorials/3D_analyst/3D_5.htm), so I was hoping to achieve the same using open source GIS. GRASS has a v.drape module, which can drape vector layers over 3D surfaces, NVIZ seems to do the same. It would be great if someone could tell me that GRASS (or any other open source GIS package) is capable of draping raster images over either other rasters or vectors containing elevation information. Many thanks and all the best, Anna __ Anna Kathrin Hodgkinson MA MSt (Oxon) AIFA EES Supervisor Geomatics Oxford Archaeology: Exploring the Human Journey Mobile: +44 (0)7906 638308 (personal) Direct Dial: +44 (0)1865 980 855 http://thehumanjourney.net Files attached to this email may be in ISO 26300 format (OASIS Open Document Format). If you have difficulty opening them, please visit http://iso26300.info for more information. This email has been processed by SmoothZap - www.smoothwall.net ___ 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] Draping of Raster images
If you are using spearfish dataset demo (http://grass.osgeo.org/download/data6.php) g.region rast=elevation.10m nviz elev=elevation.10m Cheers Saber Thanks, but there doesn't seem to be a way of elevating a 2D raster surface. Do you have any experience with this, and would know what command to use? Many thanks, Anna - Original Message - From: Saber Razmjooei razmjoo...@faunalia.co.uk To: Anna Hodgkinson anna.hodgkin...@thehumanjourney.net Cc: grass-user@lists.osgeo.org Sent: Thursday, 3 November, 2011 1:52:36 PM Subject: Re: [GRASS-user] Draping of Raster images This page might be helpful: http://grass.osgeo.org/wiki/Help_with_3D Cheers Saber Dear all, I was wondering whether it is possible to use GRASS to drape raster images over existing 3D surfaces or vector layers. I know ArcGIS is capable of this, is seems a fairly straight-forward procedure there (http://webhelp.esri.com/arcgisdesktop/9.3/tutorials/3D_analyst/3D_5.htm), so I was hoping to achieve the same using open source GIS. GRASS has a v.drape module, which can drape vector layers over 3D surfaces, NVIZ seems to do the same. It would be great if someone could tell me that GRASS (or any other open source GIS package) is capable of draping raster images over either other rasters or vectors containing elevation information. Many thanks and all the best, Anna __ Anna Kathrin Hodgkinson MA MSt (Oxon) AIFA EES Supervisor Geomatics Oxford Archaeology: Exploring the Human Journey Mobile: +44 (0)7906 638308 (personal) Direct Dial: +44 (0)1865 980 855 http://thehumanjourney.net Files attached to this email may be in ISO 26300 format (OASIS Open Document Format). If you have difficulty opening them, please visit http://iso26300.info for more information. This email has been processed by SmoothZap - www.smoothwall.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user Files attached to this email may be in ISO 26300 format (OASIS Open Document Format). If you have difficulty opening them, please visit http://iso26300.info for more information. This email has been processed by SmoothZap - www.smoothwall.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Draping of Raster images
Alternatively, you can use: r.shaded.relied That will be 2D only. Thanks, but there doesn't seem to be a way of elevating a 2D raster surface. Do you have any experience with this, and would know what command to use? Many thanks, Anna - Original Message - From: Saber Razmjooei razmjoo...@faunalia.co.uk To: Anna Hodgkinson anna.hodgkin...@thehumanjourney.net Cc: grass-user@lists.osgeo.org Sent: Thursday, 3 November, 2011 1:52:36 PM Subject: Re: [GRASS-user] Draping of Raster images This page might be helpful: http://grass.osgeo.org/wiki/Help_with_3D Cheers Saber Dear all, I was wondering whether it is possible to use GRASS to drape raster images over existing 3D surfaces or vector layers. I know ArcGIS is capable of this, is seems a fairly straight-forward procedure there (http://webhelp.esri.com/arcgisdesktop/9.3/tutorials/3D_analyst/3D_5.htm), so I was hoping to achieve the same using open source GIS. GRASS has a v.drape module, which can drape vector layers over 3D surfaces, NVIZ seems to do the same. It would be great if someone could tell me that GRASS (or any other open source GIS package) is capable of draping raster images over either other rasters or vectors containing elevation information. Many thanks and all the best, Anna __ Anna Kathrin Hodgkinson MA MSt (Oxon) AIFA EES Supervisor Geomatics Oxford Archaeology: Exploring the Human Journey Mobile: +44 (0)7906 638308 (personal) Direct Dial: +44 (0)1865 980 855 http://thehumanjourney.net Files attached to this email may be in ISO 26300 format (OASIS Open Document Format). If you have difficulty opening them, please visit http://iso26300.info for more information. This email has been processed by SmoothZap - www.smoothwall.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user Files attached to this email may be in ISO 26300 format (OASIS Open Document Format). If you have difficulty opening them, please visit http://iso26300.info for more information. This email has been processed by SmoothZap - www.smoothwall.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] WPS with GRASS and PyWPS
Dear list We have developed a Web Processing service for calculating catchment outlines in UK using GRASS and PyWPS: http://www.faunalia.co.uk/catchmentfinder Any feedback is appreciated. Regards Saber ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] WPS with GRASS and PyWPS
Valentin Many thanks for your feed back. The DTM is quite old (over 20 years) and very coarse (50 metre..but grass region setting is I think 100 metre). The reason you don't see the point is that it doesn't use the point and it uses another one snapped to the stream (r.watershed output). But we definitely add something to show the point, either user input or the snapped one. In terms of validation, we have Flood Estimation Handbook which unfortunately is not free and using the same DTM in our WPS. For copyright reasons there is no vector output from FEH either. So you need to create screenshots. I don't have a license for it. There will be more bells and whistles to it (hydrological parameters: urban extent, main drain slope, max min elevation, etc). Cheers Saber Nice tool, I like it! Three things came up in my mind when I used it. 1) After I selected a point and press calculate I couldn't see any more the point I selected... 2) It seams like there are some wrong calculation. For example for the randomly chosen point (527714, 328638). And the catchment border could be a little smoother ;) 3) I can only select a visual point, it would be nice to select the point by coordinates... it would be more useful. Did you try to validate your data with a independent dataset? I know from France there exists free hydrological surface catchments shapes over the whole country maybe it possible to make a validation with public DEM-data like SRTM. regards Valentin On Oct 26, 2011 9:27 AM, Saber Razmjooei razmjoo...@faunalia.co.uk mailto:razmjoo...@faunalia.co.uk wrote: Dear list We have developed a Web Processing service for calculating catchment outlines in UK using GRASS and PyWPS: http://www.faunalia.co.uk/catchmentfinder Any feedback is appreciated. Regards Saber ___ grass-user mailing list grass-user@lists.osgeo.org mailto: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] WPS with GRASS and PyWPS
Dear Valentin Thanks for your through inspection. That part of UK is flat with lots of artificial drains and canals. There is a loop river and also 2 parallel channel. The DTM doesn't look great either and I think streams have not been extracted properly in r.watershed. Cheers Saber On 26/10/11 21:07, Valentin Wittich wrote: Hi Saber, could you recreate the wrong calculation I've get with the point (527714, 328638). The point is a little above the connection to the tributary, could it be that the snap function doesn't work proper? V. On 10/26/2011 12:08 PM, Saber Razmjooei wrote: Valentin Many thanks for your feed back. The DTM is quite old (over 20 years) and very coarse (50 metre..but grass region setting is I think 100 metre). The reason you don't see the point is that it doesn't use the point and it uses another one snapped to the stream (r.watershed output). But we definitely add something to show the point, either user input or the snapped one. In terms of validation, we have Flood Estimation Handbook which unfortunately is not free and using the same DTM in our WPS. For copyright reasons there is no vector output from FEH either. So you need to create screenshots. I don't have a license for it. There will be more bells and whistles to it (hydrological parameters: urban extent, main drain slope, max min elevation, etc). Cheers Saber Nice tool, I like it! Three things came up in my mind when I used it. 1) After I selected a point and press calculate I couldn't see any more the point I selected... 2) It seams like there are some wrong calculation. For example for the randomly chosen point (527714, 328638). And the catchment border could be a little smoother ;) 3) I can only select a visual point, it would be nice to select the point by coordinates... it would be more useful. Did you try to validate your data with a independent dataset? I know from France there exists free hydrological surface catchments shapes over the whole country maybe it possible to make a validation with public DEM-data like SRTM. regards Valentin On Oct 26, 2011 9:27 AM, Saber Razmjooeirazmjoo...@faunalia.co.uk mailto:razmjoo...@faunalia.co.uk wrote: Dear list We have developed a Web Processing service for calculating catchment outlines in UK using GRASS and PyWPS: http://www.faunalia.co.uk/catchmentfinder Any feedback is appreciated. Regards Saber ___ grass-user mailing list grass-user@lists.osgeo.orgmailto: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] PyWPS and GRASS - Example?
http://pywps.wald.intevation.org/documentation/process.html We are working on some demos. Should be up and running soon(ish). Cheers Saber Hi everyone, i'm looking for an examplary PyWPS process that makes use of GRASS GIS. So far I found a lot of tutorials concerning the setup of PyWPS but what I need is a working process - if possible with a fitting OpenLayers (or similar) frontend. Can anybody help me out? Thanks in advance, Georg ___ 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] Location of minimum value in raster
First make a 2nd copy of your raster r.info to get the value of minimum. r.null to null out all the values more than minValue + 0.01 (depends on the decimal values in your raster) through to maxValue. Cheers Saber Hi, 2011/8/3 Hamish hamis...@yahoo.com: Joshua wrote: [snip] two steps, eval `r.info -r data_map` to get min= r.mapcalc binary_map = if(data_map == $min, 1, 0) or perhaps min from 'r.info -r' and coord(s) from grep'ing output of 'r.stats -1g'. perhaps you need to do abs(data_map - $min) 0.001 for the equivalence test in r.mapcalc, or some similar awk test from the printf() formatted output of r.stats. perhaps you need to do r.buffer from the binary map to make the single pixel(s) easier to locate. d.rast with the cat or value list option might quickly do it too, but a single cell might not visibly render. Maybe you can convert the resulting raster map into vector points for visualization and further processing? Soeren ___ 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] Mapcalc to replace a range of values
Or maybe you can use r.null: r.null map=YourRASTER setnull=-0.1-0.1 null=0 On Wed, 2011-07-13 at 14:10 +0100, Rebecca Bennett wrote: Hello Grass users, I have a (hopefully) quick question - I would like to replace all values in a raster that fall in the range -0.01 to 0.01 (i.e. 0.01-0.01) with 0 but can't quite work out how to express this in mapcalc. Can anyone set me on the right track? Thanks for reading, Rebecca ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to draw the earthquake-magtintude map...
It's a bit off topic and an item in my wishlist, but will be interesting to have AnuGA integrated in GRASS. http://sourceforge.net/projects/anuga/ http://en.wikipedia.org/wiki/ANUGA_Hydro and a tsunami risk assessment in New Zealand: http://www.aon.com/attachments/reinsurance/200812_ab_analytics_nz_tsunami.pdf We did a plugin a while back in QGIS for preprocessing the gis layers for AnuGA. Cheers Saber On Wed, 2011-06-15 at 12:51 +0300, Micha Silver wrote: On 15/06/2011 10:46, sgw00...@nifty.com wrote: Dear grass users I'm from Japan and begginer user of grass GIS. I tried to draw earthquake map of tohoku Japan. First, I summarized of some earthquake data, for example lat, lon, depth and m agnitude. - lon latdepth Magnitude - 139.9, 37.7, 10, 3.1 141.0, 35.8, 10, 2.6 142.7, 40.8, 20, 4.4 ... - Second, this text data was imported to grass data by using v.in.ascii, as foll ows; cat data.txt | v.in.ascii fs=',' cat=0 x=1 y=2 out=quake columns='x double pre cision, y double precision, depth integer, mag double precision' So, I typed to 'v.vect quake' These points where the earthquake were displayed by cross (x) . Next, I tried to display new map as circle size of earthquake-magnitude. You probably want here: d.vect.thematic map=quake col=mag type=point themetype=graduated_points icon=basic/circle There are more options for coloring, etc. v.vect quake size=magnitude However, unsuccessfulWhat should I do? Yasuo shimada ( Japan weathere association ) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Hevra Logo Micha Silver Arava Development Co 052-3665918 http://www.surfaces.co.il ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: how to only fill small sinks with r.fill.dir
You probably need a full 2D hydraulic model to do work based on water level gradient. Maybe something like like ANUGA: http://sourceforge.net/projects/anuga/ On Tue, 2011-04-26 at 11:54 -0600, Bryan Keith wrote: Great ideas, thank you. I did something like this: r.fill.dir input=wl200310_krig_g elevation=wl200310_fill_g direction=wl200310_dir_g r.mapcalc wl200310diff_g=wl200310_fill_g - wl200310_krig_g r.mapcalc wl200310diffsel_g = if(wl200310diff_g 5.0, 0, wl200310diff_g) r.mapcalc wl200310sel_g = wl200310diffsel_g + wl200310_krig_g r.fill.dir input=wl200310sel_g elevation=wl200310sel_fill_g direction=wl200310sel_dir_g r.flow elevin=wl200310sel_fill_g flout=test_ln I think my problem is now with r.flow. I would like to see flows continue through the filled sink and out the other end (as if it were a lake), but that doesn't happen. I don't think it can happen with r.flow because r.flow encounters a flat area. Not sure what to try next. This was a problem even before messing around with the sink tolerance. Bryan On Tue, Apr 26, 2011 at 10:31, Saber razmjoo...@faunalia.co.uk wrote: Sorry...I thought the sinks are as null values. Ignore my previous email. I can't think of anything apart from improving your DEM to get rid of the artifacts. One way will be to get the initial result from the r.fill.dir. Identify the areas where the depth is more than say 10 metre (r.nulls map=initial_result setnull 0-9.99) You can then add the new depth raster to the DEM to create the filled hole. (first r.mask to match initial_result and then r.mapcalc filled_holes=initial_result+DEM) Then remove the mast and r.patch r.patch input=filled_holes,DEM output=patched_DEM You can then use the patched_DEM for r.fill.dir Hope that helps. Hmmm, I'm not following all your steps, but I don't want to manually identify valid sinks vs. invalid sinks. I have quite a lot of surfaces to process, and I'm happy to fill all sinks except those deeper than 10m. On Tue, Apr 26, 2011 at 10:04, Saber razmjoo...@faunalia.co.uk wrote: Here is a (not a very quick) solution I can think of: 1- Create an inverse mask 2- convert the mask to vector (r.to.vect) 3- clean the vector file to get rid of the areas not required (v.clean) 4- convert the vector back to raster (v.to.rast) 5- Use the new raster as a mask for filling holes (probably you need step 4.5 to slightly buffer the raster) Hope that helps Any ideas on this? Bryan On Mon, Apr 18, 2011 at 14:49, Bryan Keith bke...@itascadenver.com wrote: Hello, Is there a way to only fill sinks in a DEM if they are smaller than a certain tolerance? I have a couple areas of internal drainage that really are correct (and quite large), and I have a quite few small sinks that are incorrect and artifacts of the interpolation. Any ideas how to this? r.fill.dir attempt to fill all the sinks. I'm trying to improve the results from r.flow Too many tracks end in small holes where I'd like them to continue downstream. Thanks. Bryan ___ 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.where in GRASS 7.0
Hi all I am trying to pipe d.where output to a new vector using example 6 of the manual: http://grass.fbk.eu/grass70/manuals/html70_user/v.in.ascii.html In grass 6.x, it was possible to set the active display by d.mon But in grass 7.0 this has been removed. Therefore, d.where sits there without outputing x and y. Is there a way to set the active display in grass 7? I tried to follow this...but no luck: http://grass.fbk.eu/grass70/manuals/html70_user/cairodriver.html Thanks in advance. Cheers Saber ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Vector map text legend
How about using Grass Plugin in QGIS to load your vectors and create legend there? Not sure how it ties up with your python script. Cheers Saber On Thu, 2011-03-31 at 11:39 +0200, Martin Album Ytre-Eide wrote: Hello! I am looking for a way to make a legend to some vector area maps. The maps contain four different types of areas and are given different names. I have made a python script running d.mon, d.vect etc. with the png driver producing about 14 png files, which then is used to create a gif file. A way to create a vector legend is to do v.to.rast and use the d.legend from the dummy rast - however, I don't see how you can get a legend with text in this way? I don't think you can since raster should have numerical values. You could make a legend without category numbers and try to add the text with d.text, but I don't know if that's a good way of doing things. Then there is ps.map which solves this in a nice way - the problem is - I would like to keep the png format. I could just make the legend with ps.map - and add it to my png files with some image tools in python, but I am unsure of the result is worth the effort. Any other ideas to do this? Thanks Martin Ytre-Eide ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] New User
If you are intending to use GIS for visualization and data conversion, maybe better to use something easy to start off, like QGIS: http://www.qgis.org/wiki/Download#Release_2 You have GRASS plugin there if needed. But usually you should be able to do export the files to other OGR/GDAL supported formats from within QGIS. Regards Sab On Sun, 2011-02-27 at 23:25 -0800, Hamish wrote: Thomas wrote: I am a new user, just discovered grass and OSGeo. My current interest is accessing the files from NaturalEarthData, converting them into matlab format, and, using octave to manipulate and/or plot. I use Linux and Unix. Octave 3.4.0, Gnuplot, Grass 6.4.0. This workstation has Ubuntu 10.04. I downloaded the 10m-admin-0-countries.zip and 10m-admin-0-boundary-lines-land.zip files from NaturalEarthData. I want to use grass to display some portion of the data. Then, use r.out.mat to generate a octave compatible binary file. I started looking at the manual, but, quickly became confused. that would be my fault. Is there a simple how to for what I want to do? Could you explain your workflow a bit more? there are many choices. In general it's not so bad, and you can do matlab/oct to and from grass without much hassle. natural earth is polyline data, so vector thus v.* modules. 2D array is raster data, so r.* modules. see v.to.rast and r.to.vect modules. do you have 2D arrays in matlab/octave from contourf() or similar and you want to fill(), area(), or patch() a polygon coastline over the top? see http://grass.osgeo.org/wiki/Addons#v.out.ascii.mat are you using the (free/gratis) m_map toolbox? http://www.eos.ubc.ca/~rich/map.html or matlab's mapping toolbox? (which is using proj4, same as us) do you wish to rasterize the natural earth polygons into a 2D array? see v.to.rast + r.out.mat do you want to bring your matlab data into grass and plot it there? see r.in.mat, or v.in.lines, v.in.mapgen. the first step is to bring the data into grass with v.in.ogr. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: How to merge rasters using feathering/blending ?
If you are working with Images, maybe try to use GIMP. It has stitching and smoothing/blending options. I use them quite often for my photos. Regards Sab On Mon, 2011-02-28 at 04:05 -0800, kaipi wrote: You can create a buffer zone between your 2 dataset where they merge and set the buffer to null values. Then using r.fillnulls Thanks. This works fine for elevation data but not for satellite imagery. g.region rast=map1,map2 r.series in=map1,map2 out=both_maps method=average Thanks. But this is not really a smooth blending/fading/feathering ? Any alternatives ? Kaipi ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] calculate raster statistics for vector polygons
http://grass.fbk.eu/grass64/manuals/html64_user/v.rast.stats If you have lots of polygons, you probably need to look into v.rast.stats2 , an addon Cheers Sab On Mon, 2011-02-28 at 21:32 +, Kim Besson wrote: Greetings I have a vector file with polygons and I want to obtain a raster map with statistics of an input raster for each polygon. What can I use to do this? Thanks Kim ___ 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 .shp
You can just type in terminal where grass is running: v.in.ogr it will bring up the right window. Sab On Thu, 2011-02-24 at 15:31 -0800, Rich Shepard wrote: I've lost the GUI menu item for importing .shp files in -6.5svn. If it's on the File-Import Vector menu please point out which item I should use. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to merge rasters using feathering/blending ?
I guess you can do it with a combination of the existing tools: You can create a buffer zone between your 2 dataset where they merge and set the buffer to null values. Then using r.fillnulls (http://grass.osgeo.org/gdp/html_grass64/r.fillnulls.html) and you can fill the gap using both datasets. Hope that helps. Cheers Sab On Thu, 2011-02-24 at 10:53 -0800, kaipi wrote: I need to patch rasters containing NULL areas with alternative rasters. I can do this with r.patch but due to different quality/resolution of the rasters this produces prominent seamlines. Is there a way to feather/blend between merged rasters using Grass GIS? Thanks, Kaipi ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Multithreading in GRASS 6.4
I haven't done it for 6.4 but tried it in 7.0 You can follow the instruction on compiling grass here: http://grass.osgeo.org/wiki/Compile_and_Install Once you have all the prerequisite, you can set the environment variables CFLAGS=-O3 -Wall -Werror-implicit-function-declaration -fno-common -fopenmp LDFLAGS=-lgomp and then compile grass. Cheers Sab On Wed, 2010-12-29 at 16:03 -0800, LeeDaniel wrote: Hello, I'm trying to use multithreading in GRASS 6.4 in OpenSUSE 64-bit installed from the community repository and am wondering what the best way would be to set my script. What I'm wanting to do is several r.sun analyses parallel to one another using a Python script. I've read the wiki page about OpenMP at http://grass.osgeo.org/wiki/OpenMP#Multithreaded_jobs_in_GRASS but don't quite understand it - it seems that I can activate multithreading capabilities in GRASS by modifying /lib/gpde/Makefile but I don't have that file. I also wasn't able to find a whole lot of information about it that I could understand - I meet the requirement listed at http://grass.osgeo.org/programming6/gpdelib.html of having a compiler that supports OpenMP (I've got gcc 4.5), but am unsure of how I should proceed. Does anyone have any tips? Happy New Year's and thanks a bunch in advance! Daniel ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] batch conversion of ESRI to MapInfo
Mik As suggested, better to use ogr2ogr. if you use linux: 1- Create an empty text file in the folder where your shapefiles are located 2- Open the text file and copy paste these: mkdir ./mapinfo for file in *.shp do ogr2ogr -f MapInfo File mapinfo/$file.tab $file done 3-save the file as convert.sh 4- from the terminal, run convert.sh $ cd /PATH/TO/SHAPEFILES/ $ sh convert.sh 5- all your mapinfo files should appear under mapinfo folder. Cheers Sab On Sat, 2011-01-29 at 09:15 +, Leonidas Liakos wrote: What about using ogr2ogr? ogr2ogr -f MapInfo File mapinfo.tab input.shp You have to write a script for batch conversion. --- Στις Σάβ., 29/01/11, ο/η Mick bare...@tpg.com.au έγραψε: Από: Mick bare...@tpg.com.au Θέμα: [GRASS-user] batch conversion of ESRI to MapInfo Προς: grass user grass-user@lists.osgeo.org Ημερομηνία: Σάββατο, 29 Ιανουάριος 2011, 9:43 I have downloaded several sheets (approx 1600 files) from UK Ordnance Survey open data web site in ESRI shapefile format and OS grid references and now I need to convert them to MapInfo format and real world co-ordinates. I currently have available GRASS and Quantum GIS in linux and Mapinfo 4.1 5 running in windows emulation. Don't laugh, poverty dictates I can't afford to update. Does anyone know of a better solution than loading each shapefile into QGIS then save as Mapinfo and keeping the OS grid ref.? mik ___ 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] Re: batch conversion of ESRI to MapInfo
On Sat, 2011-01-29 at 12:48 +0100, Hermann Peifer wrote: ogr2ogr -f MapInfo File mapinfo/$file.tab $file You could cut away the shp suffix like this: ogr2ogr -f MapInfo File mapinfo/${file%shp}tab $file indeed! cheers Sab I am not fully sure what real world co-ordinates means, but this additional switch might do what you expect: -t_srs epsg:4326 Hermann On 29/01/2011 10:53, Saber Razmjooei wrote: Mik As suggested, better to use ogr2ogr. if you use linux: 1- Create an empty text file in the folder where your shapefiles are located 2- Open the text file and copy paste these: mkdir ./mapinfo for file in *.shp do ogr2ogr -f MapInfo File mapinfo/$file.tab $file done 3-save the file as convert.sh 4- from the terminal, run convert.sh $ cd /PATH/TO/SHAPEFILES/ $ sh convert.sh 5- all your mapinfo files should appear under mapinfo folder. Cheers Sab ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.patch: number of rasters
I am having about 700 LiDAR tiles (~4 GB). Did it a couple of hours ago and all worked well...the DTM looks as I expected! The problem I had last year when I did it last, there were 1000+, values for one tile were changed. Never investigated the problem. Simply changed the region and things worked smoothly. Many thanks Sab On Tue, 2010-11-02 at 16:54 +0100, Markus Neteler wrote: On Tue, Nov 2, 2010 at 4:41 PM, razmjoo...@faunalia.co.uk wrote: Hi Markus Many thanks for your response. In the past versions, the patched raster used to have an odd values. No, this is rather impossible (perhaps you refer to the odd number of pixels in a moving raster window). I use Grass 6.4 from Ubuntu 10.10 64 bit repo. Master question: how many files do you plan to patch? Markus Cheers Sab On Tue, Nov 2, 2010 at 2:03 PM, razmjoo...@faunalia.co.uk wrote: Dear all Is there a limit on the number of rasters for r.patch? How funny, we had the same question at the same time... :) From GRASS 6.4+ there is no more limit despite the limit of open files per process (which depends on the operating system; for Linux it is usually 1024 files). To increase that: The limit of open files per process was reached (64bit Linux: 1024 files; ulimit -n) To overcome this limit, the nofile parameter needs to be increased in the /etc/security/limits.conf file. Here an increase to 1500 files (use * instead of neteler to aply it to all users): #domain type item value neteler hardnofile 1500 But myself, I got troubles when having registered the maps with r.external: GRASS 6.4.1svn (patUTM32):~ r.patch in=`g.mlist type=rast pattern=dtm00* sep=comma mapset=dtm_1m_external` out=dtm_trentino_2m ERROR 1: TIFFOpen:/geodata/base_cartography/pat_DTM_LIDAR_2007_UTM_WGS84/geotiffs/dtm000501.tif: Too many open files WARNING: Unable to open '/grassdata/patUTM32/dtm_1m_external/cell_misc/dtm000501/f_format' WARNING: quantization file [dtm000501] in mapset [dtm_1m_external] missing WARNING: Unable to open raster map dtm000...@dtm_1m_external ERROR: One or more input raster maps not found So far I could not find out the reason. Here I am using 500 files registered with r.external. Help welcome, Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GRASS GIS
On Sat, 2010-10-09 at 17:28 +1100, azadeh farjadnejad wrote: Hi, I am a software developer and I am working on a forestry application involving remote sensing, 3D spatial modeling and image analysis for processing full-waveform airborne laser scanner (or lidar) data. I am looking for a high performance open source library for that. I want to know 1- Can I use GRASS GIS libraries to develop my own application? You can, but like all other software you need to comply with its license. GRASS is released under GPL (http://grass.fbk.eu/intro/general.php) here is more on GPL: http://en.wikipedia.org/wiki/GNU_General_Public_License 2- Do I have any limitation to choose any programming languages if I only want to use GIS libraries to develop my application? The core of GRASS is in C, C++. For scripting you can choose Python. 3- What if I use Windows platform. Do I still have the capability of using libraries with the same performance and speed? I guess yes, but I haven't used it under windows 4- Is there any other open source libraries you may recommend specifically for remote sensing, 3D spatial modeling and lidar data? http://en.wikipedia.org/wiki/ILWIS Thanks. I look forward hearing from 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] GRASS GIS
Hope this helps On Sat, 2010-10-09 at 21:24 +1100, azadeh farjadnejad wrote: Hi, Thanks for the proper answer. one more question: Assume my application is written in Java. Is it a way to call functions of GRASS libraries in my Java application to use its functionality? http://jgrass.wiki.software.bz.it/jgrass/JGrass_Wiki or I have to use C, C++ and Python to develop the whole application? On the whole is it possible to use other languages for development using GRASS libraries rather than these three languages? Regards On Sat, Oct 9, 2010 at 6:59 PM, Saber Razmjooei razmjoo...@faunalia.co.uk wrote: On Sat, 2010-10-09 at 17:28 +1100, azadeh farjadnejad wrote: Hi, I am a software developer and I am working on a forestry application involving remote sensing, 3D spatial modeling and image analysis for processing full-waveform airborne laser scanner (or lidar) data. I am looking for a high performance open source library for that. I want to know 1- Can I use GRASS GIS libraries to develop my own application? You can, but like all other software you need to comply with its license. GRASS is released under GPL (http://grass.fbk.eu/intro/general.php) here is more on GPL: http://en.wikipedia.org/wiki/GNU_General_Public_License 2- Do I have any limitation to choose any programming languages if I only want to use GIS libraries to develop my application? The core of GRASS is in C, C++. For scripting you can choose Python. 3- What if I use Windows platform. Do I still have the capability of using libraries with the same performance and speed? I guess yes, but I haven't used it under windows 4- Is there any other open source libraries you may recommend specifically for remote sensing, 3D spatial modeling and lidar data? http://en.wikipedia.org/wiki/ILWIS Thanks. I look forward hearing from 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] QGIS GRASS plugin shadedrelief issue
Have you tried QGIS plugin? It should be under 3rd party plugins. Cheers Sab On Sun, 2010-09-19 at 04:03 -0700, nunosousa84 wrote: Hi, Im new to the forum and to GIS (both GRASS and QGIS). Im using QGIS GRASS Plugin. Im having some issues using GRASS shadedrelief functionality. When i manipulate the values of the altitute of the sun above the horizon (in degrees) i catch some black spots/dots in the final result. This maybe normal since we are working with different altitudes so shadow projections vary and some areas might be out of reach? Anyway if you could help me with this (try to diminish the problem or try another method) i would apreciate. Here are 3 examples: http://osgeo-org.1803224.n2.nabble.com/file/n5547443/shaded30270.jpeg http://osgeo-org.1803224.n2.nabble.com/file/n5547443/shaded50270.jpeg http://osgeo-org.1803224.n2.nabble.com/file/n5547443/shaded80270.jpeg ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: QGIS GRASS plugin shadedrelief issue
There is a native QGIS plugins: http://spatialserver.net/pyqgis_1.0/contributed/shadedrelief.zip To activate it in QGIS: Fetch Plugins, a new window will pop up: Click on Repositories tab Add 3rd Party repository Click on Plugins Install DEM relief Shader You should be able to create relief shade as TIFF format from your dem. Cheers Sab On Sun, 2010-09-19 at 04:37 -0700, nunosousa84 wrote: Sorry didnt understand. You mean the GRASS plugin? available in QGIS. Plugins - Manage Plugins - Grass. Thats what im using. Working inside QGIS but through GRASS plugin. Sorry for any misunderstanding. Nuno ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing multiple files with r.in.xyz
Hanile Try this: 1- Install Paths: http://www.textpad.com/add-ons/files/utilities/paths.zip It's a useful tool to copy the path and filenames in windows. It adds a function to rightclick: Pathcopy. You can select all your rasters and rightclick to copy the paths: Paste them in a text file c:\example\raster1.xyz c:\example\raster2.xyz ... 2- change them to something like this (use find/replace..macros, etc:) r.in.xyz -s -g input=c:\example\raster1.xyz output=raster1 tmpRegion SET /p myregion= tmpRegion DEL tmpRegion g.region %myregion% r.in.xyz --overwrite input=c:\example\raster1.xyz output=raster1 fs=, 3- Save that as a import.bat 4- From the Start menu, run Grass in Txt mode, change directory to the folder where import.bat is 5- import.bat If you google around you should be able to make the script a bit smarter using DOS and loop (http://www.robvanderwoude.com/ntfor.php). In linux for example you can shrink it to something like this: cat list_of_files.txt | while read line; do echo $line+Being processed # or whaterver you want to do with the $line variable raster=$line r.in.xyz -s -g input=$line output=$line tmpRegion myregion= `head -n 1 tmpRegion` g.region $myregion r.in.xyz --overwrite input=$line output=$line fs=, done and your list_of_files.txt is raster1 raster2 . Hope that helps Saber On Tue, 2010-08-10 at 14:59 +0200, Hanlie Pretorius wrote: Hi, I'm working with WinGRASS 6.4RC6 on Win XP. I'm following the steps at http://grass.osgeo.org/wiki/Import_XYZ to import hundreds files containing TRMM rainfall data. When I try to follow the instruction: - Just amend above procedure to use wildcards. Change in above example all occurencies of cat VTL2733.XYZ | ... to cat *.XYZ | ... and use a more reasonable output name of course. That's all to import even thousands of files (tiled DEM) easily. - Not sure what a 'reasonable output name is', so I run it as shown below and get the following error: - cat F:\Hanlie\UCT\M.Sc.\Data\TRMM\2000\02_Februrarie\*.txt| r.in.xyz -s in=- fs=, out=test cat: F:\Hanlie\UCT\M.Sc.\Data\TRMM\2000\02_Februrarie\*.txt: No such file or directory - If I try it with one file only, it works: - cat F:\Hanlie\UCT\M.Sc.\Data\TRMM\2000\02_Februrarie\3B42.000201.0.6.nc.lieb.txt| r.in.xyz -s in=- fs=, out=test Range: x:-28.625000-27.375000 y: 28.125000 28.625000 z: 0.00 0.00 - How do I get this to work? And how would I get the output filenames to be related to the input filenames? For example, if the input filename is 3B42.000201.0.6.nc.lieb.txt, then I want the output raster to be named '3B42.000201.0.6'. Thanks Hanlie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user