继续手搓DICOM文件解析 (2)

news2024/10/6 1:41:27

继续手搓DICOM文件解析 (2)
书接上回,我们使用了简单的函数进行了DICOM文件的解析,那么我们就继续吧。

// to be continued (part 2)
    std::cout<<std::string(25,'%')<<" Part 2 continued "<<std::string(25,'%')<<std::endl;
    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));```
    以上,就是重复使用读取TAG函数,读取VR函数,然后使用getShort函数来读取数据的长度。还好很多数据都是字符串形式,可以使用readString函数来读取,下面是读取的结果。
    
```cpp
%%%%%%%%%%%%%%%%%%%%%%%%% Part 2 continued %%%%%%%%%%%%%%%%%%%%%%%%%
[0008,0005] :  CS 10 # ISO_IR 100
[0008,0008] :  CS 24 # DERIVED\SECONDARY\AXIAL 
[0008,0012] :  DA 8 # 20180828
[0008,0013] :  TM 14 # 145109.969000 
[0008,0016] :  UI 26 # 1.2.840.10008.5.1.4.1.1.2
[0008,0018] :  UI 54 # 1.2.276.0.20.1.4.5.964244510648.1152.1535439069.968750
[0008,0020] :  DA 8 # 20160113
[0008,0023] :  DA 8 # 20180828
[0008,0030] :  TM 0 # 
[0008,0033] :  TM 14 # 145109.969000 
[0008,0050] :  SH 0 # 
[0008,0060] :  CS 2 # CT
[0008,0070] :  LO 8 # BrainLAB
[0008,ff90] :  PN 0 # 
[0008,1010] :  SH 14 # RTWORKSTATION2
[0008,1030] :  LO 20 # RT^1_RT_Head (Adult)
[0009,0010] :  LO 20 # BrainLAB_Conversion 
[0009,1001] :  LO 24 # ARIA/Varis/OBI-Platform 
[0010,0010] :  PN 14 # Zheng^Mingtang
[0010,0020] :  LO 10 # 0000246783
[0010,0030] :  DA 8 # 19410113
[0010,0040] :  CS 2 # F 
[0018,0050] :  DS 16 # 2.9925389620816 
[0018,0060] :  DS 0 # 
[0018,1020] :  LO 6 # 4.1.4 
[0018,5100] :  CS 4 # HFS 
[0020,000d] :  UI 54 # 1.2.276.0.20.1.2.5.964244510648.3588.1451978243.765625
[0020,000e] :  UI 54 # 1.2.276.0.20.1.3.5.964244510648.1152.1535439069.765625
[0020,0010] :  SH 2 # 1 
[0020,0011] :  IS 2 # 2 
[0020,0012] :  IS 0 # 
[0020,0013] :  IS 2 # 50
[0020,0032] :  DS 50 # -259.15500957601\-184.74776528391\87.274973824935 
[0020,0037] :  DS 16 # 1\-0\-0\-0\1\-0 
[0020,0052] :  UI 56 # 1.2.276.0.20.1.30.5.964244510648.1152.1535439069.781250
[0020,1040] :  LO 0 # 

接下来我们会发现VR的数值变成US类型,此时就需要使用getShort函数来读取数值了,因为该数值是unsigned short 类型。需要注意到是小端字节序,读取之后自己进行拼接。代码如下:

std::cout<<"-------------[0028,0002]------------ "<<std::endl;
    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    unsigned short usvar=getShort(file);
    std::cout<<usvar<<std::endl;
    //std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    usvar=getShort(file);
    std::cout<<usvar<<std::endl;
    //std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    usvar=getShort(file);
    std::cout<<usvar<<std::endl;
    //std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    std::cout<<"-----------[0028,0100]----------"<<std::endl;
    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    usvar=getShort(file);
    std::cout<<usvar<<std::endl;
    //std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    usvar=getShort(file);
    std::cout<<usvar<<std::endl;
    //std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    usvar=getShort(file);
    std::cout<<usvar<<std::endl;
    //std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    usvar=getShort(file);
    std::cout<<usvar<<std::endl;
    //std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    length=getShort(file);
    std::cout<<length<<" # ";
    readString(file,str,length);
    std::cout<<std::string(str)<<std::endl;
    std::memset(str,0,sizeof(str));

    readTag(file);
    read2Byte(file);
    long len=getLong(file);
    std::cout<<len<<" # ";

然后输出的结果如下:

