Search Results

Add a new field to the point layer, same type as the polygon field. Go to Layer Properties - Attributes Form - Drag and drop designer, drag the newly added field from Available Widgets to Form Layout …

Find the total bounds (xmin, ymin, xmax, ymax), create two lists of x and y coordinates within them. List all combinations of the x and y coordinates, create points. import geopandas as gpd import num …

Use spatial join import geopandas as gpd import numpy as np poly = gpd.read_file(r'D:\temp_diverse\polygons.shp') point = gpd.read_file(r'D:\temp_diverse\points.shp') df = gpd.sjoin(poly, point[['ge …

I get the same count in geopandas and QGIS. Your code is a bit complicated, try using groupby: import geopandas as gpd dfpoint = gpd.read_file(filename=r'/home/bera/Desktop/GIStest/1000points.shp') d …

With SQL and a Virtual Layer you can try: select town, type, count(*) pointcount, geometry from pointlayer group by town, type

Intersect points with polygons Dissolve with statistics field unique: Unique—The number of unique values of the specified field will be counted

Try: Join attributes by location Takes an input vector layer and creates a new vector layer that is an extended version of the input one, with additional attributes in its attribute table.Takes an …

Reproject so the dataframes are in the same crs: Transform all geometries in an active geometry column to a different coordinate reference system I want both in EPSG:4326: print(poly.head(2)) # 0 P …

This is now very easy with .sample_points (see the Changelog Version 0.13 (May 6, 2023)) import geopandas as gpd df = gpd.read_file(r"/home/bera/Desktop/gistest/district.gpkg") ax = df.plot(figsize=( …

You can use the expression: array_majority( array:=overlay_intersects( layer:='points', expression:="rasterval" ) )[0] It will create an array of all values in the co …

You can use "Join attributes by location" tool: Takes an input vector layer and creates a new vector layer that is an extended version of the input one, with additional attributes in its attribute ta …

You can create a new field using Field Calculator to fetch the region name from the region layer. Then symbolize the point layer using the new field. Example: array_to_string(overlay_intersects('distr …

Try buffering your polygons with 0 distance, sometimes it fixes geometry errors. If you want to ignore the invalid geometries, you can set it Options - Processing - General:

I use this pyqgis code to create a fixed number of sample plots in forest stands/polygons. It will create a point layer in memory. Click/highlight your polygon layer, adjust the number of plots below …

You can create a virtual layer. The query will calculate the azimuth from each polygons centroid to all intersecting points. Replace "centroids" and "grid" with the names of your layers: select pnt.*, …