Re: [Qgis-user] Advance Digitising Tools- compass bearings
Hi Calvin: An interesting question. Another consideration is the coordinate system for your project. If it is a projected system, then you should consider convergence which is the angle between true north and grid north. In a common example, UTM zone meridians are the center of each zone which are given an arbitrary easting of 500,000m. If you are converting from magnetic to true north close to the meridian, then true north will be very close to grid north. If you are east of the meridian, perhaps close to the adjacent zone boundary, then grid north and true north will vary noticeably. We developed a script for calculating convergence in UTM as I am unaware of a plugin provides that level of functionality. Ultimately, it will depend on the coordinate system you are working in. If you reach out to me privately, I can share the script with you. Kirk Schmidt k...@nortekresources.com On 5/17/2024 9:15 AM, Jeremy Jackson via QGIS-User wrote: On 2024-05-17 07:00, C Hamilton via QGIS-User wrote: Phil, The Shape Tools plugin has an "Azimuth distance digitizer" and for what you are talking about an "Azimuth distance sequence digitizer". With the latter function you click on a starting location and then give a list of azimuth, distance sequences. I have used this for old surveys like this. There is also an declination offset just in case your surveys were using magnetic north and not true north. In that case you have to look up what the A related question, what about grid to ground conversions and vice-versa, are there plugins which support this? Surveys here in Canada always specify a "combined scale factor", which is tedious to apply manually to each vertex when digitizing. declination offset was for your 1880s date. Regards, Calvin On Thu, May 16, 2024 at 9:50 PM Philip Ryan via QGIS-User wrote: Hello, I am trying to digitise a series of old gold lease "portion plans" produced in the 1880s in Australia by a mining surveyor using a compass and chain. The data for each gold lease includes a table of 1) compass bearings in degrees and 2) distance measurements in links (0.2 m) for each vertex of the lease polygon (mostly a rectangle or trapezoid but sometimes more complex). The Advance Digitizing Tools would seem to be ideal for this data except for the angle measure (a) which I can't input as a compass bearing (0-360 deg). The current settings with 90 deg E as the zero makes conversion of compass bearings very tedious. Is it possible to change the angle variable (a) to a compass bearing? Regards, Phil Ryan ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] gridding with a polygon
Hi Byron: If I understand this correctly, you are concerned with how the contours will be generated at the edge of your raster data. If I have this correct, I would grid the data that you have based on the attribute you want to contour. Create a vector clip layer for your study area that is at least 1 grid cell smaller than the the boundary size of the original vector data. Then create the contours from this new raster and you should be good to go. Kirk Schmidt On 5/2/2024 2:07 PM, Byron Veilleux via QGIS-User wrote: i have a large geologic dataset that i want to grid and contour. There is an irregular zero edge that bounds these data on one side and a straight bounding area edge on the other sides. If i create an irregular polygon that encloses these data, is there a way to grid and contour such that these data are truncated against the polygon and do not extend beyond? -- Byron Veilleux, MSc. P.Geo Conjugate Geologic Services Limited Calgary, Alberta Canada by...@conjugategeo.com Cell:4037108414 ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Problem with contour algorithm
Hi Grant: Based on the other comments, I tested a shapefile and the following GDAL command works gdal_contour -b 1 -a ELEV -i 10.0 -f "ESRI Shapefile" "K:/Work2/McCallum Environmental/Wedgeport/Land Clearing/CHM Clearing.sdat" K:/Temp/contour_test.shp The only difference is that in your example you have /vsimem/xxx.tif. There are no double quotes so I do not know if you are using a path variable or similar. I would try an actual path to the tif dem file. Kirk Schmidt On 3/13/2024 12:37 AM, Grant Boxer via QGIS-User wrote: I am trying to contour a grid file using the QGIS raster > extract>contour and I get getting errors saying “Warning 1: Model file not specified” or “ERROR 4: No such file or directory”. I am using 3.34 and 3.36 and get the same error. Any thoughts? Regards Grant Boxer, Perth WA ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Problem with contour algorithm
Hi Grant.Check to make sure your destination directory actually exists and you have write permission. I seem to remember running into a similar issue but I am not close to a computer to test. Kirk SchmidtSent from my Galaxy Original message From: Grant Boxer via QGIS-User Date: 2024-03-13 12:55 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Problem with contour algorithm I am trying to contour a grid file using the QGIS raster > extract>contour and I get getting errors saying “Warning 1: Model file not specified” or “ERROR 4: No such file or directory”. I am using 3.34 and 3.36 and get the same error. Any thoughts? Regards Grant Boxer, Perth WA ___ 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] Issues Extracting LiDAR From a Tif file to a Golden 7 Software Binary Grid (.GRD)
Hi Hudson: I just tried to export to golden 7 binary with QGIS version 3.20 and 3.28 running on windows 8 and it worked fine. Kirk On 2/13/2024 12:27 PM, Hudson Mealy via QGIS-User wrote: Hello, I just purchased a new computer. ROG GEOFORCE RTX STRYX and have downloaded QGIS 3.34.3. I was previously using a microsoft surface with GQIS 3.26.3. We use QGIS to extract lidar from a large area to individual fields (160-700 acres) to be used for drainage applications. The large lidar files are uploaded into QGIS in a tif file format. From there I use the raster extraction tool and clip raster by extent (The field boundary). From there I use the raster calculator tool, I choose the file destination which is into the ADMS drainage software. I then choose the golden software 7 binary grid file type for an output, and in the bottom box I multiply the clipped extent by 3.2808 to convert from metres to feet. From here I hit save, however QGIS does not respond and no matter how long I leave it it will not calculate the operation, I have to force quit the application. However, when going back to my old microsft surface everything works as normal. I am confused on to what has changed other than the newer QGIS download, why is my 5 year old laptop running this application without a hitch but my new desktop computer will not handle it. If anyone has any solutions that would be very helpful! Thanks -- Hudson Mealy Impact Tile Drainage Box 172 Minto, Manitoba R0K 1M0 204-305-0877 impacttiledrain...@gmail.com ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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 Change Layer coordinate fields
Hi Alana.You cannot edit the csv file, however, you can save it as a shapefile or other vector format. You can choose the coordinate system you want the data to reside in when you save it.Load the newly minted vector file snd you can recalculate or add new coordinate fields if required, however. If you imported the original csv and correctly set the coordinate system, then a recalculation of the coordinate fields will not be required to display the data correctly.Kirk SchmidtSent from my Galaxy Original message From: Alanna Mitchell via QGIS-User Date: 2024-01-24 5:55 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] How to Change Layer coordinate fields Dear group, I have a qgis project I have been working on. It was point data uploaded as a delimited text layer, and in the excel file (.csv) I have the data in lat/long and easting/northing (irish transverse mercator) and I haven't opened it in a while, but when I did it changed one of the fields to longtidue, so now one is projecting in easting and one in longitude. I wish to alter this from the actual file, I don't want to reupload it again as I have a lot of project layouts saved that I only need to alter the symbology on and resave. I haven't been able to figure this out. In my excel file I have the coordinates, so I wish to change the column the field is derived from. Is this possible? When I hover over the layer on QGIS, I can see the file path so even if I had a way to access this and change it, I think I would manage. Thank you for any help. Kind regards. ___ 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] building gdal for mrsid
Hi Don..Did you download the SDK from lizardtech?Kirk Schmidt Sent from my Galaxy Original message From: Don Harter via QGIS-User Date: 2024-01-20 10:10 p.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] building gdal for mrsid I have been following this procedure to build gdal with mrsid capabilities. I had done this a while back previously but the different dependencies between packages cause a mess. https://trac.osgeo.org/ubuntugis/wiki/TutorialMrSid I tried to run this command: sudo apt-get install libgdal-mrsid-src But that package is not available in the unbuntugis repository. ___ 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] Calculate upslope catchment areas for multiple points?
Hi Ben.The python script should work in your version since it is written in python 3. It us,a relatively simple routine which iterates over each sample point and rund the saga tool. It nay take a while to run 100s of points but it should finish.Saga is a third party, open source software suite that is developed independantly of QGIS. It is very powerful and fast for most raster processing. Saga also has native batch or scripting capability which may be an option for you. Just a note that the module syntax can change with versions. So and older script may not complete on a new release.Kirk SchmidtSent from my Galaxy Original message From: Ben Jones via QGIS-User Date: 2024-01-14 2:23 a.m. (GMT-04:00) To: Qgis Subject: [Qgis-user] Calculate upslope catchment areas for multiple points? Hello All,I have used SAGA to generate channel networks and r.watershed to generate catchment sub-basins from Copernicus 30 m DEM data to assist with planning a stream sediment sampling program. Once samples are collected I would like to generate polygons that define the upslope area for each sample point.I understand that this can be done using SAGA Terrain Analysis ‘Upslope area’ tool that requires manually copying and pasting the X and Y coordinates for each sample point (one point at a time which works but is not practical for 100’s of samples, or in batches which is still not ideal and doesn’t appear to work?).Is there an easy way to import a point file with multiple points (100’s of points) to SAGA or another tool to automatically calculate upslope catchment areas for each of the points?I have come across two potential solutions:1) ‘Qgis SAGA Upslope Area for several points automation’:https://gis.stackexchange.com/questions/342333/qgis-saga-upslope-area-for-several-points-automation2) ‘Snap pour points and delineate multiple catchments in QGIS’:https://www.youtube.com/watch?v=Hqc_1CMadhAHowever, both are three years old and it appears there are complications using both these methods in newer versions of QGIS3. Are there alternative methods?I am using QGIS 3.28.11a on macOS Monterey.Thanks,Ben___QGIS-User mailing listQGIS-User@lists.osgeo.orgList info: https://lists.osgeo.org/mailman/listinfo/qgis-userUnsubscribe: 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] Pansharpening Problem
Hi Grant: If you are ending up with null values and zero's, check the values in the raster you are pan sharpening against the training or higher resolution raster you are using to sharpen. Is there null data in the high res raster. Review the input requirements of the routine to ensure the tool will work with your data. Perhaps pan sharpening a 3 channel raster with a single channel raster creates the problem. Kirk Schmidt On 11/15/2023 9:52 PM, Grant Boxer via QGIS-User wrote: I am attempting to pansharpen a PRISMA hyperspectral RGB image and I am getting areas within the image that are clipped (no data). These appear to bands with zero values in two of the three bands. Anyone have a work around or suggest how this can be avoided. I am using the gdal pansharpen tool and latest release QGIS 3.34. Cheers Grant Perth WA ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] intersection question
Hi Matt.The trick is to use the attributes to id each layer before you run the intersection. So create a new field for each layer and give each polygon the same number. So layer 1, maybe 1. Do the same for the second or even third layer, use 2, or something different from the first layer. Complete your intersect and you can determine the origin of each resulting polygon.If you use 1 and 2 to identify the layers. New polygons with the attribute of 1 or 2 did not intersect. Polygons with 1 and 2 did.Kirk SchmidtSent from my Galaxy Original message From: Matthew Carroll via QGIS-User Date: 2023-11-11 3:45 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] intersection question Hi there,I'm fairly new to qgis, and hoping someone can point me in the right direction for a problem I'm trying to solve.I need to intersect two polygon layers, and identify which polygons from the original layer are, essentially, the ones that were "split" by the intersection process.In this example I'd want to identify polygons 1,2,4,5 since 3,6 are the same as in the original layer.Ideally those polygons would be identified somehow in the attribute table for the intersection layer.In addition, is there any way to build some tolerance into the process to account for spots where the polygon boundaries on layer 2 don't align 100% perfectly with layer one.Thanks so much for the help,Matthew ___ 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] filtering data from highest to lowest in QGIS
Hi Simon; Why yes of course you can. Please see the Filtering in Layer Properties section at the following link which provides the basics. https://mapscaping.com/filtering-in-qgis-a-comprehensive-guide/ You can use the array functions to select and sort attributes as required. Kirk Schmidt On 8/17/2023 12:09 PM, Simon via QGIS-User wrote: Hi all, Just had a quick question. Is it possible to filter and display the 10 highest values in a field using the Query Builder or Field calculator? Thanks, Simon ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Clipping ecw to mask not working
Hi Michael: Enhance Compression Wavelet (ECW) is a proprietary compression format. I would try exporting your image as a geotiff to see if it is being read correctly and does not have internal write protection on it. Kirk Schmidt On 7/26/2023 2:02 AM, Michael Dufty via QGIS-User wrote: I'm having issues trying to clip ecw images to a mask layer. Using raster-extraction-clip raster to mask layer. Default parameters. The tool runs extremely quickly (successfullyin 0.12 seconds), but produces a raster entirely black (0,0,0 for all pixels). With the correct extent of the mask layer. If I tick the box to create an alpha output band it produces an entirely transparent result. Source and mask layers both in EPSG 28350 Anyone else having this issue? I cannot see any bug reports. Tried in 328.8, 3.28.9, 3.32.0 and 3.32.1 with same result. Seems to work fine if input raster is a tif. ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] YOLO Algorithm on Oil Palm Tree Counting
Hi Loke; I know that You Only Look Once (YOLO) is a neural network approach that is not very good at detecting small objects in an image. This may be problematic depending on the scale of your plantation image mosiacs. The image.darknet package in R is a YOLO solver. It can be implemented in QGIS. Kirk Schmidt On 7/21/2023 1:17 AM, Loke Jian Munn via QGIS-User wrote: Hi Everyone, I have been recommended to try YOLO Algorithm on oil palm tree recognition and counting. (DSM and Orthomosaic ) Just to ask if someone have any insights of this algorithm and provide any reference or guidance ? Thank you. *Loke Jian Munn* Information Technology (IT) Department *KUALA LUMPUR KEPONG BERHAD* Wisma Taiko, No. 1, Jalan S.P. Seenivasagam, 3 Ipoh,Perak Darul Ridzuan,Malaysia *T : *+605 240 8000 (ext. 2443) A close up of a sign Description automatically generated <https://clicktime.symantec.com/15tpDKjABDCCmu4drBEZD?h=4-LzQeEUOcKmp7sCwswGafke694l_dZCDQJ4Rb8N8bc==http://www.klk.com.my/>*A close up of a sign Description automatically generated* <https://clicktime.symantec.com/15tpJ9vSdpsoBqtZPjdhq?h=bX-LFOj0HrSJV_t-3L_WhMP3i5PRx5AtNBaHlI02C7k==http://www.linkedin.com/company/klkberhad> DISCLAIMER : This email and any attachments are confidential and intended solely for the addressee. If you are not the intended recipient, any disclosure, copying, distribution either whole or partial, is prohibited and may be unlawful. In such case, please delete this message and kindly notify the sender. *Kuala Lumpur Kepong Berhad * ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Editing expression in QGIS
The short answer is yes. You will need all the formula variables within the feature table. Your syntax will change with ield names replaceing cell references. You can copy your final expression and save as a text file for future use. Alternatively, you can write a simple python script or use the graphical modeler to save your expressionKirk SchmidtSent from my Galaxy Original message From: Simon via QGIS-User Date: 2023-07-10 11:58 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Editing expression in QGIS Hi all, Is it possible to add an Excel equation into QGIS via the field calculator in the attribute table window. The expression is fairly large, and the references will have to be changed but is it even do-able? If not, is there a better tool or plugin I can use? The goal is to add the equation, but also to be able to edit it in the future if the data changes. Thanks, Simon ___ 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] Reclassify field values from txt rule file in QGIS
Hi Fabio: This is GRASS function and you can find an explanation at: https://www.ibiblio.org/pub/packages/gis/grass/grass63/manuals/html63_user/v.reclass.html Instead of a simple lookup table, some sql sytax is required. 9110 - (Hainsimsen-Buchenwald) 9110 you would require: 9110 where feature = '9110 - (Hainsimsen-Buchenwald)'or '9110 - (Hainsimsen-Buchenwald)(100.0%)' Note the single quotations for text. An easier approach would be top calculate a new field using the field calculator and use the string left expression - left(/landclass field/, 4) Kirk Schmidt On 6/28/2023 9:52 AM, fabio castelli via QGIS-User wrote: Dear Qgis users, I kindly ask for your support with an issue related to the reclassification of attribute tables with the v.reclass function. Here's the link to the question posted on Stackexchenge: https://gis.stackexchange.com/questions/462318/reclassify-field-values-from-txt-rule-file-in-qgis Thank you and best regards, Fabio Castelli <https://www.avast.com/sig-email?utm_medium=email_source=link_campaign=sig-email_content=webmail> Privo di virus.www.avast.com <https://www.avast.com/sig-email?utm_medium=email_source=link_campaign=sig-email_content=webmail> ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Palm Counting using QGIS
Hi LokeRichard suggested using a NDVI approach which would require data in the Near Infrared (NIR) band. An orthomisaic would contain 3 colour bands, but not the NIR. NDVI provides an understanding of plant greenness and health and may not be suitable for detecting oil palm's. If you are getting 90 percent plus with your current approach , you may find that you will spend 4 times longer capturing the additional 10 percent. NIR data,are available from various satellite sensors however, their resolution would not work at the individual plant level.Some drone sensors record the NIR band.Kirk Schmidt Sent from my Galaxy Original message From: Loke Jian Munn Date: 2023-06-06 6:18 a.m. (GMT-04:00) To: Richard McDonnell , kirk , qgis-user@lists.osgeo.org Cc: Moke Kwai Cheong , Cheong Chee Phang , Leong Kah Choon , Mulwanto , Endang Sudita Prayoga , M Rizki Ramdani , Fira Yunita , Qurata Ayun Wijanarko Subject: RE: [Qgis-user] Palm Counting using QGIS Hi Kirk & Richard, Thanks for the great inputs and we have a clearer direction now. Just to have a quick check with you, can the Imagery be Orthomosaic map ? Thank you. Best Regards, Loke Jian Munn Information Technology (IT) Department KUALA LUMPUR KEPONG BERHAD Wisma Taiko, No. 1, Jalan S.P. Seenivasagam, 3 Ipoh,Perak Darul Ridzuan,Malaysia T : +605 240 8000 (ext. 2443) From: Richard McDonnell Sent: Thursday, 1 June, 2023 10:55 PM To: kirk ; Loke Jian Munn ; qgis-user@lists.osgeo.org Cc: Moke Kwai Cheong ; Cheong Chee Phang ; Leong Kah Choon ; Mulwanto ; Endang Sudita Prayoga ; M Rizki Ramdani ; Fira Yunita ; Qurata Ayun Wijanarko Subject: RE: [Qgis-user] Palm Counting using QGIS Hi Loke, Just to add to Kirk’s comments, have you considered using Imagery along with the DSM? More specifically a NDVI of the plantation. You will need to determine the signature range for Palm Tree Leaves, but where you have both a Palm tree Signature and your Local Minima/Maxima it could improve your detection. Kind Regards, Richard —— Richard McDonnell MSc GIS, FME Certified Professional FRM Data Management —— Oifig na nOibreacha Poiblí Office of Public Works Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36 Jonathan Swift Street, Trim, Co Meath, C15 NX36 —— M +353 87 688 5964 T +353 46 942 2409 https://gov.ie/opw —— To send me files larger than 30MB, please use the link below https://clicktime.symantec.com/15uBY3eXTuVgDRHFkBtiB?h=f_sNL5U9oSfDKaodrCLg96KXSTHDaiQO0T0VIEDipgY==https://filetransfer.opw.ie/filedrop/richard.mcdonnell%40opw.ie Email Disclaimer: https://clicktime.symantec.com/15uBhi36P8rs3Jw6qJh1R?h=7BRo4J_IELubeqOZ5acMrrBwdwGb5BtTEWtxK9bd_8w==https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/ —— MSc GIS, FME Certified Professional —— Oifig na nOibreacha Poiblí Office of Public Works Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36 Jonathan Swift Street, Trim, Co Meath, C15 NX36 —— M +353 87 688 5964 T +353 46 942 2409 https://https://gov.ie/opw —— Email Disclaimer: https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/ From: QGIS-User On Behalf Of kirk via QGIS-User Sent: 01 June 2023 10:58 To: Loke Jian Munn ; qgis-user@lists.osgeo.org Cc: Moke Kwai Cheong ; Cheong Chee Phang ; Leong Kah Choon ; Mulwanto ; Endang Sudita Prayoga ; M Rizki Ramdani ; Fira Yunita ; Qurata Ayun Wijanarko Subject: Re: [Qgis-user] Palm Counting using QGIS Hi Loke. A couple of comments. 1) The 95 % success rate is very good. If you were to have a close look at the missed trees you may realize that the highest point on a particular palm tree is lower than the neighboring tree or palm frond. What ever algorithm you use, you are limited by the physical data, accuracy and resolution being the key factors. The physical condition of the trees you are trying to identify also plays a role,as mentioned previously. A DSM is an approximation as it is derived from some form of point data, either lidar or photo based and then generated by extracting the maximum, 95th percentile, maximum or perhaps the average height if all the points that land in a cell. Your issue is a limitation of the data and the basic logic of a max/minima algorithm. It cannot detect what is not reflected in the data. You can interpret the colour imagery but have a look at the dsm at the missed maxima. You could try developing a routine by processing the,raw point cloud data or try generating a dsm yourself. You can also try changing the algorithm parameters. Perhaps in areas if missed maxima. Although, this may not be effective as the maxima you are looking for probably is not reflected the dem. 2) Plug-ins,are created by third parties and not by the Qgis team directly. You can reach out to the author and discuss with him/her. Kirk Schmidt Sent
[Qgis-user] (no subject)
Hi LokeRichard suggested using a NDVI approach which would require data in the Near Infrared (NIR) band. An orthomisaic would contain 3 colour bands, but not the NIR. NDVI provides an understanding of plant greenness and health and may not be suitable for detecting oil palm's. If you are getting 90 percent plus with your current approach , you may find that you will spend 4 times longer capturing the additional 10 percent. NIR data,are available from various satellite sensors however, their resolution would not work at the individual plant level.Some drone sensors record the NIR band.Kirk SchmidtSent from my Galaxy___ 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 calculate the visible area of a viewshed?
Hi David.Here is a link to a tutorial on viewshed analysis;https://dges.carleton.ca/CUOSGwiki/index.php/Conducting_a_Viewshed_Analysis_in_QGISKirk SchmidtSent from my Galaxy Original message From: david cebrian via QGIS-User Date: 2023-06-03 12:51 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] How to calculate the visible area of a viewshed? Good morning guys, I need to calculate the visible area of a viewshed in order to do some statistics for a paper. However, I can't find the way to do it. Any help would be highly appreciated. Many thanks. All the best, David ___ 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] Palm Counting using QGIS
Hi Loke.A couple of comments. 1) The 95 % success rate is very good. If you were to have a close look at the missed trees you may realize that the highest point on a particular palm tree is lower than the neighboring tree or palm frond. What ever algorithm you use, you are limited by the physical data, accuracy and resolution being the key factors. The physical condition of the trees you are trying to identify also plays a role,as mentioned previously. A DSM is an approximation as it is derived from some form of point data, either lidar or photo based and then generated by extracting the maximum, 95th percentile, maximum or perhaps the average height if all the points that land in a cell. Your issue is a limitation of the data and the basic logic of a max/minima algorithm. It cannot detect what is not reflected in the data. You can interpret the colour imagery but have a look at the dsm at the missed maxima.You could try developing a routine by processing the,raw point cloud data or try generating a dsm yourself. You can also try changing the algorithm parameters. Perhaps in areas if missed maxima. Although, this may not be effective as the maxima you are looking for probably is not reflected the dem.2) Plug-ins,are created by third parties and not by the Qgis team directly. You can reach out to the author and discuss with him/her. Kirk SchmidtSent from my Galaxy Original message From: Loke Jian Munn via QGIS-User Date: 2023-06-01 6:28 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Cc: Moke Kwai Cheong , Cheong Chee Phang , Leong Kah Choon , Mulwanto , Endang Sudita Prayoga , M Rizki Ramdani , Fira Yunita , Qurata Ayun Wijanarko Subject: [Qgis-user] Palm Counting using QGIS Hi Team, We have been using QGIS to do Palm Counting. Palm Counting here mean oil palm, we use QGIS processing tool, Features-Raster Tools (Local Minima and Maxima) to do the counting on a DSM layer. Example as below: Based on our industrial experiences, the accuracy of the Local Maxima algorithm is approximately 90%-95%. However, we still need to do toggle editing such as we utilize the “ Add Point Feature “ function when there are points missed, the “Vertex Tool” is used to adjust the position of points and the “Select Feature” for removal of any unnecessary points. We would like to have an improvement plan on the accuracy for more than 95% to reduce the extra manual works ( toggle edit ). Appreciate if the team could share with us any latest QGIS palm count plugin or any other solutions to achieve the accuracy more than 95% when doing palm count. Your input is highly appreciated. Thank you. Regards, Loke Jian Munn Information Technology (IT) Department KUALA LUMPUR KEPONG BERHAD Wisma Taiko, No. 1, Jalan S.P. Seenivasagam, 3 Ipoh,Perak Darul Ridzuan,Malaysia T : +605 240 8000 (ext. 2443) DISCLAIMER : This email and any attachments are confidential and intended solely for the addressee. If you are not the intended recipient, any disclosure, copying, distribution either whole or partial, is prohibited and may be unlawful. In such case, please delete this message and kindly notify the sender. Kuala Lumpur Kepong Berhad DISCLAIMER : This email and any attachments are confidential and intended solely for the addressee. If you are not the intended recipient, any disclosure, copying, distribution either whole or partial, is prohibited and may be unlawful. In such case, please delete this message and kindly notify the sender. Kuala Lumpur Kepong Berhad ___ 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] 3D Moon models with QGIS
Hi MeganThanks for the interesting question and link. Yes QGIS gas the tools to complete this analysis.First generate a slope raster from your dem. You can use the slope tool in the Processing Toolbox>raster terrain analysis folder.The result will be a raster with the slopes of your study area. You can then use the raster calculator to convert any pixels with a slope greater than your threshold value to either a large number or the raster null value. The remaining cells will allow you to visualize a route between 2 points. See this post for more details.https://gis.stackexchange.com/questions/166104/using-dtm-slope-map-to-find-shortest-route-between-two-locations-using-qgisTo actually generste a route between 2 points , use a least cost path. There is a least cost path plug-in which will generate a vector lines to connect 2 points once you have the slope raster generated.KirkSent from my Galaxy Original message From: Megan Eskey via QGIS-User Date: 2023-05-05 7:16 p.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] 3D Moon models with QGIS I am curious about the capabilities of QGIS. I came across this interesting project to make a 3D map of the Moon, but now I would like to generate routes (rover traverses) based on slope data, specifically the low slope route between two point locations.Is that possible? Any tips or tricks would be helpful.https://www.ursa.fi/fileadmin/ursa2010/Harrastus/Ryhmatapaamiset/aktapaaminen/2021/Kuukuvat_QGIS_ropponen.pdfMeganᐧ ___ 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] Qchainage plugin
Hi Stewart: You could write a simple script using python or modeler to do what you are looking for. I would reproject to a UTM or similar projection and then the math for each segment gets easier. Reach out to me privately as I have some ideas on how to implement it. Kirk Schmidt k...@nortekresources.com On 3/30/2023 4:56 PM, Stewart Holt via QGIS-User wrote: Tim, Depending on how the GPS data was collected for the trail, your line shape may be longer than the actual trail. Each GPS point has some error and if the points were collected at one per second, the additional length caused by points being off the centerline, could be significant. GPS height is more noisy with error than the horizontal distance. That can contribute a lot of extra length. I use a sub-meter GNSS unit and filter the points to be no less than 10 feet apart. I add elevation from a USGS layer using the QGIS drape function because that is more accurate than the GPS height. By doing this, the 3D length is within about 1% as compared to a measuring wheel. I made the earlier post for work I do for the Appalachian Trail. What I want is a way to calculate the 3D distance from one point to the next and put the running total in the M value of each point. In my case, the layers are in geo coordinates and would require the horizontal distance be calculated with lat/lon. I read the descriptions of every plugin and determined that nothing currently does what I wanted. I am a programmer but have very little Python experience. I did get one response which was helpful in writing plugins but have not yet pursued it. An example of a plugin which takes a vector layer input, does something (or nothing) with it and creates an output layer would be a good starting point. If I can achieve what I want to do, what you need would be a simple change. Stewart Holt Georgia Appalachian Trail Club On Thu, Mar 30, 2023 at 3:10 PM Tim Fearn via QGIS-User wrote: Hi everyone, I’m trying to use the qchainage plugin to create ½ mile markers along a vector line that represents a trail. The line is a single feature with no multi-parts or dangles. It is about 318 miles long so I would expect right around 636 points if there is one every half mile. My project default is set to miles for length measurements as is the qgis main preference. I set the unit measurement to miles and input .5 as the interval. I continually get 771 points instead of the expected 636. Also when I use the automatic labeling feature of the plugin it appears to be in degrees not an aggregate mileage. Several weeks ago a fellow had queried the list on the same subject. He was working with the Appalachian Trail I believe and wanted to place points at given distances. No one seemed to really understand so I’m hoping that if anyone remembers the original poster, and has come up with any thoughts or other ideas it would be much appreciated. This is something that can be done for sure. All of the National Scenic Trail maps in this country are demarcated with some form of mile or half mile makers which when labeled show the accumulated distance. Hoping for some traction on this. I’ve only just started learning python so I wouldn’t describe myself as a programer and am looking for a relatively simple solution. With much gratitude, Tim Fearn ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Shape Tools New Capabilities - Compass cardinal directions
Nice Update Calvin. What is your 0 or 360 degree reference? I assume it is grid north if you use a projected coordinate system. Kirk Schmidt On 3/10/2023 8:04 AM, C Hamilton via QGIS-User wrote: I thought I would highlight a couple of updates to the Shape Tools plugin. For those who are not aware, Shape Tools has a measuring tool that not only includes the distance from point A to B, but also displays the azimuth from A to B and from B to A. It shows each segment with these features and the total distance. It also has a measurement layer algorithm that creates a new layer from a polygon or line layer with already styled measurements of the line segments. The recent addition is to also include the compass cardinal directions such as N, S, E, W or NE, WSW, or NEbE depending on the number of compass points that are selected. This has to be enabled from Vector->Shape Tools->Settings and then look for "Points of the compass" and you can select None, 4, 8, 16, or 32. It defaults to None. The other addition is a related expression function with some examples st_compass( azimuth[, number_of_points, mode] ) st_compass(33) → 'NNE' st_compass(33, 8) → 'NE' st_compass(33, 32) → 'NEbN' st_compass(33, 16,'full') → 'north-northeast' st_compass(33, 32,'full') → 'northeast by north' st_compass(33, 32,'trad') → 'Quarto di Greco verso Tramontana' Enjoy, Calvin ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] LIDAR
Hi Allan: One way is to combine the files into a super file which may not be an option due to size. The other would be to find the total range of z values in your project area and set the min and max symbology for the Z attribute. while using the attributes by ramp option. You can save the style and use it for the remaining point clouds. Kirk Schmidt On 2/13/2023 11:49 PM, Allan Younger via QGIS-User wrote: Hi everyone, I have multiple Lidar downloaded las files, on opening find each has its individual range of levels. How can I shade uniform levels? 200m plus points in several of blocks. Chrs AY ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Azimuth and orientation of 2 features
Hi Alex.You could use the coordinate of the centrists to determine which polygon is east or north. For example, if your coordinate system is utm. the centroid with the highe Easton will he the most easterly, the one with the smallest northing will be the most southerly.Kirk Schmidt Sent from my Galaxy Original message From: Alexandre Rodier via QGIS-User Date: 2023-02-12 1:57 p.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Azimuth and orientation of 2 features Dear list.I have a question I am not able to resolve:I need to calculate the orientation (North, East, South or North) of a feature (Polygon, a parcel) respect of another one (Line, a road).I calculated the azimuth of the two features using the centroids of each but the formulas I found don't really work, or I am doing something wrong (it's probably that :-) ).I am sure that's a simple problem but I am quite blocked.Thanks in advance for your help.Alex___ 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] CSV layer, Time column null expressie to test for NULL/empty?
You can also useCASEWhen field is null then 123else 321endKirk SchmidtSent from
my Galaxy
Original message From: Richard Duivenvoorde via QGIS-User
Date: 2023-01-18 1:44 p.m. (GMT-04:00) To:
qgis-user@lists.osgeo.org Subject: Re: [Qgis-user] CSV layer,
Time column null expressie to test for NULL/empty? Pff, ok found a
solution:"bezoek.begintijd" IS NULLworksSo if I want to test for empty strings,
or NULL (like strings), AND real NULLS, I use now:( "bezoek.begintijd" IS NULL
OR "bezoek.begintijd" in ('', ' ', 'NULL') )Regards,Richard DuivenvoordeOn
1/18/23 18:31, Richard Duivenvoorde via QGIS-User wrote:> We have csv data,
which is loaded as layer via the Delimited Text Layer provider.> > The data
contains both a date and a time column, which we have to format as isodatetime
using an QGIS Expression.> > BUT sometimes time column is empty/no value (in
which case we have to format it as isoTIME only).> > So we have to check IF the
time column is empty, then only format the Date, else use Date AND Time.> > My
problem is that I fail to create an expression for the test IF the column (see
screenshot) is empty.> > I thought to use (our column is called
"bezoek.begintijd"):> > "bezoek.begintijd" in ('', ' ', NULL)> > BUT this does
NOT return True or False, but NULL !!> > But:> > > ("bezoek.begintijd" in ('',
' ', NULL)) = NULL> does ALSO NOT resolve to True... (also NULL)> > When the
record has no data (the csv holds no data for that record), it shows> >
if it has data:> > > Anybody an idea what test I can use to
test if the time column (which is typed as Time) is empty/null/none?> > Thanks
for any pointers.> > Regards,> > Richard Duivenvoorde> >
___> 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 listQGIS-User@lists.osgeo.orgList info:
https://lists.osgeo.org/mailman/listinfo/qgis-userUnsubscribe:
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] Mapa de visibilidad
Hi Sara: It depends what you mean by automatically. You can create a static viewshed of one or several points and click on a location to see if any or all are visible. If you want to select a point and do the calculation on the fly, that is a different process and would take a bit of thinking on how to implement it. The Visibility Analysis plugin will give you the functionality I described first. Kirk Schmidt On 1/11/2023 10:06 AM, Sara Gómez Zaragozà via QGIS-User wrote: Buenos días, Mi nombre es Sara Gómez, y querría preguntarles sobre una cuestión relacionada con la aplicación del QGIS (versión 3.8 para MAC OS X Capitan 10.11.6). Para generar un mapa de cuencas visuales de determinados elementos del paisaje (en este caso, arte rupestre o yacimientos arqueológicos) es posible que el programa marcando un punto en el mapa nos pudiera marcar los elementos que se visualizan desde este punto de forma automática? Gracias, un saludo -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Airborne camera rectification
Hi Scott. A couple of technical points to consider:What is the accuracy of the imu on the camera? 0.1 degrees at 30 km results in a lot of error. Do you have the interior orientation of the camera? This is the lens distortion which makes up a critical part if the orthorectification process.The grass library has a number of ortho tools that may suite your purposes (i.ortho).Simple manual georeferencing may also meet your requirements.Feel free to reach out to me privately if you have any questions.Kirk Schmidt Nortek Resource Solutions Inc.kirk@ nortekresources.comSent from my Galaxy Original message From: "Madry, Scott via QGIS-User" Date: 2022-12-30 3:48 p.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Airborne camera rectification Hello list. First of all, thanks to all of you for helping make the QGIS community the very special thing that it is. I will be involved in a stratospheric balloon launch in January (altitude ~30 km), and it will have a nadir-pointing RGB camera that will have it’s images telemetered to the ground. It is nadir-pointing, but will have random swaying motion in all axes. Does anyone have experience in reading the 6 parameters of exterior orientation (X, Y, Z, w, f, k (Rotations about x, y, z axes respectively)) from the EXIF header attached to the JPEG and using it to georectify the images as well as can be done in QGIS? We are hoping to compare the various images at different altitudes with satellite imagery acquired at about the same time. Thank you so much, and happy new year to everyone! Scott Scott Madry, Ph.D. Research Associate Professor of Archaeology The University of North Carolina at Chapel Hill Tel 1-919-448-4493 Email: mad...@email.unc.edu https://scottmadry.web.unc.edu Skype: scott madry ___ 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] SQL query in a QGIS model
Good Day.Try dropping the @ as you are using field values and not variables. If the were variables, then would show up with @ in the list.Kirk SchmidtSent from my Galaxy Original message From: ALT SHN via Qgis-user Date: 2022-12-03 5:26 a.m. (GMT-04:00) To: "qgis-user-requ...@lists.osgeo.org" Subject: [Qgis-user] SQL query in a QGIS model Hello everyone,Does someone knows how to solve this:https://gis.stackexchange.com/questions/446636/how-to-use-an-input-parameter-as-part-of-a-sql-query-in-a-qgis-modelThank you-- ---Sociedade de História NaturalDepartamento de Informação GeográficaPolígono Industrial do Alto do AmialPav.H02 e h06i.geograf...@alt-shn.orgwww.shn.ptwww.alt-shn.blogspot.com FacebookTLM: 964138188 ___ 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] large DEM not loading tiles correctly
Perhaps you can consolidate small groups of tiles into single rasters. Kirk SchmidtSent from my Galaxy Original message From: "Michael.Dodd via Qgis-user" Date: 2022-09-28 1:32 p.m. (GMT+01:00) To: "'qgis-user@lists.osgeo.org'" Subject: [Qgis-user] large DEM not loading tiles correctly I am trying to put together the 50m uk DEM e.g. https://www.data.gov.uk/dataset/835cf20a-8feb-4394-8b30-dcfe840ac13d/os-terrain-50-dtm but it often does not work. The data are supplied as thousands of small tiles that takes ages to unpack and put into a folder to get into qgis, however when you try to load the tiles some work and others don’t and even the ones that do work and display fine sometimes don’t merge to form large tiles and sometimes will form virtual rasters but not merged files. Basically it is a mess. I suspect others have done this successfully so any suggestions? It is particularly odd that some of the tiles will merge and others won’t and that some won’t display if you ask it to load hundreds at a time but if you just add them in a few at a time then they will load and display (in addition to the hundreds already loaded). This suggests various bugs assuming the original data are not in a whole variety of slightly different formats (should all be .asc) ___ 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] Symbology
If the data is being read as a character versus a integer or decimal, you cannot use graduated symbols. Double left click on the layer in the legend and select properties. Select the Fields tab and see what format the attributes are being read as. I have run into this when generating a csv from a spreadsheet from different language. I would convert the format in excel or livre office and regenerate the csv file. Alternatively you can use text editor to replace the quotes on each value in the existing csv.Kirk SchmidtSent from my Galaxy Original message From: Allan Younger via Qgis-user Date: 2022-06-18 3:13 a.m. (GMT-04:00) To: "Kobben, Barend (UT-ITC)" , chris hermansen Cc: QGIS User Subject: Re: [Qgis-user] Symbology csv file -- Original Message -- From "Kobben, Barend (UT-ITC)" To "chris hermansen" ; "Allan Younger" Cc "QGIS User" Date 18/06/2022 12:32:07 PM Subject Re: [Qgis-user] Symbology Usually this happens if you have numbers (quantitatvie data) that are stored in your file/db as Strings (text) -- Barend Köbben On 18/06/2022, 04:25, "Qgis-user on behalf of chris hermansen via Qgis-user" wrote: Allan and list, On Fri, Jun 17, 2022 at 7:00 PM Allan Younger via Qgis-user wrote: Hi, I'm trying to show assay results as simple dot plots as graduated colour and size scaled points based on a selected element. Problem is that in the layer properties/symbology tab, I cannot select anything. Had the problem in Madeira 3.4 & assumed was just a glitch that had developed in that copy over time. Installed 3.22.7, same problem. Csv file has 162000 records & 30 column with co-ords for every point; with locations plotting a-ok. At the top of the Layer Properties > Symbology dialogue window, have you selected one of the options other than "No Symbols"? For example, "Graduated"? Then in the field Value one of your layer's attributes? Then pressed the "Classify" button? -- 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
Re: [Qgis-user] How to make joined fields editable?
Hi Roland I have not tried that checkbox bug just save your joined layer and you can edit the resultant file. Kirk SchmidtSent from my Galaxy Original message From: Roland Berger via Qgis-user Date: 2022-05-15 5:25 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] How to make joined fields editable? Hi, I have a point vector layer A and joined some fields from a table B. Now I would like to edit values in B by looking at the attribute table of A. In the 'Edit Vector Join' Dialog' of A I have ticked the checkbox 'Editable join layer'. If I now try to edit the attribute table of A the joined fields of B are shown but still not editable. What does this checkbox 'Editable join layer' in the 'Edit Vector Join' Dialog really do? What do I have to do to make joined fields editable? Thanks for any help. Roland Remarks: Layer A is in a Postgis database. Table B is in a Geopackage File (gpkg) -- ___ 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 convert depth belw ground level to depth below mean sea level in QGIS
Hi NavneetThis is classic raster analysis. You will need the corresponding Digital Elevation File (DEM). If your water table data are in vector fomat, you will have to rasterize the data, preferably to the same resolution as the dem. Then you can use the grid calculator to subtract the water tabke depth from the dem. The resulting raster will be your water table height above sea level or the zero reference for the dem. You can contour the raster after that if required.Kirk SchmidtSent from my Galaxy Original message From: nr via Qgis-user Date: 2022-05-03 5:08 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] How to convert depth belw ground level to depth below mean sea level in QGIS Dear friends,The data provided for depth of water table in my country by the CGWB is depth below ground level , but for drawing the water table contours in modflow software the depth of water below mean sea level is needed.So how the depth below ground level can be converted to depth with reference to mean sea level in QGISRegardsNavneetQGIS 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] Issues with SAGA's tools
Hi Lucas; Try installing SAGA as a stand alone and run the tool. I use this tool regularly and it works fine. Kirk Schmidt On 4/19/2022 8:58 AM, Lucas Pomares via Qgis-user wrote: Hello, I’m trying to use fill sinks tool (Wang & Liu) of SAGA in Qgis Toolbox. However, when the process ends the program fails and doesn´t generate the result layers. I´ve tried it using different versions of Qgis (3.16 , 3.20 and 3.24) but I couldn´t solve the problem. SAGA version is 7.8.2. I also tried to uninstall and reinstall the program and I didn´t achieve. I tried to solve it by myself but I don´t know where is the problem and would want to ask if anyone knows how to solve it. I attach a picture of the fail and the log file in txt. In English says that: The following layers doesn´t generate correctly. • C:/Users/Lucas/AppData/Local/Temp/processing_vRdlpU/1ee9a17ddce54df9967efcee956d7b95/FDIR.sdat • C:/Users/Lucas/AppData/Local/Temp/processing_vRdlpU/83bae0d547e5462290b9d6530fb03bf3/WSHED.sdat • C:/Users/Lucas/AppData/Local/Temp/processing_vRdlpU/a7ad2897f43e4c7bbfff4898a4a2cc07/FILLED.sdat You can look the << register message panel>> in main window of Qgis for more information about algorithm. Thank for the help!! *Lucas Pomares Zambrano* *Dpto. Ingeniería* GEOLÉN Advertencia legal: *GEOLEN INGENIERIA S.L.,*CON CIF B41954124 y domicilio social en C/ Séneca, 1 - Local 7 (Valencina de la Concepción) - Sevilla, le informa que es responsable del Tratamiento de datos de carácter personal, donde se incluyen sus datos, cuya finalidad es la Gestión de sus servicios, la relación administrativo-comercial o la relación con el personal a su cargo, que incluye o puede incluir el envío de documentación. Le informamos que podrá ejercer los derechos de acceso, rectificación, limitación del tratamiento y portabilidad de sus datos personales en la dirección indicada mediante escrito, concretando su solicitud y al que acompañe fotocopia de su Documento Nacional de Identidad. El envío de este email responde con la totalidad de la legislación vigente en materia de protección de datos de carácter personal y Ley de Servicios de la Sociedad de la Información. ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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 vector layer from slope map with few classes
Hi Maria.The symbology does not reclassify the underlying raster, it just changes the way it is displayed. Use the raster calculator to generate a new raster with the Classes you require. The native qgis calculator works well. You can find it in the upper word menu under raster. Once you have your raster reclassified, you can then convert it to a vector format Kirk SchmidtSent from my Galaxy Original message From: "Shinoto, PD Dr. Maria via Qgis-user" Date: 2022-03-29 6:32 p.m. (GMT-04:00) To: qgis forum Subject: [Qgis-user] creating vector layer from slope map with few classes Hi, I have created a slope map from a DTM. I have limited the symbology to three meaningful classes and wish to create a vector layer with polygons of three classes from that. The conditions are One-Channel-Pseudocolor, Channel1 (Gray), interpolation discrete, modus "fortlaufend" -- not identical interval or quantile, 3 classes with values <=20, 20 - 40, > 40. The classes show up in the layer design panel symbology tab, but the layer browser still (thankfully!) shows the greyscale ramp. I wish to create a vector layer with polygons from the 3 classes. All solutions I have looked up work under slightly different condition. After all I do not know how to get the data of the DTM into three classes like "Flat-Gentle-Steep"; I think the vectorization should work easily then. Can anybody hint me to a solution?Thanks and best wishes, Maria___Qgis-user mailing listQgis-user@lists.osgeo.orgList info: https://lists.osgeo.org/mailman/listinfo/qgis-userUnsubscribe: 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] shapefile import error
Hi David Just to Clarify Richards comments. A shape file is a collection of files, all with the same prefix. At a minimum, you should have 3 files, *.shp, *.shx and *.dbf. There may also be a few more files as well depending on how the file was created. You may just have the *.shp file which is creating the error. Kirk Schmidt On 3/3/2022 3:24 PM, Richard Greenwood via Qgis-user wrote: Do you have the "side car" files .shx and .dbf and possibly others? On Thu, Mar 3, 2022 at 9:31 AM David Witton wrote: Hello I'm pointing the Data source Manager|Vector at a .shp file and getting the error "Invalid Data Source". The file exists on an external drive, and the path is valid - any ideas on what could be causing this? -- David Witton 307 413 2072 ___ 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 -- Richard W. Greenwood www.greenwoodmap.com <http://www.greenwoodmap.com> ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Labeling Lines with Bearings
Hi Springfield: You can place the python code in the C:\Users\User\AppData\Roaming\QGIS\QGIS3\profiles\default\processing\scripts folder on your local machine and then it will show up under scripts in the processing toolbox. Kirk Schmidt On 2/10/2022 4:03 AM, Springfield Harrison wrote: Hello Everyone, Thanks very much for the suggestions. The main_angle option seem to do the trick. I didn't try the Python code. Not familiar with that side of things although it is probably simple enough. Could not find the recipe for adding it to the plug-in list, which I presume is where it is supposed to go. Thanks again, very helpful… Cheers, Spring On 04/Feb/2022 07:10, Kirk Schmidt wrote: Hi Springfield: Sure it would , I have attached my little convergence tool. Just create a point in the area where you want to know the convergence at let her rip. The vary last column is the convergence in minutes. So, depending on what side of the central meridian of your coordinate system you are, it will either be negative or positive. SO, I just used it yesterday for and area that is east of the meridian, UTM easting > 500,000, therefore the resulting convergence was 22 minutes, which need to be subtracted from the grid north value to obtain true north. It is not a sophisticated workflow, but I only require it intermittently. I usually use NRCAN magnetic declination calculator (online) to verify the magnetic declination for the day. But it seems to be offline. I convert everything to magnetic for field work and then there are no mistakes due to declination calculations and settings. Have fun Kirk On 2/4/2022 6:34 AM, Springfield Harrison wrote: Aha, thanks Kirk. That's why my lines are not as I intended. This probably affects making the lines with the Shape Tool. Yes, it seems to draw in grid north, angle_main reports in True DD. The map is 26907, UTM 7N. I calculated the difference between as drawn and as reported and used that to change the input figure. It works OK, not exact as the convergence varies with longitude, I believe. Thanks Kirk . . . . Cheers, Spring On 03/Feb/2022 02:59, kirk wrote: Sent from my Galaxy Original message From: kirk Date: 2022-02-03 5:55 a.m. (GMT-04:00) To: Springfield Harrison , qgis-user@lists.osgeo.org Subject: RE: [Qgis-user] Labeling Lines with Bearings Hi Spingfield. I use a little routine and the core if it is: main_angle($geometry) which you can use in field calculator. Note this will give you the azimuth in relation to grid north for the coordinate system you are using. For projections, you will need to lookup or calculate convergence to get the adjustment to true north. Kirk Schmidt Sent from my Galaxy Original message From: Springfield Harrison Date: 2022-02-02 7:30 p.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Labeling Lines with Bearings Hello, I need to put the line direction (azimuth) into a line label. Some searching has not come up with anything. Is there a way? Thanks very much . . . . Cheers, Spring ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] FW: RE: Labeling Lines with Bearings
Sent from my Galaxy Original message From: kirk Date: 2022-02-03 5:55 a.m. (GMT-04:00) To: Springfield Harrison , qgis-user@lists.osgeo.org Subject: RE: [Qgis-user] Labeling Lines with Bearings Hi Spingfield.I use a little routine and the core if it is:main_angle($geometry) which you can use in field calculator.Note this will give you the azimuth in relation to grid north for the coordinate system you are using. For projections, you will need to lookup or calculate convergence to get the adjustment to true north. Kirk SchmidtSent from my Galaxy Original message From: Springfield Harrison Date: 2022-02-02 7:30 p.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Labeling Lines with Bearings Hello,I need to put the line direction (azimuth) into a line label. Some searching has not come up with anything.Is there a way? Thanks very much . . . .Cheers, Spring___Qgis-user mailing listQgis-user@lists.osgeo.orgList info: https://lists.osgeo.org/mailman/listinfo/qgis-userUnsubscribe: 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] Labeling Lines with Bearings
Hi Spingfield.I use a little routine and the core if it is:main_angle($geometry) which you can use in field calculator.Note this will give you the azimuth in relation to grid north for the coordinate system you are using. For projections, you will need to lookup or calculate convergence to get the adjustment to true north. Kirk SchmidtSent from my Galaxy Original message From: Springfield Harrison Date: 2022-02-02 7:30 p.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Labeling Lines with Bearings Hello,I need to put the line direction (azimuth) into a line label. Some searching has not come up with anything.Is there a way? Thanks very much . . . .Cheers, Spring___Qgis-user mailing listQgis-user@lists.osgeo.orgList info: https://lists.osgeo.org/mailman/listinfo/qgis-userUnsubscribe: 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] Adding alpha band to raster
Hi Nicolas: I just had an idea, why not use gdal directly? I know gdalwarp has a -dstalpha switch to code null data Assuming this is what you want the alpha band for) Kirk On 1/21/2022 6:06 PM, Nicolas Cadieux wrote: Hi, I don't see it in the current processing algorithms. I see lot's or r.out... but not that one:( I guess I could do it in Grass... Thanks for the idea! Nicolas On 2022-01-21 1:37 p.m., Kirk Schmidt wrote: Hi Nicholas: Try r.out.gdal https://grass.osgeo.org/grass78/manuals/r.out.gdal.html Kirk Schmidt On 1/21/2022 2:25 PM, Nicolas Cadieux wrote: Hi, What is the easiest way, in QGIS3, to add an alpha band to a raster? Alpha bands are handy with the georeferencer. Currently, I am just adding it in GIMP before I start the georeferencer. Thanks -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Adding alpha band to raster
Hi Nicholas: Try r.out.gdal https://grass.osgeo.org/grass78/manuals/r.out.gdal.html Kirk Schmidt On 1/21/2022 2:25 PM, Nicolas Cadieux wrote: Hi, What is the easiest way, in QGIS3, to add an alpha band to a raster? Alpha bands are handy with the georeferencer. Currently, I am just adding it in GIMP before I start the georeferencer. Thanks -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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 Modeling, using an model input layer name
Hi Nicolas The field calculator will not process raster data. Try using the raster calculator Kirk Schmidt On 1/7/2022 2:56 PM, Nicolas Cadieux wrote: Hi, I can figure this on out!!! I want to use the layer name of the input layer (a raster) as value for the field calculator. layer_property(@Inputraster , 'name') does not work. Any suggestions? Thanks -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] SAGA Issues
Hi All: My short term work around is to run simple batch files (windows environment) to handle the saga processing. Not very elegant, but it gets the job done. Be aware that different SAGA versions may have different syntax. Kirk Schmidt On 1/4/2022 3:48 PM, i...@jcis.net.au wrote: I am certainly not an expert, but the SAGA issue has been discussed by the redoubtable Hans van der Kwast in a video in the QGIS Open Day session on Processing Provider plugins: https://youtu.be/9fWQZqPjPVk . I believe the issue lies with more recent versions of QGIS using a more recent version of SAGA but SAGA not being integrated into QGIS in the way it previously was. All sorts of QGIS Gurus are discussing this and at some point the clouds will open and a stone tablet will be issued. It would be really great to have a fix for tis issue as the SAGA tools are really useful. Cheers *Dr Iain Stuart* */JCIS Consultants /* P.O. Box 2397 Burwood North NSW, 2134 (02) 9701 0191 (0413) 380116 (m) ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Create new coverage
Hi Paddy: The easiest approach would be to generate a new grid using the extents of your existing grid an defined at the 5 degree resolution. Select Vector> Research Tools>Create Grid and follow the prompts. Kirk Schmidt On 11/3/2021 10:54 AM, Pat Brown wrote: Hi, I have a shapefile which consists of a grid (polygons), each of which is 1 quarter degree in extent. I would like to create a second coverage where each polygon is again subdivided into smaller polygons, each being 5 degrees by 5 degrees in extent. Essentially each degree square would consist of 144 polygons. Appreciate any suggestions, thanks, Paddy ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] Batch Aerial Imagery GeoReferencing
Hi Richard: It sounds like you are in the world of manually georeferencing or ortho creation. This would entail identifying ground control points for each image and processing. Assuming you do not need all 6000 images at one time, you could tag each image to its central location and then ortho or geo it if it is needed. Kirk On 11/3/2021 7:45 AM, John Layt wrote: Hi Richard, The Historic England aerial photography team have been working on a similar sort of project which they recently released (see https://historicengland.org.uk/whats-new/news/new-virtual-aerial-map-allows-everyone-to-explore-englands-archaeology-from-the-air/). Perhaps they might be able to provide some advice on how they did it? There's an email address at https://historicengland.org.uk/research/methods/airborne-remote-sensing/. Cheers! John. On Wed, 3 Nov 2021 at 10:27, Richard McDonnell wrote: Hi Kirk and Charles, Thanks for your reply. This is what we know… CAMERA: Wild RC8 and RC10 152mm Focal Length. ALTITUDE: 4,560m with Block Adjustments for Mountainous areas IMAGE OVERLAP: 60% along track 25% side lap and all flight line Orientations East – West except for offshore. Im afraid though that GPS just wasn’t a thing when these images where captured. We have been looking at reverse geotagging the Images, that is we know there location…ish and we can convert them to jpegs with EXIF values assigned in theory. Though the biggest issue is, in talking on the guys on the Open Drone Mapper forum, our overlaps aren’t sufficient for SFM, so an automated process will be out of the question. Once again, thanks for getting back, if you can think of anything else, please let me know. Regards, Richard —— Richard McDonnell MSc GIS, FME Certified Professional /FRM Data Management/ —— Oifig na nOibreacha Poiblí Office of Public Works Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36 Jonathan Swift Street, Trim, Co Meath, C15 NX36 —— M +353 87 688 5964 T +353 46 942 2409 https://gov.ie/opw —— To send me files larger than 30MB, please use the link below https://filetransfer.opw.ie/filedrop/richard.mcdonn...@opw.ie Email Disclaimer: https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/ *From:*Qgis-user *On Behalf Of *Kirk Schmidt *Sent:* 01 November 2021 12:59 *To:* qgis-user@lists.osgeo.org *Subject:* Re: [Qgis-user] Batch Aerial Imagery GeoReferencing Hi Richard: You would require a couple of things: 1) Camera interior orientation which would include the lens calibration 2) exterior orientation which would include precise coordinate at time of exposure as well as the pitch, roll and yaw of the camera. 3) the offsets from the camera film or sensor center to the center of the aircraft gps antenna. 4) an accurate DEM or your target area If you have the above data, you can autoreference your images. If the images have a high endlap and sidelap, you could use commercial software such as Pix4D or Global Mapper to orthorectify the images based on structure from motion. Check into drone orthorectification software. Kirk On 10/29/2021 9:13 AM, Richard McDonnell wrote: Afternoon all, I have access to over 6000 Aerial Images covering Ireland, we also have scans of the Flight Paths and Image Orientation and Extents (Scans from Paper Maps) What I was wondering is if there is a way to Batch Georeference all of these Images? Any Ideas, help, guidance welcome. I have thought about maybe even reverse engineering the Flight Plan, then using OpenDroneMap to georeference the Images. Regards & Thanks, Richard McDonnell * ——* *Richard McDonnell MSc GIS, FME Certified Professional* /FRM Data Management/ —— *Oifig na nOibreacha Poiblí* Office of Public Works *Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36 * Jonathan Swift Street, Trim, Co Meath, C15 NX36 —— M +353 87 688 5964 T +353 46 942 2409 https://gov.ie/opw <https://gov.ie/opw> —— To send me files larger than 30MB, please use the link below https://filetransfer.opw.ie/filedrop/richard.mcdonn...@opw.ie <https://filetransfer.opw.ie/filedrop/richard.mcdonn...@opw.ie> Email Disclaimer: https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/ <https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/> ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info:https://lists.osgeo.org/mail
Re: [Qgis-user] Batch Aerial Imagery GeoReferencing
Hi Richard: You would require a couple of things: 1) Camera interior orientation which would include the lens calibration 2) exterior orientation which would include precise coordinate at time of exposure as well as the pitch, roll and yaw of the camera. 3) the offsets from the camera film or sensor center to the center of the aircraft gps antenna. 4) an accurate DEM or your target area If you have the above data, you can autoreference your images. If the images have a high endlap and sidelap, you could use commercial software such as Pix4D or Global Mapper to orthorectify the images based on structure from motion. Check into drone orthorectification software. Kirk On 10/29/2021 9:13 AM, Richard McDonnell wrote: Afternoon all, I have access to over 6000 Aerial Images covering Ireland, we also have scans of the Flight Paths and Image Orientation and Extents (Scans from Paper Maps) What I was wondering is if there is a way to Batch Georeference all of these Images? Any Ideas, help, guidance welcome. I have thought about maybe even reverse engineering the Flight Plan, then using OpenDroneMap to georeference the Images. Regards & Thanks, Richard McDonnell —— Richard McDonnell MSc GIS, FME Certified Professional /FRM Data Management/ —— Oifig na nOibreacha Poiblí Office of Public Works Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36 Jonathan Swift Street, Trim, Co Meath, C15 NX36 —— M +353 87 688 5964 T +353 46 942 2409 https://gov.ie/opw —— To send me files larger than 30MB, please use the link below https://filetransfer.opw.ie/filedrop/richard.mcdonn...@opw.ie Email Disclaimer: https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/ ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ 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] georeferencing and displaying lat long.
Hi Byron: So if you want to great a grid/graticule , the simplest way is to set the map view to your desired coordinate reference system (EPSG:4326 for example). If you have some data that you can add to the map canvas that covers you can add it to use as the extent of the new graticule. You can set the extent of the grid from either; calcutate from an existing layer, using the current map canvas extent, or you can draw a box in the canvas. To add the graticule, select Vector>Research Tools> Create Grid and fill out the parameters. Kirk Schmidt On 10/1/2021 2:30 PM, Firstname Lastname wrote: i have several maps to georeference that have reference grid in lat/long. can someone please tell me how to display the coordinate grid in lat long on the map so that i get some nice grid lines for reference. actually, knowing how to post the reference lines in any coordinate system would be helpful. -- Byron Veilleux, MSc. P.Geo Conjugate Geologic Services Limited Calgary, Alberta Canada by...@conjugategeo.com <mailto:by...@conjugategeo.com> Cell:4037108414 ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] ordinary kriging parameters
Hi Maria: Just from the cuff, the units are pixels or meters. Try reprojecting your source data to a projected coordinate system (UTM, for example) and try kriging again. Set the raster resolution to a course value that makes sense for your data to run a test raster and then run at the final resolution you require if you are concerned with processing time. By using geographic corrdinates, you probably were generating a hugh raster which caused your computer to hang. Kirk Schmidt On 9/30/2021 3:03 PM, Maria Eugenia Gomez wrote: Hi, I am using QGIS 3.16.9 Hannover. I want to create a surface (with kriging) from a list of points whose SRC is WGS84. They are in decimal degrees covering a surface of 10km x 9km. I can see the points correctly- My problem comes with the "Ordinary Kriging" setup. I unchecked "block kriging" and used the following options: quality measure: stand. dev. Logarithmic transf. variogram model: linear fit:nodes search range: local search directions: all directions number of points: max. number of points But when I have to fill "cell size", "maximum search distance, block size.. etc", I do not understand the units and I think this is the problem. Which type of units should I use? I tried to use decimal degrees units but then the number of iterations was so high that my computer hangs. Can somebody help me with this? I searched on the QGIS documentation but it does not say anything about units. Thanks in advance, Maria ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. 26 Church Road Sutherlands River, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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 on a mobile platform
Josh: Another option is to use Avenza Maps (commercial) to view your georeferenced pdf's directly on a phone or tablet. This makes a simple to use mobil solution that does not require an intenet connection Kirk On 9/10/2021 2:02 PM, Greg Troxel wrote: If you just want something that will show a georeferenced and your current look at viking and gpsd. I am not sure if it does georerferenced pdf, but am 99% sure it will do geotiff. I use it with openstreetmap TMS for looking at tracks. I don't know if any of that works on windows, but with cygwin/X I would expect it would. Maybe someone knows of a less-posixy approach. You probably could just run qgis and use live gps, but that may be too heavy for your win7-32 box. Your subject said "mobile" which usually means phone/tablet. Not what you asked, but there are various qgis-world apps: input, qfield, SMASH that might do you what want, and Trail Sense on android, which can just about do georeferenced pdf, if the pdf is in web mercator, but I'm having a hard time with qgis-produced maps. Also not what you asked, but if anyone knows of a command-line tool to just print out the georeferencing information from a pdf and check it for conformance, please let me know, and it could be useful to Josh in debugging next steps too. Greg ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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 can I export Vector & Raster layers together to import into mapshaper
Hi Kris; I am not familiar with MapShaper, however, the vector data and the raster layer (Openstreetmap) are fundamentally different animals. You could: 1) generate a geopackage which can store both raster and vector data types, but I do not know if you can read such a database with mapshaper? 2) create a map using print composer and export the vector and raster components as a georeferenced image or pdf file. This export a pretty map but you would lose all the feature attributes of your vector layer. Kirk Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com On 8/12/2021 10:40 AM, Kris Lagocki wrote: My issue seems to be with the raster layer, can't seem to figure out how to have it export together with the vector layer. On Thu, Aug 12, 2021 at 5:27 AM Richard McDonnell mailto:richard.mcdonn...@opw.ie>> wrote: Hi, There is a TopoJSON Writer Plugin available for QGIS, you might give that a try. It should write the vector dataset out as TopoKSON for you. Regards, Richard —— Richard McDonnell MSc GIS, FME Certified Professional /FRM Data Management/ —— Oifig na nOibreacha Poiblí Office of Public Works Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36 Jonathan Swift Street, Trim, Co Meath, C15 NX36 —— M +353 87 688 5964 T +353 46 942 2409 https://gov.ie/opw <https://gov.ie/opw> —— To send me files larger than 30MB, please use the link below https://filetransfer.opw.ie/filedrop/richard.mcdonn...@opw.ie <https://filetransfer.opw.ie/filedrop/richard.mcdonn...@opw.ie> Email Disclaimer: https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/ <https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/> *From:*Qgis-user mailto:qgis-user-boun...@lists.osgeo.org>> *On Behalf Of *Kris Lagocki *Sent:* 11 August 2021 14:32 *To:* qgis-user@lists.osgeo.org <mailto:qgis-user@lists.osgeo.org> *Subject:* [Qgis-user] How can I export Vector & Raster layers together to import into mapshaper Hi - I created a shape file using QGIS with a US Postal Code shape file that I found on the Census website and joined it with a data file listing the Postal Codes and which Sales Territory they belong to. I added in a raster layer OpenStreetMap. I'd like to export both layers together and import them into mapshaper so that I can export as Topojson. It has to be in this format to use within Power BI. Has anyone accomplished this? Any help is appreciated. 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Clipping WMS Layer for Export to Geo PDF
Hi StewartI assume you are using a rater as a background layer and therefore you will always end up with a rectangle as your final extent. If you clip to an irregular polygon, you will still have a rectangle but the outer area that had been clipped out will be transparent. This is how raters work as they are regular grids. The disadvantage of clipping to an irregular polygon is that you will create an additional alpha channel which defines the transparency but adds to the file size.I would consider making a number of smaller maps with some overlap to use with avenza.Kirk SchmidtSent from my Galaxy Original message From: Stewart Holt Date: 2021-07-30 4:38 p.m. (GMT-04:00) To: User List QGIS Subject: [Qgis-user] Clipping WMS Layer for Export to Geo PDF I am working on a GeoPDF map to use in the Avenza map app. It is an 80 mile hiking trail which stair steps from SW to NE. I am using a WMTS layer as a background. I want to avoid having the PDF as large as a rectangle containing this trail. I created a polygon and used it to select the extent for the map, but it did not clip the map but instead it created a rectangular extent containing the polygon.Next I created a larger rectangular polygon covering the area and cut out a ring for the trail and included this layer in the map. The result was the inclusion of all of the WMTS map layer as before with my mask layer able to block unwanted areas. This does not reduce file size, ofcourse.I think what that did is the correct behavior. What I need is a way to clip the WMTS layer on the fly. I think having the export to PDF being able to do this would be the best approach.Can anyone think of how to accomplish this without exporting the WMTS layer and then clipping it with a polygon?Thanks,Stewart ___ 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] Building Historical Maps
Hi JamesIt may be more efficient to digitize directly from scanned and georeferenced historical maps to segment the current data. For example lakeshore can be copied over from the current digital data do you do not have to digitize Georgian Bay and its many islands. As you go back in time, maps tend to be less accurate so you may run into problems. If you use township lines of baselines (ie. the Niven Line in Northeastern Ontario) as you key features to georeference your historical maps Yolo, you might have better luck getting everything to fit on the modern database. Georeferncing is a straight forward procedure. Kirk SchmidtSent from my Galaxy Original message From: James Cobban Date: 2021-07-22 10:16 p.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Building Historical Maps I am running a project to present historical information about Ontario, Canada. To support this I have imported GIS information on administrative boundaries from the Ontario Government. The only problem is this information only describes the current administrative organization of the province, not the organization it had for most of its history up until the early 1990s when the system was "rationalized". I am looking for suggestions of how to extract as much information from the current GIS data as possible.For example: From 1850 to the early 1990s most of Ontario was administered as counties with lower level townships. In the 1990s the more urbanized portions of the province were reorganized as "Regional Municipalities" which largely ignored the old county boundaries either by including only the urban portion of a former county, or by incorporating urban portions of multiple counties. So, for example, if I wish to reconstruct the boundaries of the former county of Welland I have to concatenate the boundaries of the individual townships that made up the former county, and then remove all of the internal, and redundant, boundaries. I have manually done that for one county but it was a pain in part because the polygons of the individual townships did not begin at a "corner" shared with an adjacent township within the county. Just finding exactly where the polygon began and which direction it went took time.The second issue is that boundaries changed multiple times over the period between 1763 when Britain assumed control of the area of Ontario, and the present. In particular city boundaries change as areas are transferred from the surrounding rural municipalities to the city as a result of urban growth. Where do I get the old boundary information aside from tracing them from old maps?-- James Cobban911-500 Springbank DrLondon, ON, CanadaN6J 4G6mailto: webmas...@jamescobban.net+1-226-504-7603___Qgis-user mailing listQgis-user@lists.osgeo.orgList info: https://lists.osgeo.org/mailman/listinfo/qgis-userUnsubscribe: 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] Did scale change outputting to PDF?
Hi All: I always add a scale bar to a printed map so that if there is variaition between the pdf/print configuration, It can always be confirmed on the final printed map. Kirk Schmidt On 5/27/2021 4:27 AM, j.hu...@post-ist-da.de wrote: Hi John, as I understand it, you created the base map in QGIS. If you use the measure tool in QGIS to get the distance between two distinct features in the map (e.g. road intersections) and then measure the same distance on the printed base map with a ruler, it should be possible to calculate the scale. Maybe use two distances, one aligned more or less horizontally and one vertically, to check if the scaling is proportional. As Andreas pointed out, it is probably a problem with the settings when the PDF was printed. In my experience it is a good idea to go to print shops usually working for architects and engineers since they are familiar with the importance of scaling (for advertising etc. it is more important that the whole content is printed, so that scaling might be used to fit the output to the printable area without potential cropping). You can print directly to a plotter in QGIS if you have access to the device, avoiding the PDF detour. EPSG 2264 should be fine. Units should be US feet. Regards Jochen Am 27.05.21 um 07:15 schrieb John Antkowiak: Hi. This plan was too simple to fail - but it failed. The charity whose project this is needed a large (that is... massive) paper wall map on which to plot and rethink its delivery driver assignments. Both drivers and delivery addresses are subject to change from week to week but it's not a pizza delivery; this is a regular run to supply people in a bad way. So the plan was to print the base map (roads and road names and county boundaries only) and then print 8.5 x 11 address maps with parcel data and orthos. That way, the base maps don't change but the physical parcel layer is flexible. (On top of that is a third paper layer indicating which drivers go where so someone can stand back and take in the whole picture graphically. Not a cutting-edge state of the digital art solution, but not everyone is cut out for that. It is what it is.) In order for this to work, the parcel maps have to be the same scale as the base map. Which they were... in QGIS. We have to convert all the maps to PDF to print them, and we had to send the base map PDFs to FedEx/Kinkos to print the 9 map grid panels at 42" by 62" each. When we got the big base maps up on the wall, we discovered the scale did not match the 8.5" x 11" parcel maps output to PDF and printed from home. It's not off by a lot, but it's enough to be painfully obvious from a single standard size sheet of paper. I don't know how to reverse engineer the big map scale precisely enough to enter a new scale number in the QGIS Print Layout. I didn't foresee it because this never would've been a conceivable scenario at the engineering firm where I picked up my meager GIS skills. (ArcMap sent a map directly to the plotter without interim steps.) There was no scale bar on the map. It shouldn't have been needed for this. Did something happen to the map scale when QGIS output the map to PDF? Could the size of the image on the pdf page have been adjusted manually or otherwise when being sent to a plotter with 42" paper? Could the image have been distorted horizontally differently from vertically? For the life of me, I cannot trial-and-error guess at a scale to enter. I've gone through dozens of new 8.5" x 11" test maps trying to guess the correct scale. Any ideas? Thank you all - John A. ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Aerial photo dates in legend without RGB patches?
Hi Michael: Just to be clear, these are your own georeferenced aerial photo's or professionally produced ortho photo's (as opposed to satellite imagery). If so, where would the date field come from, the name of the image files or perhaps embedded in the geotiff header (if using this file type)? Kirk Schmidt On 5/25/2021 10:22 AM, Michael Dufty wrote: Has anyone got a good solution to display aerial photo dates in the legend? I can change the name of the layer, but then it displays a huge useless R G B legend. Is there any way to turn this off? Setting the patch size to 0 doesn’t seem to be possible for local images, although it does work for WMS layers. Alternatively is there a way to just put a manual line of text in the legend, other than with dummy layers/groups? *Michael Dufty* ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Model Designer - Updating Field in Place
Hi Ryan; Yes, you would run a selection on then table first, apply the updates and save. If you require interactive selection criteria, you can use a number or string in input in your model. The user will then define them when it is run. Look under the vector selection Algorithms for select and/or extract functions in the process modeler Kirk Schmidt On 5/17/2021 11:27 AM, ryanpeel wrote: Perfect! It worked! And would you know how to set up model designer to have the model input to choose only selected rows versus whole table? -- Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-User-f4125267.html ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Model Designer - Updating Field in Place
Hi Ryan; There is a network plugin which may solve your issue. see https://gis.stackexchange.com/questions/279650/processing-modeler-to-update-table Kirk Schmidt On 5/17/2021 10:50 AM, Ryan Peel wrote: Hello, I have created a new Model in Model Designer and I'm using a Point Layer as Input and a Field Calculator that calls a custom function. I would like the output/result of the field calculator to update an existing field in the attribute table. I cannot get this to work. I can only get it to output a new temporary scratch table. Is my idea possible? -- Sincerely, Ryan ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] validatore topologico Qgis3
Hi Elisa: This may be best answered by the plugin developer. I am not familiar with that plugin but you can use the grass module v.clean or in process modeler to fix topological errors. Kirk Schmidt On 5/12/2021 12:45 PM, Elisa Franchi wrote: Buongiorno a tutti, ho notato che il plugin validatore topologico Qgis3 non rileva correttamente gli errori: segnala errori che sulla versione Qgis2 non ci sono e non ne segnala altri che in Qgis2 vengono rilevati. Ci sono impostazioni particolari da fare rispetto alla versione di Qgis2? Elisa -- ** *Elisa Franchi* *Direzione Tecnica* Email: *e.fran...@nemeasistemi.com* Mobile: +39 377 4314224 Skype: e.franchi1 *NeMea Sistemi Srl* *www.nemeasistemi.com <https://www.nemeasistemi.com>* *Sede Legale:* Via Mazzini 78, 07041, Alghero (SS) - Sardegna Tel: 079 4801861 *Sede Operativa:* Via Palazzo 96, 18038, Sanremo (IM) - Liguria Tel: 0184 872105 P.IVA: 02897110041 PEC: nemea.sistemi@pcert.it CUF: 66OZKW1 *NOTE LEGALI:* Le informazioni contenute in questo messaggio di posta elettronica e relativi allegati sono riservati e confidenziali. Qualora Lei non fosse la persona cui il presente messaggio è destinato, La invitiamo ad eliminarlo e a darcene gentile e tempestiva comunicazione. La diffusione e riproduzione in qualunque modo eseguita è severamente vietata ai sensi della normativa in vigore per la tutela della privacy. *EN*: This mail and/or attachments are confidential and may also be legally privileged. If you are not the intended recipient, you are hereby notified that any review, dissemination, distribution or copying of this e-mail and/or attachments is strictly prohibited. Please notify us immediately by e-mail and delete this message and all its attachments. *Rispetta l'ambiente. Non stampare questa mail se non è necessario.* ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Kriging using Qgis 3.18.2
Hi Jean:Export your excel file as a csv file. Then you can import the file and define the x and y fields. Then open processing toolbox and search for kriging. Select the tool and file out ghd parameters.Kirk SchmidtSent from my Galaxy Original message From: jean Lukusa Date: 2021-04-29 7:56 p.m. (GMT-04:00) To: Qgis Users List Subject: [Qgis-user] Kriging using Qgis 3.18.2 Hi all.I have an excel sheet with xy coordinates and Au grades.So I want to create Au grade contour using Kriging method.How can I do it?Jean LukusaB.Sc. GeologyTrained on Jorc CodeGIS Operator +243812156128 ___ 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] (no subject) (elevation question)
Hi Nicolas.You can set the barometer my using a known point and elevation or by setting the current atmospheric pressure. As with an anologue altimeter, your readings will change throughout the day as the pressure changes. This is less of a problem if you are out on a sunny day with a high pressure cell with a small pressure gradient. but becomes bit more problematic if a low pressure system is moving in.Sent from my Galaxy Original message From: Nicolas Cadieux Date: 2021-04-21 9:14 a.m. (GMT-04:00) To: David Strip Cc: Qgis Users List Subject: Re: [Qgis-user] (no subject) (elevation question) Hi,Interesting but I very much doubt they are not applying a geoid on top of the ellipsoid. If they can put maps of the earth on the device, they can put a geoid with a posting of 15´x15´! I have been wondering about those barometric altimeters GPS for back packing. I have a none electronic altimeter with me normally but it only really good if you consider the temperature changes. Nicolas Cadieuxhttps://gitlab.com/njacadieuxLe 20 avr. 2021 à 23:03, David Strip a écrit : In the course of trying to find what geoid model Garmin uses, I found this on a Garmin support page: . It is not uncommon for satellite heights to be off from map elevations by +/- 400 ft. Use these values with caution when navigating. Apparently some, but not all, models of the Garmin64 have a barometric altimeter. Depending on your actual model, this may or may not apply. On 4/20/2021 8:23 PM, Nicolas Cadieux wrote: Hi, David gave you a very good answer. Your looking at a 11 meter difference. That is not much. The first thing would be to look at the gps measurement. (Any trees? Good satellite constellation? Any obstacles?). At best, your looking at a 6 to 10 meter error vertically. SRTMs are affected by trees. That would make the SRTM higher (not your case). Last thing to look at the the differences between the gps geiod model (a model of the average sea level). SRTM use EGM96. Your GPS probably uses something a bit more modern. You can easily get a 5 meter error between two Geiods. ___ 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] (no subject)
Hi jean: Another issue is that you can set your elevation reference point in most Garmin units to either use the height above the ellipsoid or and approximation mean sea level. This setting will also have an impact on your observed elevation. Kirk Schmidt On 4/20/2021 8:45 PM, jean Lukusa wrote: Hi all . I hope you are doing fine. I have experienced an issue: The elevation taken in the field with gps garmin64 is 623 m but the one given by SRTM download from Qgis is 612m. How can I interprete this gap knowing that SRTM resolution is 30 meters. Jean Lukusa B.Sc. Geology Trained on Jorc Code GIS Operator +243812156128 ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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-3 single line/point layers from hand-written co-ordinates
Hi John: Yes, life is too short to waste time. I would add the lat/long into a spreadsheet or simple ascii text file. This way, you can double check your entries and once you are ready, you can bring them into QGIS as a *.txt or *.csv file. A simple text format would be : 1,45.3741, -62.1548 where 1 is the point id and the next is the longitude and latitude. These values are separated by a comma, however a semi colon or tab would work as well. You can also reverse the lat/long if it is easier for you. Instructions on how to import the file and convert it to a spatila layer are found at: (I know, another tutorial) https://www.qgistutorials.com/en/docs/3/importing_spreadsheets_csv.html If your coordinates are in degrees/minutes, you can use them as well. Reach out to me and I can give you the correct syntax if you need it. Kirk Schmidt On 4/19/2021 11:36 AM, John Moyle wrote: I am a retired Chartered Engineer, Physician and Anaesthetist!. In retirement I have completed a PhD "The Maintenance of the Victorian Submarine Telegraph Cable System". I want to correlate N Atlantic cable interruptions with geological events of that era. I have found a raster-file of the sea bed and a shape file of the tectonic plate boundaries. Fantastic! Now I want to make separate point or line layers for each interruption, volcano, earthquake and landslide. I have all the data in lat/long, hand written in a notebook whilst searching the archives at the National Maritime Museum. My problem is that all QGIS-3 courses take many hours turning one into a town and country planner before the simple layers I want to make! At 75, life is too short. I want to input the co-ordinates, by hand, layer by layer. Can anyone give me step by step instructions please. Thank you John Moyle ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Help! QGIS crashed and has deleted my QGZ file
Hi John: I am sorry to hear about your crash. It is important to note that QGIS provides a means of querying, updating, analysing and displaying data that is stored outside of the software (in a variety of formats). SO, I would suggest looking for the actual spatial files(s) that you were editing when it crashed. I assume they still exist, however, depending on the data format, you probably would have lost the attributes you were editing at the time of the crash. The QGIS project file may also have been corrupted at the time of the crash. You may have to use a backup to restore your last saved version. I hope this helps Kirk Schmidt On 4/19/2021 9:56 AM, j...@john-firth-editor.co.uk wrote: I’m using QGIS 3.16.0 on Windows 10, standalone PC. I have taken some training, but am a pretty inexpert user. The map that I’ve spent several months developing (for a book I’m writing) crashed on Saturday, and when I tried to load it just now it is greyed out in the Open Recent list (with an exclamation mark next to it) and File Explorer no longer finds it on my hard disk. I attempted to do what your website advised (downloaded DeBug View, which did not appear to run on my computer). When the file crashed, I had just added a couple of entries to an attribute table (it’s not a large one, though there are others in the file). I’ve no doubt someone more expert than I am could tell me how to extract a crash report, if that would be helpful. Is there any chance of saving my map? John Firth Advanced Professional Member of the Chartered Institute of Editing and Proofreading On 1 March 2020 a Royal Charter was given to the Society for Editors and Proofreaders (SfEP), which is now the *Chartered Institute of Editing and Proofreading*. I would like to keep your contact details in my records, so I can stay in touch. If you are happy for me to do this, please reply to this email, copying the following text into your reply: ‘I am happy for John Firth to record my name and email address, to keep me informed about his services and availability’. Avast logo <https://www.avast.com/antivirus> This email has been checked for viruses by Avast antivirus software. www.avast.com <https://www.avast.com/antivirus> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2> ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Generate the bearing in the attribute table
Hi Carlos; It seems you have the calculation of azimuth figured out. I would generate a new filed in the attribute table and use the field calculator to calculate your bearings. If you are comfortable with python, you can create a script as well. If you are using a projected coordinate system, take note that North would be actually be grid north and your calculate azimuth would be in relation to it. If you want to reference actual true north, you would have to calculate the convergence between true and grid north for your project area. If you want the coordinates referencing magnetic north (perhaps for field work), then you would also have to calculate the magnetic declination for your project. Kirk Schmidt On 4/8/2021 6:30 PM, Sebastian Gutwein wrote: Carlos- Interesting problem- You might get better responses on stack exchange. I have the same confusion as Ujaval. Each side of your polygon will have a different bearing. If all you care about is the average angle then you can use this expression main_angle( $geometry) This will return azimuth of the polygon in degrees from north. I don't know how to convert this to bearing in QGIS. On Thu, Apr 8, 2021 at 3:28 PM carlos h <mailto:carloslda...@gmail.com>> wrote: Hello, for example have this in the attribute table: https://ibb.co/wWn4PfM <https://ibb.co/wWn4PfM> What I need is to have the bearing in the polygon's attribute table. Thanks CA image.gif El jue, 8 abr 2021 a las 10:49, Ujaval Gandhi (mailto:uja...@spatialthoughts.com>>) escribió: Can you elaborate with an example? It's not clear what the bearing would be for a polygon. Logo <https://spatialthoughts.com/> Ujaval Gandhi Spatial Thoughts mobile: +91-8095684687 email: uja...@spatialthoughts.com <mailto:uja...@spatialthoughts.com> LinkedIn icon <https://www.linkedin.com/in/spatialthoughts/> Twitter icon <https://twitter.com/spatialthoughts> On Thu, Apr 8, 2021 at 8:50 PM carlos h mailto:carloslda...@gmail.com>> wrote: Hello, nice to write here. How can I generate the bearing in the attribute table of a polygon? the bearing not in azimuth, but for example N30E etc Thanks CA image.gif ___ 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> ___ 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> ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] WMS Image Projection
Hi Springer. My mistake, you will require a local grid shift file, BC_93_05.gsb. The NTv2 is a national NAD27to NAD83 grid shift file. The BC_93_95 provides localized conversion from Nad83 (original) to Nad83 (csrs). Sent from my Galaxy Original message From: Springfield Harrison Date: 2021-04-08 5:24 a.m. (GMT-04:00) To: Andreas Neumann Cc: qgis-user Subject: Re: [Qgis-user] WMS Image Projection Hi Andreas, Thanks for the reply, I'll try to clarify: Version is 3.16.4 The two municipal sources are being viewed in their own web windows and are displaying the same data They are the source of the property lines for my GIS (downloaded shapefile) and the orthophoto (WMS link, not resident TIFF or other) The property lines are EPSG:3157 - NAD83(CSRS) / UTM zone 10N - Projected, and from the municipal source. I've attached the shapefile if you have the time. In QGIS, the ortho is loaded from the WMS source https://mapservices.crd.bc.ca/arcgis/services/OrthoImagery2019/MapServer/WMSserver EPSG:26910 - NAD83 / UTM zone 10N - Projected The project CRS is EPSG: 26910 but it doesn't seem to change anything by changing that or the source projections (at least within versions of NAD83 UTM Zone 10N) Thus, the same information displays differently online versus in QGIS. The shift seems to be 1.5 m +/-, quite a lot for what I'm trying to do (building permit with tight setbacks). These are commonly used projections for our neck of the woods The two images I sent earlier depict the differences; in QGIS, earlier orthophotos show a similar or worse displacement. Not so on the municipal web viewer. I've read up on projections in the tutorial and Help file, no revelations there, at least that I could see. Thanks very much for any thoughts you may have. - Cheers, Spring On 08/Apr/2021 00:47, Andreas Neumann wrote: Hi, Are there two different CRS involved here? Are the vectors (parcel boundaries) in a different CRS than the orthoimages? This might be a problem of a missing grid-shift file in QGIS. It is also unclear for me if you load your Orthoimages locally (from tiff-files as well) or from WMS (or both) and only one of them is shifted? Can you also add information what QGIS version with what proj version you use? You can find this information in the about dialogue in the help menu. Greetings, Andreas On 2021-04-08 02:37, Springfield Harrison wrote: Hello All, Re obtaining orthophotos from a WMS server. Locally we have orthophotos that can be brought into QGIS using the WMS service. They are also viewable on the local government web map service. In both cases I can compare them to the property lines which are from the same source. On the government web pages (2), the orthophotos and the lot lines appear to coincide very well. In QGIS, the same two layers are quite misaligned. One would suspect a CRS problem but changing the CRS of the Project and the WMS image in QGIS has no effect. Even followed the tutorial to no effect. I'm attaching screen shots of the two results. Ant thoughts would be appreciated, thank you . . . . - Cheers, Springfield Harrison ___ 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] WMS Image Projection
Hi Springer.I believe you will need to apply a grid shift file (NTv2). It is available from Natural Resouces Canada and I will send you a link when I get to the office.The installation is a simple procedure. Kirk SchmidtSent from my Galaxy Original message From: Andreas Neumann Date: 2021-04-08 6:52 a.m. (GMT-04:00) To: Springfield Harrison Cc: qgis-user Subject: Re: [Qgis-user] WMS Image Projection Hi, The approx 1.5m offset indicates that QGIS doesn't take a necessary grid-shift file into account. The grid-shift file is a separate file that defines the displacement vectors of coordinate systems with different datums. As a European GIS user I am not comfortable with North-American CRS, but looking the two CRS up at epsg.io I can see that http://epsg.io/?q=26910 vs http://epsg.io/?q=3157 EPSG 26910 claims to have an accuracy of 4 m while EPSG 3157 claims to have an accuracy of 1m. There might be a grid shift file necessary between the two CRS (I am personally unsure if that is the case here). Since both, the project and the WMS are in the same CRS, I would assume that the offset problems occur with the parcel shapefile. Do you know a local GIS expert in your region who knows if the transformation between EPSG 26910 and 3157 requires a grid-shift file? I would assume yes, because at epsg.io it indicates that the datums of the two CRS are slightly different: NAD83(CSRS) (http://epsg.io/4617) vs NAD83 (http://epsg.io/4269). From there - and if we have a proper corresponding grid-shift file - we can help you set up the transformation in QGIS. Greetings, Andreas On 2021-04-08 10:24, Springfield Harrison wrote: Hi Andreas, Thanks for the reply, I'll try to clarify: Version is 3.16.4 The two municipal sources are being viewed in their own web windows and are displaying the same data They are the source of the property lines for my GIS (downloaded shapefile) and the orthophoto (WMS link, not resident TIFF or other) The property lines are EPSG:3157 - NAD83(CSRS) / UTM zone 10N - Projected, and from the municipal source. I've attached the shapefile if you have the time. In QGIS, the ortho is loaded from the WMS source https://mapservices.crd.bc.ca/arcgis/services/OrthoImagery2019/MapServer/WMSserver EPSG:26910 - NAD83 / UTM zone 10N - Projected The project CRS is EPSG: 26910 but it doesn't seem to change anything by changing that or the source projections (at least within versions of NAD83 UTM Zone 10N) Thus, the same information displays differently online versus in QGIS. The shift seems to be 1.5 m +/-, quite a lot for what I'm trying to do (building permit with tight setbacks). These are commonly used projections for our neck of the woods The two images I sent earlier depict the differences; in QGIS, earlier orthophotos show a similar or worse displacement. Not so on the municipal web viewer. I've read up on projections in the tutorial and Help file, no revelations there, at least that I could see. Thanks very much for any thoughts you may have. - Cheers, Spring On 08/Apr/2021 00:47, Andreas Neumann wrote: Hi, Are there two different CRS involved here? Are the vectors (parcel boundaries) in a different CRS than the orthoimages? This might be a problem of a missing grid-shift file in QGIS. It is also unclear for me if you load your Orthoimages locally (from tiff-files as well) or from WMS (or both) and only one of them is shifted? Can you also add information what QGIS version with what proj version you use? You can find this information in the about dialogue in the help menu. Greetings, Andreas On 2021-04-08 02:37, Springfield Harrison wrote: Hello All,Re obtaining orthophotos from a WMS server.Locally we have orthophotos that can be brought into QGIS using the WMS service. They are also viewable on the local government web map service. In both cases I can compare them to the property lines which are from the same source.On the government web pages (2), the orthophotos and the lot lines appear to coincide very well. In QGIS, the same two layers are quite misaligned.One would suspect a CRS problem but changing the CRS of the Project and the WMS image in QGIS has no effect. Even followed the tutorial to no effect.I'm attaching screen shots of the two results. Ant thoughts would be appreciated, thank you . . . .-Cheers, Springfield Harrison ___Qgis-user mailing listQgis-user@lists.osgeo.orgList info: https://lists.osgeo.org/mailman/listinfo/qgis-userUnsubscribe: 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] off topic Deimos 2 images
Hi Carlos: I just looked at the technical specifications for DEIMOS-2 data and data are available at 2 processing levels, L1B and L1C, both of which are calibrated and radiometrically corrected. The L1C level data are orthorectified. So your data is most likely already corrected if you purhased it from a vendor. Kirk Schmidt On 4/4/2021 7:13 PM, carlos h wrote: Hello, a pleasure to greet you, has someone carried out atmospheric correction of images of deimos 2? Thanks CA ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] off topic Deimos 2 images
Hi Carlos: Do you have the calibration files? These are usually available from the data provider. Kirk Schmidt On 4/4/2021 7:13 PM, carlos h wrote: deimos 2 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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 access r.series command
Hi Jeff: Under the top right word menu, click on Processing>Toolbox to activate the toolbox which will appear on the right hand side of the canvas. Add r.series to the search window, or navigate to grass> raster>r.series If grass is not activated, left click the tool icon in the processing toolbox. Look under providers and activate Grass. Kirk Schmidt On 3/30/2021 1:19 PM, gafna jeff wrote: Hello everyone, I would like to sum up pixel values of several raster layers using the r.series command but I don't seem to know where to find it on QGIS after loading all my rasters. I would appreciate any help Jeff ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] First steps with 3D Scan Data?
Hi Bernd: Just a clarification on LAStools. The library has open as well as licensed algorithms. Typically, the basic functionality is available for free and implemented in QGIS. The advanced tools are available through licensing. Rapidlasso GmbH also developed the laz format which is an efficient compression standard to make large datasets/point clouds more manageable. PDAL and R can be used to process the data but if your data is compressed, you would have to use the unlicensed Laszip routine available in QGIS 3.x to convert the data to an uncompressed las format. Kirk On 3/26/2021 6:23 AM, Karl Magnus Jönsson wrote: Hi! Lots of suggestions here. But if I understand you correct it is lidar data and all about photogrammetry is not in question. Why not pdal? Open source, powerful and what QGIS is using in 3.18. https://pdal.io/ Karl-Magnus -Ursprungligt meddelande- Från: Qgis-user För Bernd Vogelgesang Skickat: den 26 mars 2021 09:35 Till: Benoit de Cabissole ; qgis-user Ämne: Re: [Qgis-user] First steps with 3D Scan Data? Hi Benoit, ok, thanks. las then. But when I read about the lastools plugin, the description says that I need to install the lastools software, and the website there tells me that the license is about 1500€ which is a little too much for some learning efforts. Any other options? Cheers, Bernd On 26.03.21 09:28, Benoit de Cabissole wrote: Hi Bernd, OK, that more info! The files are likely to be in Timble RealWorks format. So you will have to get your contractor to provide you with .las files, otherwise your stuck buying RealWorks! Or, perhaps, the viewer is able to export in .las?? When you get las files, you can use lastools plugin for QGIS to process them. Cheers, Benoit On 26/03/2021 10:01, Bernd Vogelgesang wrote: Hi Benoit, as I was not present when the scans were made, I can only assume that this was a laser scanner on a tripod. With the sent data, there is a Trimble Scan Explorer Viewer software where I can make a "virtual walk" through the scene, so there have to be photos but also laser scan data as I can see from different view modes. As I said, this is the very first time I get in touch with this stuff, and Google was of no big help so far. Besides of tons of other files, the bigger ones with around 50 to 80 MB are .rwcx and .tzf files. I hoped that someone would know what these file explicitly hold, and what I have to do to get them e.g. into QGIS or other open source software. Cheers, Bernd On 26.03.21 08:40, info wrote: Hi Bernd, Well, could you give more info? E.g. what kind of scan: laser/lidar/photo/other? This will influence your workflow. +/-Open source laser/lidar could be based on lastools [1] Open source photogrammetry could be based on AliceVision meshroom [2,3], MeshLab [4] and blender [5] In photogrammetry: Agisoft metashape [6] and pix4D [7] as good commercial products. BTW: 32 scans, if photos, is quite a very low number of photos for reconstructing a 3D scene! Hope this helps. Cheers, Benoit de Cabissole [1] http://lastools.org/ [2] https://alicevision.org/#photogrammetry [3] https://github.com/alicevision/meshroom [4] https://www.meshlab.net/ [5] https://www.blender.org/ [6] https://www.agisoft.com/ [7] https://www.pix4d.com/ On 25/03/2021 23:07, Bernd Vogelgesang wrote: Hi folks, in a project on a former dumpside now to be transformed into a sand lizard habitat, I saw the chance to get my hands dirty for the first time on 3D scan data (terrestrial). I received 10GB of stuff, and now I'm completely lost. The data consists of .rwcx files, .tzf files and lots of other stuff. Is there any chance that I can process any of this data with an open source stack? Or what should be done, to make them digestible? Second step: There seem to be 32 different scans of the area form all angles. In case I succeed to somehow load the data, how do I combine them ? For the first steps, I actually "only" would like to generate a decent DEM from the data to be able to place some features on a map for the guys in charge to do the habitat stuff. Any kind of hints or links to reads about this topic would be appreciated Cheers, Bernd ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1
Re: [Qgis-user] First steps with 3D Scan Data?
Hi Bernd.It sounds like the data were acquired and processed with Trimble hardware and software. I am not familiar with their products however they tend to write proprietary files unless they are exported to a open data format within the Trimble Procesding software. I will have a look at ghd standards when I have a few moments. As the data are kased scans, Agrisofglt or Pix4d software are not suitable as these aed I are based, structure from motion solutions.Kirk SchmidtSent from my Galaxy Original message From: Bernd Vogelgesang Date: 2021-03-26 7:34 a.m. (GMT-04:00) To: Richard McDonnell , info , qgis-user Subject: Re: [Qgis-user] First steps with 3D Scan Data? Hi Richard,wow, this is amazing stuff in the video!For the type of things I am doing, this solution would be quit sufficient.I think the data I received is way too detailed, but I'll try to get intouch and ask them to convert it to .las if possible.Thanx a lot for all the hints!BerndOn 26.03.21 10:38, Richard McDonnell wrote:> Hi,> If you are looking for an Opensource Pointcloud Stack, have a look at this video on youtube by Klas Karlsson>> https://www.youtube.com/watch?v=pe04m-gHM-M=1180s>> To make the above work for you , you would have needed to receive the scan data in LAS format, it would appear that the rcwx files are proprietary format from Trimble, and are Trimble RealWorks Project files. The Pointcloud data is in the .tzf files Can you ask the provider of the data to export the scan data as a LAS Format File? You should then be able absorb the data into an Open Source Stack i.e.Cloud Compare and QGIS 3.18>> Hope that’s of some guidance and help.> Regards,>> Richard.>>> ——> Richard McDonnell MSc GIS, FME Certified Professional> FRM Data Management>> ——> Oifig na nOibreacha Poiblí> Office of Public Works>> Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36> Jonathan Swift Street, Trim, Co Meath, C15 NX36> ——> M +353 87 688 5964 T +353 1 647 6543> https://gov.ie/opw>> ——> To send me files larger than 30MB, please use the link below https://filetransfer.opw.ie/filedrop/richard.mcdonn...@opw.ie>> Email Disclaimer: https://www.gov.ie/en/organisation-information/439daf-email-disclaimer/ ——> Richard McDonnell MSc GIS, FME Certified Professional> FRM Data Management>> ——> Oifig na nOibreacha Poiblí> Office of Public Works>> Sráid Jonathan Swift, Baile Átha Troim, Co na Mí, C15 NX36> Jonathan Swift Street, Trim, Co Meath, C15 NX36> ——> M +353 87 688 5964 T +353 1 647 6543> https://www.opw.ie>> ——> To send me files larger than 30MB, please use the link below https://filetransfer.opw.ie/filedrop/richard.mcdonn...@opw.ie>> Email Disclaimer: https://www.opw.ie/en/disclaimer/> -Original Message-> From: Qgis-user On Behalf Of Bernd Vogelgesang> Sent: 26 March 2021 08:01> To: info ; qgis-user > Subject: Re: [Qgis-user] First steps with 3D Scan Data?>> Hi Benoit,>> as I was not present when the scans were made, I can only assume that this was a laser scanner on a tripod.>> With the sent data, there is a Trimble Scan Explorer Viewer software where I can make a "virtual walk" through the scene, so there have to be photos but also laser scan data as I can see from different view modes.>> As I said, this is the very first time I get in touch with this stuff, and Google was of no big help so far.>> Besides of tons of other files, the bigger ones with around 50 to 80 MB are .rwcx and .tzf files.> I hoped that someone would know what these file explicitly hold, and what I have to do to get them e.g. into QGIS or other open source software.>> Cheers,>> Bernd>> On 26.03.21 08:40, info wrote:>> Hi Bernd, Well, could you give more info? E.g. what kind of scan:>> laser/lidar/photo/other? This will influence your workflow. +/-Open source laser/lidar could be based on lastools [1]>> Open source photogrammetry could be based on AliceVision meshroom>> [2,3], MeshLab [4] and blender [5] In photogrammetry: Agisoft metashape [6] and pix4D [7] as good>> commercial products. BTW: 32 scans, if photos, is quite a very low number of photos for>> reconstructing a 3D scene! Hope this helps.>> Cheers,>> Benoit de Cabissole [1] http://lastools.org/>> [2] https://alicevision.org/#photogrammetry>> [3] https://github.com/alicevision/meshroom>> [4] https://www.meshlab.net/>> [5] https://www.blender.org/>> [6] https://www.agisoft.com/>> [7] https://www.pix4d.com/ On 25/03/2021 23:07, Bernd Vogelgesang wrote:>>> Hi folks,>> in a project on a former dumpside now to be transformed into a sand>>> lizard habitat, I saw the chance to get my hands dirty for the first>>> time on 3D scan data (terrestrial). I received 10GB of stuff, and now>>> I'm completely lost.>> The data consists of .rwcx files, .tzf files and lots of other stuff.>> Is there any chance that I can process any of this data with an open>>> source stack?>>> Or
Re: [Qgis-user] Map styles changing when I export as image
Hi. Jeremy. Is the issue that the satellite image is not being printed? If so, it is probably an issue with the image link and or copyright limitations if it originates from a internet based mapping service.Kirk SchmidtSent from my Galaxy Original message From: jeremy benwell Date: 2021-03-21 5:44 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Map styles changing when I export as image Hi all,I'm having some trouble with exporting my map as an image. When I try to export my map as an image all the styles change (see images below, I had to shrink size to make email within 100kb limit so they are a bit blurry). The dpi I choose is 100. I'm guessing the problem has something to do with what layers are selected in the project tab but can't seem to figure out how exactly. Is there some way to stop the project tab from feeding information to the print layout tab so i can export the print layout tab as an image exactly as is.Any help would be much appreciated.Cheers,JeremyHow I want my map to be exported as an imageMy project layersHow my map is actually exported as an image ___ 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] Antivirus objecting to QGIS component
Hi NeilThis has come up on the past. https://github.com/qgis/QGIS/issues/32247Kirk SchmidtSent from my Galaxy Original message From: Neil Mander Date: 2021-03-18 7:54 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Antivirus objecting to QGIS component In the last week, my AV software has been flagging a file from within my QGIS package as suspicious. It is both my QGIS2 and QGIS3 packages Can anyone please confirm to me that nircmdc.exe is a legitimate part of QGIS. I see that it is something that edits registry entries. Many thanks to anyone who can help. Neil ___ 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] Do GPX files contain CRS information?
Thanks for the research Nicolas. I think I will save this for a rainy day or when an issue comes up.Kirk SchmidtSent from my Galaxy Original message From: Nicolas Cadieux Date: 2021-03-17 7:14 p.m. (GMT-04:00) To: Greg Troxel , stewartbh...@gmail.com Cc: qgis-user Subject: Re: [Qgis-user] Do GPX files contain CRS information? Hi, This last question brought me down a rabbit hole that took me a while to find (at least partially). Of course, this is far away from the original question of "Do GPX files contain CRS information?". The answer to that was no. The datum is WGS84 (in it's most current iteration, or WGS84(G1762). Transformation parameters to and from ITRF have been published and can be found https://confluence.qps.nl/qinsy/latest/en/world-geodetic-system-1984-wgs84-182618391.html#id-.WorldGeodeticSystem1984(WGS84)v9.1-WGS84realizations This file contain multiple transformations with sources https://confluence.qps.nl/qinsy/files/latest/en/182618383/182618384/1/1579182881000/ITRF_Transformation_Parameters.xlsx This contains info on WGS84 and NAD83 https://mcraymer.github.io/geodesy/pubs/nad83_agu2007spr.pdf https://www.unoosa.org/documents/pdf/icg/2018/icg13/wgd/wgd_12.pdf Other question that were raised: What are the other reference frames used by the various GNSS services: https://www.sciencedirect.com/science/article/pii/S0273117720308292#! "The TRFs realized by the GPS, GLONASS, Galileo, and BeiDou-2 and BeiDou-3 broadcast ephemerides are the orbital realizations of WGS 84 (G1762′), PZ90.11, GTRF19v01, and BDCS respectively." More info here. https://gssc.esa.int/navipedia/index.php/Reference_Frames_in_GNSS This document made a comparison comparison between broadcast reference frames and ITRF: http://www.epncb.eu/_newseventslinks/workshops/EPNLACWS_2017/pdf/06_Open_Session/04_Broadcast-Precise.pdf "Reference frames: •GPS and Galileo broadcast reference frames are aligned with ITRF: translations are less than 0.10 m and rotations are less than 2 milli-second of arc •GLONASS M broadcast reference frame is offset to ITRF by at most 0.27 ±0.04 m in Y and maximum rotation is 4 ±2 milli-second of arc about Y. •GLONASS K broadcast reference frame is offset to ITRF by at most 1.06 ±0.17 m in Y and maximum rotation is 19 ±2 milli-second of arc about X. •BeiDoubroadcast reference frame is offset to ITRF by at most 0.26 ±0.18 in X and Y, and maximum rotation is 2 ±1 milli-second of arc about Z." Then there are questions pertaining to SBAS around the world: https://www.faa.gov/about/office_org/headquarters_offices/ato/service_units/techops/navservices/gnss/library/factsheets/media/SBAS_Worldwide_QFact.pdf It state the various services Wide Area Augmentation System (WAAS) European Geostationary Navigation Overlay Service (EGNOS) Multi-functional Transport Satellite (MTSAT) Satellite Based Augmentation System (MSAS) (Japan) GPS Aided Geostationary Earth Orbit (GEO) Augmented Navigation (GAGAN) (India) System of Differential Correction and Monitoring (SDCM) Russia Korean Augmentation Satellite System (KASS) BeiDou Satellite Based Augmentation System (BDSBAS) Also found this https://www.gps.gov/technical/ps/2008-WAAS-performance-standard.pdf And this https://www.nstb.tc.faa.gov/ Finally, there was the question of the introduction of new reference frames when a SBAS services was used. I was skeptical of this but you guys are right. Using a WAAS service does introduce a new reference frame as the "... GEO satellites do not belong to any satellite positioning service (e.g. GPS, GLONASS), ephemeris for those satellites are not externally available. Therefore, it is the SBAS that is in charge of providing the user with the GEO ephemeris. Keep in mind that all components are expressed in ECEF reference coordinates and the time offset is with respect to SBAS Network time (SNT)." (https://gssc.esa.int/navipedia/index.php/The_EGNOS_SBAS_Message_Format_Explained#SBAS_broadcast_data) Same information is stated here (https://www.sciencedirect.com/science/article/pii/S0273117720308292#!) "Also, note that the methodology and results reported here are only valid for direct, real-time unaugmented GNSS applications. Any form of augmentation including all types of differential positioning, Assisted GPS, or specialized augmentation, such as the US
Re: [Qgis-user] Learning resources/books for QGIS
Hi Rod: I would suggest a video to get a feel for how things work. Things are a lot easier without command level inputs. Try: https://www.youtube.com/watch?v=kCnNWyl9qSE There are lots of resources on youtube. Kirk Schmidt On 3/12/2021 1:39 PM, Roderick Anderson wrote: Good day. I looking for some suggestions of books on learning QGIS. I have some background (very long ago - late 80's / early 90's) with GIS but at that time is was more about mapping not analysis. In fact, I am looking to use QGIS for building maps for personal use. Something along the lines of the Washington state (US) "Map Your Neighborhood" project. The ability to easily load personal data is important. For me the QGIS Documentation assumes a more current knowledge of the technology and terminology than I have. I've looked at the QGIS book selections but not quite sure which would the best for my needs. Any and all suggestions are welcome. TIA, Rod -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Combing Vector Layers - With Attributes
Hi Springer You could try joining the databases using a common field assuming the second one does not have new point features. However, I expect you are adding new point features. In that case, I usually write a little python script to generate a new set of features in the update layer that matches those of the master layer. I then recalculate the values of the new features before dropping the initial set of fields. You can use Graphical Modeler to generate a model and access through the Processing Toolbox> Models Group. If you send me to sample shapefiles, I can put it together if you are not familiar with the Graphical Modeler. Kirk Schmidt On 3/12/2021 4:35 AM, Springfield Harrison wrote: Hello All, I have a master database (point shapefile) of trees to which I need to add new records frequently. The two files have different (but similar) attribute lists. I need to combine them and add the new attributes into the appropriate SHP database columns. Copying and pasting records from the new file into the master list does work but if the fields are not identical, no attributes are brought into the master file. Is there a way to import the new records into their corresponding fields in the master file? I've searched for help on this to no avail. Manifold GIS handles this very well, setting up a field matching dialogue but I would prefer to keep the process within QGIS. Thanks very much for any ideas . . . . - Cheers, Springfield Harrison ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Georeferencer GDAL - Tool - What is the meaning of the scaling factor in the Info
Hallo Michael: I assume you are referring to the parameters in the lower part of the georeferencing window? I qset up a test georeferencing and the screen shot is below. In this case, the values are the transformation parameters based on the 4 GCP's I used. The value after the | are the coordinates of the cursor measured from the top left corner. Is this what you are referring too? If the values are in your main map view, something has gone wrong in your georefrencing and perhaps you do not have enough GCP or one of your coorinates is wrong. Try overlaying some additional vector or raster data to confirm your image is in the right spot. Kirk Schmidt On 3/11/2021 10:47 AM, Michael Werner Maur wrote: Hello, I use the Georeferencer in the QGis Raster dialog för a long time. But yet I need further and detailed information about the quality parameters of the transformation. For example I did a Helmert Transformation and find in the lower part of my desktop the parameters. There are two scaling factors for each coordinate axis, but they are as the following value (0,05325...) ...I expected something in relation to the original scale of the scanned raster map (f.e. 1 : 1250 -> 0,0008). It will be very helpful, if there is a explanation to it ! Thanx a lot ! Mit freundlichen Grüßen Im Auftrag Michael Maur -- Michael Maur Fachgruppenleiter 2.1 VERMESSUNGS- UND KATASTERAMT OSTEIFEL-HUNSRÜCK Am Wasserturm 5a 56727 Mayen Telefon 02651 9582-120 Telefax 02651 9582-400 michael.m...@vermkv.rlp.de www.vermka-osteifel-hunsrueck.rlp.de Das Thema Datenschutz ist für uns eine wichtige Angelegenheit. Bitte beachten Sie unsere geltenden Datenschutzbestimmungen <https://lvermgeo.rlp.de/de/wichtige-Informationen/datenschutz> und unsere datenschutzrechtlichen Vorgaben bei Bewerbungsverfahren. <https://lvermgeo.rlp.de/fileadmin/lvermgeo/pdf/wir_ueber_uns/Datenverarbeitung_in_Bewerberverfahren.pdf> ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Trimble GeoXT 2005 Accuracy
Hi Springer.Have fun and your questions stimulated a good discussion and ideas.If you are collecting gps data on the clock, I would seriously consider installing a ground plane on your range pole. This will eliminate zingers and you can reduce your observation time. A electrical junction box cover works well.Waas currently is good to 0.9 m horizontal and 1.3 m vertical 95% of the time. It helps but is not a stand alone solution for high precision work.Kirk SchmidtSent from my Galaxy Original message From: Springfield Harrison Date: 2021-03-10 4:30 a.m. (GMT-04:00) To: Kirk Schmidt , qgis-user@lists.osgeo.org Subject: Re: [Qgis-user] Trimble GeoXT 2005 Accuracy Hi Kirk,OK, I'll rule out the PPP for now, although it doesn't seem to be practical for production mapping anyway.I think my antennas are probably OK, none of them have a ground plane, not practical for field mapping. I've only seen them in geodetic applications. Of course, the aircraft installation is free of multipath except while on the ground. Longer observing periods should minimize the impact of multipath zingers.The tight clustering of the GeoXT points in Map A doesn't smell of multipath to me.Earlier you noted the possibility of antenna shielding by the operator. With that in mind, I use a tripod or a range pole with the antenna above head height. WAAS is actually quite accurate, good for certain (lower order?) precision approaches.I'm still digesting all this helpful feedback and will work on a summary shortly. Thanks again . . . .-Cheers, SpringOn 09/Mar/2021 12:14, Kirk Schmidt wrote:> Hi Springfield:>> PPP requires long observation times in the order of hours and > therefore a new field survey. I would try post processing from a high > precision site with your existing data.>> Another issue may be the aviation antennae you mentioned. These are > usually mounted on the top of a fuselage and therefore the aluminum > acts as a ground plane to effectively filter out multipath. If you > use a metal disk or plate and a ground plane on your antennae and > mount it on a tripod, you may find you end up with better results. > You will notice survey grade antennae's have a large diameter plastic > case since the contain an internal 10 cm metal disk.>> Kirk Schmidt>> On 3/9/2021 1:36 PM, Springfield Harrison wrote:>> Hello Garth,>>>> Thanks very much for that reference. I have used the service in the >> past.>>>> I'm in the middle of other work but may try to re-process some Map A >> files using PPP to see what changes.>>>> Your work sounds more like surveying than mapping, interesting.>>>> I do use the GPS reception planning tools but SV availability is >> certainly better than in the early days.>>>> Thanks again . . . .>>>> ->> Cheers, Spring>>>>>>>> On 09/Mar/2021 09:01, Garth Fletcher wrote:>>> Nicolas Cadieux wrote:>>>> Kirk could be on the right track. You could try PPP using this >>>> site. You will need to create a user name and password. Then, you >>>> upload the Rinex file. I am 99% sure you can process point from >>>> outside of Canada. >>>>>> I routinely use CSRS PPP to post process RINEX files from Southern New>>> Hampshire (43°N, 72°W).>>>>>> Using an iGage iG3s (L1/L2 dual frequency, GPS + GLONASS) recording at>>> 5 second epochs in forested areas, from post processing I get>>> 30 minute recordings, 95% error ellipses < 1 meter>>> 10 hour recordings, 95% error ellipses < 2 cm>>>>>> CSRS-PPP site is>>> <https://webapp.geod.nrcan.gc.ca/geod/tools-outils/ppp.php?locale=en>>>> enter your email address>>> select Static and which NAD 83 epoch you want, or ITRF>>> Submit a RINEX file (.zip compressing the file is encouraged)>>> they will email back a detailed report>>>>>> You can submit for 1 of 3 levels of post processing>>> ultra-rapid : submit > 90 minutes after end of recording>>> rapid: submit > 24 hours after end of recording>>> final: submit ~ 3 weeks after end of recording>>>>>> Major part of corrections are in ultra-rapid, the higher levels provide>>> further refinements.>>>>>>>>> FYI: Trimble provides occupation planning aids at>>> <https://www.gnssplanning.com/#/charts>>>> just enter your location and date and it shows you the DOP as>>> a function of the hour of day.>>>>>> Cordially,>> ___>> 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] Trimble GeoXT 2005 Accuracy
Hi Springfield: PPP requires long observation times in the order of hours and therefore a new field survey. I would try post processing from a high precision site with your existing data. Another issue may be the aviation antennae you mentioned. These are usually mounted on the top of a fuselage and therefore the aluminum acts as a ground plane to effectively filter out multipath. If you use a metal disk or plate and a ground plane on your antennae and mount it on a tripod, you may find you end up with better results. You will notice survey grade antennae's have a large diameter plastic case since the contain an internal 10 cm metal disk. Kirk Schmidt On 3/9/2021 1:36 PM, Springfield Harrison wrote: Hello Garth, Thanks very much for that reference. I have used the service in the past. I'm in the middle of other work but may try to re-process some Map A files using PPP to see what changes. Your work sounds more like surveying than mapping, interesting. I do use the GPS reception planning tools but SV availability is certainly better than in the early days. Thanks again . . . . - Cheers, Spring On 09/Mar/2021 09:01, Garth Fletcher wrote: Nicolas Cadieux wrote: Kirk could be on the right track. You could try PPP using this site. You will need to create a user name and password. Then, you upload the Rinex file. I am 99% sure you can process point from outside of Canada. I routinely use CSRS PPP to post process RINEX files from Southern New Hampshire (43°N, 72°W). Using an iGage iG3s (L1/L2 dual frequency, GPS + GLONASS) recording at 5 second epochs in forested areas, from post processing I get 30 minute recordings, 95% error ellipses < 1 meter 10 hour recordings, 95% error ellipses < 2 cm CSRS-PPP site is <https://webapp.geod.nrcan.gc.ca/geod/tools-outils/ppp.php?locale=en> enter your email address select Static and which NAD 83 epoch you want, or ITRF Submit a RINEX file (.zip compressing the file is encouraged) they will email back a detailed report You can submit for 1 of 3 levels of post processing ultra-rapid : submit > 90 minutes after end of recording rapid: submit > 24 hours after end of recording final: submit ~ 3 weeks after end of recording Major part of corrections are in ultra-rapid, the higher levels provide further refinements. FYI: Trimble provides occupation planning aids at <https://www.gnssplanning.com/#/charts> just enter your location and date and it shows you the DOP as a function of the hour of day. Cordially, ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Trimble GeoXT 2005 Accuracy
Hi SpringerI just looked up the marine beacon specs and they are designed for
10m accuracy. This may explain why your non dgps data is closer to your survey
corner.Sent from my Galaxy
Original message From: Springfield Harrison
Date: 2021-03-08 3:17 p.m. (GMT-04:00) To: Nicolas
Cadieux Cc: kirk ,
Jorge Gustavo Rocha , qgis-user@lists.osgeo.org, Greg Troxel
, Dan <19dm...@gmail.com> Subject: Re: [Qgis-user] Trimble
GeoXT 2005 Accuracy
I'm resending this without the map as there is a size limit. The moderator
may let it through, I hope . . . .
Hi Nicolas, thanks for your observations. I'll try to answer your
questions, please see the attached map, especially Map A:
Note that my previous email contained information for Map B; Map A is based
on the Municipal Cadastre (NAD83 UTM zone 10N) and illustrates the problem as
well. Other locations based on Provincial Monuments and/or the Municipal
cadastre (not illustrated here) have yielded similar results.
How many “known” points have you tested?
2 in this case, Maps A and B
Also several other locations with similar results
How where those point position calculated.
From the Municipal cadastre, visible in Map A
Make sure the coordinates are in the right CRS
NAD 83 UTM 10N used throught. See workflow in previous email
When converting from the monument’s CRS to NAD83 UTM zone 10N, are you
using the correct grid files?
[No monuments in this example] These were brought into QGIS from the
Municipal GCM database CSV (NAD83(CSRS) 3.0.0.BC.1.CRD) and reprojected by QGIS
to EPSG:26910 - NAD83 / UTM zone 10N
Presumably QGIS would choose the correct grid files
Municipal Cadastre is NAD83 UTM zone 10Np, li { white-space: pre-wrap;
}
Find a geodesic point that is in the middle of a field or on the side of
a highway with no obstacles.
Map A is open sky
Make sur your observations will be done when the constellation is well
distributed in the sky
As you probably know, TerraSync provides for PDOP, HDOP, SNR and
Horizon masks to preclude collecting poor quality positions. These were set
towards the "Precision" end of the scale
What post processing techniques are you using? How far is the base
station from your unit?
Real time was SBAS or RTCM; Post processing using the Pathfinder Office
differential correction engine, baseline about 30 km
How long are the observations? Have you tried other methods of post
processing like PPP?
Logging interval is 5 sec; 33 to 2037 positions per point
Did not use PPP. This is a test of mapping best practices, not geodesy
Have you contacted Trimble?
Yes, no response
Have you looked on there site to see if there is a software update
(firmware) for the unit or the post processing software?
Yes, receiver firmware is the latest, PFO and Terrasync are older but
compatible
Thanks Nicolas. If I have missed something, I hope someone can point it
out, I've tried to cover all the bases based on my training and experience.
-
Cheers, Spring
On 08/Mar/2021 06:40, Nicolas Cadieux
wrote:
Hi Harrison,
How many “known” points have you tested? How where those
point position calculated. They could be off. If you are using
state geodesic monuments, try to find the documented precision
of the monument. States have different types of monuments, some
are very old and have different standards. Make sure the
geodesic point is not the problem. Make sure the coordinates
are in the right CRS. As an example, if the coordinates are
published in NAD83 original but you are assuming NAD83(CSRS),
then you have a problem. When converting from the monument’s CRS
to NAD83 UTM zone 10N, are you using the correct grid files?
What is the published precision for this reprojection?
You say you have houses and trees. This could be the
problem. Find a geodesic point that is in the middle of a field
or on the side of a highway with no obstacles. Make sur your
observations will be done when the constellation is well
distributed in the sky. I believe Trimble has a observation
planing software that can help you figure out the best time for
observation. This could explain why the GEoTX are to the east
unless the observations where made at the same time and same
conditions (ex leaf off).
What post processing techniques are you using? How far is the
base station from your unit? If you a
Re: [Qgis-user] Trimble GeoXT 2005 Accuracy
Hi Springer.Thanks for providing the high level of detail.I assume the
comparison of field units occurred at separate times or even days?In map 1,
where were you standing in relation to the internal antenae. Your body will
block and or distort radio signals. It is good practice to set the field unit
on a tripod or rock and move away from the unit while observing. Additionally,
how did you set up the external antenae when you used it (tripod, pack pack?).I
am not sure what the marine beacon accuracy is as we used them in eastern
Canada during the days of selective availability. WAAS may be totally useless
in this application as the system is designed for coarse aircraft navigation to
get down to 200' of a runway.Have a look at the Canadian Active Control system.
Perhaps you are close to one if the federal or western deformation array
sites. You can download the historical data for these sites which can include
broadcast emphemerifes and clock corrections. I do not know if these sites are
implemented in PFO. I assume you can import into PFO and DGPS. Kirk SchmidtSent
from my Galaxy
Original message From: Springfield Harrison
Date: 2021-03-08 3:17 p.m. (GMT-04:00) To: Nicolas
Cadieux Cc: kirk ,
Jorge Gustavo Rocha , qgis-user@lists.osgeo.org, Greg Troxel
, Dan <19dm...@gmail.com> Subject: Re: [Qgis-user] Trimble
GeoXT 2005 Accuracy
I'm resending this without the map as there is a size limit. The moderator
may let it through, I hope . . . .
Hi Nicolas, thanks for your observations. I'll try to answer your
questions, please see the attached map, especially Map A:
Note that my previous email contained information for Map B; Map A is based
on the Municipal Cadastre (NAD83 UTM zone 10N) and illustrates the problem as
well. Other locations based on Provincial Monuments and/or the Municipal
cadastre (not illustrated here) have yielded similar results.
How many “known” points have you tested?
2 in this case, Maps A and B
Also several other locations with similar results
How where those point position calculated.
From the Municipal cadastre, visible in Map A
Make sure the coordinates are in the right CRS
NAD 83 UTM 10N used throught. See workflow in previous email
When converting from the monument’s CRS to NAD83 UTM zone 10N, are you
using the correct grid files?
[No monuments in this example] These were brought into QGIS from the
Municipal GCM database CSV (NAD83(CSRS) 3.0.0.BC.1.CRD) and reprojected by QGIS
to EPSG:26910 - NAD83 / UTM zone 10N
Presumably QGIS would choose the correct grid files
Municipal Cadastre is NAD83 UTM zone 10Np, li { white-space: pre-wrap;
}
Find a geodesic point that is in the middle of a field or on the side of
a highway with no obstacles.
Map A is open sky
Make sur your observations will be done when the constellation is well
distributed in the sky
As you probably know, TerraSync provides for PDOP, HDOP, SNR and
Horizon masks to preclude collecting poor quality positions. These were set
towards the "Precision" end of the scale
What post processing techniques are you using? How far is the base
station from your unit?
Real time was SBAS or RTCM; Post processing using the Pathfinder Office
differential correction engine, baseline about 30 km
How long are the observations? Have you tried other methods of post
processing like PPP?
Logging interval is 5 sec; 33 to 2037 positions per point
Did not use PPP. This is a test of mapping best practices, not geodesy
Have you contacted Trimble?
Yes, no response
Have you looked on there site to see if there is a software update
(firmware) for the unit or the post processing software?
Yes, receiver firmware is the latest, PFO and Terrasync are older but
compatible
Thanks Nicolas. If I have missed something, I hope someone can point it
out, I've tried to cover all the bases based on my training and experience.
-
Cheers, Spring
On 08/Mar/2021 06:40, Nicolas Cadieux
wrote:
Hi Harrison,
How many “known” points have you tested? How where those
point position calculated. They could be off. If you are using
state geodesic monuments, try to find the documented precision
of the monument. States have different types of monuments, some
are very old and have different standards. Make sure the
geodesic point is not the problem. Make sure the coordinates
are in the right CRS. As an example, if the coordinates are
published in NAD83 original but you are assuming NAD83(CSRS),
Re: [Qgis-user] Trimble GeoXT 2005 Accuracy
Hi Springer.I assume your raw data files are being converted to gpx on a computer since the raw terrasync files are proprietary binary files. If you are using trimble pathfinder, you can post process differentialy correct the data if you have access to base station logged at the same time you captured your field data. Having a base station 100 miles away will not improve your results as the baseline is too long.I do not know if you can write a gpx file directly from pathfinder but I would not bother. I would write a shapefile which will contain the coordinate system you specify. Simply open in qgis and you should be good to go. If your older unit works better, I would expect it may be an issue with the setup within pathfinder or perhaps the software version.I think your consistent offset is a direct result of how you are converting your data from trimble to gpx.As I mentioned in my previous comments, there are many issues which affect accuracy. Just because the box says it is accurate you will rarely replicate that in the field. In terms of WAAS dataframes, these are processed internally on your field unit.Kirk Schmidt Sent from my Galaxy Original message From: Springfield Harrison Date: 2021-03-07 5:57 a.m. (GMT-04:00) To: kirk , Jorge Gustavo Rocha , qgis-user@lists.osgeo.org, Greg Troxel , Dan <19dm...@gmail.com>, Nicolas Cadieux Subject: Re: [Qgis-user] Trimble GeoXT 2005 Accuracy Hello All, Thanks for the comments, I'll reply more fully tomorrow. The receiver is Trimble mapping grade: "The GeoExplorer 2005 series consists of: • The GeoXH™ handheld, providing subfoot (30 cm) accuracy, or even 8-inch (20 cm) accuracy with the optional Zephyr™ antenna. • The GeoXT™ handheld offering submeter accuracy for GIS data collection and data maintenance. • The GeoXM™ handheld with 1–3 meter GPS accuracy for mobile GIS applications." "Post processed carrier accuracy: 1-30cm". This receiver was probably $5-8000 (?) new. Data collection was stationary, open sky, good satellite coverage, several minutes of 5 sec observations, good PDOP SBAS and/or post processed The concern is not the accuracy as such, but the systematic NW shift. This has been observed over several months, consistently. My old Trimble ProXR (1994?, $20K new!) is actually better in this regard than the GeoXT! The Trimble manuals make no mention of the SBAS CRS, implying "turn it on and go, the receiver will integrate the SBAS into the rover file." More tomorrow, thanks . . . . . - Cheers, Spring On 06/Mar/2021 15:56, kirk wrote: A few notes. sbas which is waas in north America is based on equatorial satellites which will get you in the 1 m range in southern Canada. you can achieve sub decimeter accuracy consistently using rtk ,either through a ntrip caster (base station) broadcasting over the intenet or with your own base station and a radio link. there are a few chip sets and break out boards that you acquire and assemble your own system. This is a very inexpensive option. Another option in Canada is to use precise point positioning (PPP) which requires 6 to 12 hours of observation data using L1, L2 or L1 and L2 data and rinex log files. This comes in handy if you need to establish a remote base station. A proper antennae with a metal ground plane is also critical to getting quality results. Observing under a forested canopy is difficult especially in summer under leaf on conditions, after a rain which creates multiparth mayhem. There is a reason survey grade equipment is relatively expensive. If you require repeatably accurate results in a variety of conditions this is an option. Kirk Schmidt Sent from my Galaxy Original message From: Jorge Gustavo Rocha Date: 2021-03-06 6:41 p.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: Re: [Qgis-user] Trimble GeoXT 2005 Accuracy Hi, I just jump in this thread to say I'm really impressed with Ardusimple. I have a RTK Handheld Surveyor Kit [1] for +- 400 € and it works really well. I use the national NTRIP service and I have consistently precisions around 10 cm with just one receiver
Re: [Qgis-user] Trimble GeoXT 2005 Accuracy
A few notes.sbas which is waas in north America is based on equatorial satellites which will get you in the 1 m range in southern Canada. you can achieve sub decimeter accuracy consistently using rtk ,either through a ntrip caster (base station) broadcasting over the intenet or with your own base station and a radio link. there are a few chip sets and break out boards that you acquire and assemble your own system. This is a very inexpensive option.Another option in Canada is to use precise point positioning (PPP) which requires 6 to 12 hours of observation data using L1, L2 or L1 and L2 data and rinex log files. This comes in handy if you need to establish a remote base station.A proper antennae with a metal ground plane is also critical to getting quality results.Observing under a forested canopy is difficult especially in summer under leaf on conditions, after a rain which creates multiparth mayhem.There is a reason survey grade equipment is relatively expensive. If you require repeatably accurate results in a variety of conditions this is an option.Kirk Schmidt Sent from my Galaxy Original message From: Jorge Gustavo Rocha Date: 2021-03-06 6:41 p.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: Re: [Qgis-user] Trimble GeoXT 2005 Accuracy Hi, I just jump in this thread to say I'm really impressed with Ardusimple. I have a RTK Handheld Surveyor Kit [1] for +- 400 € and it works really well. I use the national NTRIP service and I have consistently precisions around 10 cm with just one receiver. I use a free Android application called SW Maps [2]. My survey points, tracks and photos are collect in a geopackage that I can read in QGIS. I use it mostly to collect ground control points for my drone flights. Regards, Jorge Gustavo [1] https://www.ardusimple.com/product/rtk-handheld-surveyor-kit/ [2] https://play.google.com/store/apps/details?id=np.com.softwel.swmaps Às 20:12 de 06/03/21, Greg Troxel escreveu: Springfield Harrison writes: Thanks Dan. See my relies to Kirk and Greg. The Emlid sounds interesting, will have a look. I have an earlier Emlid Reach (not RS or RS2), which has L1 only, and I never got it to work well. Also look at the Ardusimple unit -- but it's more a parts kit than a system. You need a way to get RTK reference data in, and a good antenna. One approach is Vespucci (OSM editor for Android) as a datalogger, and the Ardusimple WiFi NTRIP master to get corrections over the phone's hotspot. https://www.ardusimple.com/product/simplertk2b/ ___ 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 -- Email Signature Geomaster Jorge Gustavo Rocha | Software Engineer e: j...@geomaster.pt | m: +351 910 333 888 g: 41.54094,-8.40490 | v: 510 906 109 a: Rua António Cândido Pinto, 67, 4715-400 Braga ___ 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] Trimble GeoXT 2005 Accuracy
That is the limitation of the unit you have. Without using differential correction or real time kinematic techniques, that is about as good as it gets. You can improve your position a bit by keeping the unit static and averaging the 1 hz position data for at least 2 minutes. This will give you an average of >120 observations.. A clear view of the sky we also help improve the accuracy. For post processing work, you will need unit that records raw rinex format data which ghe garmin units typically do not.Kirk Schmidt Sent from my Galaxy Original message From: Springfield Harrison Date: 2021-03-06 4:35 a.m. (GMT-04:00) To: qgis-user Subject: [Qgis-user] Trimble GeoXT 2005 Accuracy Hello All, I recently acquired a Trimble GeoXT 2005 Series and am puzzled by the results it produces: Compared to a variety of "known" points, it consistently records positions that appear to be in error by 1.2 - 1.5 m NW from the known point. Points are collected and then mapped in QGIS as NAD83, UTM Zone 10 N. The known points include property survey pins, Government control survey monuments, Total Station survey points derived from the above, other GPS results (Trimble ProXRS) and identifiable points on orthophotos. I'm using SBAS correction in the GeoXT. It appears to be adding a consistent offset to the GPS result although no offset has been set in TerraSync. Many thanks for any thoughts on this situation . . . . . - Cheers, Springfield Harrison ___ 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] Lidar download
Hi Jean: I am not sure what you are looking for? Do you want to download a lidar derived image (such as a surface model) or do you want to work with raw point cloud data? Kirk Schmidt On 3/4/2021 8:41 AM, jean Lukusa wrote: Hi all , how Can I download Lidar image from Qgis 3.18.0 ? Jean Lukusa B.Sc. Geology Trained on Jorc Code GIS Operator +243812156128 ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] contours and symbology
Hi Saber: You can do this directly in symbology. See link below, Use Rule Based Symbology and "Elev" % 100= 0 to separate all the contours that are multiples of 100. Set your line symbols and you should be good to go. https://gis.stackexchange.com/questions/259875/styling-contours-at-specific-numbers Kirk Schmidt On 3/3/2021 2:15 PM, Saber Razmjooei wrote: Hi Azzurra, You can use Raster styling to display it as contours with minor and major contours: https://www.qgis.org/en/site/forusers/visualchangelog314/index.html#feature-raster-layer-contour-renderer <https://www.qgis.org/en/site/forusers/visualchangelog314/index.html#feature-raster-layer-contour-renderer> Kind regards Saber On Wed, 3 Mar 2021 at 18:01, Azzurra Lentini <mailto:azzurralent...@gmail.com>> wrote: Hi to all, I have extracted contours from a dtm raster 20 m resolution file and the result is a vector file of contours with 20 m of interval between the contours. I would like to have two different symbols for these contours: for example black colour and thickness 1 for the contours with interval of 100 m (0, 100, 200, 300 etc.) and grey colour and thickness 0.5 for the contours 20, 40, 60, 80, 120, 140 .. etc.). What can I do? Normally I create another contours file only for the contours equidistant 100 m and I put a different symbol for the 2 layers. Do you know other ways without the creation of another file? maybe rule-based symbology? Thank you, Azzurra -- Lecturer GIS University "Roma Tre" Consultant Environmental Risk Prevention and Hydrogeology AZZURRA LENTINI ++ Italy Mobile Tel.: **(39) 338 24 40 676 ++ SKYPE azzurrahydro ++ *_azzurralent...@gmail.com <mailto:azzurralent...@gmail.com>_* *_azzurra.lent...@uniroma3.it <mailto:azzurra.lent...@uniroma3.it> _**__* ++ *Par respect pour l'environnement,* *n'imprimez ce mail qu'en cas d'absolue nécessité* ___ 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> -- Saber Razmjooei www.lutraconsulting.co.uk <http://www.lutraconsulting.co.uk> ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] raster with distance values from a point
Hi Nicholas: If you activate the various providers under the Processing toolbox options, you can have direct access to the various algorithms I mentioned. I often generate a quick general model, export as a python script and then refine the script as required. I am running 3.16.2-Hannover. Kirk On 2/25/2021 9:26 PM, Nicolas Cadieux wrote: Interesting Kirk, The QGIS vector to raster conversion is not available in the graphical modeler or do you have another reason behind your solution? I don’t use much the graphical modeler so I appreciate your incite. Cheers! Nicolas Cadieux https://gitlab.com/njacadieux <https://gitlab.com/njacadieux> Le 25 févr. 2021 à 15:44, Kirk Schmidt a écrit : To add to Nicolas" comment: 1) Create your point layer from a csv (easting, northing id) as a minimum 2) If you want to stay with process modeler, use the SAGA>Raster creation tools> Rasterize Algorithm to convert the point layer to your required resolution and size 3) Use the GDAL>Raster Analysis>Proximity or SAGA>Raster tools> Proximity Raster to create your output. Kirk Schmidt On 2/25/2021 4:24 PM, Nicolas Cadieux wrote: Hi, You can cut down you use of several plugins by doing the following, 1) Create a vector point file with a single point having the numerical field containing the value 1. This file could be a csv containing WKT. 2) Use the vector to raster conversion. (This will take your vector point file, an create a raster of x,y dimensions. You point will become 1 and all other values will be 0. 3) Use the Proximity raster distance gdal function. This requires no plugins so it should be easy to create a graphical modeler. My favorite way would be to go to strait Python and to call the gdal proximity function from my Python script. If you had hundred or thousands of these to do, I could develop a Python script that could automate the entire thing. Contact me if you have questions. Nicolas Cadieux On 2021-02-25 12:58 p.m., Fabio Falchi wrote: Dear all, not having received any answer, I repeat the request. I need to speed up the process of obtaining rasters with distances from a single point (or a single lat-lon coordinate) in it. Manually I use this procedure: With a georeferenced raster in a UTM projection I use raster calculator to obtain a raster with all pixels set to zero. Then, with Serval plug-in, I manually set to 1 the value of the pixel in the wanted position (either by using a satellite overlay or by using a point vector layer with the appropriate coordinates). Then I use the proximity function under Raster>Analisys menu. I want to make this via the Graphical Modeler, in order to be able to make the process fast and skip the manual introduction with Serval plug-in. At the very least, I need a mode to put the value 1 in the pixel of the selected null raster in the position givel by lat-lon coordinates. Thank you! Best, Fabio Falchi /USC - Universidade de Santiago de Compostela/ /ISTIL - Light Pollution Science and Technology Institute/ Il sabato 20 febbraio 2021, 18:04:14 CET, Fabio Falchi ha scritto: Dear all, I need to speed up the process of obtaining rasters with distances from a single point in it. Manually I use this procedure: With a georeferenced raster in a UTM projection I use raster calculator to obtain a raster with all pixels set to zero. Then, with Serval plug-in, I manually set to 1 the value of the pixel in the wanted position (either by using a satellite overlay or by using a point vector layer with the appropriate coordinates). Then I use the proximity function under Raster>Analisys menu. I want to make this via the Graphical Modeler, in order to be able to make the process fast and skip the manual introduction with Serval plug-in. At the very least, I need a mode to put the value 1 in the pixel of the selected null raster in the position givel by lat-lon coordinates. Thank you! Best, Fabio Falchi /USC - Universidade de Santiago de Compostela/ /ISTIL - Light Pollution Science and Technology Institute/ ___ 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 -- Nicolas Cadieux https://gitlab.com/njacadieux ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email:k...@nortekresources.com Web:www.nortekresources.com ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.o
Re: [Qgis-user] raster with distance values from a point
To add to Nicolas" comment: 1) Create your point layer from a csv (easting, northing id) as a minimum 2) If you want to stay with process modeler, use the SAGA>Raster creation tools> Rasterize Algorithm to convert the point layer to your required resolution and size 3) Use the GDAL>Raster Analysis>Proximity or SAGA>Raster tools> Proximity Raster to create your output. Kirk Schmidt On 2/25/2021 4:24 PM, Nicolas Cadieux wrote: Hi, You can cut down you use of several plugins by doing the following, 1) Create a vector point file with a single point having the numerical field containing the value 1. This file could be a csv containing WKT. 2) Use the vector to raster conversion. (This will take your vector point file, an create a raster of x,y dimensions. You point will become 1 and all other values will be 0. 3) Use the Proximity raster distance gdal function. This requires no plugins so it should be easy to create a graphical modeler. My favorite way would be to go to strait Python and to call the gdal proximity function from my Python script. If you had hundred or thousands of these to do, I could develop a Python script that could automate the entire thing. Contact me if you have questions. Nicolas Cadieux On 2021-02-25 12:58 p.m., Fabio Falchi wrote: Dear all, not having received any answer, I repeat the request. I need to speed up the process of obtaining rasters with distances from a single point (or a single lat-lon coordinate) in it. Manually I use this procedure: With a georeferenced raster in a UTM projection I use raster calculator to obtain a raster with all pixels set to zero. Then, with Serval plug-in, I manually set to 1 the value of the pixel in the wanted position (either by using a satellite overlay or by using a point vector layer with the appropriate coordinates). Then I use the proximity function under Raster>Analisys menu. I want to make this via the Graphical Modeler, in order to be able to make the process fast and skip the manual introduction with Serval plug-in. At the very least, I need a mode to put the value 1 in the pixel of the selected null raster in the position givel by lat-lon coordinates. Thank you! Best, Fabio Falchi /USC - Universidade de Santiago de Compostela/ /ISTIL - Light Pollution Science and Technology Institute/ Il sabato 20 febbraio 2021, 18:04:14 CET, Fabio Falchi ha scritto: Dear all, I need to speed up the process of obtaining rasters with distances from a single point in it. Manually I use this procedure: With a georeferenced raster in a UTM projection I use raster calculator to obtain a raster with all pixels set to zero. Then, with Serval plug-in, I manually set to 1 the value of the pixel in the wanted position (either by using a satellite overlay or by using a point vector layer with the appropriate coordinates). Then I use the proximity function under Raster>Analisys menu. I want to make this via the Graphical Modeler, in order to be able to make the process fast and skip the manual introduction with Serval plug-in. At the very least, I need a mode to put the value 1 in the pixel of the selected null raster in the position givel by lat-lon coordinates. Thank you! Best, Fabio Falchi /USC - Universidade de Santiago de Compostela/ /ISTIL - Light Pollution Science and Technology Institute/ ___ 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 -- Nicolas Cadieux https://gitlab.com/njacadieux ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] raster with distance values from a point
Hallo Fabio: What version of QGIS are you using? Kirk Schmidt On 2/25/2021 1:58 PM, Fabio Falchi wrote: Dear all, not having received any answer, I repeat the request. I need to speed up the process of obtaining rasters with distances from a single point (or a single lat-lon coordinate) in it. Manually I use this procedure: With a georeferenced raster in a UTM projection I use raster calculator to obtain a raster with all pixels set to zero. Then, with Serval plug-in, I manually set to 1 the value of the pixel in the wanted position (either by using a satellite overlay or by using a point vector layer with the appropriate coordinates). Then I use the proximity function under Raster>Analisys menu. I want to make this via the Graphical Modeler, in order to be able to make the process fast and skip the manual introduction with Serval plug-in. At the very least, I need a mode to put the value 1 in the pixel of the selected null raster in the position givel by lat-lon coordinates. Thank you! Best, Fabio Falchi /USC - Universidade de Santiago de Compostela/ /ISTIL - Light Pollution Science and Technology Institute/ Il sabato 20 febbraio 2021, 18:04:14 CET, Fabio Falchi ha scritto: Dear all, I need to speed up the process of obtaining rasters with distances from a single point in it. Manually I use this procedure: With a georeferenced raster in a UTM projection I use raster calculator to obtain a raster with all pixels set to zero. Then, with Serval plug-in, I manually set to 1 the value of the pixel in the wanted position (either by using a satellite overlay or by using a point vector layer with the appropriate coordinates). Then I use the proximity function under Raster>Analisys menu. I want to make this via the Graphical Modeler, in order to be able to make the process fast and skip the manual introduction with Serval plug-in. At the very least, I need a mode to put the value 1 in the pixel of the selected null raster in the position givel by lat-lon coordinates. Thank you! Best, Fabio Falchi /USC - Universidade de Santiago de Compostela/ /ISTIL - Light Pollution Science and Technology Institute/ ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Symbology blank
Do your data show up spatially at all on the screen?Sent from my Galaxy Original message From: Ruby Jardine Date: 2021-02-17 7:44 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] Symbology blank Hi, I am trying to create a standard global map on the prevalence of maternal health issues. I have added my datasets and joined them, however when I go to the ‘Symbology’ section to colour code the map, it is blank. Do you know why this is and how I can change this to be able to create my map. Thank you ___Qgis-user mailing listQgis-user@lists.osgeo.orgList info: https://lists.osgeo.org/mailman/listinfo/qgis-userUnsubscribe: 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] distance between rectangle polygon in Qgis
Another option is to convert the polygons to lines and the calculate the length if each line in the attribute table.Sent from my Galaxy Original message From: Marin Mirosevic Date: 2021-02-15 6:41 p.m. (GMT-04:00) To: Nicolas Cadieux Cc: qgis-user@lists.osgeo.org Subject: Re: [Qgis-user] distance between rectangle polygon in Qgis Hi Nicolas Yes that is what i mean Dohvati Outlook za iOS Šalje: Nicolas Cadieux Poslano: Monday, February 15, 2021 9:45:08 PM Prima: Marin Mirosevic Kopija: qgis-user@lists.osgeo.org Predmet: Re: [Qgis-user] distance between rectangle polygon in Qgis Hi, I think you will need to delete duplicate nodes, then export the polygon nodes (vertices) to a new point layer. Then, I would apply a distance matrix. Write back if you want more details. This is for the distances between the nodes within a polygon right? Nicolas Cadieux https://gitlab.com/njacadieux Le 15 févr. 2021 à 12:35, Marin Mirosevic a écrit : Hello, Is there a way to calculate distance between points of the polygon ? I tried to find in “filed calculator “ under geometry but I found just $perimeter and that is whole bounding line of the polygon. With regards/ Vennlig hilsen Marin Mirosevic Consultant Mobil +47 909 67 511 marin.mirose...@bouvet.no Bouvet Norge AS bouvet.no ___ 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] MrSIS Import Reference File Problem
Hi Baz: Do not click on the sdw file but the actual MrSID image file. It will look for the world file and load appropriately. Kirk On 2/1/2021 8:08 PM, bazza_05.excite wrote: I am an absolute newbie to QGIS but have used other GIS programs. I have tried to import a MrSID rater file with a correlating world file (.sdw) by using the menu path: Layer - Add Layer - Add Raster Layer. I navigate to the folder where the file resides that also contains the .sdw and when I click on add I get a CRS error and the image does not go to its proper coordinates and defaults to near 0,0 The project reference is NAD 83 UTM Z12 N (ESPG 26912) which is the same as the image world file. Is there a correct method to insert the image file? I am on version 3.16.2 Cheers, Baz ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] simple conversion question
Hi Cole.You need to add the coordinates to your attribute table prior to exporting to csv format. First save the data in the correct coordinates system. Open the attribute table and add 2 new columns for the x and y coordinates. Update each volume using the $x abs $y formula and save the results prior to exporting to csv.Kirk SchmidtSent from my Samsung Galaxy smartphone. Original message From: Cole Burger Date: 2020-12-28 3:00 a.m. (GMT-04:00) To: qgis-user@lists.osgeo.org Subject: [Qgis-user] simple conversion question Hello,First time QGIS user. I'm using 3.16.1-Hannover on macOS Mojave 10.14.6. I'm trying to convert latitude/longitude coordinates to EPSG: 2163 US National Atlas Equal Area. The coordinates are here: https://raw.githubusercontent.com/cpb831/webscraper1/master/institutions%20to%20convert.csv. I would like to create a csv file with these coordinates in new projection. I followed the steps here: https://3drobotics.com/convert-lat-lon-coordinates-into-state-plane-with-qgis/. The new csv file just shows the original latitude/longitude coordinates. Any additional suggestions?Thank you!Sent with ProtonMail Secure Email.___ 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] Volume calculation - Reg
Hi Mike: I usually convert the contours or breaklines into a point layer and then generate a dem from that. This is not a precise methodology but it would depend on what you are trying to do and the precision of the underlying data. Kirk Schmidt On 12/18/2020 5:28 AM, Michael Dufty wrote: The links for Nicolas are quite useful, but I think to create a raster DEM from vector lines you probably need to use the interpolate triangulation with breaklines or structure lines. The dialogue takes a bit of getting used to, you need to choose the layer the data will come from, click + then in the box below choose whether the layer is being used for points, breaklines or structure lines. I’ve had major issues with the interpolation crashing QGIS if I try to used it with breaklines for complex imported data, but it seems to work OK for simple test geometries, and may be good for a simple excavation. I’d be interested to know if anyone has had success with interpolation using breaklines for more complex data, I’ve had to resort to other software trying to import DEM supplied as dxf 3d faces, but maybe I’m missing a basic step like fixing geometry or removing duplicates. *From:*Qgis-user *On Behalf Of *Nicolas Cadieux *Sent:* Friday, 18 December 2020 12:04 AM *To:* Abhi Shek *Cc:* qgis-user@lists.osgeo.org *Subject:* Re: [Qgis-user] Volume calculation - Reg This should help https://youtu.be/EerXfUf-aHo <https://youtu.be/EerXfUf-aHo> https://www.pointsnorthgis.ca/blog/how-to-calculate-surfacecut-volume-qgis-advanced/ <https://www.pointsnorthgis.ca/blog/how-to-calculate-surfacecut-volume-qgis-advanced/> Nicolas Cadieux https://gitlab.com/njacadieux <https://gitlab.com/njacadieux> Le 17 déc. 2020 à 04:57, Abhi Shek mailto:psn.abhis...@gmail.com>> a écrit : Sir, I invite your kind attention to the above subject. I am a beginner in qgis. I want to create a raster DEM using Vector Lines. In a nutshell, I want to calculate the planned excavated volume. please guide me in this matter. -- with warm regards, P S N Abhishek, ___ 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> ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Import mapinfo tabfile
Hi Jean: The tab file is what you need to open. The other files make up the rest of the components of the vector file similar to a shapefile. You can either drag and drop the tab file into and open QGIS canvas or use the add vector layer button. If you use the later, you can select the file type to Mapinfo and the suitable files will be displayed (*.TAB files). Kirk Schmidt On 12/17/2020 7:34 AM, Weijers, Jean (Kerkrade) wrote: Hello I got mapinfo files (.DAT, .ID, .MAP, .TAB). I want to import this into QGIS, but I don’t get a vectorfile. I see only TAB-file. How can I convert the tab-files into vectorfiles. Does anaone have a solution. Thanks in advance, Jean Weijers ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] [X] Re: Trying to use the v.split to break down a Polyline into 2 point segment to add feet label.
Hi Mark: The issue appears to be with the group permissions within your network that is creating the issue. The routine starts, run into a write issue, then shuts down. I suspect it is write permission to the local\temp folder. "This program is blocked by group policy. For more information, contact your system administrator." Kirk On 12/11/2020 3:29 PM, Wolf, Mark wrote: -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] R: Importing a raster into a Geopackage
Hi Nicolas:
I reached out to Alberto for the original data as the one on google
drive did not contain georeferencing and therefore would not write to a
geopackage. I tried to write the geopackafge from Alberto's data and I
generated 3 bands. I just got back from a meeting and I will play with
it some more.
Kirk
On 12/2/2020 12:10 PM, Nicolas Cadieux wrote:
Hi Alberto,
Can I share the file and your question in the gdal dev group? Perhaps
they will have an answer to your problem. You can also do that
yourself if you prefer.
Kirt, if you have time, can you do the same run with the data in my
google link? If you can’t duplicate the problem, then we will have a
better idea if you share your QGIS version and OS. Perhaps there is
change in the gdal version between both computers.
Nicolas Cadieux
https://gitlab.com/njacadieux <https://gitlab.com/njacadieux>
Le 2 déc. 2020 à 10:11, Alberto Vavassori
a écrit :
Hi Kirk,
I tried with another raster file that represents a temperature
distribution (continuous values). In this case, when I import it into
the Geopackage, I see the same issue that you have just described:
minimum and maximum values in the legend change, but it seems that
the actual values remain unchanged, thus correct. Unfortunately this
is not the case of the land cover raster (discrete set of values)
that I have shared with you: in this case, the raster characteristics
change, as we have all noticed.
Just another question, hoping that you might have an answer: is it
possible to store the symbology of a raster file inside the
Geopackage? I know that it is possible for vector data, by clicking
on “Style > Save style” in the “Symbology” tab. But I cannot find a
way to save the symbology in case of raster files.
Thanks a lot.
Alberto
*Da:* Qgis-user *Per conto di
*Kirk Schmidt
*Inviato:* mercoledì 2 dicembre 2020 16:00
*A:* Nicolas Cadieux ;
qgis-user@lists.osgeo.org
*Oggetto:* Re: [Qgis-user] R: Importing a raster into a Geopackage
Hi all:
I ran the r.stats to get a list of each pixel value on my test data
(A small DEM) before and after conversion to a geopackage (see second
image below). The values and summary stats are identical for the two
rasters but the default legend entry has a different min and maximum
value. So, I can conclude that the data is being written correctly
but the default display parameters for the geopackage are different
than a geotiff.
Kirk
On 12/1/2020 3:24 PM, Nicolas Cadieux wrote:
Hi,
I changed the link so you have access to more files.
https://drive.google.com/drive/folders/1fMl9odowUFFRQFK61rTxgHeom_D8ymN6?usp=sharing
<https://drive.google.com/drive/folders/1fMl9odowUFFRQFK61rTxgHeom_D8ymN6?usp=sharing>
Nicolas
On 2020-12-01 1:58 p.m., Kirk Schmidt wrote:
Hi Nicolas
I did the same and the legend indicates a change magnitude
but if I subtract the geopackage raster from the original,
all pixels equal zero.
Kirk
On 12/1/2020 2:36 PM, Nicolas Cadieux wrote:
Hi,
Tried this:
gdal_translate -of GPKG -b 1
C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif
C:/Users/Nicolas/Desktop/test.gpkg
I read 4 band and still get data modifications.
Nicolas
On 2020-12-01 1:28 p.m., Kirk Schmidt wrote:
Hi Nicolas:
try the -b 1 parameter. I just tried this and wrote
a single band raster from a 3 channel tiff.
Kirk
On 12/1/2020 2:03 PM, Nicolas Cadieux wrote:
I tried
gdal_translate -of GPKG -co "BAND_COUNT=1"
C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif
C:/Users/Nicolas/Desktop/TEST.gpkg for the
creation and
C:\>gdal_translate -of GPKG -co "TILE_FORMAT =
TIFF"
C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif
C:/Users/Nicolas/Desktop/TEST.gpkg
Same problem. I wonder if we need to change the
reading options...
https://gdal.org/drivers/raster/gpkg.html
<https://gdal.org/drivers/raster/gpkg.html>
"Opening options
By default, the driver will expose a GeoPackage
dataset as a four band (Red,Green, Blue,Alpha)
dataset, which gives the maximum compatibility
with the various encodings of tiles that can be
stored. It is possible to specify an explicit
number of bands with the BAND_COUNT opening option."
The problem could be when QGIS Reads the file,
not when it writes
Re: [Qgis-user] R: Importing a raster into a Geopackage
Hi all: I ran the r.stats to get a list of each pixel value on my test data (A small DEM) before and after conversion to a geopackage (see second image below). The values and summary stats are identical for the two rasters but the default legend entry has a different min and maximum value. So, I can conclude that the data is being written correctly but the default display parameters for the geopackage are different than a geotiff. Kirk On 12/1/2020 3:24 PM, Nicolas Cadieux wrote: Hi, I changed the link so you have access to more files. https://drive.google.com/drive/folders/1fMl9odowUFFRQFK61rTxgHeom_D8ymN6?usp=sharing Nicolas On 2020-12-01 1:58 p.m., Kirk Schmidt wrote: Hi Nicolas I did the same and the legend indicates a change magnitude but if I subtract the geopackage raster from the original, all pixels equal zero. Kirk On 12/1/2020 2:36 PM, Nicolas Cadieux wrote: Hi, Tried this: gdal_translate -of GPKG -b 1 C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif C:/Users/Nicolas/Desktop/test.gpkg I read 4 band and still get data modifications. Nicolas On 2020-12-01 1:28 p.m., Kirk Schmidt wrote: Hi Nicolas: try the -b 1 parameter. I just tried this and wrote a single band raster from a 3 channel tiff. Kirk On 12/1/2020 2:03 PM, Nicolas Cadieux wrote: I tried gdal_translate -of GPKG -co "BAND_COUNT=1" C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif C:/Users/Nicolas/Desktop/TEST.gpkg for the creation and C:\>gdal_translate -of GPKG -co "TILE_FORMAT = TIFF" C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif C:/Users/Nicolas/Desktop/TEST.gpkg Same problem. I wonder if we need to change the reading options... https://gdal.org/drivers/raster/gpkg.html "Opening options By default, the driver will expose a GeoPackage dataset as a four band (Red,Green, Blue,Alpha) dataset, which gives the maximum compatibility with the various encodings of tiles that can be stored. It is possible to specify an explicit number of bands with the BAND_COUNT opening option." The problem could be when QGIS Reads the file, not when it writes the file. I explains why we see 4 bands but not why values are changing. I would make a bug report or ask the question on the gdal dev group to get Even Rouault's attention. Nicolas On 2020-12-01 12:17 p.m., Sebastian Gutwein wrote: This looks like the same issue you are having. github.com/qgis/QGIS/issues/38227 <http://github.com/qgis/QGIS/issues/38227> According to the documentation linked in the issue you need to set the band count to 1 otherwise it defaults to 4. gdal.org/drivers/raster/gpkg.html <http://gdal.org/drivers/raster/gpkg.html> I have not tried this to see if it works. On Tue, Dec 1, 2020 at 11:52 AM Nicolas Cadieux <mailto:njacadieux.git...@gmail.com>> wrote: Hi, Looking at it now. I see what you see. I will try a few things. Nicolas On 2020-12-01 11:42 a.m., Alberto Vavassori wrote: Hi Nicolas, thank you for the answer. I send the raster that I am interested in. Alberto *Da:*Nicolas Cadieux <mailto:njacadieux.git...@gmail.com> *Inviato:* martedì 1 dicembre 2020 17:40 *A:* Alberto Vavassori <mailto:alberto.vavass...@polimi.it> *Cc:* qgis-user@lists.osgeo.org <mailto:qgis-user@lists.osgeo.org> *Oggetto:* Re: [Qgis-user] Importing a raster into a Geopackage Hi, Send a sample data set if you can't figure it out. Nicolas On 2020-12-01 10:35 a.m., Nicolas Cadieux wrote: Hi, Perhaps this can help: https://gdal.orgdrivers/raster/gpkg.html <https://gdal.orgdrivers/raster/gpkg.html> “Various kind of input datasets can be converted to GeoPackage raster : ·Single band grey level ·Single band with R,G,B or R,G,B,A color table ·Two bands: first band with grey level, second band with alpha channel ·Three bands: Red, Green, Blue ·Four band: Red, Green, Blue, Alpha GeoPackage rasters only support Byte data type.” My guess is that if the original raster does not correspond to one of the categories above, it will be converted to the closest format. Nicolas Cadieux https://gitlab.com/njacadieux <https://gitlab.com/njacadieux> Le 1 déc. 2020 à 09:58, Alberto Vavassori <mailto:alberto.vavass...@polimi.it> a écrit : Hello everybody. I have an issue regarding the Geopackage format. I have tried to import a singleband raster file (it represents the land cover, thus every pixel has a single integer value corresponding to the land cover class), but this operation seems to modify the raster characteristics: the raste
Re: [Qgis-user] R: Importing a raster into a Geopackage
Hi Nicolas I did the same and the legend indicates a change magnitude but if I subtract the geopackage raster from the original, all pixels equal zero. Kirk On 12/1/2020 2:36 PM, Nicolas Cadieux wrote: Hi, Tried this: gdal_translate -of GPKG -b 1 C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif C:/Users/Nicolas/Desktop/test.gpkg I read 4 band and still get data modifications. Nicolas On 2020-12-01 1:28 p.m., Kirk Schmidt wrote: Hi Nicolas: try the -b 1 parameter. I just tried this and wrote a single band raster from a 3 channel tiff. Kirk On 12/1/2020 2:03 PM, Nicolas Cadieux wrote: I tried gdal_translate -of GPKG -co "BAND_COUNT=1" C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif C:/Users/Nicolas/Desktop/TEST.gpkg for the creation and C:\>gdal_translate -of GPKG -co "TILE_FORMAT = TIFF" C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif C:/Users/Nicolas/Desktop/TEST.gpkg Same problem. I wonder if we need to change the reading options... https://gdal.org/drivers/raster/gpkg.html "Opening options By default, the driver will expose a GeoPackage dataset as a four band (Red,Green, Blue,Alpha) dataset, which gives the maximum compatibility with the various encodings of tiles that can be stored. It is possible to specify an explicit number of bands with the BAND_COUNT opening option." The problem could be when QGIS Reads the file, not when it writes the file. I explains why we see 4 bands but not why values are changing. I would make a bug report or ask the question on the gdal dev group to get Even Rouault's attention. Nicolas On 2020-12-01 12:17 p.m., Sebastian Gutwein wrote: This looks like the same issue you are having. github.com/qgis/QGIS/issues/38227 <http://github.com/qgis/QGIS/issues/38227> According to the documentation linked in the issue you need to set the band count to 1 otherwise it defaults to 4. gdal.org/drivers/raster/gpkg.html <http://gdal.org/drivers/raster/gpkg.html> I have not tried this to see if it works. On Tue, Dec 1, 2020 at 11:52 AM Nicolas Cadieux mailto:njacadieux.git...@gmail.com>> wrote: Hi, Looking at it now. I see what you see. I will try a few things. Nicolas On 2020-12-01 11:42 a.m., Alberto Vavassori wrote: Hi Nicolas, thank you for the answer. I send the raster that I am interested in. Alberto *Da:*Nicolas Cadieux <mailto:njacadieux.git...@gmail.com> *Inviato:* martedì 1 dicembre 2020 17:40 *A:* Alberto Vavassori <mailto:alberto.vavass...@polimi.it> *Cc:* qgis-user@lists.osgeo.org <mailto:qgis-user@lists.osgeo.org> *Oggetto:* Re: [Qgis-user] Importing a raster into a Geopackage Hi, Send a sample data set if you can't figure it out. Nicolas On 2020-12-01 10:35 a.m., Nicolas Cadieux wrote: Hi, Perhaps this can help: https://gdal.orgdrivers/raster/gpkg.html <https://gdal.orgdrivers/raster/gpkg.html> “Various kind of input datasets can be converted to GeoPackage raster : ·Single band grey level ·Single band with R,G,B or R,G,B,A color table ·Two bands: first band with grey level, second band with alpha channel ·Three bands: Red, Green, Blue ·Four band: Red, Green, Blue, Alpha GeoPackage rasters only support Byte data type.” My guess is that if the original raster does not correspond to one of the categories above, it will be converted to the closest format. Nicolas Cadieux https://gitlab.com/njacadieux <https://gitlab.com/njacadieux> Le 1 déc. 2020 à 09:58, Alberto Vavassori <mailto:alberto.vavass...@polimi.it> a écrit : Hello everybody. I have an issue regarding the Geopackage format. I have tried to import a singleband raster file (it represents the land cover, thus every pixel has a single integer value corresponding to the land cover class), but this operation seems to modify the raster characteristics: the raster becomes a 4-bands raster and the pixels values change accordingly (basically, the original raster is saved in the Geopackage as an RGB image). It seems to be a bug of QGIS. I have also tried to use GDAL translate through the processing toolbox of QGIS to convert the TIFF file to a Geopackage file, but the result is the same. Has anybody faced this problem? May anybody suggest a way to work this problem out? Thank you. Alberto ___ Qgis-user mailing list Qgis-user@lists.osgeo.org <mailto:Qgis-user@lists.osgeo.
Re: [Qgis-user] R: Importing a raster into a Geopackage
Hi Nicolas: try the -b 1 parameter. I just tried this and wrote a single band raster from a 3 channel tiff. Kirk On 12/1/2020 2:03 PM, Nicolas Cadieux wrote: I tried gdal_translate -of GPKG -co "BAND_COUNT=1" C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif C:/Users/Nicolas/Desktop/TEST.gpkg for the creation and C:\>gdal_translate -of GPKG -co "TILE_FORMAT = TIFF" C:/Users/Nicolas/Desktop/Land_cover_Corine2018.tif C:/Users/Nicolas/Desktop/TEST.gpkg Same problem. I wonder if we need to change the reading options... https://gdal.org/drivers/raster/gpkg.html "Opening options By default, the driver will expose a GeoPackage dataset as a four band (Red,Green, Blue,Alpha) dataset, which gives the maximum compatibility with the various encodings of tiles that can be stored. It is possible to specify an explicit number of bands with the BAND_COUNT opening option." The problem could be when QGIS Reads the file, not when it writes the file. I explains why we see 4 bands but not why values are changing. I would make a bug report or ask the question on the gdal dev group to get Even Rouault's attention. Nicolas On 2020-12-01 12:17 p.m., Sebastian Gutwein wrote: This looks like the same issue you are having. github.com/qgis/QGIS/issues/38227 <http://github.com/qgis/QGIS/issues/38227> According to the documentation linked in the issue you need to set the band count to 1 otherwise it defaults to 4. gdal.org/drivers/raster/gpkg.html <http://gdal.org/drivers/raster/gpkg.html> I have not tried this to see if it works. On Tue, Dec 1, 2020 at 11:52 AM Nicolas Cadieux mailto:njacadieux.git...@gmail.com>> wrote: Hi, Looking at it now. I see what you see. I will try a few things. Nicolas On 2020-12-01 11:42 a.m., Alberto Vavassori wrote: Hi Nicolas, thank you for the answer. I send the raster that I am interested in. Alberto *Da:*Nicolas Cadieux <mailto:njacadieux.git...@gmail.com> *Inviato:* martedì 1 dicembre 2020 17:40 *A:* Alberto Vavassori <mailto:alberto.vavass...@polimi.it> *Cc:* qgis-user@lists.osgeo.org <mailto:qgis-user@lists.osgeo.org> *Oggetto:* Re: [Qgis-user] Importing a raster into a Geopackage Hi, Send a sample data set if you can't figure it out. Nicolas On 2020-12-01 10:35 a.m., Nicolas Cadieux wrote: Hi, Perhaps this can help: https://gdal.orgdrivers/raster/gpkg.html <https://gdal.orgdrivers/raster/gpkg.html> “Various kind of input datasets can be converted to GeoPackage raster : ·Single band grey level ·Single band with R,G,B or R,G,B,A color table ·Two bands: first band with grey level, second band with alpha channel ·Three bands: Red, Green, Blue ·Four band: Red, Green, Blue, Alpha GeoPackage rasters only support Byte data type.” My guess is that if the original raster does not correspond to one of the categories above, it will be converted to the closest format. Nicolas Cadieux https://gitlab.com/njacadieux <https://gitlab.com/njacadieux> Le 1 déc. 2020 à 09:58, Alberto Vavassori <mailto:alberto.vavass...@polimi.it> a écrit : Hello everybody. I have an issue regarding the Geopackage format. I have tried to import a singleband raster file (it represents the land cover, thus every pixel has a single integer value corresponding to the land cover class), but this operation seems to modify the raster characteristics: the raster becomes a 4-bands raster and the pixels values change accordingly (basically, the original raster is saved in the Geopackage as an RGB image). It seems to be a bug of QGIS. I have also tried to use GDAL translate through the processing toolbox of QGIS to convert the TIFF file to a Geopackage file, but the result is the same. Has anybody faced this problem? May anybody suggest a way to work this problem out? Thank you. Alberto ___ 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> -- Nicolas Cadieux https://gitlab.com/njacadieux <https://gitlab.com/njacadieux> -- Nicolas Cadieux https://gitlab.com
Re: [Qgis-user] Importing a raster into a Geopackage
Hello Alberto: I just tried saving a dem (single band raster) as a geopackage using QGIS 3.10.1. The only issue I can think of is if you select rendered vs raw data in the export dialog (see top of screen shot below). Can you share the version you are using? Kirk On 12/1/2020 10:56 AM, Alberto Vavassori wrote: Hello everybody. I have an issue regarding the Geopackage format. I have tried to import a singleband raster file (it represents the land cover, thus every pixel has a single integer value corresponding to the land cover class), but this operation seems to modify the raster characteristics: the raster becomes a 4-bands raster and the pixels values change accordingly (basically, the original raster is saved in the Geopackage as an RGB image). It seems to be a bug of QGIS. I have also tried to use GDAL translate through the processing toolbox of QGIS to convert the TIFF file to a Geopackage file, but the result is the same. Has anybody faced this problem? May anybody suggest a way to work this problem out? Thank you. Alberto ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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] Does QGIS3 freezes more often than QGIS2 ?
Hi Idan: I usually go through the "add raster Layer" button and it seems to work seamlessly with small and large rasters in 3.10.1 working on windows 8.1x 64. I wonder if it is an issue with the drag and drop function and or related to raster files size? Kirk On 11/26/2020 1:18 PM, Idan Miara wrote: Hi, I was wondering if I'm the only one who experiences sudden freezes in QGIS3 (on Win10x64) It's hard to pinpoint exactly what triggers these freezes but usually it happens more often when I drag and drop raster files from the file explorer onto the QGIS canvas. Usually after a few seconds it unfreezes and I can continue my work, otherwise I just force close and reopen QGIS. This almost never happens to me with QGIS 2.17 and keeps happening on every QGIS3 version I tested (including the latest). Is there any easy way to debug this? Idan ___ 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 -- Kirk Schmidt, MScF, BScF, RPF General Manager Nortek Resource Solutions Inc. RR # 1 Thorburn, NS B0K 1W0 Tel (902) 922.3607 Email: k...@nortekresources.com Web: www.nortekresources.com ___ 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
