About Site Map Submit Contact Us Log in | Create an account
Create an account Log In
Average Rating
User Rating:
Visitors Rating:
My rating:

Write review
See full specifications

scripts default iconCT Projection Simulator (Scripts) Publisher's description

Specifying an object as a shape results in fewer edge artifacts in the reconstructed image than with projections from radon.m

Specifying an object as a shape results in fewer edge artifacts in the reconstructed image than with projections from radon.m, which specifies the object as an image.

A comparison of reconstructed images with CTProjSim and radon.m is given at

The object starts as empty and then shapes with specified dimensions and density can be added.
The projections are the product of the length through the instances of all shape times their density.
Different instances of each shape can have different density.
The density can be a vector quantity so the projections are a
(ndetectors, nangles, ndensity) matrix.

The shapes are specified by optional parameters of one or more pairs
'shape',spec where shape can be one of 'ellipse', 'gaussian', or 'polygon'
the spec parameter depends on the shape--see the function help
The lines are equally spaced around the origin and start parallel to the y-axis.

Here is an example of a Shepp-Logan phantom with a square in the middle

% projections for Shepp-Logan phantom
% the spec is for phantom in air .. not water
s = [ ...
0,0,0.92,0.69,90,2; ...
0,-0.0184,0.874,0.6624,90,-0.98; ...
0.22,0,0.31,0.11,72,-0.02; ...
-0.22,0,0.41,0.16,108,-0.02; ...
0,0.35,0.25,0.21,90,0.01; ...
0,0.1,0.046,0.046,0,0.01; ...
0,-0.1,0.046,0.046,0,0.01; ...
-0.08,-0.605,0.046,0.023,0,0.01; ...
0,-0.605,0.023,0.023,0,0.01; ...
0.06,-0.605,0.046,0.023,90,0.01; ...

% add a small square at the center
% polygon vertexes are specified as complex vectors: x + iy
% exp(1i*pi/4) rotates the square by 45 degrees
pgon.vertex = 0.1*exp(1i*linspace(0,2*pi,5))*exp(1i*pi/4);
pgon.density = -0.1;

pjn = CTProjSim(0.0025,240,'ellipse',s,'polygon',pgon);

img = CTrecon(pjn); % see http://aprendtech.com/wordpress/?p=74 for CTrecon code
% can also use iradon but it will be slow
% nangles = 240;
% angles = linspace(0,180,nangles+1);
% img = iradon(pjn,angles(1:end-1));
imshow(img,[0.8 1.2])
for more information see

System Requirements:

MATLAB 7.5 (R2007b)
Program Release Status: New Release
Program Install Support: Install and Uninstall

CT Projection Simulator (Scripts) Tags:

Click on a tag to find related softwares

Is CT Projection Simulator (Scripts) your software?

Manage your software

Most Popular

scripts default icon ASK, OOK, FSK, BPSK, QPSK, 8PSK modulation 1.1
ASK, OOK, FSK, BPSK, QPSK, 8PSK modulation contain several functions for digital modulation simulation
scripts default icon Simulink Communication Labs 1.1
Simulink Communication Labs allows you to learn communication systems in greater depth.
scripts default icon M-QAM modulation and demodulation 1.1
M-QAM modulation and demodulation is the QAM modulation and demodulation tech.
scripts default icon LZW Compression/Decompression 1.1
LZW Compression/Decompression - Updated LZW compressor and decompressor with reasonable performance
scripts default icon InSPIRE utility to plot a 2D displacement field (Scripts) 1.0
This program plots the deformation field (displace vectors) contained in vector.txt.