4.33333

4.3 | 3 ratings Rate this file 85 Downloads (last 30 days) File Size: 16.7 KB File ID: #35548
image thumbnail

3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students

by Kyung Sang Kim

 

10 Mar 2012 (Updated 24 Apr 2013)

3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes

| Watch this File

File Information
Description

Hello,

My name is Kyung-Sang Kim.

This program is with respect to the mesh-grid based 3D cone-beam CT.

This is made for students who learn the medical imaging.

Please let me know if the program has problems.
(E-mail: kssigari@gmail)

Thank you

Required Products MATLAB
MATLAB release MATLAB 7.8 (R2009a)
Tags for This File  
Everyone's Tags
backprojection(3), cbct, cone beam, conebeam, ct, fdk, mlem, projection, reconstruction
Tags I've Applied
Add New Tags Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (7)
04 Apr 2013 Kyung Sang Kim

Hi, Masoud Hashemi,
This program is for flat detectors,
I think you should modify for the arc detector geometry,
Thanks for using my software.

30 Mar 2013 Masoud Hashemi

I am wondering if the projection function is for Arc detectors (3rd Generation) or Flat detectors?

02 Mar 2013 Tai Chieh

Great!

16 Dec 2012 ljpzgx

The codes can run but the result is not very good,can you upload the referrences you refer for the codes,so tha we can understand clearly. Thank you!

13 Mar 2012 Kyung Sang Kim

This program is similar to other programs but
much faster.

To accelerate speed,
(1) define the reconstruction object and projection images as a "single" precision,
and then (2) use "matlabpool".

12 Mar 2012 Kyung Sang Kim

You can use "matlabpool",

From

for i=1:n (angle)
projection or backprojection
end

To

parfor i=1:n (angle)
projection or backprojection
!!please delete figure plot functions
end

This increases the performance more than twice in quad cores.

12 Mar 2012 Kyung Sang Kim

Please change the function code "projection_pixel.m":

From

dist = (SAD+SDD)./sqrt((SAD+SDD)^2 + uu.^2 + vv.^2)*abs(ys(2)-ys(1));

To

dist = sqrt((SAD+SDD)^2 + uu.^2 + vv.^2)./(SAD+SDD)*abs(ys(2)-ys(1));

Thank you, :)

Updates
20 Mar 2012

I optimized the source codes

16 Nov 2012

Fixed a bug

13 Dec 2012

I changed title, tag

14 Jan 2013

Change description

24 Apr 2013

I modified codes simply.

Contact us