Hot answers tagged

3

Are you using the -n parameter? If not, nodata values will be included in your mosaic. GDAL will ignore nodata values if you specify this and will mosaic data without overlapping nodata areas included. Note (update from OP): np.nan may not work as a nodata value For example: gdal_merge.py -n 0 -o output.tif input1.tif input2.tif You will need to know the ...


3

I figured this out. First you need to create a new field merged_ID in the lot layer, giving a unique ID to each group of polygons to be merged. The attribute table would look like this: Then type the following query in the virtual layer creation window: select st_union(geometry), l.* from "lot" as l group by l."merged_ID"


3

Creating layers with hidden fields is the way to go here, using field info objects. I have a function I regularly use for this. Then it's a matter of creating your two layers and using them as inputs in your merge. Try: import arcpy fc1 = r"your\first\featureclass" fc2 = r"your\second\featureclass" mergeFc = r"your\output\merge\...


2

If you simply want to know the total area, you don't need to merge anything! Assuming you have a numeric field with the area this can be very easily summed. Projected data held as a file geodatabase featureclass would have a shape_area field which you can simply right click on and choose Statistics as shown below, you can obviously do this on any numeric ...


2

So the big issue is what "algorithm" to use. I'm doing pretty much the same as @Fezter with an overlay and checking overlaps. I take anything with an overlaps=1 and keep it as a final result. I take one feature at random with the highest number of overlaps, and drop it from the translation. So: Do Overlay Data If Overlaps = 1 then ...


2

I think I understand your question now, and I have a potential solution. I have created a simple workspace which buffers some points by 100m and then passes the points and buffers into a PointOnAreaOverlayer. Then, I test for how many overlaps there are. The workspace looks like this: The points and buffers look like this: Tester Parameters: Resulting ...


2

This is potentially less performant than Emil's answer, but much more concise: just merge everything indiscriminately, create 2 sets with the original field names, then delete the symmetric difference of those field names (belonging to either fc_a or fc_b, but not both). merged = arcpy.Merge_management([fc_a, fc_b], 'merged') fields_a = {field.name for ...


2

I think that as @MrXsquared tells you, you must first solve the problem of non data or null. This is very common, it depends on the type of data and how the raster was generated. You can search on this website to find out how to solve it. If your null value is recognized by QGis, for example -9999 (used by SAGA). You could perform the procedure with the ...


1

Solution 1: r.patch solved the problem if DSM and DEM files are inputted correctly into the algorithm. Unfortunately, changing the order of the files in the UI is not clear. You need to drag them in order in the selection window before clicking OK. Solution 2: As @MrXsquared described in a comment in the initial post, this problem occurs when "nodata ...


1

The both coordinate systems are essentially the same despite the name and the code. The software you are using obviously selects the first coordinate system that matches the proj string and shows that. projinfo -o proj epsg:8687 PROJ.4 string: +proj=utm +zone=33 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs +type=crs projinfo -o proj epsg:25833 PROJ....


1

For polygons in polygons I would suggest using either within or have their center in. Within will only select precincts that are within the county. Have their center in works well if you have messy polygon edges that overlap a lot, however, if you have extremely irregular polygons whos centroid may fall outside the polygon boundaries. Just something to ...


1

Another alternative to merge all the polygons into a single geometry is to use the Dissolve tool, which: Creates a new coverage by merging adjacent polygons, lines, or regions that have the same value for a specified item. As long as all the features share the same attribute (it seems that ID_SINAP would work in your case), it will merge all of them into a ...


1

Start editing, select all of them - Merge: And select one of them to take attributes from.


Only top voted, non community-wiki answers of a minimum length are eligible