4 DICOM成像协议编码实现-元数据组解析

2024-04-04 10:32

本文主要是介绍4 DICOM成像协议编码实现-元数据组解析,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

以下链接是本系列文章,不足之处,可在评论区讨论:
系列文章

以下链接中的代码是完整的且可运行的,链接如下,可按需下载:
dicom成像程序

  本篇文章对应 专栏 从零讲解DICOM协议-成像协议中的文章DICOM成像协议剖析和DICOM成像协议实现思路,建议先看以上两篇文章以了解DICOM底层协议,有助于理解代码实现。

上篇文章DICOM成像协议编码实现-文件头解析讲解了DICOM解析引擎的实现思路和整体代码框架,并完成

  1. 读取DICOM文件至内存中
  2. 读取文件头到文件头对象中

本篇文章将继续进行以下几部分的代码思路讲解和实现:

  1. 按照元数据组特性读取元数据组中的各个DataElement
  2. 按照数据组特性读取数据组中的各个DataElement
  3. 如果PixelData是压缩格式,则用相应的解压算法解压
  4. CT值转BMP,保存BMP图像

依赖:

class DcmRead : public FileRead
{
public:DcmRead(bool iscompress = false);DcmRead(string path, bool iscompress = false);DcmRead(char *buffer, int len, bool iscompress = false);~DcmRead();
public:virtual bool ReadFile(string path) override;vector<ElementData> ParseFile(TagName tagname);ElementFormat GetTransferSyntax();//获得除了pixeldata以外的元素ElementValue GetElement(TagName tagname);ElementData GetPixelData();DcmToBmpTag GetBMPInfo();void SaveBMP(string path);virtual void Clear() override;bool IsCompress();
public:DcmTags * GetDcmTags();DcmFile *GetDcmFile();void SetClassify(DcmFilePeriod *dcmclassify);
private://dicom文件前128个0和4个"DICM"virtual void GetHead() override;//dicom文件0x0002组,元数据组virtual void GetMetadata() override;//dicom文件信息组virtual void GetInfodata(void *dcmelementset) override;//dicom所有元素//virtual void GetElements(void *dcmelementset) override;private://压缩图像数据ElementData UnCompressPixel(vector<ElementData> pixelset);vector<ElementData> GetUnCompressPixel();void TransferUnCompress();
private:void MemCopy(int len, int &datalen, char **data);void TagClassify(int startindex,vector<DcmElement> *dcmelementset);//判断是否是元数据组0x0002(显式小端)bool IsMetaData(char *data);//判断是否是SQbool IsSQ(char * data);int GetDataLen(char *data);//是否完整获取SQ序列bool IsSQEnd();//函数是否递归操作bool IsEmbed();void SetVrDataLen(char *data);void DeletePartPointer(DcmElement dcmelement);void DeletePointer(DcmElement dcmelement);void Init(bool iscompress);
private:void SetBinding();bool(DcmRead::*CompareModel)(char *source, char *des);int(DcmRead::*ElementIntValue)(char *buffer, int len);bool BigCompareModel(char *source, char *des);bool LittleCompareModel(char *source, char *des);int BigElementLen(char *buffer, int len);int LittleElementLen(char *buffer, int len);
private:void GetDcmClassifyInfo();void FileClassify();
private:bool status;bool isbigmodel;bool isin;bool iscompress;DcmFile dcmFile;int offset;//SQ序列嵌套层数int embednum;int vrlen;int datalen;ElementFormat ef;CompressFormat cf;ElementValueType et;vector<DcmElement> sqelementset;//Pipe *jpeg2k;SocketClient *sc;private:DcmTags *dcmtags;DcmFilePeriod *dcmclassify;
};

  前面文章讲过元数据组的特点是:
  元数据组是小端显式格式
  元数据组和数据组都是由DICOM协议中的DataElement结构组成
  DataElement由Tag,VR,Length,Value组成
  Tag由group和elment组成,VR有显式和隐式之分,Length所在字节由显隐式和VR类型决定,Value长度是偶数等
  元数据组中最重要的元素是传输语法(0002,0010) Transfer Syntax UID
根据以上特点,从DcmRead类中实现GetMetadata()函数

void DcmRead::GetMetadata()
{//元数据都是显式小端ef = ElementFormat::ExplicitLittle;int startindex = sizeof(dcmFile.Head) + sizeof(dcmFile.DcmFlag);TagClassify(startindex, &dcmFile.DcmElementSet);GetTransferSyntax();
}

  startindex为从元数据组在buffer内存指针的偏移地址
  TagClassify函数的主要内容为:
  for循环内判断是否是元数据组0002,循环执行

