stir.AcquisitionData.get_subset() fails on get_output() of ListmodeToSinograms, even with stir.AcquisitionData.set_storage_scheme("memory")

[epapoutsellis]

I would like to perform stochastic optimisation using CIL for the NEMA dataset. I use this notebook.

If I run acq_data.get_subset(i) then I get this error

ERROR: write_basic_interfile_PDFS_header: Error casting the projdata to one of its geometries: Cylindrical/BlocksOnCylindrical/Genreic

[KrisThielemans]

I think this is the known problem that we cannot write subset-data yet. UCL/STIR#1003

I also think in the past we solved this by setting stir.AcquisitionData.set_storage_scheme("memory") (before getting the subset)?

[KrisThielemans]

Cause:

SIRF/src/xSTIR/cSTIR/stir_data_containers.cpp

Lines 382 to 386 in e34ae29

	std::unique_ptr<STIRAcquisitionData>
	STIRAcquisitionDataInFile::get_subset(const std::vector<int>& views) const
	{
	auto ptr_ad = new STIRAcquisitionDataInFile(std::move(_data->get_subset(views)));
	// auto ptr_ad = new STIRAcquisitionDataInMemory(std::move(_data->get_subset(views)));

but this currently fails as STIR cannot write subset data yet.

[evgueni-ovtchinnikov]

@epapoutsellis have you tried storing acquisitions in memory as Kris suggested:

I also think in the past we solved this by setting stir.AcquisitionData.set_storage_scheme("memory") (before getting the subset)?

[KrisThielemans]

@epapoutsellis I think this was resolved as above? I cannot quite remember...

[KrisThielemans]

Setting to memory does not help when using ListmodeToSinograms.get_output(). Current work-around is to create stir.AcquisitionData from the corresponding file.