-------------[0028,0002]------------ 
[0028,0002] :  US 2 # 1
[0028,0004] :  CS 12 # MONOCHROME2 
[0028,0010] :  US 2 # 512
[0028,0011] :  US 2 # 512
[0028,0030] :  DS 20 # 0.9765625\0.9765625 
-----------[0028,0100]----------
[0028,0100] :  US 2 # 16
[0028,0101] :  US 2 # 16
[0028,0102] :  US 2 # 15
[0028,0103] :  US 2 # 0
[0028,1050] :  DS 2 # 35
[0028,1051] :  DS 2 # 80
[0028,1052] :  DS 6 # -1024 
[0028,1053] :  DS 2 # 1 
[ffe0,0010] :  OW 8 #

是不是差不多了,我们可以看到数据是的行和列都是512,数据是16位数据,但是大家看到没有,最后的图像数据,TAG是【ffe0,0010】,实际上,使用DCMTK解析应该是7fe0,0010我也不知道哪里出错了,而且我也没有解析出数据的长度。OW是对的,到nema主页去查询,这个与传输与法有关联。我还在摸索中,等待下一期的结果吧,谢谢大家!如果有哪位专业人士指导原因,也可以评论区或者私信指点一下,非常感谢!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/2190905.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

408笔记|随笔记录|自用|2

文章目录 cache和TLB查找总结 接上篇选择题 打开同一个文件共享同一个内存索引结点 cache和TLB cache由SRAM组成&#xff0c;TLB通常由相联存储器组成&#xff0c;可以由SRAM组成。 DRAM需要不断刷新&#xff0c;性能比较低。 TLB缺失可以由软件或者硬件处理&#xff0c;cac…

JC系列CAN通信说明

目录 一、CAN协议二、指令格式三、通信接线3.1、一对一通信3.2、组网通信 四、寄存器定义五、指令说明4、读取电源电压5、读取母线电流6、读取实时速度8、读取实时位置10、读取驱动器温度11、读取电机温度12、读取错误信息32、设定电流33、设定速度35、设定绝对位置37、设定相对…

Java-进阶二

ArrayList的源代码分析&#xff08;扩容原理&#xff09; 1 使用空参构造的集合&#xff0c;在底层创建一个容量为0的数组。2 添加第一个元素时&#xff0c;底层会扩容创建一个容量为10的数组。3 存满时会扩容1.5倍。4 如果一次添加多个元素&#xff0c;1.5倍还放不下&#xff…

机器学习/数据分析--用通俗语言讲解时间序列自回归(AR)模型,并用其预测天气,拟合度98%+

时间序列在回归预测的领域的重要性&#xff0c;不言而喻&#xff0c;在数学建模中使用及其频繁&#xff0c;但是你真的了解ARIMA、AR、MA么&#xff1f;ACF图你会看么&#xff1f;&#xff1f; 时间序列数据如何构造&#xff1f;&#xff1f;&#xff1f;&#xff0c;我打过不少…

提升客户满意度,数字化售后工单系统的实施策略

售后工单系统助企业优化运营&#xff0c;提升客户满意度。涵盖技术支持、维修服务、客户反馈管理等应用场景。ZohoDesk凭高效定制成企业首选。未来趋势包括AI、物联网、AR/VR技术集成&#xff0c;提升服务智能化水平。 一、售后工单系统的重要性 售后工单系统是企业用来管理和…

基于Java的GeoTools对Shapefile文件属性信息深度解析

目录 前言 一、Shapefile的属性列表信息 1、属性表格信息 2、属性表格包含的要素 二、GeoTools对属性表格的解析 1、常规解析方法 2、基于dbf文件的属性信息读取 三、总结 前言 ESRI Shapefile&#xff08;shp&#xff09;&#xff0c;或简称shapefile&#xff0c;是美…

SQL优化 - 排序

文章目录 排序和索引降序索引 FilesortORDER BY 顺序问题ORDER BY LIMIT 排序和索引 如果ORDER BY操作使用了索引&#xff0c;那么就可以避免排序操作&#xff0c;因为索引本身就是按索引 key 排好序的。那什么情况下&#xff0c;ORDER BY会走索引呢&#xff1f; 例如&#…

专题十_穷举vs暴搜vs深搜vs回溯vs剪枝_二叉树的深度优先搜索_算法专题详细总结

目录 搜索 vs 深度优先遍历 vs 深度优先搜索 vs 宽度优先遍历 vs 宽度优先搜索 vs 暴搜 1.深度优先遍历 vs 深度优先搜索(dfs) 2.宽度优先遍历 vs 宽度优先搜索(bfs) 2.关系图暴力枚举一遍所有的情况 3.拓展搜索问题全排列 决策树 1. 计算布尔⼆叉树的值&#xff08;medi…

CF2013E Prefix GCD

