itkImageをDICOMとして保存すると画像の位置と向きが書き込まれない

Question

3DボリュームがitkImage<unsigned char, 3>で、DICOMシリーズとして保存したい。シリーズを保存することはこれまでのところは機能しますが、メタデータタグ "画像位置患者"と "画像方向患者"はdicomファイルに保存されません。その他のタグはすべて正しく保存されます。

メタデータ辞書をコンソールに印刷すると、これらのタグが正しく印刷されます。 ITK-SNAPでファイルを開くと、これらの2つのタグはミスインです。

0010 | 0010患者

0010 | 0020 12345

0020 | 0032 0 \ 0 \ 0

0020 | 0037 1 \ 0 \ 0 \ 0 \ 1 \ 0

これは、これらのタグがメタデータに格納されるコードの一部です。私は実際の例でgistを作成しました。

UC3ImageType::PointType position; 
UC3ImageType::IndexType index; 
index[0] = 0; 
index[1] = 0; 
index[2] = f; 
image->TransformIndexToPhysicalPoint(index, position); 
value.str(""); 
value << position[0] << "\\" << position[1] << "\\" << position[2]; 
itk::EncapsulateMetaData<std::string>(*dictionary, "0020|0032", value.str()); 
value.str(""); 
value << position[2]; 
itk::EncapsulateMetaData<std::string>(*dictionary, "0020|1041", value.str()); 

itk::EncapsulateMetaData<std::string>(*dictionary, "0020|0037", "1\\0\\0\\0\\1\\0"); 

このコードで何が問題になっていますか？これらのタグも保存されるように、必要な他のタグはありますか？

Answer 1

異なるSOPクラスUIDを指定することで問題を解決できました。デフォルトではSecondary Capture IODが使用されていて、それをUltrasound Multi Frameに変更しました。これは自分の画像にも適しています。これで、画像の向き/位置が正しく保存されます。

const std::string SOP_CLASS_UID = "0008|0016"; 
const std::string US_MULTIFRAME_UID = "1.2.840.10008.5.1.4.1.1.3.1"; 
itk::EncapsulateMetaData<std::string>(*dictionary, SOP_CLASS_UID, US_MULTIFRAME_UID); 

しかし、私にとっては、何らかの理由でメタデータタグが保存されていないときに警告または例外が表示されないというのはまだ奇妙です。