memcpy(dcmelement.tag, buffer + offset, sizeof(dcmelement.tag));
SetVrDataLen(dcmelement.tag);
MemCopy(vrlen, dcmelement.vr.len, &dcmelement.vr.data);
MemCopy(datalen, dcmelement.datalen.len, &dcmelement.datalen.data);
int len = (this->*ElementIntValue)(dcmelement.datalen.data, dcmelement.datalen.len);
MemCopy(len, dcmelement.data.len, &dcmelement.data.data);

  显式,VR长度根据27种类型判断

class ElementBytes
{
public:static const int ImplicitLenBytes = 4;//隐式VR,vr长度为0字节,len长度为4字节static const int ExplicitLenComBytes = 2;//显式VR,普通类型vr长度为2字节,len长度为2字节static const int ExplicitVrComBytes = 2;//显式VR,OB,OW,OF,SQ,UT,UN的vr长度4字节,len长度为4字节static const int ExplicitLenOtherBytes = 4;//显式VR,OB,OW,OF,SQ,UT,UN的vr长度4字节,len长度为4字节static const int ExplicitVrOtherBytes = 4;//显式VR,OB,OW,OF,SQ,UT,UN的vr长度4字节,len长度为4字节};
void DcmRead::SetVrDataLen(char *data)
{if (!IsSQ(data)){switch (ef){case ElementFormat::ImplicitLittle:vrlen = 0;datalen = ElementBytes::ImplicitLenBytes;break;case ElementFormat::ExplicitLittle:case ElementFormat::ExplicitBig:case ElementFormat::CompressPixel:char vrtype[2];memcpy(vrtype, buffer + offset, sizeof(vrtype));for (int j = 0; j < TagType::VrType.size(); j++){if (memcmp(vrtype, TagType::VrType[j], sizeof(vrtype)) == 0){vrlen = ElementBytes::ExplicitVrOtherBytes;datalen = ElementBytes::ExplicitLenOtherBytes;break;}if (j == TagType::VrType.size() - 1){vrlen = ElementBytes::ExplicitVrComBytes;datalen = ElementBytes::ExplicitLenComBytes;}}break;default:break;}}
}

  得到VR的长度和Length的长度后,从buffer内存指针中获取对应的数据,偏移地址累加

void DcmRead::MemCopy(int len, int &datalen, char **data)
{datalen = len;*data = new char[datalen];memcpy(*data, buffer + offset, datalen);offset += datalen;
}

