Re: [Qgis-user] Python code to select polygon at location - SOLVED
Hi, Just by way of finality: The CRS of the polygon SHP was not what it was purported to be. I found that the .prj file was invalid when opening it in QGIS - maybe next time I will first look at the maps I am given, rather than just dive into python scripting. All working as expected, now. Thanks to those who gave ideas to try. Kind regards, Zoltan On 2024/08/25 19:44, Zoltan Szecsei via QGIS-User wrote: So as updated plan-of-attack as below. lyr_poly = QgsVectorLayer(poly_shp, '', 'Ogr') lyr_poly_si = QgsSpatialIndex(lyr_poly.getFeatures(), flags=QgsSpatialIndex.FlagStoreFeatureGeometries) # EPSG:4326 point = QgsPointXY(lon, lat) point_rect = QgsGeometry.fromPointXY(point).buffer(0.001,5).boundingBox() # feature_ids = lyr_poly_si.intersects(point_rect)# intersects wants rectangle requests = QgsFeatureRequest().setFilterFids(feature_ids) # feature_idx = lyr_poly_si.nearestNeighbor(point,10,maxDistance=0.001) requestx = QgsFeatureRequest().setFilterFids(feature_idx) For the above code the lists feature_ids and feature_idx are always empty. The Longitude and Latitude coords are correct, also for the box created from the Point. I'd be grateful for any insight on why nothing is found. Thanks in advance, Zoltan On 2024/08/25 16:19, Zoltan Szecsei via QGIS-User wrote: Hi, I'm doing my "once a year need to do something in python with QGIS" :-/ Using QGIS 3.38 and PyCharm. I have a SHP file with 28000 polygons. *A SPREADSHEET file with 22000 Points (for which I am making a SHP file and wanting attribute vaues from underlying polygons)* Both in EPSG:4326 Without using QGIS processing, what is the quickest way to, for each point 1 by 1, find the underlying polygon and read the attribute field values of that polygon? The code snippet below is far too inefficient - but I cannot find code that selects at XY and benefits from any spatial index on the polygons. lyr_poly = QgsVectorLayer(polygon_shp, '', 'Ogr') point = QgsPointXY(lon, lat) request = QgsFeatureRequest().setFilterRect(lyr_poly.extent()) for feature in lyr_poly.getFeatures(request): geom = feature.geometry() if geom.contains(QgsGeometry.fromPointXY(point)): lyr_poly.select(feature.id()) print(f"Feature with ID {feature.id()} is selected.") break Thanks in advance, Zoltan ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Locating Polygon at Point
Hi Mike, Create the rectangle AT the POI. You have not said what your environment is (Python or GUI), but you can try this from the Processing Toolbox: Rectangles, ovals, diamonds Creates rectangle, oval or diamond-shaped polygons from the input point layer using specified width, height and (optional) rotation values. Multipart inputs should be promoted to singleparts first. HTH, Zoltan On 2024/08/26 00:51, Michael Jabot via QGIS-User wrote: I am trying to create a series of study sites for a citizen science project but I am stuck. I would like to create a 3km X 3km rectangle and then be able to move this so that the point where the two diagonals of the rectangle meet is over a point of interest. I have tried everything I can think of so appreciate any suggestions anyone has. Thank you Mike * Michael Jabot, Ph.D. SUNY Distinguished Professor, Science Education Director, Institute for Research in Science Teaching State University of New York at Fredonia US Partner - NASA GLOBE Program NASA Earth Ambassador Fulbright Specialist - Geospatial Technologies Chancellor's Award Recipient for Excellence in Research 716.320-0189 (Google Voice) ja...@fredonia.edu * ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Python code to select polygon at location
So as updated plan-of-attack as below. lyr_poly = QgsVectorLayer(poly_shp, '', 'Ogr') lyr_poly_si = QgsSpatialIndex(lyr_poly.getFeatures(), flags=QgsSpatialIndex.FlagStoreFeatureGeometries) # EPSG:4326 point = QgsPointXY(lon, lat) point_rect = QgsGeometry.fromPointXY(point).buffer(0.001,5).boundingBox() # feature_ids = lyr_poly_si.intersects(point_rect)# intersects wants rectangle requests = QgsFeatureRequest().setFilterFids(feature_ids) # feature_idx = lyr_poly_si.nearestNeighbor(point,10,maxDistance=0.001) requestx = QgsFeatureRequest().setFilterFids(feature_idx) For the above code the lists feature_ids and feature_idx are always empty. The Longitude and Latitude coords are correct, also for the box created from the Point. I'd be grateful for any insight on why nothing is found. Thanks in advance, Zoltan On 2024/08/25 16:19, Zoltan Szecsei via QGIS-User wrote: Hi, I'm doing my "once a year need to do something in python with QGIS" :-/ Using QGIS 3.38 and PyCharm. I have a SHP file with 28000 polygons. *A SPREADSHEET file with 22000 Points (for which I am making a SHP file and wanting attribute vaues from underlying polygons)* Both in EPSG:4326 Without using QGIS processing, what is the quickest way to, for each point 1 by 1, find the underlying polygon and read the attribute field values of that polygon? The code snippet below is far too inefficient - but I cannot find code that selects at XY and benefits from any spatial index on the polygons. lyr_poly = QgsVectorLayer(polygon_shp, '', 'Ogr') point = QgsPointXY(lon, lat) request = QgsFeatureRequest().setFilterRect(lyr_poly.extent()) for feature in lyr_poly.getFeatures(request): geom = feature.geometry() if geom.contains(QgsGeometry.fromPointXY(point)): lyr_poly.select(feature.id()) print(f"Feature with ID {feature.id()} is selected.") break Thanks in advance, Zoltan ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Python code to select polygon at location
Ah - bad explanation on my part. [apologies] I have the coordinates from a spreadsheet and I am busy making the Point SHP (but wanting some attribute values from underlying polygons) On 2024/08/25 18:47, oisin.kelly.wild--- via QGIS-User wrote: For the point layer create a calculated field using overlay_within() function. Oisin 07570 977449 Get BlueMail for Android On 25 Aug 2024, at 16:07, Greg Troxel via QGIS-User wrote: Zoltan Szecsei via QGIS-User writes: Hi, I'm doing my "once a year need to do something in python with QGIS" :-/ Using QGIS 3.38 and PyCharm. I have a SHP file with 28000 polygons. Another SHP *spreadsheet* file with 22000 Points Both in EPSG:4326 Without using QGIS processing, what is the quickest way to, for each point 1 by 1, find the underlying polygon and read the attribute field values of that polygon? Load it into postgis and write a select statement? QGIS-User mailing list QGIS-User@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- ===== Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Python code to select polygon at location
It's part of a much longer python script doing a bunch of other stuff, so to now change the flow to incorporate postgis is not really an option. (but thanks) On 2024/08/25 18:07, Greg Troxel via QGIS-User wrote: Zoltan Szecsei via QGIS-User writes: Hi, I'm doing my "once a year need to do something in python with QGIS" :-/ Using QGIS 3.38 and PyCharm. I have a SHP file with 28000 polygons. Another SHP file with 22000 Points Both in EPSG:4326 Without using QGIS processing, what is the quickest way to, for each point 1 by 1, find the underlying polygon and read the attribute field values of that polygon? Load it into postgis and write a select statement? ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user -- ===== Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] Python code to select polygon at location
Hi, I'm doing my "once a year need to do something in python with QGIS" :-/ Using QGIS 3.38 and PyCharm. I have a SHP file with 28000 polygons. Another SHP file with 22000 Points Both in EPSG:4326 Without using QGIS processing, what is the quickest way to, for each point 1 by 1, find the underlying polygon and read the attribute field values of that polygon? The code snippet below is far too inefficient - but I cannot find code that selects at XY and benefits from any spatial index on the polygons. lyr_poly = QgsVectorLayer(polygon_shp, '', 'Ogr') point = QgsPointXY(lon, lat) request = QgsFeatureRequest().setFilterRect(lyr_poly.extent()) for feature in lyr_poly.getFeatures(request): geom = feature.geometry() if geom.contains(QgsGeometry.fromPointXY(point)): lyr_poly.select(feature.id()) print(f"Feature with ID {feature.id()} is selected.") break Thanks in advance, Zoltan ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Rendering shared exterior boundaries of adjacent polygon features
Hi, Also for clarification: Being a County layer, I would not expect overlapping polygons. My take is that Tim means "concurrent exterior boundaries" and not "overlapping exterior boundaries" *Tim: *What is the reason for wanting only one of the concurrent(?) boundaries rendered? Maybe there are other ways to do what you want to. HTH, Zoltan On 2024/06/20 07:46, Emma Hain via QGIS-User wrote: Hi Tim So let's see if I understand you. - There are 2 polygon layers (green and black) - These polygon layers sometime overlap - You wish to show only on the pink line in the below image and not the green overlapping part? image.png On Thu, 20 Jun 2024 at 11:26, Tim Malloy via QGIS-User wrote: Fellow QGIS Users – I’m trying to determine if its possible to have QGIS render the exterior boundary of polygon features only once when that boundary is shared with an adjacent feature (don’t render overlapping exterior boundaries for both polygon features).The data I’m working with is a county parcels layer and the symbology I’m using for this layer is a cross hatched line fill with a separate component for the outline of the polygons.I’ve tried a couple of expressions (examples below) including one geometry generator and one data defined override for the stroke width of the boundary line but none of them achieves my objective.I’m using QGIS 3.34.7 and my coding skills are limited.I feel like I’ve seen this issue covered somewhere but have been unable to find a solution after extensive online exploration. *Method #1:“*boundary( difference( $geometry, aggregate( 'union', 'intersects', $geometry, other_layer := @layer)))”This expression was used on the outline component of the symbology described above. *Method #2:*“CASE WHEN intersects( $geometry, geom_to_wkt( boundary( $geometry ))) THEN 0 ELSE 1 END”Here I switched to using an expression on the stroke width of a simple fill, but I got an error message saying “Cannot convert to geometry”. Any input would be greatly appreciated. Tim ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user -- Emma Hain — Product Manager/Senior GIS Analyst e...@north-road.com https://north-road.com *North Road* Cartography • Development • Spatial Analysis *north-road.com* <http://north-road.com> <https://twitter.com/northroadgeo><https://www.linkedin.com/company/north-road-studios><https://www.facebook.com/North-Road-997236690392419/home> ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Please help!
Oops - yes, quite correct. Regards, Zoltan On 2024/04/23 10:01, Nigel Berjak - S3 Technologies wrote: Hi Zoltan Apologies, the CRS is ESRI 102567 and not EPSG 102567, for South Africa, Lo29° metres projection based upon the North/Easting orientation of the data. Cheers. --- Regards, Nigel Berjak S3 Technologies Geographic Information Systems & Large Format Printing specialists T: +27 33 3423681 F: +27 86 6721242 E: ni...@s3.co.za Website: http://www.S3.co.za <http://www.S3.co.za> Please consider the environment before printing this email. On 2024-04-22 10:44, Zoltan Szecsei via QGIS-User wrote: Hi, I'm guessing the '2053' might be wrong. You need to know if your DEM coords are +ve East and North, or +ve West and South If +ve E/N then try EPSG:102567 I forget which of the remaining 2 (2053 and 102487) Lo29 codes are left or right handed, but you only have these 3 to try. When your DEM lands in the correct place in South Africa, then worry about re-projecting to UTM35S HTH, Zoltan On 2024/04/22 10:32, Stefan Bogenhofer via QGIS-User wrote: Good day I am having some trouble in regards to georeferencing a DEM which I received. The CRS of the DEM is "EPSG:2053 - Hartebeesthoek94 / Lo29" but when I reproject the DEM to the desired projection "WGS84 / UTM zone 35S" it puts the DEM somewhere in Egypt although the correct position where it should be in my study area is located is in South Africa. I have tried assigning new coordinate systems to both my project, web image layer and the DEM but it simply does not want to reference it to the correct position. Any help would be much appreciated! Thank you ___ QGIS-User mailing list QGIS-User@lists.osgeo.org <mailto:QGIS-User@lists.osgeo.org> List info:https://lists.osgeo.org/mailman/listinfo/qgis-user <https://lists.osgeo.org/mailman/listinfo/qgis-user> Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user <https://lists.osgeo.org/mailman/listinfo/qgis-user> -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za <http://www.geograph.co.za> = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Please help!
Hi, I'm guessing the '2053' might be wrong. You need to know if your DEM coords are +ve East and North, or +ve West and South If +ve E/N then try EPSG:102567 I forget which of the remaining 2 (2053 and 102487) Lo29 codes are left or right handed, but you only have these 3 to try. When your DEM lands in the correct place in South Africa, then worry about re-projecting to UTM35S HTH, Zoltan On 2024/04/22 10:32, Stefan Bogenhofer via QGIS-User wrote: Good day I am having some trouble in regards to georeferencing a DEM which I received. The CRS of the DEM is “EPSG:2053 - Hartebeesthoek94 / Lo29” but when I reproject the DEM to the desired projection “WGS84 / UTM zone 35S” it puts the DEM somewhere in Egypt although the correct position where it should be in my study area is located is in South Africa. I have tried assigning new coordinate systems to both my project, web image layer and the DEM but it simply does not want to reference it to the correct position. Any help would be much appreciated! Thank you ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- ===== Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] "QgsProcessing" has no attribute 'run' (same for 'processing')
Nice attitude, Hannes :-) Yes, totally true about the AI trap. You can't write a program by handing the logic over to a "well read entity" I use it when my mind is plain out of google keywords. In maybe half a dozen attempts with AI, I have only hit gold once. Thanks for the links. Regards, Zoltan On 2024/03/18 10:31, Johannes Kröger (WhereGroup) wrote: Hi Zoltan, the most important thing first: Don't fall into the "AI can program" mindset trap. You *will* end up with nonsense like this if you do. Those tools can often hallucinate non-existent APIs and mislead you. They might work for some more basic or wide-spread code bases but so far I have seen very bad results for PyQGIS. Best search the web using "pyqgis" as additional keyword to your query. As https://qgis.org/pyqgis/3.34/core/QgsProcessing.html says, that class only contains data to be used by Processing things. It is not your entry point into using Processing via Python. For that please see https://docs.qgis.org/3.34/en/docs/user_manual/processing/console.html , written by actual, smart humans ;) And for the Batch file: Don't write your own, QGIS comes with a perfect one already, called |python-qgis[-ltr].bat| . See e.g. https://gis.stackexchange.com/a/428577/51035 Cheers and good luck! Hannes On 11.03.24 09:30, Zoltan Szecsei via QGIS-User wrote: Hi, Apologies for asking help on what hopefully might turn out to be trivial, but I am pretty much blue in the face googling and asking Gemini for help. There are 2 questions in this email: I have a Windows 11 OSGEO4W installation QGIS 3.36.0 and grass (and more). I use PyCharm 2023.3.4 (Community Edition) - but not yet configured "nicely" for QGIS as I am mostly in a non-Open Source environment. Most of the stuff I write is python is for "problem fixes" and these python scripts ideally run on a windows command-line with a BAT file. So, in this scenario below, I keep getting: . *has no attribute 'run' *This whether I am running inside PyCharm or on the command-line through below batch job. *Question #1* *What is the correct bat file structure for running python scripts?* I currently have [below] but is it optimal? SET QGIS_PATH=C:\OSGeo4W\apps\qgis call "C:\OSGeo4W\bin\o4w_env.bat" SET PYTHONPATH=%QGIS_PATH%\python;%PYTHONPATH% SET PATH=%QGIS_PATH%\bin;%PATH% C:\OSGeo4W\apps\Python39\python.exe Z:\Projects\Development\python_QGIS\qgis_tests\qgisrun.py pause *Question #2* I have tried many various import options, but I cannot get past the *has no attribute 'run'* error from python below code: Minimalst as this code is, it should run without adding anything else, should it? from qgis.core import QgsProcessing shp_inErfPoly = 'Z:\\Erven_1February2024.shp' shp_outErfFix = 'Z:\\Erven_1February2024_FixGeo.shp' QgsProcessing.run("native:fixgeometries", { 'INPUT': shp_inErfPoly, 'METHOD': 1, 'OUTPUT': shp_outErfFix}) Some pointers on what to do or where to read would be great. Thanks in advance, Zoltan -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- Johannes Kröger / GIS-Entwickler/-Berater * FOSSGIS Konferenz 20.-23. März 2024 in Hamburg https://www.fossgis-konferenz.de/2024/ WhereGroup-Beiträge auf der FOSSGIS https://wheregroup.com/unternehmen/aktuelles/ * WhereGroup GmbH c/o KK03 GmbH Lange Reihe 29 20099 Hamburg Germany Tel: +49 (0)228 / 90 90 38 - 36 Fax: +49 (0)228 / 90 90 38 - 11 johannes.kroe...@wheregroup.com www.wheregroup.com Geschäftsführer: Olaf Knopp, Peter Stamm Amtsgericht Bonn, HRB 9885 --- -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] "QgsProcessing" has no attribute 'run' (same for 'processing')
Hi, Apologies for asking help on what hopefully might turn out to be trivial, but I am pretty much blue in the face googling and asking Gemini for help. There are 2 questions in this email: I have a Windows 11 OSGEO4W installation QGIS 3.36.0 and grass (and more). I use PyCharm 2023.3.4 (Community Edition) - but not yet configured "nicely" for QGIS as I am mostly in a non-Open Source environment. Most of the stuff I write is python is for "problem fixes" and these python scripts ideally run on a windows command-line with a BAT file. So, in this scenario below, I keep getting: . *has no attribute 'run' *This whether I am running inside PyCharm or on the command-line through below batch job. *Question #1* *What is the correct bat file structure for running python scripts?* I currently have [below] but is it optimal? SET QGIS_PATH=C:\OSGeo4W\apps\qgis call "C:\OSGeo4W\bin\o4w_env.bat" SET PYTHONPATH=%QGIS_PATH%\python;%PYTHONPATH% SET PATH=%QGIS_PATH%\bin;%PATH% C:\OSGeo4W\apps\Python39\python.exe Z:\Projects\Development\python_QGIS\qgis_tests\qgisrun.py pause *Question #2* I have tried many various import options, but I cannot get past the *has no attribute 'run'* error from python below code: Minimalst as this code is, it should run without adding anything else, should it? from qgis.core import QgsProcessing shp_inErfPoly = 'Z:\\Erven_1February2024.shp' shp_outErfFix = 'Z:\\Erven_1February2024_FixGeo.shp' QgsProcessing.run("native:fixgeometries", { 'INPUT': shp_inErfPoly, 'METHOD': 1, 'OUTPUT': shp_outErfFix}) Some pointers on what to do or where to read would be great. Thanks in advance, Zoltan -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) +974-5058-8729 +36-20-4808362 www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Running QGis on a machine with Python 3.11 and PyCharm ?
Hi Toussaint, I am not a hard-core developer, but I script/code whenever I can automate a procedure, so: * It is my understanding that QGIS is developed/maintained in the QT environment - so go that way rather than PyCharm if you plan to dabble in QGIS development. Perhaps one of the devs that lurk on this user list will add to this comment of mine. * Get on to QGIS-dev as well - you might find some better guidance there, for PyCharm and other python scripting issues. I will try to find time to set up QGIS and PyCharm on my laptop in the next day or so - *_but in the meantime, google_* "pycharm with qgis" - there are about a billion articles and videos that come up. Oh!, Yes!, I am exaggerating. :-) Regards, Zoltan On 2023/06/12 10:06, Toussaint OTTAVI via QGIS-User wrote: Hi Zoltan, Thank you in advance for your help. For now, I didn't manage to find a suitable configuration by setting paths and environment variables in a startup ".bat" file. QGis still tries to load things from my system-wide Python folder (c:\python), even if no path or env variable points to it. Then, it fails. The only workaround I found is a bit rough : rename my c:\python folder to c:\python.bak, run QGis, then rename back to the original name ! Of course, this prevents from using QGis at the same time with PyCharm or other Python stuff, and is not an acceptable solution on the long term. Even if I don't have the skills to dive into QGis source code for now, I may try to play around with some simple plugins. Then, wouldn't it be a better idea to install QGis in "developer" mode, by cloning / forking the GitHub repo ? Any suggestions are welcome. Kind regards. Le 09/06/2023 à 22:25, Zoltan a écrit : Hi, I seem to remember that pycharm website provides help on using pycharm with qgis, so don't mess with your existing installations of python. The gist of it is that you add the pycharm startup to the end of the bat file that qgis provides for setting up all the env variables. Starting pycharm from the Taskbar will not set these environment for you so use the edited bat file I describe above. There might also be a way in your pycharm venv to pre set these env vars. I have a new laptop which I have not yet set up for this but if you are still stuck tomorrow, I'll set mine up and let you know what works. Regards Zoltan ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +974-5058-8729 (Signal, not WhatsApp) +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Running QGis on a machine with Python 3.11 and PyCharm ?
Hi, I seem to remember that pycharm website provides help on using pycharm with qgis, so don't mess with your existing installations of python. The gist of it is that you add the pycharm startup to the end of the bat file that qgis provides for setting up all the env variables. Starting pycharm from the Taskbar will not set these environment for you so use the edited bat file I describe above. There might also be a way in your pycharm venv to pre set these env vars. I have a new laptop which I have not yet set up for this but if you are still stuck tomorrow, I'll set mine up and let you know what works. Regards Zoltan Sent from Blue On 09 Jun 2023, 22:45, at 22:45, Toussaint OTTAVI via QGIS-User wrote: > > >Le 09/06/2023 à 18:11, Toussaint OTTAVI via QGIS-User a écrit : >> If I understand well, I'll have to write a customized startup script, > >> with correct path (not pointing to my c:\python\lib) so that QGis can > >> start using the included Python interpreter (and not my system-wide >> one). But I didn't find how to achieve that. > >For now, the only workaround I found to start QGis is to remove >completely my system-wide "c:\python" directory. >Trying to configure path and PYTHONPATH with priority to the QGis >folder >(such as the included .bat files do) does not work. Doing so, and using > >SysInternals Process Monitor, we can see it still tries to load things >from my "c:\python" directory, and it fails because there are >incompatible library versions. > >Of course, removing my system-wide Python installation is not an >acceptable workaround ;-) > >How to force QGis to use only the Python stuff in its folders and >sub-folders ? Virtualenv ? > > > >___ >QGIS-User mailing list >QGIS-User@lists.osgeo.org >List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] How to set PyQGIS imported in VS Code in Windows for a Python file?
Hi Li I am responding because I had a similar experience using PyCharm so below _might_ be useful things for you to check whilst you are waiting for a more informed response: Firstly Google using QGIS is VS - there will be a lot of environment variables that will need to be set up. I start PyCharm using the following .bat job: HTH, Zoltan REM https://haifengniu.com/post/how-to-set-environment-for-pyqgis-3-in-pycharm/ REM Setting PyCharm IDE REM Just double-click pycharm.cmd to start PyCharm. REM In the settings of PyCharm project, set interpreter as %OSGEO4W_ROOT%\bin\python3.exe REM In my computer, the full path of interpreter is D:\OSGeo4W64\bin\python3.exe. @echo off SET OSGEO4W_ROOT=C:\OSGeo4W call "%OSGEO4W_ROOT%"\bin\o4w_env.bat REM call "%OSGEO4W_ROOT%"\bin\qt5_env.bat REM call "%OSGEO4W_ROOT%"\bin\py3_env.bat call "%OSGEO4W_ROOT%"\apps\grass\grass78\etc\env.bat @echo off path %PATH%;%OSGEO4W_ROOT%\apps\qgis\bin path %PATH%;%OSGEO4W_ROOT%\apps\grass\grass78\lib path %PATH%;D:\OSGeo4W64\apps\Qt5\bin path %PATH%;D:\OSGeo4W64\apps\Python39\Scripts set PYTHONPATH=%PYTHONPATH%;%OSGEO4W_ROOT%\apps\qgis\python set PYTHONHOME=%OSGEO4W_ROOT%\apps\Python39 set QT_PLUGIN_PATH=%OSGEO4W_ROOT%\apps\qgis\qtplugins;%OSGEO4W_ROOT%\apps\qt5\plugins start "Pycharm Aware of OSGeo4W" /B "C:\\Program Files\\JetBrains\\PyCharm Community Edition 2021.2.2\\bin\\pycharm64.exe" On 2023/03/13 09:30, LI, Xin Stephen [MIT] via QGIS-User wrote: The OS is Windows 10, The QGIS version is 3.28 LTR The Python version is 3.11 The steps following the steps in stection 1.4.1 for the PyQGIS Cookbook (GQIS Testing), i.e. https://docs.qgis.org/testing/en/docs/pyqgis_developer_cookbook/intro.html#python-applications The “PAHT” and “PYTHONPYTHON” have been added. The problem is: When running the Python program file with “import qgis.core”, there is the error “ModuleNotFoundError: No module named 'qgis' ” Moreover, after reinstalling the Python 3.9.5 as the QGIS related, then run the Python file in VS Code by Python 3.9.5, the Error is the same. Any suggestions could help modify this setting? Thanks. ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +974-5058-8729 (Signal, not WhatsApp) +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ QGIS-User mailing list QGIS-User@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] Trapping Kyb/Mouse events during a QGIS geometry edit process
Hi Listers, I'm looking for some pointers as where to start, or keywords to google or recipes to follow (if any). I'm looking to trap kyb/mouse events whilst in the middle of a geometry edit. If a suitable event is trapped, then to run a subroutine and then return to the edit process by either handing the edit process the same event or generating another event for the edit routine to continue the geometry edit with. So I suppose my question is, is where do I insert such kyb/mse event filter. Any pointers would be welcome. Thanks in advance, Zoltan -- ===== Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Classify road network into dead-end links and links connected on both ends
that's why SHP files need to die: They're too simple to use, so also get used by people less "in the know" On 2022-05-02 11:05, Roland Spielhofer wrote: I totally agree with you, but the shape is all I got from the client - nothing better is available at his side... *Gesendet:* Montag, 02. Mai 2022 um 10:49 Uhr *Von:* "Zoltan via Qgis-user" *An:* qgis-user@lists.osgeo.org *Betreff:* Re: [Qgis-user] Classify road network into dead-end links and links connected on both ends Sigh, this is where SHP-file hell comes in - it's the scourge of GIS. Wouldn't it have been nice if your road network was in a data structure that allows you to go to each point and ask it "how many lines emanate from me"? Sorry, I know this is not helpful to you, but the SHP file format needed to have died 20 years ago. Regards, Zoltan On 2022-05-02 10:35, Roland Spielhofer via Qgis-user wrote: Hi, I have a road network (Shape linestring) where I would like to classify all links that are dead-end and compute the share of the dead-end links in relation to the total network length. I am pretty sure this is possible in QGIS or a plugin - any pointers would be appreciated. Regards, Roland ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- ===== Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user -- ===== Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Classify road network into dead-end links and links connected on both ends
On a more helpful note, you could add 2 columns to your line-string attribute table and then iterate through the start and end coords of each line, counting the number of lines selected at their endpoints. HTH, Zoltan On 2022-05-02 10:49, Zoltan via Qgis-user wrote: Sigh, this is where SHP-file hell comes in - it's the scourge of GIS. Wouldn't it have been nice if your road network was in a data structure that allows you to go to each point and ask it "how many lines emanate from me"? Sorry, I know this is not helpful to you, but the SHP file format needed to have died 20 years ago. Regards, Zoltan On 2022-05-02 10:35, Roland Spielhofer via Qgis-user wrote: Hi, I have a road network (Shape linestring) where I would like to classify all links that are dead-end and compute the share of the dead-end links in relation to the total network length. I am pretty sure this is possible in QGIS or a plugin - any pointers would be appreciated. Regards, Roland ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- ===== Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- ===== Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Classify road network into dead-end links and links connected on both ends
Sigh, this is where SHP-file hell comes in - it's the scourge of GIS. Wouldn't it have been nice if your road network was in a data structure that allows you to go to each point and ask it "how many lines emanate from me"? Sorry, I know this is not helpful to you, but the SHP file format needed to have died 20 years ago. Regards, Zoltan On 2022-05-02 10:35, Roland Spielhofer via Qgis-user wrote: Hi, I have a road network (Shape linestring) where I would like to classify all links that are dead-end and compute the share of the dead-end links in relation to the total network length. I am pretty sure this is possible in QGIS or a plugin - any pointers would be appreciated. Regards, Roland ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- ===== Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] RAM and processing speed requirements for QGIS
Hi Helena/List, For all the good intentions everyone who has responded to you, have had - they are truly "not that useful" So, before I get shot: Yes, more memory, HD->SSD->NVMe does help a lot - but did Mac's have M.2 slots 7 years back? If yes, were they configured as NVMe or SATA? Is "new" Mac memory different from 7-year-old Mac memory? (yes, of course it is) So let's get pragmatic: Buying any of the suggested HW to upgrade might be throwing good money after bad - because if they do not help you sufficiently, you will still have to buy a newer PC, and you may end up being stuck with components you cannot use in your newer PC. So, judging by your original post below - you have not even tried QGIS yet. It is free! Download it. Try it. If your electricity bill goes up because you are boiling your kettle more often than normal (while waiting for QGIS to respond), then buy a new PC - that is, if you are unable to_borrow _extra memory to see if that helps sufficiently. If you buy a new PC, then yes, look at all the responses you have received. Buy at least 16GB memory - preferably 32GB Do not buy any 'old' HDs (those with moving parts) unless you need >4TB of disk space. (very easy with GIS data) - but keep system & swap on the fastes SSDs you can afford. If you buy SSDs, make sure they are NVMe and not SATA-SSDs. No matter what GIS SW you choose, it will need a 'decent' PC if you are asking it to do stuff more than just play around with simple maps. Hope this helps you with your way forward. Kind regards, Zoltan On 2022-04-20 14:00, Patrick Dunford via Qgis-user wrote: As others have written it will depend on the tasks you are doing. The performance of any computer can be significantly improved by replacing HDDs with SSDs. Modern computers using the NVMe standards for SSD (M.2 interface etc) can offer a great deal of speed improvement for SSD storage over SATA which is very significant when working on intensive projects that use a lot of resources. The PC platform is significantly easier to expand if resource needs change. My seven year old system started out with 8 GB of RAM and currently has 32 GB, whilst the swap file performance (virtual memory) has been significantly improved with SATA and M.2 SSD storage (now at 200 GB). The latter provides almost limitless performance capability enhancements to the system at modest cost without having to upgrade the entire system. I use Ubuntu as the operating system for this computer. On 20/04/22 06:46, Helena Farrell via Qgis-user wrote: Hello, I am seeking input on whether my 7 year old Mac with 8GB RAM and a 2.7gigahertz processor is likely to perform poorly at running QGIS. Since these are like the bare minimum in terms of power, and I assume QGIS is a data-heavy program, I may need to invest in a new computer. Anyone have advice on RAM and processing speed specifications that I should look for and whether a PC be better than a Mac for running this program? Thank you! ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] again: Create Geospatial PDF producing faulty output
Hi Chris, I have not followed this thread in any detail, but as far as your labels comment below is concerned, depending on the output you choose, there will be slight changes in extents and scale, which could affect the labels becoming rendered or not. Maybe try playing with options that enforce displaying of all labels and then see if you labels still disappear? HTH, Zoltan On 2022-03-21 17:49, chris hermansen via Qgis-user wrote: Uwe and list, On Mon, Mar 21, 2022, 10:42 Uwe Fischer via Qgis-user wrote: Hello dear friends on the list, I’m sorry but I have to ask again about „Creating Geospatial PDF“ from the 3.22 Print Composer: does nobody have problems like missing features in the PDF output? The problem is very important for me, so any idea what the reason is would be appreciated. Not a single answer until now. The problem description: some features (in my case line features) are omitted when QGIS writes the geospatial PDF. Normal PDF output works fine, as well as TIF image output. So the problem cannot be in the data. Please help! No help I am afraid, but I do notice that labels appearing both on main screen and in print composer are sometimes dropped on export to image, whether from main screen or print composer. I have a simple demo in case anyone is interested. This is with latest version on Ubuntu. Chris ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- = Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Center line of an elongated BASIN
Hi, I'm coming in late to this methodology discussion, so apologies if I am off the mark. So by 'center line' I'm guessing you mean a line that is always equidistant from the polygon boundary. I had a similar problem a few years back where I had to create road center lines from a polygon road representation for a small country. After some experimentation on speed and efficacy, I came up with a Python script that: merged (dissolved) all poygons) - but you only have 1 polygon :-) rasterised the polygon (as finely as you choose relevant to your issue) ran r.thin vectorised the (now) single line of cells. ran a Douglas-Peuker style smoothing on the staircase lines Worked like a dream for my purpose. Maybe you could adapt this logic flow to your needs? HTH, Zoltan On 2022-02-11 00:46, Allan López wrote: Yes, it does. Thanks a lot indeed. I am testing it. Best, Allan El jue, 10 feb 2022 a la(s) 15:19, Alexis R.L. (alroyl...@gmail.com) escribió: As far as I know v.voronoi.skeleton in gdal is the closest thing we have that can do this in base QGIS. I hope this helps, Alex Le mer. 9 févr. 2022 à 17:57, Allan López a écrit : Dear All, I need to calculate and display the center line of an elongated river basin ( not the center of the main river trunk) from a SHP. Thank's in advance. Allan Lopez ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- ===== Zoltan Szecsei GPrGISc 0031 Director, Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Mosaic Datasets and Raster Catalogs in qGIS
Hi Jason, You just need to create a VRT. The vrt is just that: a virtual raster table. Fiddle with the options and notice that you can mix rasters that have differing projection. Test out on maybe a 100 or so Raster files and when you get the feel for it, dump the lot in. I've had >5000 files in one with no hassles. Note that if you output your vrt file into the same folder as your rasters, when you recreate your vrt, it will read in your old one and go into a loop. I logged this as a bug and was told unceremoniously that it was not a bug. Btw, when you first load your vrt into qgis, have rendering turned off. Then set your zoom display scales before turning rendering on again. This will stop qgis from trying to display all your tasters if you happen to be zoomed out too far. Have fun, it works well. Regards Zoltan Sent from BlueMail On 14 Jan 2022, 21:45, at 21:45, Jason Aleksi wrote: >I am trying to add imagery in qGIS. In ArcMap, I had to create a >Mosaic >Dataset, then a Raster Catalog. It tooks days to create due to the >large >files and quantity of files. I have 6 inch imagery of a county in a >folder. The GeoTIFFs are about 308GB in size and approx 6,610 files >(TIF >and TIFW). When I attempt to create a raster layer, it will only pull >in >100 images. > >Will I have to create 305 Layers with 100 images each? Or is there a >better way? I'm at a loss for terms to even search for. I feel I've >exhausted Google! If there is a better way, could you possibly point >me in >the right direction? > >FTW: I also have MrSID Tiles and Composites, in a separate directory. > > > > >___ >Qgis-user mailing list >Qgis-user@lists.osgeo.org >List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Wrong area count QGIS 3.20+?
Hi Bas, Please say which version of ArcPro - I might be heading for these type of issues in the next few months. Thanks and regards, Zoltan On 2022-01-03 16:46, Sebastian Gutwein wrote: Ok to add to the fun I tried the calculations in ArcPro and interestingly if I run it using calculate field with this expression *!shape.geodesicArea@hectares!* I get the QGIS 3.22 value of 31271112 but if I run Geoprossessing Add Geometry Attributes I get the QGIS 3.16 value of 31271108. image.png So if you want them to match I guess you should use the python expression in ArcGIS Pro! I guess we are back to the old adage of All Measuring is Guessing. -Bas On Mon, Jan 3, 2022 at 8:43 AM Sebastian Gutwein wrote: Yes that was my read as well except previously it was the Grass method and now it is Geographic Lib method because that is theoretically more accurate and is the same method that is used in Postgis and other places. The switched happened because someone noticed the difference between QGIS’s calculations and Postgis’s calculations. If I am reading the pull requests correctly even Grass has switched to the Geographic Lib method. I’m surprised that ArcGIS Pro uses the same method as the Grass method. On Mon, Jan 3, 2022 at 2:57 AM Sadowski Jarosław wrote: Sebastian, thank you very much for your precise answer and materials sent - I have read them carefully. So, as far as I can understand, QGIS has moved in version 3.20 to a new way of counting areas (other libraries? GRASS algorithm instead of Geographic lib? I have actually noticed a different nomenclature and definition of EPSG 2180 for some time in new versions…), that will be more accurate? I must admit that, for me as a GIS operator, it is surprising and sounds like bad news - so far surfaces have always been the same in various GIS software, whether or not it was open-source or commercial (ArcGIS Pro). So from this moment it will be that I will get different surfaces depending on the software it was calculated? :( with kind regards JS *_* *Jarosław Sadowski* Ekspert ds. Baz Danych i Analiz Przestrzennych | /Biuro Strategii i Planowania, Projektowania i Inżynierii Podprogramu Kolejowego/ /Database and Spatial Analysis Expert |Railway Subprogramme Strategy & Planning, Design & Engineering Department/ e: jaroslaw.sadow...@cpk.pl *From:* Sebastian Gutwein *Sent:* Saturday, January 1, 2022 7:13 AM *To:* Sadowski Jarosław *Cc:* qgis-user@lists.osgeo.org *Subject:* Re: [Qgis-user] Wrong area count QGIS 3.20+? Nie otrzymujesz często wiadomości e-mail z adresu b...@rdgland.com. Dowiedz się, dlaczego jest to ważne <http://aka.ms/LearnAboutSenderIdentification> There was some discussion on this in the list serve back in January <https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.osgeo.org%2Fpipermail%2Fqgis-user%2F2021-January%2F047638.html&data=04%7C01%7Cjaroslaw.sadowski%40cpk.pl%7Cf7cc24b07b7c42e8902a08d9ccedca61%7Cfa798250ca0b4a1bb47381cff4a1752b%7C1%7C0%7C637766145287655707%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=bUHOuswrG6tPdkebzFYem%2BVsSrgjFhbSLVcKCiu6TW0%3D&reserved=0> that resulted in this feature request <https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fqgis%2FQGIS%2Fissues%2F40888&data=04%7C01%7Cjaroslaw.sadowski%40cpk.pl%7Cf7cc24b07b7c42e8902a08d9ccedca61%7Cfa798250ca0b4a1bb47381cff4a1752b%7C1%7C0%7C637766145287655707%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=dY89NsDscjybmwQPMV9PjPOwckPlcv%2B8TEWjl69rTtM%3D&reserved=0> and this pull request <https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fqgis%2FQGIS%2Fpull%2F41726&data=04%7C01%7Cjaroslaw.sadowski%40cpk.pl%7Cf7cc24b07b7c42e8902a08d9ccedca61%7Cfa798250ca0b4a1bb47381cff4a1752b%7C1%7C0%7C637766145287655707%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=IRMFsaJYT7xpylDJCZNvXAUNzdmtTjNbI4VyQMKtedE%3D&reserved=0>. I tested your polygon, to the best of my limited ability, and the areas were indeed different in QGIS 3.16 and 3.22 as would be expected from that change . However it would be hoped that the area measurement would be more accurate in 3.22. On Fri, Dec 31, 2021 at 6:49 AM Sadowski Jarosław wrote: Dear all, I am asking for help with QGIS
Re: [Qgis-user] Ubuntu impish (21.10) release coming any time soon?
(Oops I shot from the hip :-) ) Sent from BlueMail On 23 Oct 2021, 08:17, at 08:17, Zoltan wrote: >Other than the ubuntu site, https://distrowatch.com/ has some >interesting info > >Regards, >Zoltan > >Sent from BlueMail > >On 23 Oct 2021, 00:44, at 00:44, chris hermansen > wrote: >>Good day everyone, >> >>Does anyone know when, at least approximately, >>https://qgis.org/ubuntu/dists/ will contain an impish release? >Better, >>can >>anyone tell me where to look for this kind of info? >> >>Thanks very much in advance. >> >>-- >>Chris Hermansen · clhermansen "at" gmail "dot" com >> >>C'est ma façon de parler. >> >> >> >> >>___ >>Qgis-user mailing list >>Qgis-user@lists.osgeo.org >>List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >>Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user > > > > >___ >Qgis-user mailing list >Qgis-user@lists.osgeo.org >List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Ubuntu impish (21.10) release coming any time soon?
Other than the ubuntu site, https://distrowatch.com/ has some interesting info Regards, Zoltan Sent from BlueMail On 23 Oct 2021, 00:44, at 00:44, chris hermansen wrote: >Good day everyone, > >Does anyone know when, at least approximately, >https://qgis.org/ubuntu/dists/ will contain an impish release? Better, >can >anyone tell me where to look for this kind of info? > >Thanks very much in advance. > >-- >Chris Hermansen · clhermansen "at" gmail "dot" com > >C'est ma façon de parler. > > > > >___ >Qgis-user mailing list >Qgis-user@lists.osgeo.org >List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] When running a Custom Function, the return value is different from when stepping through it and watching the 'Preview' value in Field Calculator
Oh "for Pete's sake". - I'd love someone to explain why that is needed and not auto-turned on or defaulted to True. Thank you for bringing my attention to that tidbit. Kind regards, Zoltan On 2021-10-19 19:07, DelazJ wrote: Hi Zoltan, You are playing with geometry so you need usesgeometry=True in your function decorator. See https://docs.qgis.org/3.16/en/docs/user_manual/working_with_vector/expression.html#function-editor or the help section in the Function Editor tab. Regards, Harrissou Le mar. 19 oct. 2021 à 18:33, Zoltan a écrit : Hi, Me again :-( I've written a custom function in Field Calculator and want to put its return value into a new Text field in the attribute table. When I step through the features I get the correct and expected values shown next to 'Preview' But when I click OK to run the Field Calculator, it runs for a long enough time, auto-closes on completion and in the attribute table the new field is created, but value is 'Nothing for None' for all features. *Note:* variable rd_at reflects correct contents when stepping through and observing 'Preview' but seems to revert to Nonetype when running the function. What am I missing regarding using custom functions? Is it an intended thing that the Field Calculator auto-closes after running a custom function? Thanks again, Kind regards, Zoltan *Note:* The msg variable was introduced in case placing 'return' in various lines in the script was causing problems. I have also made the script below "more pedantic" to try and debug it. from qgis.core import * from qgis.gui import * import processing @qgsfunction(args='auto', group='Custom', referenced_columns=[]) def del_cutboxes(feature, parent): geom = feature.geometry() attr = feature.attributes() rd_at = str(attr[5]) lyr_rcl = QgsVectorLayer('Z:/Projects/RCL_boxes.shp','lyrrcl','ogr') rcl_feats = lyr_rcl.getFeatures() msg = '' for rcl_feat in rcl_feats: if len(msg) == 0: rclbox = rcl_feat.geometry() if rclbox.contains(geom): if rcl_feat.attributes()[0] is not rd_at: #feature.delete() msg = 'ToDelete ' + rd_at else: msg = 'ToKeep ' + rd_at if len(msg) == 0: msg = 'Nothing for ' + rd_at return msg -- = Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za <http://www.geograph.co.za> = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user -- = Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] When running a Custom Function, the return value is different from when stepping through it and watching the 'Preview' value in Field Calculator
Hi, Me again :-( I've written a custom function in Field Calculator and want to put its return value into a new Text field in the attribute table. When I step through the features I get the correct and expected values shown next to 'Preview' But when I click OK to run the Field Calculator, it runs for a long enough time, auto-closes on completion and in the attribute table the new field is created, but value is 'Nothing for None' for all features. *Note:* variable rd_at reflects correct contents when stepping through and observing 'Preview' but seems to revert to Nonetype when running the function. What am I missing regarding using custom functions? Is it an intended thing that the Field Calculator auto-closes after running a custom function? Thanks again, Kind regards, Zoltan *Note:* The msg variable was introduced in case placing 'return' in various lines in the script was causing problems. I have also made the script below "more pedantic" to try and debug it. from qgis.core import * from qgis.gui import * import processing @qgsfunction(args='auto', group='Custom', referenced_columns=[]) def del_cutboxes(feature, parent): geom = feature.geometry() attr = feature.attributes() rd_at = str(attr[5]) lyr_rcl = QgsVectorLayer('Z:/Projects/RCL_boxes.shp','lyrrcl','ogr') rcl_feats = lyr_rcl.getFeatures() msg = '' for rcl_feat in rcl_feats: if len(msg) == 0: rclbox = rcl_feat.geometry() if rclbox.contains(geom): if rcl_feat.attributes()[0] is not rd_at: #feature.delete() msg = 'ToDelete ' + rd_at else: msg = 'ToKeep ' + rd_at if len(msg) == 0: msg = 'Nothing for ' + rd_at return msg -- = Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] extracting vertices from a multipolygon using Field Calculator Function Editor
Thank you Stefan, I knew it would be as simple as that, but could not find it (the 'Java style' in the QgsGeometry class must have made me skip experimenting with that ) In Afrikaans (I am South African/Hungarian) there is a saying: 'As jy dom is, moet jy ...xxx' - and that was applicable to me in this case :-) I wrote the script yesterday in a freestanding (PyCharm) environment to solve the issue, as I needed to 'get it done', but this info is useful for future. Yesterday was the first time I was trying out the function editor and, just btw, a week or two back was when I first started writing QGis python scripts to solve work issues. So yes, I am struggling to find 'what is where' even though I am 10's of years down the road picturing the algorithm needed to solve my issue. My current challenge is also QGis Atlas with multiple sub-maps per page - but I am 'getting there' with some pointers received from this list. Thanks again for your response. Kind regards, Zoltan On 2021-10-19 08:14, Stefan Giese (WhereGroup) wrote: as a starting point: /from qgis.core import *// //from qgis.gui import *// // //@qgsfunction(args='auto', group='Custom', referenced_columns=[])// //def zCutBox(feature, parent):// // """// // Calculates vertices and return the x coordinate of the first as string// // """// // geom = feature.geometry()// // vertices = []// // for v in geom.vertices():// // vertices.append(v)// // return str(vertices[0].x())/ Am 18.10.2021 um 13:42 schrieb Zoltan: Hi Listers, I can't for the life of me work out the syntax/commands to use in the FieldCalculator->FunctionEditor for extracting vertices. I have a Polygon layer with Multipolygon rectangles at various rotations. I need to extract the 4 corner points so that I can calculate (and save in a virtual layer) lines that I can later use with "Processing->Split with lines". I will google for opening a new virtual layer and saving lines to it, but if that is quick and easy then some help there would also be welcome. Please can someone point me to a sample functions covering these requirements, or just fill in the gaps below. No, I'm not asking for someone to just write it for me. Thanks in advance, Zoltan [Having a brain-dead day :-( ] from qgis.core import * from qgis.gui import * @qgsfunction(args='auto', group='Custom', referenced_columns=[]) def zCutBox(feature, parent): vertices = blah blah return vertices[0][0] -- = Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info:https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user -- Mit freundlichen Grüßen Stefan Giese Projektleiter/Consultant - Schon gewusst? In unserem Blog geben wir Tipps & Tricks zu Open-Source-GIS-Software und berichten aus unserem Experten-Alltag: https://wheregroup.com/blog/ - WhereGroup GmbH Schwimmbadstr. 2 79100 Freiburg Germany Fon: +49 (0)761 / 519 102 - 61 Fax: +49 (0)761 / 519 102 - 11 stefan.gi...@wheregroup.com www.wheregroup.com Geschäftsführer: Olaf Knopp, Peter Stamm Amtsgericht Bonn, HRB 9885 -- = Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] extracting vertices from a multipolygon using Field Calculator Function Editor
Hi Listers, I can't for the life of me work out the syntax/commands to use in the FieldCalculator->FunctionEditor for extracting vertices. I have a Polygon layer with Multipolygon rectangles at various rotations. I need to extract the 4 corner points so that I can calculate (and save in a virtual layer) lines that I can later use with "Processing->Split with lines". I will google for opening a new virtual layer and saving lines to it, but if that is quick and easy then some help there would also be welcome. Please can someone point me to a sample functions covering these requirements, or just fill in the gaps below. No, I'm not asking for someone to just write it for me. Thanks in advance, Zoltan [Having a brain-dead day :-( ] from qgis.core import * from qgis.gui import * @qgsfunction(args='auto', group='Custom', referenced_columns=[]) def zCutBox(feature, parent): vertices = blah blah return vertices[0][0] -- ===== Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] QGIS Atlas gymnastics
Sounds promising. Thanks Alexandre and Roland Now where did I put my 'do not disturb' sign.? Regards, Zoltan On 2021-10-15 11:13, Alexandre Neto wrote: Ah! good to see that video put to good use :-) Yes, Zoltan, all you ask seems possible to me. The secret is in the way you prepare your coverage layer to iterate over the items. Check the video comment as I forgot to explain one step in the video. I also have this old blog post that covers dynamic size pages. The expressions are a bit outdated. that is, there are now functions that make the expressions easier: https://gisunchained.wordpress.com/tag/atlas/ Alexandre Neto QGIS Support www.qcooperative.net <http://www.qcooperative.net> On Fri, Oct 15, 2021 at 9:53 AM Roland Spielhofer wrote: Hi Zoltan, I guess your first item is covered here: https://www.youtube.com/watch?v=NCsnTt6uxXo I call the video "Atlas on steroids" ;-) Regards, Roland *Gesendet:* Freitag, 15. Oktober 2021 um 10:39 Uhr *Von:* "Zoltan" *An:* "qgis-user" *Betreff:* [Qgis-user] QGIS Atlas gymnastics Hi, I wonder if someone can confirm if below is possible or not, using the Atlas facilities QGIS 3.20.3 in WIndows 10. Two issues, the second relying on the first. *#1 Variable size Map Item.* I would like to vary the size of the Map Item window, depending on the zoom window of the Atlas page. I can set the "Position and Size" using an Expression or an attribute Field, but this seems to be fixed rather than dynamic ie: Is it possible to have Expression and/or Field values dynamic per Atlas coverage page? *#2 Multiple Sub-maps per MapSheet.* Is it possible to have an Atlas page layout such that (say) in the top half of the page you have your zoomed map (as in #1 above) and in the bottom half of that sheet-page you place n-mapboxes as insets predefined as n-AOI rectangles inside the zoom-window in the top of this MapPage? *A practical example/use would be:* Say you are doing an aerial analysis of a stretch of river and you have identified areas of interest by placing oriented rectangles around them. You now want to show an overview of segments of the river, and then a zoom-in to each identified AOI rectangle to show larger scale prints below this river segment. *Method could/would be to:* Using the width of your intended page size, calculate the length of river you can fit across 1 page of your Atlas. Let's assume 1Km Now a river meanders so you cannot run down the polyline depicting the river, you need to rotate and cut the river to the width of your atlas page. This means your atlas Map-Item width will be constant, but the height, due to curvatures in the river, will vary. Once you have your oriented MBRs calculated along the river, you can use these as your Atlas Coverage layer. To show this with optimal page-area usage, you would need the Atlas facility I describe in *#1 above*. Now, along the river you will have pre-placed your oriented AOI rectangles. These will be much smaller rectangles than the 1Km river-length Atlas Coverage rectangles. Now, for each AOI rectangle inside the current 1Km coverage rectangle, I want to show a larger scale image below the current coverage rectangle. If*#2 above* was possible and flexible, one could use the calculated (current) bottom left paper coord of the coverage rectangle to calculate the top-left position of the first AOI area. And so on across then down the remaining Atlas page for subsequent AOI boxes, which will likely be of varying size. I suppose this is rather like wanting an Atlas page-set within and Atlas page. Anyone achieved something like this? Does this sound useful? I guess, unless I've missed some facility in Atlas, I will be doing this by using the Atlas facility to spit out carefully named jpeg images, and then assembling the pages using Python and LibreOffice (or something). All thoughts/directions welcome. Kind regards, Zoltan -- = Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za <http://www.geograph.co.za> = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user ___ Q
[Qgis-user] QGIS Atlas gymnastics
Hi, I wonder if someone can confirm if below is possible or not, using the Atlas facilities QGIS 3.20.3 in WIndows 10. Two issues, the second relying on the first. *#1 Variable size Map Item.* I would like to vary the size of the Map Item window, depending on the zoom window of the Atlas page. I can set the "Position and Size" using an Expression or an attribute Field, but this seems to be fixed rather than dynamic ie: Is it possible to have Expression and/or Field values dynamic per Atlas coverage page? *#2 Multiple Sub-maps per MapSheet.* Is it possible to have an Atlas page layout such that (say) in the top half of the page you have your zoomed map (as in #1 above) and in the bottom half of that sheet-page you place n-mapboxes as insets predefined as n-AOI rectangles inside the zoom-window in the top of this MapPage? *A practical example/use would be:* Say you are doing an aerial analysis of a stretch of river and you have identified areas of interest by placing oriented rectangles around them. You now want to show an overview of segments of the river, and then a zoom-in to each identified AOI rectangle to show larger scale prints below this river segment. * **Method could/would be to:* Using the width of your intended page size, calculate the length of river you can fit across 1 page of your Atlas. Let's assume 1Km Now a river meanders so you cannot run down the polyline depicting the river, you need to rotate and cut the river to the width of your atlas page. This means your atlas Map-Item width will be constant, but the height, due to curvatures in the river, will vary. Once you have your oriented MBRs calculated along the river, you can use these as your Atlas Coverage layer. To show this with optimal page-area usage, you would need the Atlas facility I describe in *#1 above*. Now, along the river you will have pre-placed your oriented AOI rectangles. These will be much smaller rectangles than the 1Km river-length Atlas Coverage rectangles. Now, for each AOI rectangle inside the current 1Km coverage rectangle, I want to show a larger scale image below the current coverage rectangle. If*#2 above* was possible and flexible, one could use the calculated (current) bottom left paper coord of the coverage rectangle to calculate the top-left position of the first AOI area. And so on across then down the remaining Atlas page for subsequent AOI boxes, which will likely be of varying size. I suppose this is rather like wanting an Atlas page-set within and Atlas page. Anyone achieved something like this? Does this sound useful? I guess, unless I've missed some facility in Atlas, I will be doing this by using the Atlas facility to spit out carefully named jpeg images, and then assembling the pages using Python and LibreOffice (or something). All thoughts/directions welcome. Kind regards, Zoltan -- ===== Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] Results from QgsGeometryUtils.lineCircleIntersection
Hi, I am struggling with using QgsGeometryUtils.lineCircleIntersection in a python script I am writing outside of QGIS 3.20.3 [495fbaec], using PyCharm. In https://qgis.org/pyqgis/3.2/core/Geometry/QgsGeometryUtils.html#qgis.core.QgsGeometryUtils.lineCircleIntersection it says that: If the intersection has two solutions (points), the closest point to the initial intersection point is returned. I am wanting to find the coordinate on the polyline, that is at a set distance from the polyline startpoint. Note that this polyline meanders so I cannot do this using point on line at distance 'd' from startpoint Unfortunately this utility extends my vertex to vertex line segments and I get results where the intersection point is not actually ON my line-segment. This is easy to check for, but the problem is that this utility only returns the first intersection point found, which is not always the one actually sitting on my line-segment. Am I approaching this problem incorrectly, or is this a known 'feature' of this utility? Any guidance welcome. Thanks and regards, Zoltan -- ===== Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] i have been digitizing polygons for a few days and QGIS crashed
What's a floppy disk? ;-) On 2021-10-08 22:55, Garth Fletcher wrote: Alexandre Neto wrote: Hi, You are missing one of the most important aspects of QGIS (and most GIS software): The data is independent of the project. The project > save action saves the project, which include what layers you have loaded, what order you gave them, their styles, and so on. When you edit vector layers, you need to use the save icon in the digitising toolbar, the one next to pencil icon that you have probably clicked to start editing. This is done layer by layer. The only time I have seen QGIS prompt the user to save an edited but not yet saved layer is when you Quit. At that moment it does seem to remind one about any unsaved layers. However, if it crashes then any unsaved new layers, or unsaved changes to pre-existing layers, are lost. The crucial step,as Alexandre says, is to save the changes in the layer by *clicking on the floppy diskette icon* next to the pencil... Alternatively, toggling the Pencil to turn of editing will also prompt you to save changes. -- ===== Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] QGIS Keybindings
Presumably you have looked at Settings->Keyboard Shortcuts ? The trick would be to extend the available list on that sub-menu to include items that are not listed there. (Let me know when you find out how) Regards, Zoltan On 2021-10-04 17:45, Hugh Kelley wrote: Hi all, I periodically need to make extensive edits to vector datasets in QGIS. It would be very helpful to have keybindings that allowed me to switch between the drag map extent, select, edit vertices, draw feature and maybe a few other cursor types. Web searches don't turn much up. I did find this: https://defkey.com/qgis-3-shortcuts <https://defkey.com/qgis-3-shortcuts>, but that doesn't seem to have what I need. thanks for any suggestions about an alternative workflow or custom keybindings or what not. -- Hugh Kelley ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user -- ===== Zoltan Szecsei GPrGISc 0031 Geograph (Pty) Ltd. GIS and Photogrammetric Services Cape Town, South Africa. Mobile: +27-83-6004028 (Signal, not WhatsApp) www.geograph.co.za = ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Vector Affine trouble
Hi Paolo, the QGIS 3.4 upgrade is a pull request what the original author of qgsAffine hasn't merged yet (2 monthes old). If he is on this list ... Regards, Zoltan On Wed, 6 Feb 2019, Paolo Cavallini wrote: Hi Suki, Would you mind asking authors why they didn't upload to the official repo? I think it would be good to have. All the best. Il 6 febbraio 2019 17:29:15 CET, Siki Zoltan ha scritto: Hi, there is a qgsAffine version for 3.4 onGitHub: https://github.com/pohmelie/qgsAffine Unfortunatelly it is not available from QGIS 3 Plugins. You can install it manually, download zip from github and copy the qgisAffine folder from zip to the QGIS plugin directory on your machine. To find the plugin folder read: https://gis.stackexchange.com/questions/274311/qgis-3-plugin-folder-location Regards, Zoltan On Wed, 6 Feb 2019, Olivier Dalang wrote: Hi, Vectorbender has just been ported to qgis 3 by Jorge Almerio (many thanks again). I just released the new version to the plugin repo a few seconds ago. Cheers, Olivier On Sat, Jan 12, 2019, 23:37 Rutger Emmelkamp wrote: Hello everybody, Hope someone can help me out this Saturday night (or Sunday or Monday) but no hurry! How to georeference an Illustrator vector image (contour line map) from a specific area? I have qgis 3.4.3 installed and qgisAffine and Vector Blender plugins don˙˙t work for this version. I managed to turn the file into a .shp file and I see it somewhere off the coast of Ghana (coordinates 0,0) but how do I drag it into place? I think I˙˙m looking for a Vector Affine Transformation plugin or menu but can˙˙t find it. Hope someone can help me in the right direction. Thanks and have a nice weekend. Rutger ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user -- Sorry for being short___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Vector Affine trouble
Hi, there is a qgsAffine version for 3.4 onGitHub: https://github.com/pohmelie/qgsAffine Unfortunatelly it is not available from QGIS 3 Plugins. You can install it manually, download zip from github and copy the qgisAffine folder from zip to the QGIS plugin directory on your machine. To find the plugin folder read: https://gis.stackexchange.com/questions/274311/qgis-3-plugin-folder-location Regards, Zoltan On Wed, 6 Feb 2019, Olivier Dalang wrote: Hi, Vectorbender has just been ported to qgis 3 by Jorge Almerio (many thanks again). I just released the new version to the plugin repo a few seconds ago. Cheers, Olivier On Sat, Jan 12, 2019, 23:37 Rutger Emmelkamp Hello everybody, Hope someone can help me out this Saturday night (or Sunday or Monday) but no hurry! How to georeference an Illustrator vector image (contour line map) from a specific area? I have qgis 3.4.3 installed and qgisAffine and Vector Blender plugins don˙˙t work for this version. I managed to turn the file into a .shp file and I see it somewhere off the coast of Ghana (coordinates 0,0) but how do I drag it into place? I think I˙˙m looking for a Vector Affine Transformation plugin or menu but can˙˙t find it. Hope someone can help me in the right direction. Thanks and have a nice weekend. Rutger ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] changing polygon to point
Dear Boaz, in QGIS 2.x you can use the realcentroid plugin. In QGIS 3.x you can use the point on surface from the processing toolbox->vector geometry->point on surface Zoltan On Sun, 16 Sep 2018, Boaz Bar Ilan wrote: hi I have a layer of polygons and I want to convert every polygon to be a point type is there a way to do that ? thanks boaz ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] NTV2 Grid Shits
Dear Tyler! proj4 projections definitions are stored on ubuntu in the file /usr/share/proj/epsg, edit this file to add +geoidgrids=your_grid. QGIS projection definitions are stored in a sqlite database, on ubuntu you can find the srs.db in the /usr/share/qgis/resources directory. You can update the record of a projection in this database. I have sofar used +nadgrids for horizontal coordinates, I don't know if QGIS allow +geogrids, too. Try to update parameters column in tbl_srs table: UPDATE tbl_srs SET parameters=parameters || ' +geoidgrids=your_gtx_file' where srid=your_epsg_code; Be careful, if you update your QGIS version you probably have to update also the srs.db. Best regards, Zoltan On Tue, 24 Oct 2017, Tyler Veinot wrote: Since you guys both know more about PROJ4 than I do I thought I would just tag this question on; Where is QGIS/PROJ4 with support for vertical datums? Is there a way I can attach a vertical datum to a CRS in PROJ4? Tyler On Mon, Oct 23, 2017 at 9:09 AM, Tyler Veinot wrote: Andre; It seems QGIS does not like my gsb file, I downloaded the NTV2 gsb from NRCan and pointed to it and it worked. For some reason the one specific for PEI does not. However, the grid shift could be the same as far as I know; cannot find any information on what mathematical principles are used to define the shift to correct the coordinates. Tyler On Mon, Oct 23, 2017 at 8:01 AM, André Joost wrote: Am 20.10.2017 um 20:19 schrieb Tyler Veinot: Andre; Again, sorry accidently sent the email mid writing. Where was I Right so I added nadgrids=PE7783V2.gsb a number of ways trying to sort it based on online information and how some other CRS's were setup within QGIS but nothing worked yet. It might be better to post the full string. When I add the directory to my gsb where do I add it in the line and does it need quotes? If you have blanks in folder names, surely. HTH, André Joost On Fri, Oct 20, 2017 at 2:55 PM, Andre Joost wrote: Am 20.10.2017 um 14:31 schrieb Tyler Veinot: Anyone knowledge in using grid shift files (gsb) with QGIS? I have a local CRS NAD27 Double Stereographic Imperial that I need to transform to NAD83 CSRS 2010 Epoch. I have the CRS properly defined using the provincial description Can you add your custom CRS? Maybe you have +datum=nad27 which uses the main nadcon transformation file. use +ellps=clrk66 +nadgrids=/path/to/your/file but QGIS using its Clark 1866 to GRS80 transformation puts the points a little off from the NTV2 coordinate transformation software we have. I have the GSB used (even though it says ATS77 the provincial geomatics guys tell me it shits NAD27 and ATS77 into CSRS 2010 NAD83). I followed some directions on putting the file in a specific folder and restarted QGIS set the settings to always ask for grid shift file for transformations yet I do not see the one I put in there in my list to pick from. Can you explain where that setting should be? I also tried to install a plugin for NTV2 transformations but cannot activate it. When I go to Activate Additional Providers in my Processing Toolbox where NTv2 Datum Transformations resides I get the "Specified Path does not exist: C:\PROGA~1\QGIS~1.18\bin" error that I have not yet managed to resolve. To my understanding has something to do with GRASS plugins and past versions of QGIS. You can workaround that GRASS error by creating that folder, or changing the folder entry for GRASS 7 to an existing folder. But to use your grid in the plugin, you would need some python code around it, like the other grids have in C:\Users\\.qgis2\pyt hon\plugins\ntv2_transformations Without that processing would not know about your grid. HTH, André Joost ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] ogr2ogr - how to use -s_srs, -t_srs and -a_srs
Dear Erik, if there is a .prj file to your shape and it is valid you need not to use -s_srs. Only one of -a_srs and -t_srs is neccessary in your case, if you want to reproject coordinates use -t_srs. If you want only change the projection but not the coordinates use -a_srs. Best regards, Zoltan On Sun, 6 Aug 2017, Erik Josefsson wrote: Hello, I have a handful of zip archives containing original shapefiles with different coordinate reference systems (e.g. SGER:3008 and SGER:3021). I want to extract and convert them all to SGER:3006 and thought that ogr2ogr would do the job, but I don't find intelligible instructions on how to use the srs-parameters: -s_srs, -t_srs and -a_srs. Do I have to know the coordinate reference system of each original shape file to use ogr2ogr successfully? If so, how do I retrieve that information? Or should I just override the original coordinate reference system like this? ogr2ogr -f 'ESRI Shapefile' -s_srs EPSG:3006 -t_srs EPSG:3006 -a_srs EPSG:3006 A.shp B.shp Best regards. //Erik ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] DXF export.
Dear Charles, Rasters cannot be involved in DXF export (as Jürgen wrote). You had better to use QuickMapServices or Openlayers plug-in, no need to manual georeferencing. These plug-ins are able to display google, osm, ... maps from the web. Best regards, Zoltan On Wed, 12 Jul 2017, Charles Barrows wrote: I'm rather new to QGIS but very impressed with its capabilities. I'm trying to georeference a raster photo from Google Earth. Using the GeoReferencer plugin I go through the procedure using data points in a decimal degree format. WGS84 survey. The procedure seems to go OK, including the "start georeferencer" and the Load to QGIS. I am using the Transformational type Polynomial 2. But, the DXF export fails. What am I missing? I have been entering the longitude into the x/box as a negative number because I am in the western hemisphere. Does the program not like negative numbers? Charley B. ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Creating Systematic sampling plot design
Dear David, you should turn on Processing plug-in. Regards, Zoltan On Thu, 23 Mar 2017, David Kabay wrote: Hi I have the QGIS manual for 2.8 but are using 2.16 version. I want to do the above. The manual says go to vector and then to Research Tools. I go to vector (along the top) of my 2.16 but there is no ""research tool" tab. Where is it. I looked in other places but cannot find it. Is 2.16 different?? If not where can I find it or something similar. Thanks Regards Dave Dave Kabay KABAY rehabilitation, environmental and biological CONSULTANTS Pty Ltd 19A Priestley Street Embleton WA 6062 Phone 08 6142 6756 Mob 0417 950 508 Work Power bb <mailto:dka...@iinet.net.au> dka...@iinet.net.au Je Suis cid:image001.png@01D00E62.B86A9140 Charlie https://au.linkedin.com/in/dave-kabay-94aa7292 http://www.the-numbers.com/video/Despicable-Me-2/Despicable-Me-2-poster.jpg ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] inquiry
Dear Ahmed, you can share your shape files (.shp, .shx, .dbf) simply giving the files to your collegues. To share a QGIS project file with ArcGIS users is not available. I have heard about QGIS plug-in to read an ArcGIS mxd. http://gis.stackexchange.com/questions/23717/converting-arcgis-mxd-file-to-qgis-qgs-file Best regards, Zoltan On Sat, 12 Nov 2016, Ahmed Kheir wrote: Dear all How can I save my qgis project or my shp.file layers to be opened by ARC-GIS? on other meaning I have prepared different layers by qgis and I want to share these layers to my friend (ARC-GIS user), because he can not use qgis but he is using ARC-GIS? Regards Ahmed Egypt ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] csv/csvt import
Dear Evan, thanks, I could get the same result using org2ogr. I would like to get point layer at the same time not "only" an attribute table. Regards, Zoltan On Sat, 29 Oct 2016, Even Rouault wrote: Le samedi 29 octobre 2016 10:18:57, Siki Zoltan a écrit : Hi, I tried to use csv/csvt files to add Delimited Text Layer. The widths of the fields are not considered by QGIS when saving into other formats. In the layer properties dialog the width and precision of fields are zeros after loading csv/csvt. I suppose the width and precision of fields should be preserved after reading of csv/csvt files. I tested it in QGIS 2.2/2.8/2.14 but I got maximal width for fields (depending on the field type, character 254, float 23.15) when I save the file as a shapefile. My csvt file: "String(20)","Real(11.2)","Real(11.2)","Real(7.2)","String(20)","String(20) " My csv file: point_id,e,n,z,pc,pt 2002,650769.92,237368.68,104.560,, 2004,650730.28,237418.54,104.520,, 2005,650706.42,237412.43,104.270,, 2006,650707.76,237449.22,104.590,, Is it a feature or bug? Any workaround? If you open through the OGR provider (Add vector layer), this will work. Note that the OGR CSV driver will warn (if you do ogrinfo for example) that the values in the z field do not respect the Real(7.2) formatting Even -- Spatialys - Geospatial professional services http://www.spatialys.com ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] csv/csvt import
Hi, I tried to use csv/csvt files to add Delimited Text Layer. The widths of the fields are not considered by QGIS when saving into other formats. In the layer properties dialog the width and precision of fields are zeros after loading csv/csvt. I suppose the width and precision of fields should be preserved after reading of csv/csvt files. I tested it in QGIS 2.2/2.8/2.14 but I got maximal width for fields (depending on the field type, character 254, float 23.15) when I save the file as a shapefile. My csvt file: "String(20)","Real(11.2)","Real(11.2)","Real(7.2)","String(20)","String(20)" My csv file: point_id,e,n,z,pc,pt 2002,650769.92,237368.68,104.560,, 2004,650730.28,237418.54,104.520,, 2005,650706.42,237412.43,104.270,, 2006,650707.76,237449.22,104.590,, Is it a feature or bug? Any workaround? Best regards, Zoltan ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Plugin Building
Dear Tyler, As far as I remember you have to generate the python file from the Qt ui file. It is done automaticly by the make file. Did you "make" your plugin? Regards, Zoltan On Fri, 22 Apr 2016, Tyler Veinot wrote: Hello All; I am using QGIS 2.14.1 and am trying to build a python plugin using the plugin builder using this tutorial; https://anitagraser.com/2014/04/26/getting-started-writing-qgis-2-x-plugins/ and the information here; http://docs.qgis.org/testing/en/docs/pyqgis_developer_cookbook/plugins.html#init-py But when I complete the first few steps and get to installing the plugin to see a blank interface I get error when calling it's classfactory method in the utils.py and __init__.py and I am not sure what is causing this error. I think I must be missing some required modules or something. These scripts are generated by the plugin builder and are suppose to call up the basic qgis and qt functions to my understanding. Has anyone ran into this before? Thanks Tyler ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] New Custom Projection
Dear Grant, instead of e (east) and s (south) use sign (+/-) e.g. +proj=aea +long_0=121 +lat_1=-17.5 +lat_2=-31.3 I suppose you should give +datum +ellps +units, too (I don't know if there are default and they are good for you) Regards, Zoltan On Mon, 11 Jan 2016, Grant Boxer wrote: Hi All, I have a jpeg image of a map in Albers Equal Area Projection (central meridian 121 deg E, St Parallels of 1.5 deg S and 31.5 deg S, Western Australia) and I need to bring it into QGIS. I have tried making a custom projection for this but have failed. The projection info I have created is "+proj=aea +long_0=121e +lat_1=17.5s +lat_2=31.3s". Any suggestions on where I have gone wrong? Regards Grant ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Problem managing data
Dear Eva, AFAIK there is no direct solution to edit ods/xls/xlsx files from QGIS. There are two possible scenarios: 1. Edit your speadseet data in Libre|Open|MS office and open it in QGIS and join it to your spatial data. If your spreadsheet is open in both QGIS and Excel it can cause trouble. 2. Open you spreadsheet in QGIS and save as a QGIS editable format (e.g. dbf, sqlite) and edit these data. If you need again these data in a spreadsheet convert it back (probably through csv format). Hope it helps you, Zoltan On Mon, 28 Dec 2015, Eva Mosquera wrote: Hello everybody!I need to use the data I have in my shpfile in qgis in excel. I need an excel spreadsheet that is able by itself to keep updated with the data changes in qgis. And I don't know how to do this. I'm trying to import the database (dbf) in excel but I have problems with the type and characteristics... Thank you very much ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Obtaining cs2cs command line tool
Dear Philip probably you have installed QGIS standalone automaticly into the "Program files" folder. Some of the utilities installed with QGI doesn't like special characters in the path. Solutions: 1. Use the OSGeo4W installer, you can use/install more utilities/software 2. Reinstall QGIS into a path whith no special character Regards, Zoltan On Wed, 23 Dec 2015, Philip Condit wrote: I am using the QGIS 2.2 version (Valmiera) on a Windows 8.1 64-bit laptop. When I tried to open the MSYS (Minimum System) icon, I got a message that it could not open because the C:\Program directory could not be found. I then created a C:\Program folder, and copied all of the files from the C:\Program Files\QGIS Valmiera\apps\msys folder to the C:\Program folder. I then was able to open MSYS, and get a terminal window to show up on my desktop. However, when I typed cs2cs into that window, I got a message that the command could not be found. How do I obtain the cs2cs command line tool? ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Measuring area of a selected row in attribute table
Dear Manoj, you can create an expression for a virtual field in the table view This way all the areas will be filled, not only the selected. See attached picture. You can get the area of a polygon in the actual layer using the "i" (identify) tool, too. In the info window look for Derived group. Regards, Zoltan On Mon, 21 Dec 2015, Manoj Soni wrote: Dear Users I am very new to QGIS. We can draw area and can measure area using measuring tools. According to attribute table, i just want to measure area of a selected row in table. How can we measure it. For example, i have a shape file having some cities name in attribute table, when i selected one row, that city area is highlighted. How can i measure that selected area only?. I do not want to make a selection manually. Many Thanks Regards ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] GPS TO QGIS
Dear Brent, I've never used gpsd on windows. There are some information on this page: https://code.google.com/p/gpsd-4-win/ Regards, Zoltan On Sun, 20 Dec 2015, Brent Wood wrote: The mobile unit of an SP80 doesn't output NMEA, just the base station, & I don't think it is supported by gpsd. The base station should be, if you switch the output to NMEA format. Does QGIS work with GPSD on Windows? Brent Wood From: Siki Zoltan To: g.corral Cc: "qgis-user@lists.osgeo.org" Sent: Monday, December 21, 2015 8:20 AM Subject: Re: [Qgis-user] GPS TO QGIS Hi, if you can connect your GPS to gpsd then you can connect QGIS to it. See QGIS Live GPS Tracking. Another chance can be to use NMEA output of your GPS... Regards, Zoltan On Fri, 18 Dec 2015, g.corral wrote: Hi all!! I want to connect a precision GPS (GNSS receiver, SP80) to qgis. We know connect qgis to integrated GPS tablet but we can't connect it to a GPS of major precision. Has someone tried this connection?? Cordialement, Guayente Corral Broto Hydrogéologue / EurGeol n°1336 "L'eau fait partie du patrimoine commun de tous. Sa protection, sa mise en valeur et le développement de cette ressource durable, dans le respect des équilibres naturels, sont d'intér˙˙t universel". Chemin d'Enrobert 32200 Gimont - Tél: +33 (0)5 62 65 67 65 - Mob: +33 (0)6 48 74 43 67 - Mail: g.cor...@topodoc.com http://www.survey-groupe.fr ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] GPS TO QGIS
Hi, if you can connect your GPS to gpsd then you can connect QGIS to it. See QGIS Live GPS Tracking. Another chance can be to use NMEA output of your GPS... Regards, Zoltan On Fri, 18 Dec 2015, g.corral wrote: Hi all!! I want to connect a precision GPS (GNSS receiver, SP80) to qgis. We know connect qgis to integrated GPS tablet but we can't connect it to a GPS of major precision. Has someone tried this connection?? Cordialement, Guayente Corral Broto Hydrogéologue / EurGeol n°1336 "L'eau fait partie du patrimoine commun de tous. Sa protection, sa mise en valeur et le développement de cette ressource durable, dans le respect des équilibres naturels, sont d'intér˙˙t universel". Chemin d'Enrobert 32200 Gimont - Tél: +33 (0)5 62 65 67 65 - Mob: +33 (0)6 48 74 43 67 - Mail: g.cor...@topodoc.com http://www.survey-groupe.fr ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] extracting elevations from DEM at certain points
Dear Szilard, there is a point sampling plugin for QGIS: https://plugins.qgis.org/plugins/pointsamplingtool/ Regards, Zoltan On Wed, 16 Dec 2015, Szilard Albert wrote: Hello, I need to extract elevation values from a DEM at certain points for which I have X and Y values. My points are contained in several profiles, each profile with hundreds of points. My input is a text file, one file per profile, with the following format: PointA X Y PointB X Y .. PointZ X Y. I would like a solution where I can feed the above text input, and get my output as another text file with similar format, Point X Y Z (or Point Z X Y). Alternatively just a text file with one column containing Z, but exactly same number of rows, so I can copy and paste to the input. One single command that solves one point only would do as well, I would just create a script file from my every text file for every line Any suggestions would be greatly appreciated. I am using QGIS 2.8 Wien in Ubuntu. Thank you. Szilard ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Photo2Shape plugin
Hi Nuno, exifread is an external python package what you must install outside QGIS. See: https://pypi.python.org/pypi/ExifRead or http://gis.stackexchange.com/questions/132697/getting-no-module-name-exifread-error-from-qgis-photo2shape-plugin Regards, Zoltan On Tue, 8 Dec 2015, Nuno de Santos Loureiro wrote: Hello all, I am trying to install the plugin Photo2Shape on QGIS 2.12.1-Lyon running on a Mac with OS X Yosemite 10.10.5. However I can't load the plugin and I receive the following message: The plugin is broken. Python said: No module named exifread Can somebody send me some tips to fix the problem? Many thanks in advance, Nuno -x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x- Nuno de Santos Loureiro (PhD) Faculdade de Cincias e Tecnologia - DCTMA - Ed. 8 Universidade do Algarve Campus de Gambelas 8005-139 FARO (Portugal) -x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x-x- This message was sent using IMP, the Internet Messaging Program. ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: http://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Action "open" to open a pdf-File
Hello Peter, this may help you: http://gis.stackexchange.com/questions/29464/how-to-add-links-to-pdfs-in-a-shapefile Zoltan On Mon, 14 Sep 2015, Peter Schürch wrote: Hi Sylvain Thanks for your hint. Yes this is one possibility - but a very poor one because: - If several users on different machines use this project file to access the data it only works if the acrobat installation path is the same everywhere... This is very unlikely... Also I thought (and read somewhere) that the action type "open" does not require to specify the program with which a file should be opened. Am I wrong?? Peter Von: PIERRE Sylvain [mailto:sylvain.pie...@bas-rhin.fr] Gesendet: Montag, 14. September 2015 13:28 An: Peter Schürch Betreff: RE: Action "open" to open a pdf-File Hi, You just have to give the path to acrobat reader first (depends on your OS): "C:\Program Files\Adobe\Reader 11.0\Reader\AcroRd32.exe " [% "PDF File" %] Sylvain De : qgis-user-boun...@lists.osgeo.org<mailto:qgis-user-boun...@lists.osgeo.org> [mailto:qgis-user-boun...@lists.osgeo.org] De la part de Peter Schürch Envoyé : lundi 14 septembre 2015 10:03 ? : Qgis Mailinglist (qgis-user@lists.osgeo.org<mailto:qgis-user@lists.osgeo.org>) Objet : [Qgis-user] Action "open" to open a pdf-File QGIS 2.8.1-Wien / Windows 8 / MSSQL-Spatial Hi List I want to open a PDF-File from the QGIS-Object-Attribute Window (which pops up when I click with the "info-Tool" on a e.g. point-geometry). What I have done so far: - I have a geometry-table in the MSSQL spatial database with point geometries which works fine (displays in map, is editable etc..) - I have an attribute which contains an absolute path to a PDF-file on our server (The path is correct and works in Explorer). - I have defined an action of type "Open", given it a name, added the action [% "PDF File"%] referring to the attribute with the file path and added it to the action list. When I open the Object Attribute Window and right click somewhere I can select my action BUT nothing happens, not even an error message... I'm clueless. Using the Action Tool from the toolbar - same result, nothing happens, no error message. Any ideas? Thanks, best regards Peter Schürch -- Meier und Partner AG Freiestrasse 26 8570 Weinfelden www.meierpartner.ch<http://www.meierpartner.ch/> T +41 71 626 5115 ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] qgis export to excel
Dear Matthias, in the popup manu of the layer (right click on the layer name in Layers list) select Save as, select Coma separated value (CSV) or DBF file for the format. Both CSV and DBF can be opened by Excel (or Libre Office Calc). Regards, Zoltan On Sat, 12 Sep 2015, i...@architekten.coop wrote: dear qgis users i use qgis as an architect, for roominformations of big building. is there a possibility to export queryresultats to excel ? thank you for any helb matthias moser ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] DXF Export issue
Dear Russel, you cannot export raster to dxf, it is available only for vector layers. If your cad software can handle geoferenced raszter (e.g. AutoCAD Map), you needn't export it only open. Check your target software what raster formats are supported. Regards, Zoltan On Mon, 17 Aug 2015, Russell Fincher wrote: Greetings! First post here. I'm having trouble successfully exporting to DXF. These are my steps: 1. In QGIS, choose Layer>Add Layer>Add Raster Layer 2. I select an .lbl file next to an associated .img file that I got from the Mola Mars data here <http://pds-geosciences.wustl.edu/mgs/mgs-m-mola-5-megdr-l3-v1/mgsl_300x/meg128/> . 3. A height map successfully appears in the viewer window. 4. I choose Project>DXF Export and leave the settings at default. 5. I get an exported DXF, but it is basically empty. 16 kb and only 2700 lines of what mostly looks like header data. I'm unable to import it into 3ds Max, getting an error on import about an incorrect sequence. These three files can be downloaded here <https://drive.google.com/file/d/0B65j6EIp4zMQVzExQy1sWFVsOWs/view?usp=sharing> if anyone wants to look at them. Any ideas what I'm doing wrong? - Russell ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] [Qgis-developer] DXF/DWG import requirements document
Dear Andreas, beside "Import of layer names as attributes" I would prefer 1. setting layer filter in the import dialog (e.g. as a regexp or a selection list from the available CAD layers), unselected layers are skipped in input. It is not clear for me if the layer name rules can be used for this purpose. 2. import each CAD layer into separate QGIS layers by feature type (point, annotation, polyline, polygon), it could be a checkbox in the import dialog. Optionally during the import polygon layers sould be cheched for islands (polygons inside other polygon) and remove the area of islands from the container polygon. Regards, Zoltan On Tue, 23 Jun 2015, Andreas Neumann wrote: Hi, To all who have an interest in contributing to a DXF/DWG import function in QGIS: I have started working on a requirements document - you can view it here: https://docs.google.com/a/qgis.org/document/d/1NJv1lwd5TnWdV-fTuW7g5ts1O8r3vb46xYeAIQ_BwGw/edit?usp=sharing If you want to contribute to the document, please send me your gmail account so I can give you write access. After a couple of days (maybe 1.5-2 weeks) we can hopefully finish this document and ask for a quote. I will then ask again for organizations or companies who can contribute financially. Thank you in advance for your support! Andreas ___ Qgis-developer mailing list qgis-develo...@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-developer ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Shapefile at wrong coordinates
Dear Barbara, it is probably a projection issue. Check the SRS for both layers on both machine and the SRS for the canvas of QGIS. Are they the same on both machines? If they match, check the definition of the used SRS in the srs.db on both machines. Are they the same? You can find srs.db in the resource subdirectory of your qgis installation. srs.db is a sqlite database. Regards, Zoltan On Fri, 22 May 2015, Estevan Barbará wrote: Hello guys At our lab, we've tried to compare a shapefile of island contours with an geotiff of a nautical chart. In some of the working stations, the positions were perfectly matched. But on one station, the shapes are dislocated nearly 6 miles. We've checked the positions in other stations and softwares, and the error is only at qgis at this machine. How can we solve the problem? Regards Estevan Barbara ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] text
Dear Charel, please always cc your mail to the qgis user list, other users can learn from this, too. If you would like to assign an event to a field in the attribute table you must use the Action tab of the layer property dialog. There are some examples in the user manual (Action Menu section). Regards, Zoltan On Thu, 21 May 2015, Charel gmail wrote: Good evening Siki, It was a stupid failor of me.I thought that the maximum of the text is 255 characters. But only the length is limited. Maybe I can ask you another question... I make a field for a link to a website. I thought that when you click the browser is opening.. But nothing is happening. Greetings charel -Oorspronkelijk bericht- Van: Siki Zoltan [mailto:s...@agt.bme.hu] Verzonden: donderdag 21 mei 2015 20:31 Aan: Charel gmail CC: qgis-user@lists.osgeo.org Onderwerp: Re: [Qgis-user] text Dear Charel, ESRI shape files uses dbf (dBase) format, which limits the width of text field to 255 chars. If you need longer text field (not usual in GIS), you must use another editable format of QGIS, e.g. spatialite or postgis. regards, Zoltan On Thu, 21 May 2015, Charel gmail wrote: Hello, This week I start with Qgis. The most things I want, I found it on internet... In the table for objects there is the possibility for a text field. But it is limited to 255 character. I can't find a possibility for more text such as a memo field. Is there a possibility to use more text ? (sorry my English is quit bad) Greetings charel ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] text
Dear Charel, ESRI shape files uses dbf (dBase) format, which limits the width of text field to 255 chars. If you need longer text field (not usual in GIS), you must use another editable format of QGIS, e.g. spatialite or postgis. regards, Zoltan On Thu, 21 May 2015, Charel gmail wrote: Hello, This week I start with Qgis. The most things I want, I found it on internet... In the table for objects there is the possibility for a text field. But it is limited to 255 character. I can't find a possibility for more text such as a memo field. Is there a possibility to use more text ? (sorry my English is quit bad) Greetings charel ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] QGIS Android, GPS
Deae All, I've installed QGIS Android on a WayteQ xTAB-7X tablet. I would like to use the position from the on board GPS. How can I solve it? Thanks, Zoltan ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] how to calculate the overlap area of several polygons of the same shape qgis 2.4
Dear Francesc, Create a unique field in the attribute table of your polygon layer, for example an ordinal number of the features ($rownum in field calculator) Add the same polygon layer twice to your QGIS project. Rename one of them not to mix them up. Use the Vector/Geoprcessing/Union (or intersection) between the two layers (refering to the same shape file). In the result you will get the columns of the two layers, names are changed to be unique (e.g. if you created an id column, then you get id and id_2). You can ignore/delete those polygons where id = id_2. Add an area column to the result using the field calculator. I hope it helps you. Best regards, Zoltan On Sat, 11 Apr 2015, Francesc Caas wrote: Hi, I'm using qGis 2.4 and I want to know the overlap area of many polygons of the same shape. How can I do it? There any tool to do it or I have to do some steps? There are too many polygons to create separate shapes for everyone and make clips. Thanks a lot. -- *Francesc Ca˙˙as <http://francescsig.info/>* Ambientoleg SIG <http://francescSIG.info> <http://aixonoesunblog.blogspot.com> ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] solar potential
Hi Szilard, did you try the nativ GRASS command without QGIS toolbox? Zoltan On Tue, 3 Mar 2015, Szabo Szilard wrote: Hi! Does anyone has some experience with the r.sun GRASS option (Processing toolbar). I am strugling with the results because they are not reasonable. Data can be calculated for a day and in January the solar income is greater than in July ... I have tried with a national CRS (EPSG 23700) and with WGS84 (4326) but all output seems a bit strange (too small and the issue of January-July). Do you have idea what am I doing wrong? DSM: 20 cm resolution objects: a building with roof planes derived from photogrammetric technique QGIS 2.6.1 Win 8.1 DSM is here (EPSG 4326): https://share.tresorit.com/?0024vtkidgwfi9ib4u02i883&GNNzY6Ig7hqNHeIFeYPP69CIIDynxO9L&1&dXVQoIiKa3juWxNm8I4UnbvuUffB6voJuoBqxyE9T2GgXjtgfGsmV8FWLKpVPkCl#VQr6AtH1cY2D1VqWMDnokthlHd1AQ5I4sLWQ7+EUlcs=&ZGVtX25lazN3Z3MudGlm Thanks in advance -- Dr. Szilard SZABO associate professor, head of department Department of Physical Geography and Geoinformation Systems University of Debrecen https://www.researchgate.net/profile/Szilard_Szabo/?ev=hdr_xprf Publications <https://vm.mtmt.hu/search/slist.php?lang=0&AuthorID=1172> ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] What languages are the QGIS Manuals availabe in
Dear Thomas, go to the documentation page of qgis (e.g. docs.qgis.org/en/docs/user_manual/) and select another language from the list at the upper right corner. Regards, Zoltan On Wed, 11 Feb 2015, Leahy, Thomas A. wrote: Hello, What languages are the QGIS Manuals available in? Where can I find a list and access to the manuals in other languages. Thank you, Tom ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Cutting lines with a polygon boundary
On 2015/01/25 20:23, Andre Joost wrote: Am 25.01.2015 um 13:51 schrieb Zoltan Szecsei: So, does anyone know of a way to efficiently cut features in multiple layers using a polygon, and then to be able to select the pieces (multiple layers) that are either inside or outside the polygon, and then delete them? Disregarding the "multiple": "Vector -> Geoprocessing -> Clip" deletes the items outside the polygon Nope, doesn't. - In my first point I said: */Vector Geometry Processing Clip, but it wants to extract the data inside the polygon and create a new file with the extracted data. (To be expected) /*but meant: */Vector-> GeoProcessingTools-> Clip, but. /*(my bad, apologies)*/ /*It does not delete, it clips them out into a new layer/- which is not what I'm hoping to achieve./*/ /* Thanks for the thought anyway, though. Regards, Zoltan*/ /* "Vector -> Geoprocessing -> Difference" deletes the items inside the polygon You would have to do it for every layer you want to be cut, the dialogue box will stay open after cutting. If the polygon layer contains just one polygon, no selecting is necessary. BTW: make sure all layers share the same CRS (i.e. turn on-the-fly-reprojection off) HTH, André Joost ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] Cutting lines with a polygon boundary
Hi, I'm battling finding an easy way to cut some lines with a polygon boundary, and then deleting the line portions that are OUTSIDE the polygon. Scenario: 4 separate SHAPE files with lines in them 1 shapefile with one polygon in it. I need to cut the lines in all 4 shapefiles at the point that their lines intersect the polygon boundary. Thereafter I want to delete the line portions that are outside of this polygon boundary. In other words, to be able to crop multiple layers according to some polygon boundary. I would have thought this to be a trivial issue, but I have tried: * Vector Geometry Processing Clip, but it wants to extract the data inside the polygon and create a new file with the extracted data. (To be expected) * Plugin 'Clipper': I select my lonely polygon, and then activate one of the line layers (which is also in edit mode). The plugin seems to do nothing when I choose 'intersection' or 'clipping' preview, and when I choose 'clip' it tells me to select at least one feature (but my polygon is still selected.) The home-page link comes up as '404 not found' * Didn't notice anything in the Advanced Digitizing plugin. * CadTools and CadDigitize - no relevant facility found, but really nice 'Help' doc - thanks guys! * Feature 'Cut with Polygon from another Layer' deletes the pieces I want to keep. * Feature: 'Split selected features with line from another layer' is to some extent what I want, but I want '... with polygon from another layer' So, I eventually gave up and converted my polygon into a line layer, and used 'Split selected features with line from another layer', but this is incredibly tedious because: * It only does one layer at a time * I then have to select the line segments that I want to delete. So, does anyone know of a way to efficiently cut features in multiple layers using a polygon, and then to be able to select the pieces (multiple layers) that are either inside or outside the polygon, and then delete them? Regards & thanks in advance, Zoltan -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Plotter out of memory
On 2015/01/20 10:52, Richard Duivenvoorde wrote: On 20-01-15 09:48, Zoltan Szecsei wrote: Hi, I'm trying to print on A4 an orthphoto with an HP500ps Plus plotter and I keep getting 'out of memory' failures from the plotter. Is there a way to force the rendering in Windows and then feeding the print to the plotter in a way to avoid memory errors (on the plotter)? Using QGIS 2.6 on Win7 Prof 64 bit. I have tried both "winprint/RAW" and "hpcpp112/NTF EMF1.008" print processors. Hi Zoltan, I know it sound like Windows-bashing, but the one time I needed to print out a huge A0 print via Windows I also had memory problems, rebooting into Linux and doing the exact same thing on the same version of QGIS 'just did it'. Not sure if it was coincedence, but it is worth a try? Regards, Richard Duivenvoorde ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user Yep. You're 100% correct. Just completed the whole plot on my Linux box. Would still be nice to know what Win driver/renderer does work Regards & thanks, Zoltan PS: I love windows bashing - my main desktop is Win7 64bit, but whenever I need to do anything 'production' related, I use my Linux setup. Now if only these Win-Twits would wean themselves of WORD/Excel/Etc and get happy with LibreOffice.... ;-) -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] Plotter out of memory
Hi, I'm trying to print on A4 an orthphoto with an HP500ps Plus plotter and I keep getting 'out of memory' failures from the plotter. Is there a way to force the rendering in Windows and then feeding the print to the plotter in a way to avoid memory errors (on the plotter)? Using QGIS 2.6 on Win7 Prof 64 bit. I have tried both "winprint/RAW" and "hpcpp112/NTF EMF1.008" print processors. Regards & TIA, Zoltan -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Creating oblique views of QGIS layers
Cool - I'll take a look at GMT. Thanks Brent. Regards, Zoltan On 2015/01/18 19:33, Brent Wood wrote: Hi, I have not done this in QGIS, but have done it often with GMT. I've also used it to generate a page of maps showing a time series of survey results, each as a 3D perspective plot. A script draws a layer, then the next, applying a Z offset to each layer as it renders them, building the stack. You can use pstext to add text for each layer (which can be on the page - does not have to be in the map extent. At the end ps2raster can convert the postscript to bitmap which you cam add to your layout in the QGIS composer. Depending on how complex your layers are, the script can get a bit complicated... the nature of the beast. QGIS (or GDAL) can save your QGIS vector layers in GMT format, rasters can also be converted to a format GMT can render. Brent Wood ---- *From:* Zoltan Szecsei *To:* "qgis-user@lists.osgeo.org" *Sent:* Monday, January 19, 2015 2:47 AM *Subject:* [Qgis-user] Creating oblique views of QGIS layers Hi, I need to create a "layer stack" view of some QGIS vector data. Is there a way to do this in QGIS? (Rather like what you get when you Google "Images of GIS layer") I want it for a printed map view, and be in place of a traditional 'Key' in the one corner of the map plot. Regards & TIA, Zoltan -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org <mailto:Qgis-user@lists.osgeo.org> http://lists.osgeo.org/mailman/listinfo/qgis-user -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] Creating oblique views of QGIS layers
Hi, I need to create a "layer stack" view of some QGIS vector data. Is there a way to do this in QGIS? (Rather like what you get when you Google "Images of GIS layer") I want it for a printed map view, and be in place of a traditional 'Key' in the one corner of the map plot. Regards & TIA, Zoltan -- ======= Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Which tool/referecence for adress geocoding ?
Hi, I'd be interested in this in a South African context. If you get replies that are off-list, please make sure to share them with this list. Thanks & regards, Zoltan On 2015/01/14 11:15, image93 wrote: Dear all, i have to geocode several corporates sites ( a big excel file). I would like geocode this file thanks to the adress. 1/ Could you tell me the best tool for this task (accuracy, performance,ease of use,etc.). On my station, i have *Qgis, postgis, mapinfo* 2/ is it better to use a web tool or a desktop GIS? 3/ I guess I must have a reference, a database for the streets network? How do I choose it? And how to get it? In advance, thank you very much. -- View this message in context: http://osgeo-org.1560.x6.nabble.com/Which-tool-referecence-for-adress-geocoding-tp5181397.html Sent from the Quantum GIS - User mailing list archive at Nabble.com. ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user -- ======= Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] importing Text from dxf file
Dear Andreas, a quick workaround to get text angle and height, I've written a gawk script, which makes csv file from the TEXT entities in a DXF file. CSV can be loaded to QGIS as delimited text layer. Angle and height columns can be used in labelling. Best Regards, Zoltan On Mon, 12 Jan 2015, Andreas Neumann wrote: Hi, There is certainly room for improvement regarding DXF. We first concentrated on the DXF export, which is in already quite good shape in QGIS master. Once this is done, we should invest in DXF reading as well. Until then, there are only workarounds, sorry. Andreas Am 2015-01-12 12:36, schrieb Siki Zoltan: Dear Giacomo, yes, that is the way what you found. There is no text/annotation support in most of the GIS formats (e.g. shape file, PostGIS, SpatiaLite, etc.). You can make a label from the Text column. What is not imported from the DXF file are the angle, size and font of the text. Best regards, Zoltan On Mon, 12 Jan 2015, Giacomo Minotti wrote: Good morning, is it possible to import in qgis a dxf file that contains text? When I import the file, i see the text like point thanks Giacomo ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user BEGIN { entity = ""; last =""; i = 0; printf "id;layer;x;y;angle;height;text\n"; } /ENTITIES/,/EOF/ { if ($0 ~ "^ *0?$" && entity == "TEXT") { i++; printf "%d;\"%s\";%.3f;%.3f;%.4f;%.1f;\"%s\"\n", i, layer, x, y, angle, height, text; entity = ""; } if (last ~ "^ *0$") { entity = $1; } if (last ~ "^ *1$") { text = $1; } if (last ~ "^ *8$") { layer = $1; } if (last ~ "^ *10$") { x = $1; } if (last ~ "^ *20$") { y = $1; } if (last ~ "^ *40$") { height = $1; } if (last ~ "^ *50$") { angle = $1; } last = $0; } ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] importing Text from dxf file
Dear Giacomo, yes, that is the way what you found. There is no text/annotation support in most of the GIS formats (e.g. shape file, PostGIS, SpatiaLite, etc.). You can make a label from the Text column. What is not imported from the DXF file are the angle, size and font of the text. Best regards, Zoltan On Mon, 12 Jan 2015, Giacomo Minotti wrote: Good morning, is it possible to import in qgis a dxf file that contains text? When I import the file, i see the text like point thanks Giacomo ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Auto Open Form
Dear Prasun, when a single feature is identified the form of attributes is opened automatically, if the layer is editable you can change attributes too. Regards, Zoltan On Sun, 28 Dec 2014, Prasun Kumar Gupta wrote: Hi All, Kind of naive question...but could someone help me understand what is "Auto Open Form" checkbox used for. There was no information about it on manual page [1] Same as been posted here also [2] Environment: Windows 7 QGIS 2.6 - 64 bit Feature - Identify Results Regards, PKG [1] http://docs.qgis.org/2.6/en/docs/user_manual/introduction/general_tools.html [2] http://gis.stackexchange.com/questions/127679/qgis-2-6-auto-open-form ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] open file .mxd
Hi Maria, there is a solution on the other side (in ArcGIS) see: https://dragons8mycat.wordpress.com/2014/04/17/mxd-to-qgis/ Regards, Zoltan On Thu, 11 Dec 2014, Alex Mandel wrote: On 12/11/2014 10:04 AM, Maria Susana Lopreiato wrote: Hi everybody, i want to know if know how to open files .mxd in qgis Thank you María No, you can't open .mxd. That format has not been reverse engineered to my knowledge by any non-ESRI applications. Thanks, Alex ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Switch on only 1 layer of a group
Click the little '+' sign next to the group name and the group will open, and then you can turn each layer on or off at will. When you save that project it should remember which maps in that group were on or off. If you change that state and you (or your user) saves the project, then yes, it's messed up, so are you asking how you can set a "default" state of which maps are open? That would be quite useful to have 'snapshots' of a project state, and be able to call them back, or to force one of these snapshots to be the default state at project QGIS fire-up time. (Perhaps this could be implemented by putting a flag in the .QGS file, that is only writeable by the user that originally set it?) Regards, Zoltan On 2014/11/14 08:33, Andrea Peri wrote: Hi, I need to set a qgis project to turn on only 1 layer from a group. When a layer is turn on the others should be turn off. I guess this is not possibile from group properties, but is possibile to set an action to turn off every layer in a group other than the actual layer turn on ? Thxs, Andrea ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Duplicating a vector layer
On 2014/11/13 18:38, Adrian Beasley wrote: I am new to QGIS, still very much learning how to use it. What I wish to do is generate a number of similar vector layers, to be able to generate slightly different maps from them. The obvious way would seem to be to duplicate an existing layer, renaming the new copy, and editing it. But this seems not to produce a new, distinct layer, as making changes to either the original or the new layer causes the same changes in both. I can find no mention of this facility in the user guide nor in the training documentation. Can anyone help? Thanks in anticipation. Adrian Beasley ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user Hi Adrian, Depends how similar you mean by similar. If the map geometry is the same but you just want to render the maps differently, take a look at the "style" options. There are some very power features for rendering according to field contents, etc. HTH, Zoltan -- ======= Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Android Camera App with compass bearing
Hi, Thanks for the thought - but see orange text below: Does not write bearing into EXIF fields. Regards, Zoltan On 2014/11/05 12:17, Lene Fischer wrote: GPS Essentials - App Regards Lene -Oprindelig meddelelse- Fra: qgis-user-boun...@lists.osgeo.org [mailto:qgis-user-boun...@lists.osgeo.org] På vegne af Zoltan Szecsei Sendt: 5. november 2014 07:02 Til: t...@wildintellect.com; qgis-user@lists.osgeo.org Emne: Re: [Qgis-user] Android Camera App with compass bearing Hi Alex/All, Thanks for the thoughts - my comments interspersed with your. Kind regards, Zoltan On 2014/11/04 23:17, Alex Mandel wrote: On 11/04/2014 04:55 AM, Zoltan Szecsei wrote: Hi, I need to do a field trip and take some photos. I've been googling around, but cannot find an app that will burn both GPS position _and compass bearing_ into either the image, or just the exif data for that image. Anyone got any experience in this need? Thanks & regards, Zoltan Maybe, http://geopaparazzi.github.io/geopaparazzi/ "georeferenced and orientated pictures" And now my long rant before I remembered the above link... *I manually downloaded the latest version (newer than that on Play Store) and it still does not record bearing* - even though the app does have a bearing tool. I checked the website and recording bearing into EXIF is not even on the developer's long wish-list I have emailed the developer, with a hope to at least get it onto the wish-list. GPS data is an on/off toggle on most Android devices. Bearing does appear possible according to the exif specs, but in practice is rare. http://www.sno.phy.queensu.ca/~phil/exiftool/TagNames/GPS.html Yes, I checked that last night too - TAGs 0x0017 and 0x0018 relate to bearing info. However I don't see any obvious apps for telling the phone you want to capture that. I'm also not sure all phones are capable of calculating bearing when not moving. So my suggested work around would be a GPS Logging application of which there are several. Then you can pull the bearing calculation out after the fact, if it was while moving. Sigh, nope. (In my kind of) Field trip data collection (& photography) often involves driving to a feature, leaning out the car window and taking a snap-shot, so recording moving direction will not help at all. Note that bearing accuracy is NOT that critical, it needs only to be reasonably indicative. If stationary, you might need to write your own app which lets users indicate which way they were facing when taking a picture. For someone who's made an Android application I don't think this is hard. http://developer.android.com/guide/topics/sensors/sensors_overview.htm l Indicates all version of Android support the Magentic sensor. This example shows how to get the bearing from the API https://www.codeofaninja.com/2013/08/android-compass-code-example.html So GPS+Orientation+Camera should be all required. Yep, true, but the purpose of this thread is to see if this wheel has already been invented :-) Kinda busts my GPS logging idea, since the bearing information comes from a different sensor than the GPS it won't be in the NMEA strings even if you captured them. What? You not 'hard core' enough to intercept the NMEA stream and insert "bearing-on-the-fly" by back-calculating the incoming coords? (Just pulling your leg - thanks for your interest in this topic :-) ) Thanks, Alex -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Android Camera App with compass bearing
Hi Alex/All, Thanks for the thoughts - my comments interspersed with your. Kind regards, Zoltan On 2014/11/04 23:17, Alex Mandel wrote: On 11/04/2014 04:55 AM, Zoltan Szecsei wrote: Hi, I need to do a field trip and take some photos. I've been googling around, but cannot find an app that will burn both GPS position _and compass bearing_ into either the image, or just the exif data for that image. Anyone got any experience in this need? Thanks & regards, Zoltan Maybe, http://geopaparazzi.github.io/geopaparazzi/ "georeferenced and orientated pictures" And now my long rant before I remembered the above link... I manually downloaded the latest version (newer than that on Play Store) and it still does not record bearing - even though the app does have a bearing tool. I checked the website and recording bearing into EXIF is not even on the developer's long wish-list I have emailed the developer, with a hope to at least get it onto the wish-list. GPS data is an on/off toggle on most Android devices. Bearing does appear possible according to the exif specs, but in practice is rare. http://www.sno.phy.queensu.ca/~phil/exiftool/TagNames/GPS.html Yes, I checked that last night too - TAGs 0x0017 and 0x0018 relate to bearing info. However I don't see any obvious apps for telling the phone you want to capture that. I'm also not sure all phones are capable of calculating bearing when not moving. So my suggested work around would be a GPS Logging application of which there are several. Then you can pull the bearing calculation out after the fact, if it was while moving. Sigh, nope. (In my kind of) Field trip data collection (& photography) often involves driving to a feature, leaning out the car window and taking a snap-shot, so recording moving direction will not help at all. Note that bearing accuracy is NOT that critical, it needs only to be reasonably indicative. If stationary, you might need to write your own app which lets users indicate which way they were facing when taking a picture. For someone who's made an Android application I don't think this is hard. http://developer.android.com/guide/topics/sensors/sensors_overview.html Indicates all version of Android support the Magentic sensor. This example shows how to get the bearing from the API https://www.codeofaninja.com/2013/08/android-compass-code-example.html So GPS+Orientation+Camera should be all required. Yep, true, but the purpose of this thread is to see if this wheel has already been invented :-) Kinda busts my GPS logging idea, since the bearing information comes from a different sensor than the GPS it won't be in the NMEA strings even if you captured them. What? You not 'hard core' enough to intercept the NMEA stream and insert "bearing-on-the-fly" by back-calculating the incoming coords? (Just pulling your leg - thanks for your interest in this topic :-) ) Thanks, Alex -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] Android Camera App with compass bearing
Hi, I need to do a field trip and take some photos. I've been googling around, but cannot find an app that will burn both GPS position _and compass bearing_ into either the image, or just the exif data for that image. Anyone got any experience in this need? Thanks & regards, Zoltan -- ======= Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Import gpx files into Postgis db (windows)?
QGIS cannot load gpx to postgis directly, but db manager plug-in can. Zoltan On Fri, 31 Oct 2014, dandrigo wrote: Hello, Ok thanks. I will have a look on ogr2ogr tool. Regarding qgis, what qgis tool could import gpx files directly in the postgis database? Perhaps you alluded to the Qgis DB manager? Thanks. -- View this message in context: http://osgeo-org.1560.x6.nabble.com/Import-gpx-files-into-Postgis-db-windows-tp5170663p5170677.html Sent from the Quantum GIS - User mailing list archive at Nabble.com. ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Import gpx files into Postgis db (windows)?
Dear Dandrigo, if you would like to automatize the import use ogr2ogr (from GDAL/OGR utilities). Another interactive solution may be QGIS. Regards, Zoltan On Fri, 31 Oct 2014, dandrigo wrote: Dear all, I'm working with windows 8 & postgis 2.1. I have several gps (gpx) files. I would like to import *directly* those gpx files into my postgis database. What would be the easiest way? (OGR tool ?) In advance, thank you to throw light for me. Regards. -- View this message in context: http://osgeo-org.1560.x6.nabble.com/Import-gpx-files-into-Postgis-db-windows-tp5170663.html Sent from the Quantum GIS - User mailing list archive at Nabble.com. ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Raster VRT SRS syntax
On 2014/10/29 14:37, Even Rouault wrote: Perhaps the following can work too (untested) provided a 'ZANGI' file exists in the PROJ.4 resource directory that contains the 'epsg' file, grids, and other files. +init=ZANGI:ZANGI:HBKNO27 Hi Even, When I let gdalbuildvrt simply add its own SRS from the image files, when I load that VRT into QGIS, it loads correctly without asking me to specify the CRS for that layer. However When I query that layer to see its CRS, it shows me a Generated CRS, even though the HBKNO27 is in the known list (see screenshot below). I want the user to see the named CRS and not the generated (more confusing) CRS. The first few lines of my VRT is as follows, yet I still get the pop-up asking me to set the CRS. (I'm still trying different permutations, but I just wanted to give feedback) Regards, Zoltan +init=ZANGI:ZANGI:HBKNO27 9.5763066019661390e+04, 5.e-01, 0.e+00, -2.53362775e+06, 0.e+00, -5.e-01 0.00E+00 Red -- ======= Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] Raster VRT SRS syntax
Hi, I'm running an older version of gdalbuildvrt that does not recognise the a_srs parameter, so I want to replace line 2 of the resulting VRT file with my own SRS entry. I have tried: HBKNO27 and various other permutations with and without the 'ZANGI:' prefix. From the GDAL tutorial, I know the EPSG code syntax style is EPSG:4326 The CRS I need is ZANGI:ZANGI:HBKNO27 or 'South African CRS: HBK_NO_27' Can someone please guide me to the correct syntax. TIA, Zoltan -- ======= Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] How to calculate the transformation parameters for the QGIS Affine Plugin
Hello Marcel, to calculate affine transformation parameters you must have at least 3 points in both coordinate systems. LV95_y = Y0 + a * local_y + b * local_x LV95_x = X0 + c * local_x + d * local_y If you set up the equations for the three points, you get to linear equation systems, one for Y0, a, b unknowns, and another fo X0, c, d unknowns. It is highly recommended to use more points and least squares estimation for the parameters. Best regards, Zoltan On Tue, 7 Oct 2014, Marcel Stadelmann wrote: Hi, I would like to make an affine transformation with a vector layer shape file, using the Affine Plugin. I know from two points the correct coordinates: local_y local_x LV95_y LV95_x 9.0650567.70714 2753308.44 1166381.26 9.065056-5.085236 2753295.68 1166391.31 How can I calculate the transformation parameters? Marcel ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] FW: Export Data
Dear Paolo, the new version is available from the QGIS offcial repository as an experimental plugin. Can the original owner of the plugin give me rights to upload newer versions? Best regards, Zoltan On Tue, 30 Sep 2014, Paolo Cavallini wrote: Il 30/09/2014 13:22, Siki Zoltan ha scritto: Dear Paolo, originally it is not my project, I cloned and upgraded it and sent back to the original developer to upload to the plugin page. It has not happend yet :( I have not right to upload to that plugin. I didn't want to start a new project under my name. Right. Could you please remind the original author? In case of no reply, I can grant you upload access to the original one. All the best, and thanks. -- Paolo Cavallini - www.faunalia.eu Corsi QGIS e PostGIS: http://www.faunalia.eu/training.html ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] FW: Export Data
Dear Paolo, originally it is not my project, I cloned and upgraded it and sent back to the original developer to upload to the plugin page. It has not happend yet :( I have not right to upload to that plugin. I didn't want to start a new project under my name. Zoltan On Tue, 30 Sep 2014, Paolo Cavallini wrote: Il 29/09/2014 21:48, Siki Zoltan ha scritto: the BulkVectorExport plugin can do it for vector data. The version on the plugins.qgis.org page works with 1.8 version. I've updated it to 2.x and added some new functionality. It is available from GitHub https://github.com/zsiki/bulkvectorexport Hi Zoltan, why not upgrading the pacikage in http://plugins.qgis.org/plugins/BulkVectorExport/ ? We can ask (again?) the original author, and add a new version in case he does not reply in a reasonable tiome frame (hey, it's GPL after all!). All the best, and thanks. -- Paolo Cavallini - www.faunalia.eu Corsi QGIS e PostGIS: http://www.faunalia.eu/training.html ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] FW: Export Data
Dear Nathan, the BulkVectorExport plugin can do it for vector data. The version on the plugins.qgis.org page works with 1.8 version. I've updated it to 2.x and added some new functionality. It is available from GitHub https://github.com/zsiki/bulkvectorexport Regards, Zoltan On Mon, 29 Sep 2014, Daniel S. wrote: Thanks for the fast respons. This Plugin looks like what I need. But after manually installation it show me a malfunction (invalid syntax). Maybe you know a nother Plugin, Way? -- View this message in context: http://osgeo-org.1560.x6.nabble.com/Export-Data-tp5164722p5164767.html Sent from the Quantum GIS - User mailing list archive at Nabble.com. ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] QGIS and Shellshock
On 2014/09/25 14:54, Steven Campbell wrote: Hi all Does anyone know if QGIS is affected at all by the bash bug? Linux is affected by this bug, not QGIS. So if you're running Linux, then yes, you may be affected. You can check by opening a terminal window, and running this (cut & paste it into the window: env x='() { :;}; echo vulnerable' bash -c 'echo hello' If you get the "Vulnerable" echo from this, then yes, your version of bash has the bug. You can fix this by updating to the latest bash version. sudo apt-get update && sudo apt-get install bash (Debian based Linux assumed, eg: Ubuntu) After running this, rerun the "env" line above to see if it is fixed. If still buggy, then you probably have an older version of linux, so you'll have to download bash and recompile it on your system. (or upgrade to a newer linux, then add the newest bash using the above apt-get..) To recompile bash, you'll need build-essentials installed. If you have an older (no longer supported linux, you may be up against a wall if you can't run gcc & make), but assuming your installed linux version already has gcc, make etc installed, run the following: mkdir src cd src wget http://ftp.gnu.org/gnu/bash/bash-4.3.tar.gz #download all patches for i in $(seq -f "%03g" 0 25); do wget http://ftp.gnu.org/gnu/bash/bash-4.3-patches/bash43-$i; done tar zxvf bash-4.3.tar.gz cd bash-4.3 #apply all patches for i in $(seq -f "%03g" 0 25);do patch -p0 < ../bash43-$i; done #build and install ./configure && make && make install cd .. cd .. rm -r src Then reboot (to be safe all previous bash shells have closed) and of course run the "env." command to see if your bug has gone away. Disclaimer: I got this methodology from http://askubuntu.com/questions/528101/what-is-the-cve-2014-6271-bash-vulnerability-and-how-do-i-fix-it My one linux box is current and the apt-get install bash fixed the problem. My server is still running an older outdated linux, and the above compile-install DID work on it. You use this methodology at your own risk - but it worked on both my linux boxes (this morning). Hope this helps. Regards, Zoltan Thanks Steve *Steve Campbell|GIS Manager* Corporate Strategy and Communications Borough of Poole | Civic Centre | Poole BH15 2RU Tel: 01202 633 362 Email: s.campb...@poole.gov.uk <mailto:s.campb...@poole.gov.uk> Website:www.boroughofpoole.com <http://www.boroughofpoole.com/> */Think Green! Please Recycle/* DISCLAIMER: This email and any files transmitted with it may be confidential, legally privileged and protected in law and are intended solely for the use of the individual to whom it is addressed. The copyright in all documentation is the property of the Borough of Poole and this email and any documentation must not be copied or used other than as strictly necessary for the purpose of this email, without prior written consent which may be subject to conditions. Any view or opinions presented are solely those of the author and do not necessarily represent those of the Borough of Poole. The Borough of Poole reserves the right to inspect incoming and outgoing emails. If you have received this email in error please contact the sender by return and confirm that its contents have been destroyed. Telephone enquiries should be directed to the Borough switchboard on 01202 633633.' ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] selecting features containing a polygon
Hi Giacomo, try Spatial Query plugin. See: https://qgis.readthedocs.org/en/latest/docs/user_manual/plugins/plugins_spatial_query.html http://www.qgistutorials.com/en/docs/performing_spatial_queries.html Zoltan On Mon, 8 Sep 2014, Giacomo Fontanelli wrote: Hello everybody I have a shapefile named e.g. "pol_A" made by many polygons. How to select the polygons containing another polygon constituted by the shape file "pol_B"?. Thank you very much. --- Questa e-mail ? priva di virus e malware perché ? attiva la protezione avast! Antivirus. http://www.avast.com ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] polygon centroids
AFAIK xmin, ymin is not available in QGIS 2.2, you must upgrade to 2.4. Zoltan On Thu, 28 Aug 2014, Alexandre Neto wrote: Ups, my bad... You can get your features centroid using the function centroid(). It will return a point geometry, then you can get the x and y by using either x_min or x_max and y_min or y_max. Like this: xmin( centroid( $geometry )) || ' , ' || ymin( centroid( $geometry )) Alexandre Neto On Thu, Aug 28, 2014 at 10:29 AM, Leo Kris Palao wrote: Hi Alexander, Thanks for your reply. In my understanding, this only applies if I have point coordinates, right? If I have polygons I think I need another syntax? I can create centroids in Vector > Geometry > Polygon Centroids tool then apply $x and $y to extract lat lon coordinates, but in some cases I want to skip this process if I am dealing with more than 300 thousand features in my polygon layer. Thanks, -Leo On Thu, Aug 28, 2014 at 3:01 PM, Alexandre Neto wrote: Hello Leo, That's quite straightforward actually. In the field calculator (or in any other instance of expression builder) there are 2 geometry functions to achieve what you need, $x and $y. You can either put them in separate fields or create a string to populate a single one. Something like this, $x || ' , ' || $y Would produce something like this: '1234.56 , 9876.43' Notice that the double pipe character || works as a string concatenator. And that the result is a string therefore not useful for any arithmetic operations after. Hope this helped. Best regards, Alexandre Neto Em 27/08/2014 22:37, "Leo Kris Palao" escreveu: Hi QGIS users, I don't know if this question has been asked before. Can you help me create a syntax in field calculator to add in the attribute table the centroids coordinates for each or selected polygon features? I am using Qgis 2.2 "Valmiera" version in Windows. Thanks in advance for help, -Leo ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Error when loading a shape file
Hi, your shape file contains an error message: "Error with /opt/atmail/tmp/nik...@stigfinnaren.nu/objekt15.shp - Please check the directory exists and the session has not timed out" It is the same for the dbf and shx files. Best regards, Zoltan On Thu, 28 Aug 2014, Niklas Groop wrote: Can anyone take a look at this file and see if there is any problem. Niklas Groop Polerarvägen 1 783 50 Gustafs 076-3411000 ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Field calculator QGIS2.4 division error
Hi Nathan, are you sure? QgsExpression('100.0 * 50 / 3').evaluate() 1666.6667 QgsExpression('50 / 3 * 100.0').evaluate() 1600.0 Zoltan On Sat, 23 Aug 2014, Nathan Woodrow wrote: Hey, There is no left to right here, / is higher order then * so both these are the same: exp = QgsExpression('100.0 * 50 / 2') exp2 = QgsExpression('50 / 2 * 100.0') exp.evaluate(), exp2.evaluate() (2500.0, 2500.0) What you are seeing is the QGIS expression engine doing the multiply part result first, which returns a float, then it will say "is the left side a int AND the right side a int" if that comes back as NO then it does a float operation. Basically this means once one of the values in the equation is a float the result will come out as float. If you want to see the details on how this works you can find it here https://github.com/qgis/QGIS/blob/master/src/core/qgsexpression.cpp#L2093 - Nathan On Sat, Aug 23, 2014 at 8:22 PM, ralfwessels wrote: Hi intereseting, didn't know the left to right rule explained by Zoltan: Siki Zoltan wrote 100.0 * "HC01_VC04" / "HC01_VC03" will give you the right result. First the multiplication is evaluated (left to right rule),... Learning never stops! ralf -- View this message in context: http://osgeo-org.1560.x6.nabble.com/Re-Field-calculator-QGIS2-4-division-error-tp5157715p5157980.html Sent from the Quantum GIS - User mailing list archive at Nabble.com. ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Field calculator QGIS2.4 division error
Dear Ivan, it is not strange, c/c++ evaluate expressions this way. 100.0 * "HC01_VC04" / "HC01_VC03" will give you the right result. First the multiplication is evaluated (left to right rule), the result is float, then float over integer gives float result. / is really two operators, integer division and float division. If both operands are integer then integer division is used. Cheers, Zoltan On Fri, 22 Aug 2014, Ivan Santiago wrote: Hello all: Thanks Ralf. The columns used are integer. The resulting column is real, lenght=10 precision=1 Writing "HC01_VC04" / "HC01_VC03" * 100.0 Does not work either. Applying convert function toreal() to one of the columns, produces correct results. I checked it against LibreOffice Calc. BTW: The calculator has a strange behavior though. I need to sort and sort again the calculated column to 'refresh' the table values. QGIS 2.4 32bit under WIN7. --- Iván Santiago GIS Specialist Information Technologies Office of Management and Budget 787.725.9420 x 2378 Calle Cruz 254 PO Box 9023228 San Juan, PR 00902-3228 http://gis.pr.gov<http://gis.pr.gov/> From: Ivan Santiago Sent: Thursday, August 21, 2014 4:04 PM To: 'qgis-user@lists.osgeo.org' Subject: RE: Field calculator QGIS2.4 division error Hello all: Does anyone have had problems with field calculations using the division operator? I've tried many times to update a new field using both the old field calculator and the new interface but they just don't work with division. Simple operations like: "HC01_VC04" / "HC01_VC03" * 100 Do not work. The only workaround I found was to put the function abs into the denominator: "HC01_VC04" / abs("HC01_VC03") * 100 Then it gives the result but this should be avoidable. I'll send the table used if anyone want to try it. Qgis version is 2.4 (32bit) under WIN7. --- Iván Santiago GIS Specialist Information Technologies Office of Management and Budget 787.725.9420 x 2378 Calle Cruz 254 PO Box 9023228 San Juan, PR 00902-3228 http://gis.pr.gov<http://gis.pr.gov/> ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Extract lines from raster
Dear Matt, there is Polygonize menu option in Raster/Convert. It uses gdal and creates vector polygons for all connected regions of pixels in the raster sharing a common pixel value. Probably this is not the best for you if you have lines in your raster. For such tasks I usually use GRASS commands 0. r.in.gdal import your raster into GRASS 1. r.thin to narrow raster lines to 1 pixel 2. r.to.vect to create line features from raster 3. v.generalize to simplify vector lines (there are a lots of methods available) 4. v.out.org export your vector (you can open and edit GRASS maps in QGIS too) Simplification is available also in QGIS see: Vector menu -> Geometry tools -> Simplify geometries Regards, Zoltan On Tue, 19 Aug 2014, Matt Boyd wrote: Hi List, I've got a scanned map I want to extract / digitise the contours from as lines. The lines are all black and if I can get everything that's black into a shapefile I'd be able toedit and assign values once there. Is there anything in qgis that coudl do this? I had a quick look at orfeo but for some reason the website is down today. Thanks Matt ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Flip
Dear Diego, here is a short python solution to flip all selected polylines in the actual layer. You can use it in the python console. layer = qgis.utils.iface.mapCanvas().currentLayer() for feature in layer.selectedFeatures(): geom = feature.geometry() nodes = geom.asPolyline() nodes.reverse() newgeom = QgsGeometry.fromPolyline(nodes) layer.changeGeometry(feature.id(),newgeom) qgis.utils.iface.mapCanvas().refresh() Zoltan On Wed, 13 Aug 2014, Diego Lawrie wrote: Dear all, how can I switch line direction in QGIS? Thanks ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Enhancing Raster Images
Hi, I suppose blending is not the best solution, because i depends on other layers. Why don't you use transparancy tab of the properties of the raster? There you can set pixel values to be transparent (0-100%). In your case from 0 to 0 value 100% transparent Then change render type to Singleband pseudocolor on the style tab and set individual color for any value. In your case add 1 to the list and select a color for it. Zoltan On Tue, 12 Aug 2014, G. Garibaldi wrote: Beginning with a black and white raster, black lines on a white background, I'd like to end up with a raster image with no white background with only the black lines showing. Additionally, I'd like the black lines to be a different color, bright red, bright green, etc. This is what I've done: Select raster layer properties>color rendering>blending mode>burn. This removes the white background and shows only the black lines. Going further and selecting properties>hue>colorize only lightly tints the black lines to the selected color, not really the desired result. How can I make the properties>hue>colorize change the black lines to some bright highly saturated color? ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Update field from other table?
Dear Tomas, you need spatial join, see http://www.qgistutorials.com/en/docs/performing_spatial_joins.html Regards, Zoltan On Tue, 5 Aug 2014, Tomas Renco wrote: Hello, I am struggling with elementary problem: I have 2 layers with attribute tables. I need to update one attribute in table A with values from table B based on geometry - objects in layer A are contained in objects from layer B. I tried to make this query in Attribute calculator, but I couldn't figure out how to work with 2 layers. I am looking forward for your advice, I am sure I am missing something elementary. TIA Tomas Renco ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Installation 2.4
On 2014/08/05 15:57, Udo Christiansen wrote: thanks, your answer was very helpful ;-) What was embarrassing (for me) was that, even though being born in Hungary, I had to use Google to translate for me :-) Good luck with solving your issue. Cheers for now, Zoltan Büro für Vegetationskunde und Landschaftsökologie Dipl.-Biol. Udo Christiansen Benediktinerstr. 22 67549 Worms Tel. 06241-594684 Mobil 0178 868 00 43 emailudo_christian...@online.de http://www.vhoe.de/bueros/bvl-christiansen/ Am 05.08.2014 15:46, schrieb Zoltan Szecsei: Elnézést, nem tudok németül, csak magyarül :-) On 2014/08/05 15:30, Udo Christiansen wrote: Hallo, ich hoffe, ich bekomme auch in Deutsch eine Antwort ... Ich habe problemlos QGis 2.4 installiert (glaube ich) und wollte dann ein recht umfangreiches Projekt aufbauen. Ab einer gewissen Anzahl von Luftbildern stürzt QGis ab mit der Bemerkung, ich würde informiert, wenn das Problem lösbar sei. Gibt es irgendwelche Einschränkungen, die ich beachten muss? Oder wo liegt der Fehler? Danke für Antworten Udo -- Büro für Vegetationskunde und Landschaftsökologie Dipl.-Biol. Udo Christiansen Benediktinerstr. 22 67549 Worms Tel. 06241-594684 Mobil 0178 868 00 43 emailudo_christian...@online.de http://www.vhoe.de/bueros/bvl-christiansen/ ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] Installation 2.4
Elnézést, nem tudok németül, csak magyarül :-) On 2014/08/05 15:30, Udo Christiansen wrote: Hallo, ich hoffe, ich bekomme auch in Deutsch eine Antwort ... Ich habe problemlos QGis 2.4 installiert (glaube ich) und wollte dann ein recht umfangreiches Projekt aufbauen. Ab einer gewissen Anzahl von Luftbildern stürzt QGis ab mit der Bemerkung, ich würde informiert, wenn das Problem lösbar sei. Gibt es irgendwelche Einschränkungen, die ich beachten muss? Oder wo liegt der Fehler? Danke für Antworten Udo -- Büro für Vegetationskunde und Landschaftsökologie Dipl.-Biol. Udo Christiansen Benediktinerstr. 22 67549 Worms Tel. 06241-594684 Mobil 0178 868 00 43 emailudo_christian...@online.de http://www.vhoe.de/bueros/bvl-christiansen/ ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user -- === Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax:+27-86-6115323 www.geograph.co.za === ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user