Hi I want to eliminate polygons with a different attribute value (gridcode=0) within the larger polygons (Grid Code=1) to create clean polygons for further delineation of polygons within the extent of the larger polygons. I cannot use the eliminate tool because it is only available in the advanced v...

Visit the post for more.

I hope I'm using the right vocabulary here but basically what I want to do is have Modelbuilder read the gridcode compare that to an index that defines what that means and the proceeds to write that definition in the class field. I figure it's an iteration but thats where I'm stuck. So for example gridcode 21 = Developed, Open Space

I have a shapefile, with multiple features inside. I am trying to dissolve the nearby/adjacent features based on the shapefile particular attribute value (gridcode), but the difference shouldn't be greater than 3. If the difference between the gridcode is greater than 3, it shouldn't get dissolve. final output should be a shapefile with the dissolve features with an attribute inside having all the gridvalues which got dissolve. Image before dissolve image after dissolve import ogr from rtree import index from shapely.ops import unary_union from shapely.wkt import loads from fiona import collection shapefile = #read shp driver = ogr.GetDriverByName('ESRI Shapefile') shp_read = driver.Open(shapefile) layer = shp_read.GetLayer() featureCount = layer.GetFeatureCount() idx = index.Index(interleaved=False) for feat in layer: fid = feat.GetFID() geom = feat.GetGeometryRef() geom_poly_envelope = geom.GetEnvelope() xmin, xmax, ymin, ymax = geom_poly_envelope idx.insert(fid,(xmin, xmax, ymin, ymax)) print(idx) append_feature =[] for fid2 in range(0,layer.GetFeatureCount()): feature2 = layer.GetFeature(fid2) poly_collection=[] geometry2 = feature2.GetGeometryRef() poly_collection.append(loads(geometry2.ExportToWkt())) grid2= feature2.GetFieldAsInteger('gridcode') xmin, xmax, ymin, ymax = geometry2.GetEnvelope() for fid in list(idx.intersection((xmin, xmax, ymin, ymax))): feature1 = layer.GetFeature(fid) grid1= feature1.GetFieldAsInteger('gridcode') geometry1 = feature1.GetGeometryRef() if fid2!= fid: if geometry2.Intersects(geometry1): print ('{} intersects {}'.format(fid2, fid)) print(grid2,grid1) if (grid2-grid1=-5): poly_collection.append(loads(geometry1.ExportToWkt())) shp_union= unary_union(poly_collection) print("featurecount"+ str(fid2)) append_feature.append(shp_union) print("appended") union2= unary_union(append_feature)

I have a shapefile, with multiple features inside. I am trying to dissolve the nearby/adjacent features based on the shapefile particular attribute value (gridcode), but the difference shouldn't be greater than 3. If the difference between the gridcode is greater than 3, it shouldn't get dissolve. final output should be a shapefile with the dissolve features with an attribute inside having all the gridvalues which got dissolve. Image before dissolve image after dissolve import ogr from rtree import index from shapely.ops import unary_union from shapely.wkt import loads from fiona import collection shapefile = #read shp driver = ogr.GetDriverByName('ESRI Shapefile') shp_read = driver.Open(shapefile) layer = shp_read.GetLayer() featureCount = layer.GetFeatureCount() idx = index.Index(interleaved=False) for feat in layer: fid = feat.GetFID() geom = feat.GetGeometryRef() geom_poly_envelope = geom.GetEnvelope() xmin, xmax, ymin, ymax = geom_poly_envelope idx.insert(fid,(xmin, xmax, ymin, ymax)) print(idx) append_feature =[] for fid2 in range(0,layer.GetFeatureCount()): feature2 = layer.GetFeature(fid2) poly_collection=[] geometry2 = feature2.GetGeometryRef() poly_collection.append(loads(geometry2.ExportToWkt())) grid2= feature2.GetFieldAsInteger('gridcode') xmin, xmax, ymin, ymax = geometry2.GetEnvelope() for fid in list(idx.intersection((xmin, xmax, ymin, ymax))): feature1 = layer.GetFeature(fid) grid1= feature1.GetFieldAsInteger('gridcode') geometry1 = feature1.GetGeometryRef() if fid2!= fid: if geometry2.Intersects(geometry1): print ('{} intersects {}'.format(fid2, fid)) print(grid2,grid1) if (grid2-grid1=-5): poly_collection.append(loads(geometry1.ExportToWkt())) shp_union= unary_union(poly_collection) print("featurecount"+ str(fid2)) append_feature.append(shp_union) print("appended") union2= unary_union(append_feature)

I have a shapefile, with multiple features inside. I am trying to dissolve the nearby/adjacent features based on the shapefile particular attribute value (gridcode), but the difference shouldn't be greater than 3. If the difference between the gridcode is…

I have a shapefile, with multiple features inside. I am trying to dissolve the nearby/adjacent features based on the shapefile particular attribute value (gridcode), but the difference shouldn't be greater than 3. If the difference between the gridcode is…