【题目大意】 给定一个长度为 n n n 的数列 a 1 … n a_{1 \dots n} a1…n​&#xff0c;你可以将 a 1 … n a_{1 \dots n} a1…n​ 按照任意顺序进行重排&#xff0c;使得&#xff1a; ∑ i 1 n gcd ⁡ { a 1 , a 2 , a 3 , … , a n } \sum\limits_{i1}^{n}\gcd\left \{…

10.5学习

1.GateWay GateWay⽬标是取代Netflflix Zuul&#xff0c;它基于Spring5.0SpringBoot2.0WebFlux等技术开发&#xff0c;提供统⼀的路由⽅式&#xff08;反向代理&#xff09;并且基于 Filter(定义过滤器对请求过滤&#xff0c;完成⼀些功能) 链的⽅式提供了⽹关基本的功能&…

探索 Python 虚拟环境的奥秘:virtualenv 的魔法世界

文章目录 探索 Python 虚拟环境的奥秘&#xff1a;virtualenv 的魔法世界背景&#xff1a;为何选择 virtualenv&#xff1f;虚拟环境的守护者&#xff1a;virtualenv 是什么&#xff1f;安装 virtualenv&#xff1a;简单几步&#xff0c;开启隔离之旅掌握 virtualenv 的基本用法…

Relu激活

ReLU&#xff08;Rectified Linear Unit&#xff09;激活函数 是卷积神经网络&#xff08;CNN&#xff09;以及许多深度学习模型中最常用的激活函数之一。它的主要作用是引入非线性&#xff0c;使模型能够学习和表达更复杂的特征。以下是对ReLU激活函数的详细解释。 1. ReLU的…

C语言进阶版第16课—自定义类型:结构体

文章目录 1. 结构体类型的声明和初始化2. 结构体自引用3. 结构体内存对齐3.1 结构体内存对齐规则3.2 修改默认对齐数 4. 结构体传参4. 结构体实现位段5. 位段使用的注意事项 1. 结构体类型的声明和初始化 结构体在使用之前都要对其类型进行声明&#xff0c;关键字是struct&…

15分钟学 Python 第36天 :Python 爬虫入门(二)

Python 爬虫入门&#xff1a;环境准备 在进行Python爬虫的学习和实践之前&#xff0c;首先需要准备好合适的开发环境。本节将详细介绍Python环境的安装、必要库的配置、以及常用工具的使用&#xff0c;为后续的爬虫编写奠定坚实的基础。 1. 环境准备概述 1.1 为什么环境准备…

mp4转gif在线转换怎么转?7个视频转动图方法不容错过!(超简单)

mp4转gif在线转换怎么转&#xff1f;如今&#xff0c;将mp4视频转换为gif动图格式&#xff0c;满足了人们对易于分享和网络传播内容的需求。与mp4视频相比&#xff0c;gif动图文件体积相对较小&#xff0c;几乎所有网络平台都支持这种格式&#xff0c;无需额外安装插件或软件。…

CSID-GAN:基于生成对抗网络的定制风格室内平面设计框架论文阅读

CSID-GAN: A Customized Style Interior Floor Plan Design Framework Based on Generative Adversarial Network 摘要前言II. CSID-GAN METHODA. Overall FrameworkB. Algorithm and Loss Function III. DATASETS AND EVALUATION METRICSA. DatasetsB. Evaluation Metrics IV.…

信息安全工程师(34)访问控制模型

前言 访问控制模型是实现访问控制的基础&#xff0c;不同的访问控制模型提供了不同的访问控制策略和机制&#xff0c;以适应不同的安全需求。 一、自主访问控制模型&#xff08;DAC&#xff09; 定义&#xff1a;指资源的所有者有权决定谁可以访问其资源以及访问的方式。资源的…

如何实现事件流操作

文章目录 1 概念介绍2 使用方法3 示例代码我们在上一章回中介绍了通道相关的内容,本章回中将介绍StreamProvider组件.闲话休提,让我们一起Talk Flutter吧。 1 概念介绍 在Flutter中Stream是经常使用的组件,对该组件的监听可void main() {///让状态栏和程序的appBar融为一体…

macos安装mongodb

文章目录 说明安装和配置安装mongodb配置PATH变量 验证日志及数据存放目录 mac启动和关闭mongodb后台启动失败问题mongodb-compass(GUI) 说明 Homebrew core 列表目前已经将 MongoDB 移除,不再为其提供支持。但是使用国内镜像的brew还是可以安装的&#xff01;这里直接从官网下…

使用 Python 代码连接 PostgreSQL

Python 是一个功能非常强大的编程语言&#xff0c;尤其在与数据库交互时&#xff0c;提供了丰富的解决方案。在实际项目中&#xff0c;我们经常需要通过 Python 连接并操作数据库。为了简化这种操作&#xff0c;ORM&#xff08;对象关系映射&#xff09;框架提供了便利。ORM 能…