KeyError: 'Name' about imctools HOT 11 CLOSED

Hey @algebio
Thanks a lot for the bug report!

This seems to be the same issue that is causing: BodenmillerGroup/ImcSegmentationPipeline#62

@plankter Could you look into this?

from imctools.

Hi,

An update: I just freshly installed everything and run tests and cannot see that anything looks suspicious.

Would you be able to share an example .ome.tiff ?

from imctools.

How did you generate the .ome.tiff? From .mcd files via imctools? imctools v1 or v2?

from imctools.

Hi @votti

Thanks for your help!

I created the ome.tiff from a .mcd. You can find everything in this folder. I used imctools v2 which I installed it yesterday from GitHub. I also updated Anaconda before running the Jupyter notebook. Something that I noticed when updating the external dependencies is that some of them never got updated and remained in the "need to be updated" list. However, imagecodecs, pandas, xmltodict, xtiff updated fine.

I used the same .mcd a few days ago with the previous version of Anaconda and imctools v1 and I got everything without any error messages.

Regards
Juan

from imctools.

Thanks @algebio!!
Having this data was extremly helpful.
Unfortunately I am not longer working in the mass cytometry field, so I so input from the community is really required to keep things runnig.

To the issue:

@plankter : The issue seems to because the ChannelLabel attribute that is used to fill the Name can be empty in the mcd .

This is a legit behavior, as the 'antibody name' is merely a label of the channel. The unique channel identifier is allways the isotope.

In imctools v1 I set the Isotope channel name as label, if the label was not in the mcd:

I think this filling-in is not longer done in imctools v2.

There are two t things to do:

1. Make sure the downstream processing also works with ome.tiff where the Name channel attribute is not set
   -> This is crucial to support also ome.tiff generated with imctools so far

Implementation: The quick fix would be to change line:

 int(chan.attrib["ID"].split(":")[2]): (chan.attrib["Name"], chan.attrib["Fluor"]) for chan in channels

to

 int(chan.attrib["ID"].split(":")[2]):(chan.attrib.get("Name",  chan.attrib["Fluor"]), chan.attrib["Fluor"]) for chan in channels

2. (Optional) Bring back the previous behavior to fill the Name attribute with the Channel name in cases that the Channel Label is left empty.

@plankter : Does this make sense? Could you implement this?

@nilseling : I heard you may be doing some acquisitions for a tiny test dataset: This issue could be reproduced by adding a channel manually to the acquisition and not setting a label. Would it be feasible to add this to your test acquisition??

Thanks team for helping to keep this running!

from imctools.

Hi @votti

I'm glad that I could help improve this already great pipeline. I can't wait to see these new changes implemented and use imctools v2 with my dataset.

Regards
Juan

from imctools.

Hey @votti, I will manually add Ag to capture the counterstain as one of the channels. This can then be used to test Label free channels. To me it would also make sense to bring back the old behaviour and fill the Name attribute with the metal isotope.

from imctools.

Hi

Thank you for addressing this issue. Could you let me know when it is sorted?

Regards
Juan

from imctools.

Hi,

I shall release an update today.

Best regards,
Anton Rau

from imctools.

Thank you @votti for the suggestions!

Version 2.1.5 with the fix is released.

Best,
Anton

from imctools.

Beautiful!

Thank you Vito, Anton and Nils. Great job!

Regards
Juan

from imctools.