使用python实现一个dicom影像解析入库程序demo

2024-03-14 01:52

本文主要是介绍使用python实现一个dicom影像解析入库程序demo,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

简介

DICOM(Digital Imaging and Communications in Medicine)是医学图像和相关信息的国际标准。它定义了医学影像的格式和通信协议,使得不同设备和系统之间可以交换和共享医学图像和相关数据,如CT扫描、MRI图像、超声波图像等。
DICOM标准包含了多个部分,其中包括数据模型、文件格式、网络通信协议、图像存储和检索、安全性、服务类等方面的内容。DICOM文件通常以.dcm为扩展名,其中包含了医学图像数据以及与之相关的患者信息、设备信息、图像采集参数等。

具体操作

首先,安装程序所用到的库文件,如下所示:
Windows安装

pip install pynetdicom pydicom mysql-connector-python

直接在cmd里面运行这个命令就行。
Ubuntu安装

sudo apt install python3-pip
#采用国内源下载相关库文件
sudo pip install pynetdicom pydicom mysql-connector-python -i https://pypi.tuna.tsinghua.edu.cn/simple
sudo apt install python-is-python3

pynetdicom: pynetdicom 是一个用于实现 DICOM(Digital Imaging and Communications in Medicine)网络通信协议的 Python 库,提供了对 DICOM 标准中定义的各种服务类的支持,如存储服务类(Storage Service Class)、查询服务类(Query/Retrieve Service Class)、检索服务类(Modality Worklist Service Class)等。
pydicom: pydicom 是一个用于处理 DICOM 文件的 Python 库。它提供了读取、写入和操作 DICOM 文件的功能,可以解析DICOM文件的内容,并提取其中的图像数据、患者信息、设备信息等。
mysql-connector-python: mysql-connector-python 是一个用于连接和操作 MySQL 数据库的官方 Python 客户端库。
代码实现

from pynetdicom import AE, evt, build_role, StoragePresentationContexts
from pynetdicom.sop_class import PatientRootQueryRetrieveInformationModelFind, CTImageStorage
from pydicom.dataset import Dataset
import logging
import os
import mysql.connector# 配置日志
logging.basicConfig(level=logging.DEBUG)
logger = logging.getLogger(__name__)# 影像存放的基本路径
base_store_path = "D:/test/images"def handle_store(event):"""处理接收到的DICOM影像并保存到本地路径。"""try:dataset = event.datasetdataset.file_meta = event.file_meta# 为每个检查创建不同的存储路径study_folder = os.path.join(base_store_path, dataset.StudyInstanceUID)os.makedirs(study_folder, exist_ok=True)filename = os.path.join(study_folder, f"{dataset.SOPInstanceUID}.dcm")dataset.save_as(filename, write_like_original=False)logger.info(f"存储DICOM图像: {filename}")# 将DICOM标签信息存储到数据库insert_into_db(dataset)return 0x0000except Exception as e:logger.error(f"存储DICOM图像时出错: {e}")return 0xC210def insert_into_db(dataset):"""将DICOM信息插入到MySQL数据库中。"""try:conn = mysql.connector.connect(host='127.0.0.1',user='root',password='123456',database='dicom')cursor = conn.cursor()# 在插入之前将PersonName类型的字段转换为字符串patient_name = str(dataset.PatientName) if 'PatientName' in dataset else Nonepatient_id = str(dataset.PatientID) if 'PatientID' in dataset else Nonestudy_instance_uid = str(dataset.StudyInstanceUID) if 'StudyInstanceUID' in dataset else Nonesop_instance_uid = str(dataset.SOPInstanceUID) if 'SOPInstanceUID' in dataset else Nonesql = """INSERT INTO dicom_images (StudyInstanceUID, SOPInstanceUID, PatientID, PatientName) VALUES (%s, %s, %s, %s)"""val = (study_instance_uid, sop_instance_uid, patient_id, patient_name)cursor.execute(sql, val)conn.commit()logger.info(f"成功插入记录到数据库: {cursor.rowcount}")except Exception as e:logger.error(f"数据库操作失败: {e}")finally:if conn.is_connected():cursor.close()conn.close()# 初始化应用实体(AE)和配置DICOM服务上下文
ae = AE(ae_title='lay')
ae.add_requested_context(PatientRootQueryRetrieveInformationModelFind)
ae.add_requested_context(CTImageStorage)
ae.supported_contexts = StoragePresentationContexts# 创建C-STORE SCP的角色并启动服务
role = build_role(CTImageStorage, scp_role=True)
handlers = [(evt.EVT_C_STORE, handle_store)]
ae.start_server(('0.0.0.0', 104), evt_handlers=handlers)
logger.info("C-STORE SCP服务启动,准备接收DICOM图像。")

创建表

CREATE TABLE IF NOT EXISTS dicom_images (id INT AUTO_INCREMENT PRIMARY KEY,StudyInstanceUID VARCHAR(255),SOPInstanceUID VARCHAR(255),PatientID VARCHAR(255),PatientName VARCHAR(255),ImagePath VARCHAR(255),InsertedAt TIMESTAMP DEFAULT CURRENT_TIMESTAMP
);

代码详解

配置日志

logging.basicConfig(level=logging.DEBUG)
logger = logging.getLogger(__name__)

这部分代码配置了Python的日志系统,设置了日志级别为DEBUG。这意味着DEBUG及以上级别(INFO, WARNING, ERROR, CRITICAL)的日志都会被记录。logger对象用于在脚本的不同部分记录日志消息。
设置基本路径

base_store_path = "D:/test/images"

定义了一个基础路径,用于保存接收到的DICOM影像文件。
C-STORE响应处理

def handle_store(event):# 省略了具体实现

