Re: [GRASS-user] help with testing the location wizard for the upcoming release
On 26/05/13 13:20, Hamish wrote: known issues: when reading from a .prj or WKT file, or creating from a geo- referenced file, with unspecified datum transform opts, you will not be prompted to enter them. the location will be created not mentioning them, but in the background PROJ.4 will assign some default value (of its choosing) silently. Trying to define a sector in Belgian Lambert (1972 or 2008) projections, I see the following oddity (linked to the above ?): Using EPSG codes, I do get prompted with a choice for either no transformation or one choice of transform parameters, but even when I chose Do not apply any datum transformations I get a sector with the towgs84 parameter defined for transformation. The same happens when I try to define these projections manually by choosing in a list. When I enter the proj4 string without the towgs84 parameter it works as expected, i.e. no towgs84 defined. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] help with testing the location wizard for the upcoming release
since there are so many projections and datums and ways to select them, it's pretty hard to test for working/not working as we usally could do. Please try throwing your local odd-ball CRS at it and see what happens. (check all the way through to 'g.proj -p' and 'g.proj -j' in the final location session. make sure it does what you'd expect it to do. EPSG 3416 ETRS89 / Austria Lambert System Info GRASS version: 6.5.svn GRASS SVN Revision: 56422 GIS Library Revision: 50936 (2012-02-25) GDAL/OGR: 1.9.2 PROJ4: Rel. 4.8.0, 6 March 2012 Python: 2.7.4 wxPython: 2.8.12.1 Platform: Windows-7-6.1.7601-SP1 (OSGeo4W) Select datum transformation 0 Do not apply any datum transformation g.region -p projection: 99 (Lambert Conformal Conic) zone: 0 datum: etrs89 ellipsoid: grs80 north: 1 south: 0 west: 0 east: 1 nsres: 1 ewres: 1 rows: 1 cols: 1 cells: 1 g.proj -p -PROJ_INFO- name : Lambert Conformal Conic proj : lcc datum : etrs89 ellps : grs80 lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33 x_0: 40 y_0: 40 no_defs: defined -PROJ_UNITS unit : metre units : metres meters : 1 g.proj -j +proj=lcc +lat_1=49 +lat_2=46 +lat_0=47.5 +lon_0=13.33 +x_0=40 +y_0=40 +no_defs +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000 +to_meter=1 1 Used in whole etrs89 region: towgs84=0.000,0.000,0.000 (Default 3-Parameter Transformation) g.region -p projection: 99 (Lambert Conformal Conic) zone: 0 datum: etrs89 ellipsoid: grs80 north: 1 south: 0 west: 0 east: 1 nsres: 1 ewres: 1 rows: 1 cols: 1 cells: 1 g.proj -p -PROJ_INFO- name : Lambert Conformal Conic proj : lcc datum : etrs89 ellps : grs80 lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33 x_0: 40 y_0: 40 no_defs: defined towgs84: 0.000,0.000,0.000 -PROJ_UNITS unit : metre units : metres meters : 1 g.proj -j +proj=lcc +lat_1=49 +lat_2=46 +lat_0=47.5 +lon_0=13.33 +x_0=40 +y_0=40 +no_defs +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000 +to_meter=1 GRASS version: 6.4.3svn GRASS SVN Revision: 56413 GIS Library Revision: 50937 (2012-02-25) GDAL/OGR: 1.9.2 PROJ4: Rel. 4.8.0, 6 March 2012 Python: 2.7.4 wxPython: 2.8.12.1 Platform: Windows-7-6.1.7601-SP1 (OSGeo4W) Select datum transformation 0 Do not apply any datum transformation g.region -p projection: 99 (Lambert Conformal Conic) zone: 0 datum: etrs89 ellipsoid: grs80 north: 1 south: 0 west: 0 east: 1 nsres: 1 ewres: 1 rows: 1 cols: 1 cells: 1 g.proj -p -PROJ_INFO- name : Lambert Conformal Conic proj : lcc datum : etrs89 ellps : grs80 lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33 x_0: 40 y_0: 40 no_defs: defined -PROJ_UNITS unit : metre units : metres meters : 1 g.proj -j +proj=lcc +lat_1=49 +lat_2=46 +lat_0=47.5 +lon_0=13.33 +x_0=40
Re: [GRASS-user] help with testing the location wizard for the upcoming release
since there are so many projections and datums and ways to select them, it's pretty hard to test for working/not working as we usally could do. Please try throwing your local odd-ball CRS at it and see what happens. (check all the way through to 'g.proj -p' and 'g.proj -j' in the final location session. make sure it does what you'd expect it to do. gdalinfo iseldem.tif Driver: GTiff/GeoTIFF Files: iseldem.tif Size is 3117, 2601 Coordinate System is: PROJCS[Lambert Azimuthal Equal Area, GEOGCS[grs80, DATUM[European_Terrestrial_Reference_System_1989, SPHEROID[Geodetic_Reference_System_1980,6378137,298.257222101, AUTHORITY[EPSG,7019]], AUTHORITY[EPSG,6258]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433]], PROJECTION[Lambert_Azimuthal_Equal_Area], PARAMETER[latitude_of_center,52], PARAMETER[longitude_of_center,10], PARAMETER[false_easting,4321000], PARAMETER[false_northing,321], UNIT[metre,1, AUTHORITY[EPSG,9001]]] Origin = (4471032.62800724990,2686202.52683254010) Pixel Size = (22.717005998822735,-22.717005998819641) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 4471032.628, 2686202.527) ( 11d58'53.33E, 47d16'16.87N) Lower Left ( 4471032.628, 2627115.594) ( 11d57'41.52E, 46d44'22.62N) Upper Right ( 4541841.536, 2686202.527) ( 12d54'57.39E, 47d15' 1.05N) Lower Right ( 4541841.536, 2627115.594) ( 12d53'11.78E, 46d43' 7.68N) Center ( 4506437.082, 2656659.061) ( 12d26'11.02E, 46d59'45.72N) location wizard = georeferenced file iseldem.tif g.region -p projection: 99 (Lambert Azimuthal Equal Area) zone: 0 datum: etrs89 ellipsoid: grs80 north: 2686202.52683254 south: 2627115.59422961 west: 4471032.62800725 east: 4541841.53570558 nsres: 22.717006 ewres: 22.717006 rows: 2601 cols: 3117 cells: 8107317 g.proj -p -PROJ_INFO- name : Lambert Azimuthal Equal Area proj : laea datum : etrs89 ellps : grs80 lat_0 : 52 lon_0 : 10 x_0: 4321000 y_0: 321 no_defs: defined -PROJ_UNITS unit : metre units : metres meters : 1 g.proj -j +proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=321 +no_defs +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000 +to_meter=1 - best regards Helmut -- View this message in context: http://osgeo-org.1560.x6.nabble.com/help-with-testing-the-location-wizard-for-the-upcoming-release-tp5055791p5055897.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] help with testing the location wizard for the upcoming release
since there are so many projections and datums and ways to select them, it's pretty hard to test for working/not working as we usally could do. Please try throwing your local odd-ball CRS at it and see what happens. (check all the way through to 'g.proj -p' and 'g.proj -j' in the final location session. make sure it does what you'd expect it to do. testgkat28.prj PROJCS[MGI_Ferro_M28, GEOGCS[GCS_MGI_Ferro, DATUM[D_MGI, SPHEROID[Bessel_1841,6377397.155,299.1528128]], PRIMEM[Ferro,-17.67], UNIT[Degree,0.0174532925199433]], PROJECTION[Transverse_Mercator], PARAMETER[False_Easting,15.0], PARAMETER[False_Northing,0.0], PARAMETER[Central_Meridian,28.0], PARAMETER[Scale_Factor,1.0], PARAMETER[Latitude_Of_Origin,0.0], UNIT[Meter,1.0]] GRASS version: 6.4.3svn GRASS SVN Revision: 56413 GIS Library Revision: 50937 (2012-02-25) GDAL/OGR: 1.9.2 PROJ4: Rel. 4.8.0, 6 March 2012 Python: 2.7.4 wxPython: 2.8.12.1 Platform: Windows-7-6.1.7601-SP1 (OSGeo4W) Select WKT or PRJ file g.region -p projection: 99 (Transverse Mercator) zone: 0 datum: hermannskogel ellipsoid: bessel north: 1 south: 0 west: 0 east: 1 nsres: 1 ewres: 1 rows: 1 cols: 1 cells: 1 g.proj -p -PROJ_INFO- name : Transverse Mercator proj : tmerc datum : hermannskogel ellps : bessel lat_0 : 0 lon_0 : 28 k : 1 x_0: 15 y_0: 0 pm : ferro no_defs: defined -PROJ_UNITS unit : Meter units : Meters meters : 1 g.proj -j +proj=tmerc +lat_0=0 +lon_0=28 +k=1 +x_0=15 +y_0=0 +pm=ferro +no_defs +a=6377397.155 +rf=299.1528128 +towgs84=653.000,-212.000,449.000 +to_meter=1 no asking which transformation should be used. Used in whole hermannskogel region: towwgs84=653.000,-212.000,449.000 is preselected, but not the best e.g. if location used for Austria. - best regards Helmut -- View this message in context: http://osgeo-org.1560.x6.nabble.com/help-with-testing-the-location-wizard-for-the-upcoming-release-tp5055791p5055901.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Getting Started with GRASS Python Scripts
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Il 24/05/2013 21:13, Martin Lacayo ha scritto: Hello Paolo, Yes, my initial thought was to just use GRASS via SEXTANTE. I tired to do this, but did not have any luck. I posted previously on 14/5/2013 to qgis-users with the subject How to run SEXTANTE algorithms outside of QGIS Python console? but did not get a solution to the problem I encountered, which I thought was related to creating a QGIS instance so that I could make calls to SEXTANTE. If you have any suggestions that would be much appreciated. Hi Martin, better ask on qgis-dev or on IRC, most developers read it. All the best. - -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.12 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/ iEYEARECAAYFAlGjMj8ACgkQ/NedwLUzIr4HIwCfSk0xqZOnfaFjc24GaHvWF9gP FvYAn3xKhrBdvuUijTwLdgQIqXw3Ly7H =3q/V -END PGP SIGNATURE- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.stream.order (GRASS 6), Hack, fails at some confluences
Hi everyone, Does anyone have experience with r.stream.order failing at some stream confluences for Hack ordering, such that the rest of stream network upstream gets assigned an order of -1? When I re-run the algorithm starting just upstream of the confluence, it works. I have seen many instances of this type of failure. Please let me know if you have any ideas on why it might be happening and what I can about it. Thank you very much, Hanan This message has been scanned for malware by Websense. www.websense.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Suggestion for creating a script to create a map of agricultural capacity
I'm working on the preparation of the statement of Agricultural Capacity in general terms the *final map* has classes ranging 1-8. To get this map I use three themes with their respective notes: - Map of Slope: with classes ranging 1 to 8; - Soil Map: with classes ranging 1 to 8; - Map of flood risk: with classes ranging 1 to 8; I - After converting all to a vector format and merge the three layers, I have multiple polygons with various combinations, eg: a) 1,3,1 b) 5,2,1 c) 4,2,1 d) 1,1,8 II - To determine the final value of the class is chosen the highest value: a) 3 b) 5 c) 4 d) 8 I need to automate this task because the last step (II) is very slow (actually manually), which would be the best way to do this? I think the whole process can be done with the data in raster format, but have not found a formula that I can apply in r.mapcalc to this problem. Thanks in advance, -- *Marcello Benigno B. de Barros Filho* Prof. do Curso Superior de Tecnologia em Geoprocessamento - IFPB Mestre em Ciências Geodésicas e Tecnologias da Geoinformação - UFPE Doutorando em Tecnologia Ambiental e Recursos Hídricos - UFPE http://profmarcello.blogspot.com http://about.me/marcello.benigno ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Suggestion for creating a script to create a map of agricultural capacity
On 27/05/2013 16:41, Marcello Benigno wrote: I'm working on the preparation of the statement of Agricultural Capacity in general terms the final map has classes ranging 1-8. To get this map I use three themes with their respective notes: Map of Slope: with classes ranging 1 to 8; Soil Map: with classes ranging 1 to 8; Map of flood risk: with classes ranging 1 to 8; I - After converting all to a vector format and "merge" the three layers, I have multiple polygons with various combinations, eg: a) 1,3,1 b) 5,2,1 c) 4,2,1 d) 1,1,8 II - To determine the final value of the class is chosen the highest value: a) 3 b) 5 c) 4 d) 8 Not sure I understand the logic of taking the highest value, but you can do this simply with the original rasters using: r.series input=slope,soils,floods out=agri_risk method=max Then convert to vectors only at the end. I need to automate this task because the last step (II) is very slow (actually manually), which would be the best way to do this? I think the whole process can be done with the data in raster format, but have not found a formula that I can apply in r.mapcalc to this problem. Thanks in advance, -- Marcello Benigno B. de Barros Filho Prof. do Curso Superior de Tecnologia em Geoprocessamento - IFPB Mestre em Cincias Geodsicas e Tecnologias da Geoinformao - UFPE Doutorando em Tecnologia Ambiental e Recursos Hdricos - UFPE http://profmarcello.blogspot.com http://about.me/marcello.benigno This mail was received via Mail-SeCure System. ___ 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. -- Micha Silver GIS Consulting 052-3665918 http://www.surfaces.co.il ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Suggestion for creating a script to create a map of agricultural capacity
Hi Micha, About: *Not sure I understand the logic of taking the highest value * This was defined by an agronomist, I also do not understand why =/ However, your answer is exactly what I need, did not know r.series, thank you! If you can help me also in this tip, I would be very grateful, I'm reclassifying the slope as follows: r.mapcalc class_1 = 'if(slope = 2.0, 1, null())' r.mapcalc class_2 = 'if(2.0 slope = 5, 2, null())' r.mapcalc class_3 = 'if(5 slope = 10, 3, null())' r.mapcalc class_4 = 'if(10 slope = 15, 4, null())' r.mapcalc class_5 = 'if(15 slope = 45, 5, null())' r.mapcalc class_6 = 'if(45 slope = 70, 6, null())' r.mapcalc class_7 = 'if(slope 70, 7, null())' r.mapcalc slope.reclass = class_1 + class_2 + class_3 + class_4 + class_5 + class_6 + class_7 How do I create a single expression for this problem? Regards -- *Marcello Benigno B. de Barros Filho* Prof. do Curso Superior de Tecnologia em Geoprocessamento - IFPB Mestre em Ciências Geodésicas e Tecnologias da Geoinformação - UFPE Doutorando em Tecnologia Ambiental e Recursos Hídricos - UFPE http://profmarcello.blogspot.com http://about.me/marcello.benigno ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user