def shepp2d(size=512, dtype='float32'):
    """
    Load test Shepp-Logan image array.

    Parameters