STIR icon indicating copy to clipboard operation
STIR copied to clipboard

SPECTUB:Sensitivity appears inverted along x axis

Open danieldeidda opened this issue 5 years ago • 11 comments

I came across a repeatable issue. The sensitivity appears inverted along the x-axis when reconstructing SPECT data. As a consequence attenuation correction can be wrong. This happens in both master and in #445 but with inverted axis.

danieldeidda avatar Feb 21 '20 11:02 danieldeidda

this is in #445?

KrisThielemans avatar Feb 21 '20 11:02 KrisThielemans

can you please clarify if this was a new bug introduced in #445, or an existing one.

KrisThielemans avatar Feb 21 '20 12:02 KrisThielemans

can you please clarify if this was a new bug introduced in #445, or an existing one.

This happens in both master and in #445 but with inverted axis.

danieldeidda avatar Feb 21 '20 12:02 danieldeidda

ouch. @RebeccaGillen worth checking with @danieldeidda what's happening.

Let's try to merge #445 soon. I think there was 1 thing outstanding, but can't remember... continue the conversation on that there.

KrisThielemans avatar Feb 21 '20 12:02 KrisThielemans

ouch. @RebeccaGillen worth checking with @danieldeidda what's happening.

Let's try to merge #445 soon. I think there was 1 thing outstanding, but can't remember... continue the conversation on that there.

I did not notice at the beginning because the phantoms I was using were quite symmetrical but as soon as I tried something slightly more complex it was visible

danieldeidda avatar Feb 21 '20 12:02 danieldeidda

I double checked with the master and I can reproduce the sensitivity flip along x axis.

image

danieldeidda avatar Mar 03 '20 11:03 danieldeidda

Unfortunately I came across the same problem again: I reconstructed data from GE OPTIMA and the attenuation appears inverted along x in all the patients

danieldeidda avatar Aug 25 '21 15:08 danieldeidda

Unfortunately I came across the same problem again: I reconstructed data from GE OPTIMA and the attenuation appears inverted along x in all the patients

I mean in the sensitivity image. but can definitely change how the reconstructed image looks like especially where there is air (lungs)

danieldeidda avatar Aug 25 '21 15:08 danieldeidda

how do you get the attenuation? what file format?

I presume this means that the SPECT image (with atten) comes out x-flipped. This is what we'd need to fix. But let's be 100% sure about the attenuation.

KrisThielemans avatar Aug 25 '21 19:08 KrisThielemans

how do you get the attenuation? what file format?

I presume this means that the SPECT image (with atten) comes out x-flipped. This is what we'd need to fix. But let's be 100% sure about the attenuation.

The attenuation is estimated from the CT (dicom) using the STIR utility. The SPECT image is not flipped but the lung area is I noticed this when comparing vendor and STIR images. You can see below that the lesions and the liver are in the same place but the lungs shape appears inverted. when looking at the umap image (.v) used for the reconstruction it is inverted compared to the sensitivity and the reconstructed image

image

danieldeidda avatar Aug 26 '21 09:08 danieldeidda

Note that this issue is still present in the current version of STIR (v5.1.0). Tested on a Windows system with an asymmetric attenuation dataset: image

The sensitivity data are flipped in the x-direction: image

RebeccaGillen avatar Feb 13 '23 10:02 RebeccaGillen