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,从而实现动态类型转换。

Nanonets
Nanonets

基于AI的自学习OCR文档处理,自动捕获文档数据

下载

核心解决方案:实现动态类型转换函数 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, `

热门AI工具

更多
DeepSeek
DeepSeek

幻方量化公司旗下的开源大模型平台

豆包大模型
豆包大模型

字节跳动自主研发的一系列大型语言模型

WorkBuddy
WorkBuddy

腾讯云推出的AI原生桌面智能体工作台

腾讯元宝
腾讯元宝

腾讯混元平台推出的AI助手

文心一言
文心一言

文心一言是百度开发的AI聊天机器人,通过对话可以生成各种形式的内容。

讯飞写作
讯飞写作

基于讯飞星火大模型的AI写作工具,可以快速生成新闻稿件、品宣文案、工作总结、心得体会等各种文文稿

即梦AI
即梦AI

一站式AI创作平台,免费AI图片和视频生成。

ChatGPT
ChatGPT

最最强大的AI聊天机器人程序,ChatGPT不单是聊天机器人,还能进行撰写邮件、视频脚本、文案、翻译、代码等任务。

相关专题

更多
Python 时间序列分析与预测
Python 时间序列分析与预测

本专题专注讲解 Python 在时间序列数据处理与预测建模中的实战技巧,涵盖时间索引处理、周期性与趋势分解、平稳性检测、ARIMA/SARIMA 模型构建、预测误差评估,以及基于实际业务场景的时间序列项目实操,帮助学习者掌握从数据预处理到模型预测的完整时序分析能力。

80

2025.12.04

Python 数据清洗与预处理实战
Python 数据清洗与预处理实战

本专题系统讲解 Python 在数据清洗与预处理中的核心技术,包括使用 Pandas 进行缺失值处理、异常值检测、数据格式化、特征工程与数据转换,结合 NumPy 高效处理大规模数据。通过实战案例,帮助学习者掌握 如何处理混乱、不完整数据,为后续数据分析与机器学习模型训练打下坚实基础。

32

2026.01.31

数据类型有哪几种
数据类型有哪几种

数据类型有整型、浮点型、字符型、字符串型、布尔型、数组、结构体和枚举等。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

338

2023.10.31

php数据类型
php数据类型

本专题整合了php数据类型相关内容,阅读专题下面的文章了解更多详细内容。

225

2025.10.31

c语言 数据类型
c语言 数据类型

本专题整合了c语言数据类型相关内容,阅读专题下面的文章了解更多详细内容。

138

2026.02.12

string转int
string转int

在编程中,我们经常会遇到需要将字符串(str)转换为整数(int)的情况。这可能是因为我们需要对字符串进行数值计算,或者需要将用户输入的字符串转换为整数进行处理。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

1031

2023.08.02

mysql标识符无效错误怎么解决
mysql标识符无效错误怎么解决

mysql标识符无效错误的解决办法:1、检查标识符是否被其他表或数据库使用;2、检查标识符是否包含特殊字符;3、使用引号包裹标识符;4、使用反引号包裹标识符;5、检查MySQL的配置文件等等。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

210

2023.12.04

Python标识符有哪些
Python标识符有哪些

Python标识符有变量标识符、函数标识符、类标识符、模块标识符、下划线开头的标识符、双下划线开头、双下划线结尾的标识符、整型标识符、浮点型标识符等等。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

323

2024.02.23

Python异步编程与Asyncio高并发应用实践
Python异步编程与Asyncio高并发应用实践

本专题围绕 Python 异步编程模型展开,深入讲解 Asyncio 框架的核心原理与应用实践。内容包括事件循环机制、协程任务调度、异步 IO 处理以及并发任务管理策略。通过构建高并发网络请求与异步数据处理案例,帮助开发者掌握 Python 在高并发场景中的高效开发方法,并提升系统资源利用率与整体运行性能。

37

2026.03.12

热门下载

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

精品课程

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

共4课时 | 22.5万人学习

Django 教程
Django 教程

共28课时 | 5万人学习

SciPy 教程
SciPy 教程

共10课时 | 1.9万人学习

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

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