好的,所以我正在做的是在现有的 DICOM 文件中添加一个新的叠加层并保存它(DICOM 文件现在有两个叠加层)。一切都保存没有错误,两个 DICOM 查看器 Sante 和 ClearCanvas-Workstation 打开文件,但只有 Sante 显示两个叠加层。现在,当我查看 DICOM 文件中的标签时,OverlayData(6000) 'VR' 是 'OW' 而 OverlayData(6002) 'VR' 是 'OB'。所以我的问题是如何创建一个“VR”为“OW”的新标签,因为这是用于 OverlayData 的正确标签。
这是我用来将新叠加层添加到 DicomFile.DataSet:: 注意的代码,在创建叠加层后,我确实将可见像素数据写入其中。
void AddOverlay()
{
int newOverlayIndex = 0;
for(int i = 0; i != 16; ++i)
{
if(!DicomFile.DataSet.Contains(GetOverlayTag(i, 0x3000)))
{
newOverlayIndex = i;
break;
}
}
//Columns
uint columnsTag = GetOverlayTag(newOverlayIndex, 0x0011);
DicomFile.DataSet[columnsTag].SetUInt16(0, (ushort)CurrentData.Width);
//Rows
uint rowTag = GetOverlayTag(newOverlayIndex, 0x0010);
DicomFile.DataSet[rowTag].SetUInt16(0, (ushort)CurrentData.Height);
//Type
uint typeTag = GetOverlayTag(newOverlayIndex, 0x0040);
DicomFile.DataSet[typeTag].SetString(0, "G");
//Origin
uint originTag = GetOverlayTag(newOverlayIndex, 0x0050);
DicomFile.DataSet[originTag].SetUInt16(0, 1);
DicomFile.DataSet[originTag].SetUInt16(1, 1);
//Bits Allocted
uint bitsAllocatedTag = GetOverlayTag(newOverlayIndex, 0x0100);
DicomFile.DataSet[bitsAllocatedTag].SetUInt16(0, 1);
//Bit Position
uint bitPositionTag = GetOverlayTag(newOverlayIndex, 0x0100);
DicomFile.DataSet[bitPositionTag].SetUInt16(0, 0);
//Data
uint dataTag = GetOverlayTag(newOverlayIndex, 0x3000);
DicomFile.DataSet[dataTag].SetNullValue();//<<< Needs to be something else
byte[] bits = new byte[(CurrentData.Width*CurrentData.Height)/8];
for(int i = 0; i != bits.Length; ++i) bits[i] = 0;
DicomFile.DataSet[dataTag].Values = bits;
}
public static uint GetOverlayTag(int overlayIndex, short element)
{
short group = (short)(0x6000 + (overlayIndex*2));
byte[] groupBits = BitConverter.GetBytes(group);
byte[] elementBtis = BitConverter.GetBytes(element);
return BitConverter.ToUInt32(new byte[]{elementBtis[0], elementBtis[1], groupBits[0], groupBits[1]}, 0);
}
所以在我看来,会有一些方法,比如 'DicomFile.DataSet[dataTag].SetNullValue();' 使用“OW”的“VR”创建标签。或者也许有一种完全不同的方式在 ClearCanvas idk 中添加叠加层......