# Why Are There Negative Voxel Values in My RTKFDK X-ray Reconstruction?

**URL:** https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498
**Category:** Algorithms
**Created:** [April 21, 2025, 2:06am UTC](https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498 "2025-04-21T02:06:47Z")
**Posts on this page:** 9
**Page:** 1

<div class="post-metadata">

### Author: ![zhangxiang](https://discourse.itk.org/letter_avatar_proxy/v4/letter/z/c77e96/32.png) [@zhangxiang](https://discourse.itk.org/u/zhangxiang)
#### Post date: [April 21, 2025, 2:06am UTC](https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498/1 "2025-04-21T02:06:47Z")

</div>

![image](https://discourse.itk.org/uploads/default/original/2X/1/1a5792856896e7d82d6893c0467bcad3661b0d14.jpeg)  
 ![image](https://discourse.itk.org/uploads/default/original/2X/e/eb2f3f03d58e6c65d1847e52385f63558f8d2c06.png)  
I’m working on 3D reconstruction using cone-beam X-ray projection images with the **FDK algorithm in RTK (Reconstruction Toolkit)**. After reconstruction, I noticed that **some voxel values in the output volume are negative**.

---

<div class="post-metadata">

### Author: ![zhangxiang](https://discourse.itk.org/letter_avatar_proxy/v4/letter/z/c77e96/32.png) [@zhangxiang](https://discourse.itk.org/u/zhangxiang)
#### Post date: [April 21, 2025, 2:07am UTC](https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498/2 "2025-04-21T02:07:43Z")

</div>

What I reconstruct is a head phantom.

---

<div class="post-metadata">

### Author: ![zhangxiang](https://discourse.itk.org/letter_avatar_proxy/v4/letter/z/c77e96/32.png) [@zhangxiang](https://discourse.itk.org/u/zhangxiang)
#### Post date: [April 21, 2025, 2:10am UTC](https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498/3 "2025-04-21T02:10:00Z")

</div>

@simon.rit

---

<div class="post-metadata">

### Author: ![zhangxiang](https://discourse.itk.org/letter_avatar_proxy/v4/letter/z/c77e96/32.png) [@zhangxiang](https://discourse.itk.org/u/zhangxiang)
#### Post date: [April 21, 2025, 2:37am UTC](https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498/4 "2025-04-21T02:37:10Z")

</div>

This is 3D reconstruction of a vertebra phantom.

 ![image](https://discourse.itk.org/uploads/default/original/2X/9/943795ca02c4d368ef71c347e1c95c1776fe5e00.png)  
 ![image](https://discourse.itk.org/uploads/default/original/2X/b/bf4628e71ba3c1e2531465dc9f5291d9d7dc65c9.jpeg)

---

<div class="post-metadata">

### Author: ![zhangxiang](https://discourse.itk.org/letter_avatar_proxy/v4/letter/z/c77e96/32.png) [@zhangxiang](https://discourse.itk.org/u/zhangxiang)
#### Post date: [April 21, 2025, 3:30am UTC](https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498/5 "2025-04-21T03:30:23Z")

</div>

I didn’t use -Hann -radius -multiplier .

my code is : .\rtkfdk.exe -g -o -p -r “.\*.dcm” -v “on” --neworigin “-214.76,214.76,0” --newdirection “1,0,0,0,-1,0,0,0,1” --hardware “cuda” --newspacing “0.56,0.56,1” --spacing 0.28 --dimension 240 --origin “-40,-80,-30”

and

.\rtkfdk.exe -g -o -p -r “.\*.dcm” -v “on” --neworigin “-214.76,214.76,0” --newdirection ‘1,0,0,0,-1,0,0,0,1’ --newspacing “0.56,0.56,1” --hardware “cuda” --spacing 0.28 --dimension 560

---

<div class="post-metadata">

### Author: ![simon.rit](https://discourse.itk.org/letter_avatar_proxy/v4/letter/s/f08c70/32.png) [@simon.rit](https://discourse.itk.org/u/simon.rit)
#### Post date: [April 21, 2025, 7:40pm UTC](https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498/6 "2025-04-21T19:40:46Z")

</div>

It’s not unusual to have negative values as FDK does not enforce positivity. For example, noisy measurements in air oscillate around zero.  
I don’t know why the head and vertebra phantoms seem have negative values without looking at the projection images, you should probably check for negative values in them or truncation artefacts.

---

<div class="post-metadata">

### Author: ![zhangxiang](https://discourse.itk.org/letter_avatar_proxy/v4/letter/z/c77e96/32.png) [@zhangxiang](https://discourse.itk.org/u/zhangxiang)
#### Post date: [April 22, 2025, 1:42am UTC](https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498/7 "2025-04-22T01:42:06Z")

</div>

thank you for your help.

---

<div class="post-metadata">

### Author: ![zhangxiang](https://discourse.itk.org/letter_avatar_proxy/v4/letter/z/c77e96/32.png) [@zhangxiang](https://discourse.itk.org/u/zhangxiang)
#### Post date: [April 22, 2025, 2:02am UTC](https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498/8 "2025-04-22T02:02:54Z")

</div>

The voxel values of my input 2D images are 0 for black (phantom) and 255 for white(nothing or air). Do I need to inverse it?

---

<div class="post-metadata">

### Author: ![simon.rit](https://discourse.itk.org/letter_avatar_proxy/v4/letter/s/f08c70/32.png) [@simon.rit](https://discourse.itk.org/u/simon.rit)
#### Post date: [April 22, 2025, 7:45am UTC](https://discourse.itk.org/t/why-are-there-negative-voxel-values-in-my-rtkfdk-x-ray-reconstruction/7498/9 "2025-04-22T07:45:47Z")

</div>

Yes, it would make sense as the linear attenuation coefficient of air is ~0 and positive for tissues. Are you using RTK to forward project your phantom?
