[CCPPETMR/SIRF] mSTIR.AcquisitionSensitivityModel constructor documentation is unclear (#95)

4 views
Skip to first unread message

Kris Thielemans

unread,
Feb 8, 2018, 7:14:00 AM2/8/18
to CCPPETMR/SIRF, Subscribed

we do not understand from the doc that the 2nd argument has to be a previously constructed mSTIR.AcquisitionSensitivityModel

This needs to be in a demo as well. I would put it in the reconstruction_from_listmode demo. It's the natural place. @rijobro is going to suggest some changes on that one.

Can you also add in the doc that the attenuation image needs to be in units cm^-1?


You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub, or mute the thread.

Kris Thielemans

unread,
Feb 8, 2018, 9:00:02 AM2/8/18
to CCPPETMR/SIRF, Subscribed

In addition to the above, I think most users (including me) will get confused by

unnormalise      Multiplies the argument (AcquisitionData) by (1/n).
forward          Returns the argument multiplied by (1/n).

These 2 functions currently do the same (they're aliases effectively). that should be clear from the doc. I'd say

forward   Identical to unnormalise

evgueni-ovtchinnikov

unread,
Feb 8, 2018, 9:10:52 AM2/8/18
to CCPPETMR/SIRF, Subscribed

added 'existing' and attenuation units

examples of constructors are in using_acquisition_model demos

'Identical to unnormalise': wrong: unnormalise changes (multiplies by 1/n) the argument, forward does not (returns new AcquisitionData object)

Kris Thielemans

unread,
Feb 8, 2018, 9:30:34 AM2/8/18
to CCPPETMR/SIRF, Subscribed

agreed for the "not identical.". sorry. what about saying

forward          Returns a new AcquisitionData equal the input argument multiplied by (1/n).

evgueni-ovtchinnikov

unread,
Feb 8, 2018, 9:39:26 AM2/8/18
to CCPPETMR/SIRF, Subscribed

done

evgueni-ovtchinnikov

unread,
Feb 8, 2018, 9:39:26 AM2/8/18
to CCPPETMR/SIRF, Subscribed

Closed #95.

Kris Thielemans

unread,
Feb 8, 2018, 9:41:51 AM2/8/18
to CCPPETMR/SIRF, Subscribed

Reopened #95.

Kris Thielemans

unread,
Feb 8, 2018, 9:41:51 AM2/8/18
to CCPPETMR/SIRF, Subscribed

using_acquisition_model demo seems to have only

asm = AcquisitionSensitivityModel(bin_eff)

which is fine for me. I would there just refer to other demos on AcquisitionSensitivityModel.

we cannot find a demo that uses both atn and norm. Needs to be somewhere. The recon demo I would definitely include both as that's what people will want to do.

evgueni-ovtchinnikov

unread,
Feb 8, 2018, 9:57:14 AM2/8/18
to CCPPETMR/SIRF, Subscribed

I was referring to using_acquisition_model.m, the python demo is asm_sinograms.py

Kris Thielemans

unread,
Feb 8, 2018, 10:03:27 AM2/8/18
to CCPPETMR/SIRF, Subscribed

ok, using_acquisition_model.py doesn't do the chaining. I now see there's an example of the chaining in asm_sinograms.py indeed. It's fine then as it is, but it still needs to be in the recon

Kris Thielemans

unread,
Feb 8, 2018, 11:35:18 AM2/8/18
to CCPPETMR/SIRF, Subscribed

Closed #95.

Reply all
Reply to author
Forward
0 new messages