def METHOD_NAME(self,idata,rdata,debug=False):
    output=np.zeros(idata.shape,idata.dtype)
    numvoxels=idata.shape[0]
    for voxel in range(0,numvoxels):
        blue=idata[voxel,:];
        uv=rdata[voxel,:];
        alpha=np.dot(uv,blue)/np.dot(uv,uv)
        residual=blue-alpha*uv
        output[voxel,:]=residual
    return output