Rank: 1264 based on 106 downloads (last 30 days) and 1 file submitted
photo

Kyung Sang Kim

E-mail
Company/University
KAIST

Personal Profile:

Kyung Sang Kim


 

Watch this Author's files

 

Files Posted by Kyung Sang Kim
Updated   File Tags Downloads
(last 30 days)
Comments Rating
25 Oct 2014 Screenshot 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: Kyung Sang Kim backprojection, conebeam, ct, projection, fdk, mlem 106 16
  • 4.66667
4.7 | 9 ratings
Comments and Ratings by Kyung Sang Kim View all
Updated File Comments Rating
11 Nov 2013 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: Kyung Sang Kim

Hi, Erick,

That's a good point.

This factor is used in the weighted backprojection for fan-beam and cone-beam geometries.
(see Jiang Hsieh 2nd edition, Figs 3.35 (reason) and 3.37(geometry))

Result is similar compared to the previous one, In real geometry, this factor is not critical factor for the image, because CT rotates 360 deg (average factor ~ 1) as well as itself is close to 1.

However, I think I omitted it for fast calculation, soon I will fix it for accuracy.

Thanks a lot.

ps) you can directly send e-mail to me, it is sometimes difficult to check comments here.

28 Oct 2013 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: Kyung Sang Kim

Hi Erick,

First Question:

If you have real projections ("RealProj"), then first calculate:

1. I0 = max(RealProj(:)); % if you don't know I0
2. Proj = -log(RealProj/I0);

Here, "I0" is the blank image (intensity of source) as mentioned in the text book.
If you think that the final image has a problem, check 1) "RealProj" has 0 values or 2) "I0" is too high than background values.

Solution of cases

1) step 3. Proj(isinf(Proj)) = 0;

2) step 1. Define I0 as a background (high intensity) average.
step 2. Proj = -log(min(RealProj/I0,1));

Then, you can execute FDK or MLEM as in Demo.m.

Second Question:

Parameter setting is very important. Image resolution is defined by detector size and distances of DSO, DSO. In your case, DSD is 700 and DSO is 560. so ratio is 560/700 = 0.8. Then detector size of a pixel is 50/1024, so maximum resolution of a voxel is 50/1024*0.8.

Now, voxel resolution "dx=dy=dz" is defined. then define the number of voxels that covers whole images (nx, ny, nz).
Thus, sx = nx*dx, sy=ny*dy, sz=nz*dz;

Last Question:
I modified just for fast calculation. But if you can see Jiang Hsieh's book (second edition,page 96), it is easy to understand.

Thanks,

04 Apr 2013 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: 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.

13 Mar 2012 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: 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 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: 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.

Comments and Ratings on Kyung Sang Kim's Files View all
Updated File Comment by Comments Rating
12 Jun 2014 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: Kyung Sang Kim mat

29 May 2014 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: Kyung Sang Kim Hamdi Mahdjoub

28 May 2014 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: Kyung Sang Kim Bella

Would it be possible to use this program to forward project a 3D image where dx=dy, but dx=dy!=dz? Right now, my parameters are:

%% Parameter setting %%

param.nx = 400;
param.ny = 400;
param.nz = 66;

%The real detector panel pixel density (number of pixels)
param.nu = 400;
param.nv = 400;

% Detector setting, according to Varian Trilogy OBI (real size)
param.su = 350; % mm
param.sv = 350; % mm

% X-ray source and detector setting
param.DSD = 1500; % Distance source to detector
param.DSO = 1000; % X-ray source to object axis distance

% This confusing stuff
res=(param.DSO/param.DSD)*(param.su/param.nu);
% (1000/1500)*(350/400)
param.sx = param.nx*res; % mm
param.sy = param.ny*res; % mm
param.sz = param.nz*res; % mm

but this doesn't create the correct forward projections, so far as I have been able to check. Am I trying to do something that is impossible in image processing?

Thank you!

24 May 2014 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: Kyung Sang Kim Rama Teja A S

Please include analytical equation for better understanding. Or please explain the algorithm in projection.m. If possible please include analytical equation references

11 Nov 2013 3D Cone beam CT (CBCT) projection backprojection FDK MLEM reconstruction MATLAB codes for students 3D Cone beam CT (CBCT) projection backprojection MLEM FDK reconstruction MATLAB source codes Author: Kyung Sang Kim Kyung Sang Kim

Hi, Erick,

That's a good point.

This factor is used in the weighted backprojection for fan-beam and cone-beam geometries.
(see Jiang Hsieh 2nd edition, Figs 3.35 (reason) and 3.37(geometry))

Result is similar compared to the previous one, In real geometry, this factor is not critical factor for the image, because CT rotates 360 deg (average factor ~ 1) as well as itself is close to 1.

However, I think I omitted it for fast calculation, soon I will fix it for accuracy.

Thanks a lot.

ps) you can directly send e-mail to me, it is sometimes difficult to check comments here.

Contact us