继续手搓DICOM文件解析 （2）

已于 2024-10-06 07:38:56 修改
阅读量363 收藏
点赞数 1
分类专栏： C++ & DICOM 文章标签： c++ DICOM
于 2024-10-05 21:54:19 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/kangdehua/article/details/142719496
版权

继续手搓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主页去查询，这个与传输与法有关联。我还在摸索中，等待下一期的结果吧，谢谢大家！如果有哪位专业人士指导原因，也可以评论区或者私信指点一下，非常感谢！

不拱地的猪
关注
专栏目录
继续手搓DICOM解析 （3）
不拱地的猪的专栏
10-06 436
但是这个TAG标签我们解析成了（ffe0，0010），其实应该是数据范围问题，读取使用的是char类型，它的范围是-127到127，如果是0到255就不会出现这个问题了。这个手搓项目基本完成了，除了最后一个Tag [7fe0,0010]这个有问题之外，其他的基本正确吧，有这方面的专业人员可以在评论区或者是私信中指导一下。这个DICOM文件呢，是BrainLab计划系统里面导出的图像，相对于CT机里面dicom出来的数据简化了很多信息，因此本文使用了此数据。里面的重点是数据类型，小端字节序，指针传递。
DICOM文件解析
qq_35769071的博客
12-10 996
dicom是一种神奇的格式，可以用dcmtk这个第三方库进行解析。 但是，dcmtk的编译实在令人头大，还要用到opencv库。 不就是把图片读出来嘛，吾辈当自强，我们自己也可以！！（注意！！本代码仅可解析非常非常非常基础的几种dicom，吾辈强度有限）！ 于是参照Dicom格式文件解析器 - assassinx - 博客园这位大佬的解析和实现，鄙人自己写了一个C++版本的Dicom解析工具。 DicomHelper类负责读取并解析dicom文件，并将图片数据存储到一个特殊的结构体里。结构体包括..
dicom文件详细解析
热门推荐
wenzhi的博客
07-14 8万+
dicom文件解析Dicom全称是医学数字图像与通讯，这里讲diocm格式文件的解读，读取本身是没啥难度的 无非就是字节码数据流处理。只不过确实比较繁琐。 dicom中文协议文档：http://download.csdn.net/download/wenzhi20102321/9897014 dicom文件结构： 整体结构先是128字节所谓的导言部分，跳过就是了，接着就是四个字节组成的字符串，
Dicom文件解析
木头的专栏
10-13 6874
Dicom全称是医学数字图像与通讯，这里讲的暂不涉及通讯那方面的问题 只讲*.dcm 也就是diocm格式文件的读取，读取本身是没啥难度的 无非就是字节码数据流处理。只不过确实比较繁琐。 好了 正题 分析 整体结构先是128字节所谓的导言部分，说俗点就是没啥意义的破数据 跳过就是了，然后是dataElement依次排列的方式 就是一个dataElement接一个dataEleme
医学图像DICOM文件解析——RT Structure篇
Joker 007的博客
11-07 3232
医学图像DICOM头文件解析——RT Structure篇
C# Dicom文件解析获取Tag数据
03-10
FO-DICOM是基于.NET的开源库，它为C#开发者提供了解析和操作DICOM文件的能力。PACS（Picture Archiving and Communication System）则是用于存储、管理和检索医学影像的系统。本教程将详细讲解如何使用C#通过FO-...
DICOM文件解析.rar_DICOM文件解析_dicom_dicom 读取_matlab dicom
07-15
MATLAB中读取DICOM文件，并解析其包含信息，并显示出来
DICOM.zip_dicom_dicom matlab_dicom 文件
07-15
在实际应用中，如果你有一个包含多个DICOM文件的压缩包"DICOM.zip"，你需要先解压，然后才能逐一处理这些文件。MATLAB的`unzip`函数可以完成解压操作： ```matlab unzip('DICOM.zip'); ``` 解压后，你可以遍历所有...
Dicom文件解析DLL
10-10
在您提到的“Dicom文件解析DLL”资源中，包含了三个重要的DLL文件： 1. Dicom.Native64.dll 和 Dicom.Native.dll：这些是原生库文件，分别针对64位和32位操作系统。它们实现了对底层操作系统的直接调用，用于处理与...
dicom.rar_DICOM格式文件_dicom
09-24
DICOM（Digital Imaging and Communications in Medicine）是一种标准格式，用于在医疗环境中存储、传输和查看医学图像。...此外，学习如何解析和利用元数据对于深度分析和应用 DICOM 文件是很有帮助的。
DicomParsing:js 解析 dicom 文件(基于 cornerstone.js 实现)
05-11
一个通过 HTTP (WADO-URI) 或 DICOMWeb (WADO-RS) 的 "DICOM P10 instances" 基于 的医学影像解析程序。 可以将 与 WADO-URI servers、DICOMWeb servers 或 任何返回 DICOM P10 instances (例如: 或 自定义 servers ...
上传文件和解析dicom文件所需的jar包
01-04
总之，"commons-fileupload-1.4.jar"和"dcm4che-core-5.22.6.jar"是开发基于JSP的DICOM文件处理系统的核心组件，它们分别解决了文件上传和DICOM数据解析的问题，为构建高效、安全的医疗影像管理平台提供了必要的工具...
编译器、IDE对C/C++新标准的支持
编程语言技巧经验分享
11-08 1131
很遗憾的是，微软VC++对C99的支持直到VS2015才基本完全支持，在VS2012之前，仅仅支持C89. 也正是2015年之前，这么多人骂微软并纷纷转移到其它C语言编译器的原因。VS2019开始才正式宣布支持C11/C17.对于C++标准，VS2013开始支持C++11，VS2017支持C++17, VS2019支持C++20, 显然，微软更偏爱C++.
STL set
hls0611的博客
11-09 241
练习题：P1059 [NOIP2006 普及组] 明明的随机数。
GitLab 提交 C++ 技巧
nbspzs的专栏
11-08 369
在 GitLab 中提交 C++ 代码时，通常应只提交 源代码文件（例如 .cpp、.h）和 CMake 配置文件（例如 CMakeLists.txt），而不要提交通过 CMake 编译后的中间文件或构建产物。这样，当其他开发者克隆你的项目后，他们可以在自己的机器上使用 cmake 构建项目，而不需要下载编译产物。2. 跨平台兼容性：不同的开发者和平台可能会生成不同的编译产物，提交这些文件会造成不必要的冲突和问题。1. 在 .gitignore 文件中忽略编译生成的文件和目录。
C++——智能指针
最新发布
Wangx_wang的博客
11-12 580
介绍了几种智能指针的不同和实现思路以及循环引用的问题
C++builder中的人工智能（15）：C++高斯误差线性单元（GELU）
caridle的专栏
11-08 708
激活函数（phi()），也称为转移函数或阈值函数，它根据净输入函数的给定值（sum）确定激活值（a = phi(sum)）。在这里，sum是它们权重中的信号之和，激活函数是这个和的新值，具有给定的函数或条件。换句话说，激活函数是将所有加权信号的和转换为该信号的新激活值的方法。有不同类型的激活函数，通常使用的是线性（恒等）、双极性和逻辑（sigmoid）函数。在C++（以及大多数编程语言）中，你可以创建自己的激活函数。注意，sum是净输入函数的结果，它计算所有加权信号的和。我们将使用这些作为输入函数的结果。
C++ 内存分布及 new , delete 分配问题( ~~~ 面试重要 ~~~)
GAOSHOUCYUYAN的博客
11-10 762
本篇文章笔者将会对 C++ 中的内存问题简单的讲解 , 同时对 new , delete 的面试题进行重点讲解.以上就是本篇的内容 , 希望读者认真领会.
2024年华为OD机试真题-查找充电设备-C++-OD统一考试（E卷）
面试高手的博客
11-08 329
某个充电站，可提供n个充电设备只，每个充电设备均有对应的输出功率任意个充电设备组合的输出功率总和，均构成功率集合P的1个元素功率集合P的最优元素，表示最接近充电站最大输出功率pmax 的元素。每一题都含有详细的解题思路和代码注释，精编c++、JAVA、Python三种语言解法。当充电设备输出功率50、20、20组合时，其输出功率总和为90，最接近充电站最大充电输出功率，因此最优元素为90。所有充电设备的输出功率组合，均大于充电站最大充电输出功率30，此时最优元素值为0。第 2 行为每个充电设备的输出功率。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值