Search Results

fiona): import fiona, geopandas, time #A line shapefile, total size 118 MB shape = r"/home/bera/Desktop/GIStest/ml_riks.shp" #Fiona open start = time.process_time() c = fiona.open(shape) if c.schema["geometry … "]=='LineString': print("Line") print(time.process_time()-start) #Line #0.0016448840000009568 #Geopandas checking geometry type of first feature start = time.process_time() df = geopandas.read_file …

Use the UpdateCursor for this with the SHAPE@X and SHAPE@Y tokens. This will change the data so backup before executing: import arcpy feature_class=r'C:\Default.gdb\Points' #Change to match your dat …

I'm not sure what you mean by true centroid but this will work on "normal" centroids. I use the SHAPE@X and SHAPE@Y tokens to get the centroid coordinates of each feature and then add these in two col …

Having large complex geometries will slow down processing. Try: Multipart To Singleparts Subdivide 3. And Create Spatial Index On both input layers

Code below will output a new feature class with the vertices intersecting the points removed. import arcpy precision = 1 #number of decimals to round coordinates when matching points to polygon vert …

polylines: arcpy.env.workspace = r'C:\GIS\data\testdata' #I have my polygon shapefile in this folder polygon_fc = "ak_riks" polygon_fields = [f for f in arcpy.ListFields(polygon_fc) if f.type not in ['Geometry …

Intersect with output_type=POINT POINT —Point intersections will be returned. If the inputs are line or polygon, the output will be a multipoint feature class.

You can use ST_MakeEnvelope with the intersect operator (&&) to select features in a rectangle and ST_AsText to make the geometries readable: SELECT *, st_astext(geom) FROM water WHERE water. …

Draw a new polygon covering the gap with some overlap with the old one Select this and your old one "Merge selected features"

format(lyr.crs().authid()), "simplified_line", "memory") provider = vl.dataProvider() for feature in lyr.getFeatures(): #For each line in the jagged line layer geom = feature.geometry() #Fetch the geometry … linestring from the thinned lines newline.addGeometry(newsubline) #Add this linestring to the empty multilinestring f = QgsFeature() #Create a new feature f.setGeometry(newline) #Set the geometry …

overlay_intersects( 'layer2' ) will return True if it does, else False which will become 1/0 for a integer field:

You can create a virtual layer where you join the layer to itself based on a distance between two points. Replace pnt with your layername select a.* from pnt a --Select everything from the point layer …

= gpd.points_from_xy(x=df["lon"], y=df["lat"], crs=4326) df = gpd.GeoDataFrame(data=data, geometry=geometry, crs=4326) df = df.to_crs(3829) #You need to have a projected coordinate system with meters … endy"])]), axis=1) df = df.set_geometry("line").drop(columns="geometry").rename_geometry("geometry").set_crs(3829) df.to_file(r"C:\GIS\data\testdata\lines.shp") …

Add an if statement: fc = r'C:\Users\tomlin01\Documents\ArcGIS\Default.gdb\Export_Output_8' #layer to update z_increase = -0,155 #change in meters with arcpy.da.UpdateCursor(fc, 'SHAPE@Z', …

This is now very easy with the .reverse method: import shapely line = shapely.wkt.loads("LINESTRING (0 0, 1 1, 3 1, -3 3)") multiline = shapely.wkt.loads("MULTILINESTRING ((120 105, 140 98),(160 103, …