From 1dceabe4e2aa8868f65ea635df0ba6f12f513bcc Mon Sep 17 00:00:00 2001 From: Cesar-Timana Date: Tue, 2 May 2023 15:33:34 -0300 Subject: [PATCH] Surface area used in interface area calculation (high accuracy) needs to be with voxel_size = 1 --- porespy/networks/_getnet.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/porespy/networks/_getnet.py b/porespy/networks/_getnet.py index 5ab27d354..cef11a6b9 100644 --- a/porespy/networks/_getnet.py +++ b/porespy/networks/_getnet.py @@ -248,7 +248,7 @@ def regions_to_network(regions, phases=None, voxel_size=1, accuracy='standard'): net['pore.volume'] = region_volumes(regions=im, mode='marching_cubes') areas = region_surface_areas(regions=im, voxel_size=voxel_size) net['pore.surface_area'] = areas - interface_area = region_interface_areas(regions=im, areas=areas, + interface_area = region_interface_areas(regions=im, areas=areas/voxel_size**2, voxel_size=voxel_size) A = interface_area.area net['throat.cross_sectional_area'] = A