0

0

pydicom教程:从外部数据源导入DICOM标签并自动匹配数据类型

霞舞

霞舞

发布时间:2025-10-17 09:22:24

|

490人浏览过

|

来源于php中文网

原创

pydicom教程:从外部数据源导入DICOM标签并自动匹配数据类型

本教程详细介绍了如何使用pydicom库从csv等外部数据源向dicom文件导入标准标签及其值,并解决数据类型不匹配的问题。通过实现一个基于dicom value representation (vr) 的动态类型转换机制,确保导入的数据符合dicom标准,避免pydicom警告,提高数据集成效率和准确性。

引言:DICOM数据集成中的类型挑战

在医疗影像数据处理中,经常需要将外部结构化数据(如患者信息、研究参数等,通常存储在CSV或数据库中)导入到DICOM文件中,以丰富其元数据。使用pydicom库进行此操作是常见的实践。然而,一个普遍的挑战是外部数据源的数据类型与DICOM标准中定义的Value Representation (VR) 所期望的数据类型不一致。

例如,CSV文件中的日期(如"20231110")或患者ID(如"42")通常会被Pandas等库读取为整数类型。但根据DICOM标准,AcquisitionDate (VR: DA) 和 PatientID (VR: LO) 都要求字符串类型。如果直接将整数值赋给这些标签,pydicom会发出UserWarning,提示“A value of type 'int' cannot be assigned to a tag with VR LO/DA”,这不仅可能导致数据存储异常,也可能在后续处理中引发错误。相反,像BitsAllocated (VR: US) 这样的标签则明确需要整数类型。因此,开发一个能够根据DICOM标签的VR自动进行数据类型转换的机制至关重要。

DICOM Value Representation (VR) 与 pydicom 的类型映射

DICOM标准通过Value Representation (VR) 来定义每个数据元素的预期数据类型和格式。例如:

  • DA (Date):表示日期,期望格式为YYYYMMDD,应为字符串。
  • LO (Long String):表示长字符串,如PatientID,应为字符串。
  • SH (Short String):表示短字符串,如StudyID,应为字符串。
  • UI (Unique Identifier):表示唯一标识符,如StudyInstanceUID,应为字符串。
  • IS (Integer String):表示整数字符串,如SeriesNumber,虽然是数字,但DICOM标准要求其为字符串格式,pydicom通常能自动处理整数到字符串的转换。
  • US (Unsigned Short):表示无符号短整数,如BitsAllocated,期望为整数。
  • AT (Attribute Tag):表示属性标签,其值本身是一个DICOM标签,需要特殊处理。

pydicom库内部维护了VR与Python数据类型之间的映射关系,并通过一系列集合(如STR_VR、INT_VR、FLOAT_VR、BYTES_VR)来分类不同的VR。我们可以利用这些内部定义,结合pydicom.datadict.dictionary_VR()函数获取指定标签的VR,从而实现动态类型转换。

MuleRun
MuleRun

全球首个AI Agent交易平台

下载

核心解决方案:实现动态类型转换函数 check_cast

为了解决数据类型不匹配问题,我们设计一个check_cast函数,它接收DICOM标签的VR类型和待转换的值,并返回经过适当类型转换后的值。

import pandas as pd
import pydicom
from pydicom.tag import Tag
from pydicom.datadict import dictionary_VR
from pydicom.valuerep import VR, STR_VR, INT_VR, FLOAT_VR, BYTES_VR
import warnings
import os

def check_cast(vr_type, value):
    """
    根据DICOM的Value Representation (VR) 类型,将值转换为适当的Python数据类型。

    Args:
        vr_type: DICOM标签的Value Representation (VR) 类型,例如 VR.DA, VR.LO。
        value: 待转换的原始值。

    Returns:
        any: 转换后的值,其类型与VR要求匹配。

    Raises:
        NotImplementedError: 如果遇到未实现的VR类型(如SQ)。
    """
    # 处理 NaN 值,将其转换为 None 或空字符串,取决于 VR 类型
    if pd.isna(value):
        if vr_type in STR_VR:
            return "" # 字符串类型VR的空值通常为空字符串
        elif vr_type in INT_VR or vr_type in FLOAT_VR:
            return None # 数字类型VR的空值可以为 None
        else:
            return None # 其他类型也返回 None

    if vr_type == VR.AT:  # Attribute Tag (属性标签)
        # 属性标签的值本身是一个DICOM标签,需要通过Tag()解析
        return Tag(value)
    elif vr_type in STR_VR: # 字符串类型VRs (如DA, LO, SH, UI, PN, IS等)
        return str(value)
    elif vr_type in INT_VR: # 整数类型VRs (如US, SS, UL, SL等)
        return int(value)
    elif vr_type in FLOAT_VR: # 浮点数类型VRs (如FL, FD, DS, AS等)
        return float(value)
    elif vr_type in BYTES_VR: # 字节类型VRs (如OB, OW, UN)
        # 字节类型需要特别注意,取决于原始值的格式。
        # 如果原始值是字符串,通常需要编码;如果是字节,则直接返回。
        # 这里的实现假设可以尝试直接转换为字节,具体场景可能需要更细致的处理。
        return bytes(value) if isinstance(value, str) else value
    elif vr_type == VR.SQ: # Sequence (SQ) VRs (序列)
        # 序列VR通常是一个包含多个Dataset的列表,其处理逻辑复杂,
        # 无法通过简单的类型转换实现,因此在此处抛出 NotImplementedError。
        raise NotImplementedError(f"Casting for Sequence (SQ) VR not implemented for value: {value}")
    else:
        # 对于其他未明确处理的VR类型,直接返回原始值,
        # 依赖pydicom的默认处理或后续的显式处理。
        return value

check_cast函数详解:

  • NaN值处理:Pandas读取CSV时,空单元格会被解析为NaN。check_cast首先处理NaN值,将其根据VR类型转换为适当的空值表示(如空字符串或None),避免后续转换错误。
  • VR.AT (Attribute Tag):这种VR的值本身是一个DICOM标签。函数会尝试使用Tag(value)将其解析为一个Tag对象。
  • STR_VR (字符串型VR):包括DA, `

相关专题

更多
python开发工具
python开发工具

php中文网为大家提供各种python开发工具,好的开发工具,可帮助开发者攻克编程学习中的基础障碍,理解每一行源代码在程序执行时在计算机中的过程。php中文网还为大家带来python相关课程以及相关文章等内容,供大家免费下载使用。

749

2023.06.15

python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

635

2023.07.20

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

758

2023.07.25

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

618

2023.07.31

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

1262

2023.08.03

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

547

2023.08.04

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

577

2023.08.04

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

706

2023.08.11

jQuery 正则表达式相关教程
jQuery 正则表达式相关教程

本专题整合了jQuery正则表达式相关教程大全,阅读专题下面的文章了解更多详细内容。

1

2026.01.13

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
最新Python教程 从入门到精通
最新Python教程 从入门到精通

共4课时 | 0.6万人学习

Django 教程
Django 教程

共28课时 | 3万人学习

SciPy 教程
SciPy 教程

共10课时 | 1.1万人学习

关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号 技术交流群
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn | 湘ICP备2023035733号