### Highlights from Three Dimensional Analytical Magnetic Resonance Imaging Phantom in the Fourier Domain

5.0
5.0 | 3 ratings Rate this file 7 Downloads (last 30 days) File Size: 2.32 MB File ID: #15697 Version: 1.1

# Three Dimensional Analytical Magnetic Resonance Imaging Phantom in the Fourier Domain

### Cheng Guan Koay (view profile)

25 Jul 2007 (Updated )

This software contains 2D and 3D Shepp-Logan phantoms in both the image and Fourier domains.

File Information
Description

This submission contains a collection of codes in different programming languages that implement the analytical framework proposed in

Koay CG, Sarlls JE, Ã–zarslan E.
* Three Dimensional Analytical Magnetic Resonance Imaging Phantom in the Fourier Domain. Magn Reson Med. 58: 430-436 (2007)

The article above can be obtained from:
http://dir2.nichd.nih.gov/nichd/stbb/publications.html

MATLAB release MATLAB 7.0.1 (R14SP1)
Other requirements None
04 Apr 2016 Pankaj Daga

### Pankaj Daga (view profile)

I think in the 3D Shepp Logan code, the rotation matrices are not generated correctly. The Rx matrix should be used instead of using the Rz matrix twice.

Comment only
27 Apr 2010 Cheng Guan Koay

### Cheng Guan Koay (view profile)

Sorry, I had problems with the comment submission. I accidentally deleted an earlier message that contained my explanations.

I shall try to explain it again here.

The first thing I should point out is that the field-of-view (FOV) in the image domain is fixed. Each FOV in each dimension is two units (arbitrary unit) from -1 to 1. Therefore, the resolution along the x-axis, denoted as del_x, in the image domain is determined by the number of samples, N, and it is given by del_x = FOV/N.

The resolution along the k-axis, denoted as del_kx, in the Fourier domain is determined by the FOV, and it is given by del_kx = 1/FOV.

The main thing is to know how wide one should sample the k-space or Fourier domain. The magnitude of k_max, i.e., |k_max| is given by del_kx * N/2 . Note that it is assumed that the negative and the positive regions are sampled equally.

So, "K<0.002" should not be the problem. Hope this helps. If not, feel free to email.

Comment only
27 Apr 2010 Cheng Guan Koay

### Cheng Guan Koay (view profile)

I forgot to say that there is no problem with K<0.002. :)

Comment only
27 Apr 2010 Giuseppe

### Giuseppe (view profile)

First of all really nice work!!! I have just a question the package work well for an image resolution equal or higher 128x128X128, for lower resolution the signal in kspace is not numerically unstable.......to avoid this problem is it sufficient to change the condition on K<0.002, or i have to do something else???

Thank you

12 Aug 2009 Cheng Guan Koay

### Cheng Guan Koay (view profile)

In the revised SampleTest.m, I used imtool and please make sure that you click on a menu that allows you to change contrast. min=0.0 and max=3.0 should be fine and these values should be applied to both imtool panels.

Comment only
11 Aug 2009 Cheng Guan Koay

### Cheng Guan Koay (view profile)

Due the inconvenience of updating and commenting on the site, I have decided to only update the site with major releases.

The MATLAB codes that illustrate the steps needed to perform FFT on the k-space signals of the 3D Shepp-Logan phantom can be found in a revised version of SampleTest.m, which is located on my google site: http://sites.google.com/site/hispeedpackets/Home/shepplogan

Comment only
01 Aug 2009 Jay Rod

### Jay Rod (view profile)

This package is very useful for testing MRI image reconstruction and simulation studies. 2D and 3D Shepp-Logan phantoms in the Fourier domain and image domain are critical to analysis in magnetic resonance imaging. Excellent extensible Java source codes with useful MATLAB, Mathematica, and IDL interfaces. Thanks a lot for sharing the source codes and your paper.

24 Oct 2008 John