region_surface_areas(regions, voxel_size=(1, 1, 1), strel=None)[source]

Extract the surface area of each region in a labeled image.

Optionally, it can also find the the interfacial area between all adjoining regions.

  • regions (ndarray) – An image of the pore space partitioned into individual pore regions. Note that zeros in the image will not be considered for area calculation.

  • voxel_size (tuple) – The resolution of the image, expressed as the length of one side of a voxel, so the volume of a voxel would be voxel_size-cubed. The default is 1.

  • strel (array_like) – The structuring element used to blur the region. If not provided, then a spherical element (or disk) with radius 1 is used. See the docstring for mesh_region for more details, as this argument is passed to there.


areas – A list containing the surface area of each region, offset by 1, such that the surface area of region 1 is stored in element 0 of the list.

Return type:



Click here to view online example.