boxcount#

A method for measuring the fractal dimension of an image

import matplotlib.pyplot as plt

import porespy as ps

The arguments and their defaults for this function are:

import inspect

inspect.signature(ps.metrics.boxcount)
<Signature (im, bins=10)>
im = ps.generators.sierpinski_foam([3**5, 3**5], n=5, mode=None)
fig, ax = plt.subplots(1, 1, figsize=[6, 6])
ax.imshow(im, interpolation='none', origin='lower')
ax.axis(False);
../../../_images/234452bff2ff4c152a7a26c936eb6416ae6bca58fe75775d7ee40fd11cd35271.png

im#

The image which is to be analzyed. Can be 2D or 3D.

b = ps.metrics.boxcount(im=im)
print(b)
――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――
Results of boxcount generated at Sat Jul 12 16:19:25 2025
――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――
size                      Array of size (10,)
count                     [544, 291, 160, 77, 35, 25, 9, 9, 4, 1]
slope                     Array of size (10,)
――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――

The returned object is like a dataclass with each computed value stored in an attribute. The result can be printed for inspection. The results can also be plotted as follows:

fig, ax = plt.subplots(1, 2, figsize=[12, 6])
ax[0].loglog(b.size, b.count)
ax[0].set_xlabel('box length')
ax[0].set_ylabel('number of partially filled boxes')

ax[1].semilogx(b.size, b.slope)
ax[1].plot([0, 1000], [1.9, 1.9])
ax[1].set_xlabel('box length')
ax[1].set_ylabel('slope')
ax[1].set_ylim([0, 3]);
../../../_images/1a1b43dc22d672df30b6d823f78df67454f70d05e54d339d4cb9b401a8227439.png

bins#

The box sizes to use. The default is 10. If an integer is given it computes the range of box sizes. If an array is given, these are used directly

b = ps.metrics.boxcount(im=im, bins=20)
fig, ax = plt.subplots(1, 2, figsize=[12, 6])
ax[0].loglog(b.size, b.count)
ax[0].set_xlabel('box length')
ax[0].set_ylabel('number of partially filled boxes')

ax[1].semilogx(b.size, b.slope)
ax[1].plot([0, 1000], [1.9, 1.9])
ax[1].set_xlabel('box length')
ax[1].set_ylabel('slope')
ax[1].set_ylim([0, 3]);
../../../_images/ca946daba87782a11466d6006b21c920c87b7855cffa423eea9ca2d32bdcc199.png