def test_backproject():
    N=50
    x=modified_shepp_logan((N,N,3))[:,int(N/2),1]
    x=numpy.array(x)
    rbp=ref.back_project(x)
    obp=ocv.back_project(x)
    assert is_equal(rbp,obp)
def test_backproject():
    N=50
    x=modified_shepp_logan((N,N,3))[:,int(N/2),1]
    x=numpy.array(x)
    rbp=ref.back_project(x)
    ibp=ispmd.back_project(x)
#    pylab.figure()
#    pylab.subplot(121)
#    pylab.imshow(ibp)
#    pylab.colorbar()
#    pylab.subplot(122)
#    pylab.imshow(rbp)
#    pylab.colorbar()
#    pylab.show()
    assert is_equal(rbp,ibp)