neurosynth/base/imageutils.py
def get_sphere(coords, r=4, vox_dims=(2, 2, 2), dims=(91, 109, 91)):
    """ # Return all points within r mm of coordinates. Generates a cube
    and then discards all points outside sphere. Only returns values that
    fall within the dimensions of the image."""
    r = float(r)
    xx, yy, zz = [slice(-r / vox_dims[i], r / vox_dims[
                        i] + 0.01, 1) for i in range(len(coords))]
    cube = np.vstack([row.ravel() for row in np.mgrid[xx, yy, zz]])
    sphere = cube[:, np.sum(np.dot(np.diag(
        vox_dims), cube) ** 2, 0) ** .5 <= r]
    sphere = np.round(sphere.T + coords)
    return sphere[(np.min(sphere, 1) >= 0) &
                  (np.max(np.subtract(sphere, dims), 1) <= -1), :].astype(int)
def get_sphere(coords, r=4, vox_dims=(2, 2, 2), dims=(91, 109, 91)):
    """ # Return all points within r mm of coordinates. Generates a cube
    and then discards all points outside sphere. Only returns values that
    fall within the dimensions of the image."""
    r = float(r)
    xx, yy, zz = [slice(-r / vox_dims[i], r / vox_dims[
                        i] + 0.01, 1) for i in range(len(coords))]
    cube = np.vstack([row.ravel() for row in np.mgrid[xx, yy, zz]])
    sphere = cube[:, np.sum(np.dot(np.diag(
        vox_dims), cube) ** 2, 0) ** .5 <= r]
    sphere = np.round(sphere.T + coords)
    return sphere[(np.min(sphere, 1) >= 0) &
                  (np.max(np.subtract(sphere, dims), 1) <= -1), :].astype(int)
def get_sphere ( coords , r = 4 , vox_dims = ( 2 , 2 , 2 ) , dims = ( 91 , 109 , 91 ) ) : r = float ( r ) xx , yy , zz = [ slice ( - r / vox_dims [ i ] , r / vox_dims [ i ] + 0.01 , 1 ) for i in range ( len ( coords ) ) ] cube = np . vstack ( [ row . ravel ( ) for row in np . mgrid [ xx , yy , zz ] ] ) sphere = cube [ : , np . sum ( np . dot ( np . diag ( vox_dims ) , cube ) ** 2 , 0 ) ** .5 <= r ] sphere = np . round ( sphere . T + coords ) return sphere [ ( np . min ( sphere , 1 ) >= 0 ) & ( np . max ( np . subtract ( sphere , dims ) , 1 ) <= - 1 ) , : ] . astype ( int )
# Return all points within r mm of coordinates. Generates a cube
    and then discards all points outside sphere. Only returns values that
    fall within the dimensions of the image.
# Return all points within r mm of coordinates . Generates a cube and then discards all points outside sphere . Only returns values that fall within the dimensions of the image .