def METHOD_NAME(self):
    pn = op.network.Cubic(shape=[5, 5, 1])
    pn.add_model_collection(
        op.models.collections.geometry.spheres_and_cylinders)
    pn.regenerate_models()
    im = op.visualization.generate_voxel_image(network=pn,
                                               pore_shape='sphere',
                                               throat_shape='cylinder',
                                               max_dim=500)
    assert im.shape[0] == 500