• dicom和dicomdir


    转载http://blog.sina.com.cn/s/blog_4bce5f4b01019ix5.html

    DICOM 文件内容在 Part 3 DICOM IOD 里定义。CT, MR, CR, DR, US, NM, PET, XA 等各有自己的内容定义,由共同的专有的部分 (IE 和 Modules) 组成。 

    DICOMDIR 由 Part 10 定义。是一种 mini 可变长度数据库文件结构。 

    过去我手下一个学计算机的非常聪明的印度人问我:图像文件格式不是有 BMP, GIF, TIF 和 JPEG 等等吗,为什么还要定义 DICOM 呢? 我向他解释了好几次最后还是放弃了。 任何图像文件格式无非是由两个部分组成:存参数的 header 和图点数据 (pixel data)。 BMP, JPEG, TIFF 之类的格式的 header 只描述图像的基本参数:如几行、几烈、每点用了几位、有没有压缩、调色板等等。Header 往往是固定长度的。 而医疗影像还要许多其它参数,如病人基本资料、检验基本资料、系列资料、位置资料等等。而且每种 Modality 和每种 image 所需要的内容不一样。因此,一般的图像格式不能用。 

    DICOM 的 4 个内容层次: 

    1. Patient (病人) 

    2. Study (检验) 

    3. Series (系列)

    4. Image (图像) 

    尽管头几层的内容在很多图像里是相同的,它们在每个图像文件里都要有。 每一层叫一个 Information Entity,或 IE (从 relational database schema 设计引用而来)。每一层又细分成 Module。每个 Module 里面的最小单元叫做一个 attribute 或 element。

    现在举个例子:

    CR 图像

    DICOM Part 3, A.2.3, Table A.2-1 

    1. Patient IE: a. Patient Module (参考 C.7.1.1) 

    2. Study IE: a. Study Module (参考 C.7.2.1) b. Patient Study Module (参考 C.7.2.2) 

    3. Series IE: a. General Series (参考 C.7.3.1) b. CR Series (参考 C.8.1..1) c. General Equipment (参考 C.7.5.1) 

    4. Image IE: a. Genrral Image (C.7.6.1) b. Image Pixel (C.7.6.3) c. Contrast/bolus (C.7.6.4) d. CR Image (C.8.1.2) 

    ... 

    i. SOP Common (C.12.1) 

    将这些 modules (tables) 里的所有 elements 都找出来就做成了一个 CR 图像的架构。 要注意的是这些 module 有些是一定要的 (modatory) 有些是用户选用的 (user)。 到了每个 module 里 attribute/element 表又有分五类 Type 1, 1C, 2, 2C 和 3。 Type 1 是一定要的,2 也是一定要的但是内容可以是空的。Type 3 则可要可不要。 所以浓缩一下,一个 CR 图像里的元素 (elements) 也不是太多。

    把这些表格展开后,这些 elements 组成一个 dataset。 那么写到一个文件里或通过网路传送又是个什么格式呢?这就要看 Part 5。 

    一个元素 (element) 的结构是: 

    1. group tag: 16-bit 

    2. element tag: 16-bit 

    3. length (or VR/length): 32-bit 

    4. data (bytes of length) 

    对应每一个用到的 element DICOM 标准 Part 6 都定义了一个 group tag 和 element tag。

    比如说: patient name: 0x0010, 0x0010 patient ID: 0x0010, 0x0020 ... 

    VR 说的是 element 格式,比如说 patinet name 的 VR 是 PN。格式是 last_name^first_name^middle_name^prefix^surfix。那么我的英文名字就是: Wang^JB^^Dr.^ 

    往外写时要做几个事情: 

    1. 要把所有元素按 group tage 和 element tag 理一遍 (sort)。从小排到大。 

    2. 如果是写 DICOM 介质的 DICOM file 还先写 128 bytes preamble (一般是空白),加 "DICM", 加 group 2 Meta header。(讲到 Part 10 时再细说) 

    3. 如果 dataset 里面含有 Sequence elements, sequence 里面每一个 Item 又是一个 dataset。

    作业: 有了这些知识,你就可以开始写一个小小的 BMP 到 DICOM 的转化程序。

    关键资料: 

    Modality (0008, 0060): SC 

    Photometric Interpretation (0028, 0004): RGB

    SOP Class UID: 1.2.840.10008.1.5.4.1.1.7 

    最简单的办法是写一个 structure,然后一个 

    array. typedef struct DicomElem(short int group_tag, short int element_tag, char VR[4], int length, char data[128]) DicomElem; 

    DicomElem CRDataSet [] = ( ( 0x0008, 0x0005, "CS", 10, "ISO_IR 100"), ( 0x0008, 0x0008, "CS", 16, "ORIGINAL\PRIMARY"), ... (0x0010, 0x0010, "PN", 16, "My^Test^Image^^ "), 

    (0x0010, 0x0020, "SH", 6, "123456"), 

    ... 

    (0,0,"",0,"") 

    ); 

    void WriteCDImage(FILE *fp)

    {

    DicomElem elem = CRDataSet[0]; 

    unsigned long int lComboTag; 

    int nCols, nRows; 

    unsigned char *pPixelData unsigned long int lPixelLength; 

    pPixelData = LoadBMPImgeData("MyImage.bmp", nCols, nRows, lPixelLength); 

    while(CRDataSet.group_tag)

    lComboTag = (CRDataSet.group_tag << 16) | CRDataSet.element_tag;

    switch(lComboTag) 

    {

    case 0x00280010: *((short int *)CRDataSet.data) = nCols; 

    break; 

    case 0x00280011: *((short int *)CRDataSet.data) = nRows; 

    break; 

    ... 

    // Write group and element tag 

    fwrite(&lComboTag, 1, sizeof(long), fp); 

    // Write VR 

    fwrite(CRDataSet.VR, 1, 2, fp); 

    if (lComboTag != 0x7fe00010) 

    fwrite(CRDataSet.length, 1, sizeof(short), fp); 

    fwrite(CRDataSet.data, 1, CRDataSet.length, fp); 

    }

    else

    {

    fwrite("", 1, 2, fp); 

    // Two blank bytes after VR 

    fwrite(&lPixelLength, 1, sizeof(long), fp); 

    // Length 

    fwrite(pPixelData, 1, lPixelLength, fp);

    }

    i++; 

    }

    }

    unsigned char *LoadBMPImgeData(char *fileName, int &nCols, int &nRows, unsigned long &lPixelLength) ( .... ) 

    细节自己去写。 

    更上一层楼: DICOM 文件读写最难的是两个东西: 

    o DICOM Sequence 

    o DICOM Pixel Data 

    刚刚让大家做的作业里没有 Sequence,pixel data 也很简单。 Sequence 在 C 里的类比是一个 structure 的 array。 是结构套结构。 所以读起来难。更讨厌的是 Sequence 还可以不定义长度, 即长度是 -1。要靠你自己去找 (FFFE, E0DD) 来决定 Sequence 是否结束。 Array 里面的每个 structure 单元就是 DICOM Sequence 里的一个 Item。Item 的开头是一个特定 element (FFFE, E000)。 如果 Item 的长度是 -1, 要靠找到 (FFFE, E00D) 来决定 Item 的结束。 Pixel Data (7fe0, 0010) 是一个特殊的 DICOM element。总是在所有元素的最后面。它个 Sequence 有两个相似的地方: o 长度区总是 32-位 (即 explicit VR 的情况下要在 VR 区后面填两个 bytes,然后再加 4-bytes 的长度。 o 如图像是压缩的,每幅图用一个 item 来存。 第一个 item 是个 offset table。每幅图的 offset 是一个 dword (4 bytes),第一幅图的 offset 是 0。

    DICOM 介质与 DICOMDIR: 

    DICOM 图像通过 Store SCU/SCP 来传的时候最小的组是 (0008, xxxx),一般头一个 element 是 (0008, 0000) 或 (0008, 0005)。 而存在介质上 (如光盘,硬盘和 MOD),则要加两套东西: 

    1. Preamble 和 DICOM signature: 128 bytes (一般是空的) 后面跟 DICM。 

    2. Meta header: group (0002, xxxx) 的十几个 elements. 看 DICOM Part 10。 DICOMDIR 是一个可变长度 迷你 database 文件。由 group (0002, xxxx) 和 group (0004, xxxx) 为主题。描述的是一个 4 层的树状结构 (tree structure)。 1. Patient 2. Study

    3. Series 4. Image 除了用 DICOM 惯用的 Sequene 外还用了 offset 来做 linked list。所以读写起来比较麻烦。请看 DICOM Part 10。

    DICOM 的中文支持: DICOM 本身可以支持任何文字。应用程序则可能无法破译。 DICOM 图像用中文是可以通过下列方法实现: (0008,0005)设 "GB2312" 或 “Big5"。 所有文字用双字节 GB 或 Big 5 码。 大部分的 DICOM 控件用起来都会有问题。国标用的是 8-bit bytes 而非 7-bit printable ASCII。 ANSI 字串函数如 sprintf, sscanf 会将最高一个 bit 削掉。要改写成支持 wide character 的才行。

  • 相关阅读:
    ASP.NET 4.0 与 Entity Framework 4第三篇使用Entity Framework调用存储过程
    雕虫无小技 JavaScript初学者的10个迷你技巧
    IE6下button随着文字的增多两边的内容边框也会增加的bug
    bigint ,int ,smallint ,tinyint 数据类型
    分享7个不错的jQuery游戏( 转)
    IE CSS Bug及解决方案参考手册
    利用CSS样式打印
    SQL 2005 弹出不允许对系统目录进行即席更新解决方法
    VS2010快捷键
    SQL2K,DTC错误:"该伙伴事务管理器已经禁止了它对远程/网络事务的支持"的解决办法
  • 原文地址:https://www.cnblogs.com/bayzhang/p/5707764.html
Copyright © 2020-2023  润新知