pc_to_satn(pc, im)[source]#

Converts an image of capillary entry pressures to saturation values

  • pc (ndarray) – A Numpy array with the value in each voxel indicating the capillary pressure at which it was invaded. In order to accomodateh the possibility of positive and negative capillary pressure values, uninvaded voxels should be indicated by +inf and residual phase by -inf. Solid vs void phase is defined by im which is mandatory.

  • im (ndarray) – A Numpy array with True values indicating the void space


satn – A Numpy array with each voxel value indicating the global saturation at which it was invaded.

Return type



If any -inf values are present the minimum saturation will start at a value greater than 0 since residual was present. If any +inf values are present the maximum saturation will be less than 1.0 since not all wetting phase was displaced.


Click here to view online example.