Debugging different scenarios
@@ -5,7 +5,7 @@ from porespy.tools import extend_slice
from tqdm import tqdm
-def extract_pore_network(im, pore_regions=None, dt=None, voxel_size=1,
+def extract_phase_network(im, pore_regions=None, dt=None, voxel_size=1,
dual=False, solid_regions=None):
r"""
Analyzes an image that has been partitioned into pore regions and extracts
@@ -68,6 +68,7 @@ def extract_pore_network(im, pore_regions=None, dt=None, voxel_size=1,
raise Exception('Please provide solid regions map for' +
' dual newtwork extraction')
+ if solid_regions is not None:
merge = sp.amax(pore_regions)
solid_regions = solid_regions + merge
p_on_s = pore_regions*(~im) # Expose pores labels on solid
