contour_raster_rasterio

contour_raster_rasterio.py

# From https://stackoverflow.com/questions/41487642/how-to-export-contours-created-in-scikit-image-find-contours-to-shapefile-or-geo

import rasterio as rio

flpth = '/Users/seb/Documents/Codes/VolcGIS/volcanoes/Agung/_hazard/Tephra/Agung_VEI3_P5.tif'
schema = {"geometry": "Polygon", "properties": {"value": "int"}}

maskVal = 1000

with rio.open(flpth) as src:
    image = src.read()
    # use your function to generate mask
    mask = image >= maskVal
    # and convert to uint8 for rio.features.shapes
    mask = mask.astype('uint8')
    shapes = rio.features.shapes(mask, transform=src.transform)
    # select the records from shapes where the value is 1,
    # or where the mask was True
    records = [{"geometry": geometry, "properties": {"value": value}}
               for (geometry, value) in shapes if value == 1]
    
    # Write to shapefile
    # with fiona.open('test.shp', "w", "ESRI Shapefile",
    #                 crs=src.crs.data, schema=schema) as out_file:
    #     out_file.writerecords(records)
        
    # Create GeoDataFrame
    clipMask = gpd.GeoDataFrame(records[0]['geometry'])