这个函数处理从DICOM服务器接收到的DICOM影像。它首先创建一个与影像StudyInstanceUID对应的文件夹,然后将DICOM文件保存到这个文件夹中。保存后,它调用insert_into_db函数将影像的部分信息插入到MySQL数据库中。
插入数据库

def insert_into_db(dataset):# 连接到MySQL数据库# 使用转换后的DICOM数据字段执行INSERT操作

这个函数负责将提取的DICOM信息插入到MySQL数据库中。为了避免类型不匹配问题,PatientName等字段在插入之前被显式转换为字符串类型。
初始化应用实体(AE)和启动服务

ae = AE(ae_title='lay')
ae.add_requested_context(PatientRootQueryRetrieveInformationModelFind)
ae.add_requested_context(CTImageStorage)
ae.supported_contexts = StoragePresentationContextsrole = build_role(CTImageStorage, scp_role=True)
handlers = [(evt.EVT_C_STORE, handle_store)]
ae.start_server(('0.0.0.0', 104), evt_handlers=handlers)

初始化一个AE对象,并为其设置一个AE标题(lay)。
添加所需的DICOM服务上下文,使其能够处理C-FIND和C-STORE请求。
通过build_role和start_server启动一个C-STORE SCP服务,该服务在本地机器上监听指定的端口(104),准备接收DICOM影像。
handlers列表定义了事件处理器,用于处理C-STORE相关事件。
运行截图
在这里插入图片描述

今日推荐

歌曲:《不眠之夜》
演唱:张杰
崩坏星穹铁道——匹诺康尼主题曲
在这里插入图片描述

这篇关于使用python实现一个dicom影像解析入库程序demo的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/806869

相关文章

Java中流式并行操作parallelStream的原理和使用方法

《Java中流式并行操作parallelStream的原理和使用方法》本文详细介绍了Java中的并行流(parallelStream)的原理、正确使用方法以及在实际业务中的应用案例,并指出在使用并行流... 目录Java中流式并行操作parallelStream0. 问题的产生1. 什么是parallelS

C++中unordered_set哈希集合的实现

《C++中unordered_set哈希集合的实现》std::unordered_set是C++标准库中的无序关联容器,基于哈希表实现,具有元素唯一性和无序性特点,本文就来详细的介绍一下unorder... 目录一、概述二、头文件与命名空间三、常用方法与示例1. 构造与析构2. 迭代器与遍历3. 容量相关4

Linux join命令的使用及说明

《Linuxjoin命令的使用及说明》`join`命令用于在Linux中按字段将两个文件进行连接,类似于SQL的JOIN,它需要两个文件按用于匹配的字段排序,并且第一个文件的换行符必须是LF,`jo... 目录一. 基本语法二. 数据准备三. 指定文件的连接key四.-a输出指定文件的所有行五.-o指定输出

Java中Redisson 的原理深度解析

《Java中Redisson的原理深度解析》Redisson是一个高性能的Redis客户端,它通过将Redis数据结构映射为Java对象和分布式对象,实现了在Java应用中方便地使用Redis,本文... 目录前言一、核心设计理念二、核心架构与通信层1. 基于 Netty 的异步非阻塞通信2. 编解码器三、

Linux jq命令的使用解读

《Linuxjq命令的使用解读》jq是一个强大的命令行工具,用于处理JSON数据,它可以用来查看、过滤、修改、格式化JSON数据,通过使用各种选项和过滤器,可以实现复杂的JSON处理任务... 目录一. 简介二. 选项2.1.2.2-c2.3-r2.4-R三. 字段提取3.1 普通字段3.2 数组字段四.

C++中悬垂引用(Dangling Reference) 的实现

《C++中悬垂引用(DanglingReference)的实现》C++中的悬垂引用指引用绑定的对象被销毁后引用仍存在的情况,会导致访问无效内存,下面就来详细的介绍一下产生的原因以及如何避免,感兴趣... 目录悬垂引用的产生原因1. 引用绑定到局部变量,变量超出作用域后销毁2. 引用绑定到动态分配的对象,对象

Linux kill正在执行的后台任务 kill进程组使用详解

《Linuxkill正在执行的后台任务kill进程组使用详解》文章介绍了两个脚本的功能和区别,以及执行这些脚本时遇到的进程管理问题,通过查看进程树、使用`kill`命令和`lsof`命令,分析了子... 目录零. 用到的命令一. 待执行的脚本二. 执行含子进程的脚本,并kill2.1 进程查看2.2 遇到的

SpringBoot基于注解实现数据库字段回填的完整方案

《SpringBoot基于注解实现数据库字段回填的完整方案》这篇文章主要为大家详细介绍了SpringBoot如何基于注解实现数据库字段回填的相关方法,文中的示例代码讲解详细,感兴趣的小伙伴可以了解... 目录数据库表pom.XMLRelationFieldRelationFieldMapping基础的一些代

Java HashMap的底层实现原理深度解析

《JavaHashMap的底层实现原理深度解析》HashMap基于数组+链表+红黑树结构,通过哈希算法和扩容机制优化性能,负载因子与树化阈值平衡效率,是Java开发必备的高效数据结构,本文给大家介绍... 目录一、概述:HashMap的宏观结构二、核心数据结构解析1. 数组(桶数组)2. 链表节点(Node

Java AOP面向切面编程的概念和实现方式

《JavaAOP面向切面编程的概念和实现方式》AOP是面向切面编程,通过动态代理将横切关注点(如日志、事务)与核心业务逻辑分离,提升代码复用性和可维护性,本文给大家介绍JavaAOP面向切面编程的概... 目录一、AOP 是什么?二、AOP 的核心概念与实现方式核心概念实现方式三、Spring AOP 的关