New answers tagged

0

For example use ModelBuilder with Iterate Feature Selection: Each extraction is named by an attribute in grid feature class


0

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.


-1

Well, I finally figured that out: https://github.com/piLeoni/fishNet2Layouts


2

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, ...


1

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 ...


2

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 ...


2

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 = "-...


0

You can: 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')[0] #Change Polygons to match your layer name prio_fieldname = 'prio' #1 is kept ...


0

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 ...


Top 50 recent answers are included