https://github.com/jyhmiinlin/pynufft
Tip revision: 562fe915b1338e8f4b3211a607b99c6883866c3b authored by Jyh-Miin Lin on 05 September 2022, 04:23:49 UTC
experiment torch, tf?
experiment torch, tf?
Tip revision: 562fe91
test_tf.py
# import torch
import scipy.sparse
import numpy as np
def test_tf():
import pynufft
import pkg_resources
DATA_PATH = pkg_resources.resource_filename('pynufft', 'src/data/')
# PHANTOM_FILE = pkg_resources.resource_filename('pynufft', 'data/phantom_256_256.txt')
import numpy
# import matplotlib.pyplot
import scipy
import scipy.misc
image = scipy.misc.ascent()[::2,::2]
image=image.astype(float)/numpy.max(image[...])
Nd = (256, 256) # image space size
Kd = (512, 512) # k-space size
Jd = (6,6) # interpolation size
# load k-space points
om = numpy.load(DATA_PATH+'om2D.npz')['arr_0']
A = pynufft.NUFFT_tf()
A.plan(om, Nd, Kd, Jd)
y_tf = A.forward(image)
x_tf = A.adjoint(y_tf)
# print(y_tf.shape)
# print(x_tf.shape)
nfft = pynufft.NUFFT() # CPU
# print(nfft.processor)
nfft.plan(om, Nd, Kd, Jd)
y = nfft.forward(image)
x2 = nfft.adjoint(y)
print('Forward error between tf and numpy',np.linalg.norm(y_tf - y)/np.linalg.norm(y))
print('Adjoint Error between tf and numpy', np.linalg.norm(x2 - x_tf)/np.linalg.norm(x2))
# test_torch()
# test_tensorflow()
# test_torch_class()
if __name__ == '__main__':
test_tf()
# test_random_sp()