  小端,ElementIntValue绑定小端转换函数

void DcmRead::SetBinding()
{if (ef == ElementFormat::ExplicitLittle || ef == ElementFormat::ImplicitLittle || ef == ElementFormat::CompressPixel){CompareModel = &DcmRead::LittleCompareModel;ElementIntValue = &DcmRead::LittleElementLen;}else if (ef == ElementFormat::ExplicitBig){CompareModel = &DcmRead::BigCompareModel;ElementIntValue = &DcmRead::BigElementLen;}
}
int DcmRead::LittleElementLen(char *buffer, int len)
{int data = 0;int flag = 0x000000FF;for (int i = 0; i < len; i++){data |= ((buffer[i] << i * 8) & (flag << i * 8));}return data;
}

传输语法获取

ElementFormat DcmRead::GetTransferSyntax()
{for (int i = 0; i < dcmFile.DcmElementSet.size(); i++){if (memcmp(dcmFile.DcmElementSet[i].tag, TagType::TransferSyntaxUID, sizeof(dcmFile.DcmElementSet[i].tag)) == 0){if (memcmp(dcmFile.DcmElementSet[i].data.data, TransferSyntaxType::ImplicitLittle.c_str(), dcmFile.DcmElementSet[i].data.len) == 0){ef = ElementFormat::ImplicitLittle;isbigmodel = false;}else if (memcmp(dcmFile.DcmElementSet[i].data.data, TransferSyntaxType::ExplicitLittle.c_str(), dcmFile.DcmElementSet[i].data.len) == 0){ef = ElementFormat::ExplicitLittle;isbigmodel = false;}else if (memcmp(dcmFile.DcmElementSet[i].data.data, TransferSyntaxType::ExplicitBig.c_str(), dcmFile.DcmElementSet[i].data.len) == 0){ef = ElementFormat::ExplicitBig;isbigmodel = true;}else //压缩格式都是显式小端{isbigmodel = false;ef = ElementFormat::CompressPixel;if (memcmp(dcmFile.DcmElementSet[i].data.data, TransferSyntaxType::Jpeg2000.c_str(), dcmFile.DcmElementSet[i].data.len) == 0){cf = CompressFormat::Jpeg2000;}else if (memcmp(dcmFile.DcmElementSet[i].data.data, TransferSyntaxType::Jpeglossless.c_str(), dcmFile.DcmElementSet[i].data.len) == 0){cf = CompressFormat::Jpeglossless;}}break;}}return ef;
}

  至此,元数据组解析框架已基本完成,为突出重点和整体处理流程,部分函数没有列出完整内容。
  下篇文章将进行DICOM数据组解析框架讲解。

这篇关于4 DICOM成像协议编码实现-元数据组解析的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/875543

相关文章

pandas中位数填充空值的实现示例

《pandas中位数填充空值的实现示例》中位数填充是一种简单而有效的方法,用于填充数据集中缺失的值,本文就来介绍一下pandas中位数填充空值的实现,具有一定的参考价值,感兴趣的可以了解一下... 目录什么是中位数填充?为什么选择中位数填充?示例数据结果分析完整代码总结在数据分析和机器学习过程中,处理缺失数

Golang HashMap实现原理解析

《GolangHashMap实现原理解析》HashMap是一种基于哈希表实现的键值对存储结构,它通过哈希函数将键映射到数组的索引位置,支持高效的插入、查找和删除操作,:本文主要介绍GolangH... 目录HashMap是一种基于哈希表实现的键值对存储结构,它通过哈希函数将键映射到数组的索引位置,支持

Pandas使用AdaBoost进行分类的实现

《Pandas使用AdaBoost进行分类的实现》Pandas和AdaBoost分类算法,可以高效地进行数据预处理和分类任务,本文主要介绍了Pandas使用AdaBoost进行分类的实现,具有一定的参... 目录什么是 AdaBoost?使用 AdaBoost 的步骤安装必要的库步骤一:数据准备步骤二:模型

Pandas统计每行数据中的空值的方法示例

《Pandas统计每行数据中的空值的方法示例》处理缺失数据(NaN值)是一个非常常见的问题,本文主要介绍了Pandas统计每行数据中的空值的方法示例,具有一定的参考价值,感兴趣的可以了解一下... 目录什么是空值?为什么要统计空值?准备工作创建示例数据统计每行空值数量进一步分析www.chinasem.cn处

使用Pandas进行均值填充的实现

《使用Pandas进行均值填充的实现》缺失数据(NaN值)是一个常见的问题,我们可以通过多种方法来处理缺失数据,其中一种常用的方法是均值填充,本文主要介绍了使用Pandas进行均值填充的实现,感兴趣的... 目录什么是均值填充?为什么选择均值填充?均值填充的步骤实际代码示例总结在数据分析和处理过程中,缺失数

Java对象转换的实现方式汇总

《Java对象转换的实现方式汇总》:本文主要介绍Java对象转换的多种实现方式,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录Java对象转换的多种实现方式1. 手动映射(Manual Mapping)2. Builder模式3. 工具类辅助映

如何使用 Python 读取 Excel 数据

《如何使用Python读取Excel数据》:本文主要介绍使用Python读取Excel数据的详细教程,通过pandas和openpyxl,你可以轻松读取Excel文件,并进行各种数据处理操... 目录使用 python 读取 Excel 数据的详细教程1. 安装必要的依赖2. 读取 Excel 文件3. 读

Go语言开发实现查询IP信息的MCP服务器

《Go语言开发实现查询IP信息的MCP服务器》随着MCP的快速普及和广泛应用,MCP服务器也层出不穷,本文将详细介绍如何在Go语言中使用go-mcp库来开发一个查询IP信息的MCP... 目录前言mcp-ip-geo 服务器目录结构说明查询 IP 信息功能实现工具实现工具管理查询单个 IP 信息工具的实现服

SpringBoot基于配置实现短信服务策略的动态切换

《SpringBoot基于配置实现短信服务策略的动态切换》这篇文章主要为大家详细介绍了SpringBoot在接入多个短信服务商(如阿里云、腾讯云、华为云)后,如何根据配置或环境切换使用不同的服务商,需... 目录目标功能示例配置(application.yml)配置类绑定短信发送策略接口示例:阿里云 & 腾

Spring 请求之传递 JSON 数据的操作方法

《Spring请求之传递JSON数据的操作方法》JSON就是一种数据格式,有自己的格式和语法,使用文本表示一个对象或数组的信息,因此JSON本质是字符串,主要负责在不同的语言中数据传递和交换,这... 目录jsON 概念JSON 语法JSON 的语法JSON 的两种结构JSON 字符串和 Java 对象互转