I had the same issue and I discovered my clipped raster sampling was set to 'nearest neighbor' by default. I went to Layer Properties -> Display -> Resample and changed it to 'Bilinear Interpolation (for continuous data)' and it looks much better now.
The reason you're getting this error is because ListFeatureClasses does not accept a path as an argument, the wildcard parameter is used to filter the feature classes returned to only those matching a particular file name pattern, for example 'Edit*' will return only the feature classes that start with edit.
You need to set your arcpy.env.workspace first, ...
Your where clause is using double quotes for the string literal.
The OGR SQL documentation specifies single quotes are required for strings. Double quotes are used for field names (only required if they contain special characters or reserved words):
String literals and identifiers quoting
Strict SQL92 rules are applied regarding string literals and
As @user30184 has mentioned, gdalwarp can perform exactly what you are looking for and can be called from the command line with the appropriate operations flags.
If you are working from python (as has been tagged) and all your data is in the proper coordinate systems, you might find the rasterio and fiona command line tools a bit simpler to handle. The ...
clipped[i] returns a list with one element. You meant clipped[[i]] with double brackets. Also clipped[i] within paste0() is meaningless.
Actually your use case is already covered by the package. Assuming that your shapefile has an attribute PLOTID (replace with something else if needed)
ctg <- readLAScatalog("folder/plot_a3.las", filter = "-...
1: Union to generate all possible geometries
2: Compare all geometries using .equals and keep the one with highest priority
3: Dissolve and Multipart To Singelpart
from itertools import combinations
layer_to_process = QgsProject.instance().mapLayersByName('Polygons') #Change Polygons to match your layer name
prio_fieldname = 'prio' #1 is kept ...
I try to give you a solution using a sample data set of three polygons. Polygons have a unique ID and a PRIORITY columns.
(polygon fill is 50% transparent to see the overlapping parts).
Convert the polygons to lines (we'll use these lines to split the polygon layer)
using Vector/Geometry Tools/Polygons to Lines... from the menu. You'll get layer Lines ...