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)