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

news2024/11/19 21:20:36

简介

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.dataset
        dataset.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 0x0000
    except Exception as e:
        logger.error(f"存储DICOM图像时出错: {e}")
        return 0xC210

def 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 None
        patient_id = str(dataset.PatientID) if 'PatientID' in dataset else None
        study_instance_uid = str(dataset.StudyInstanceUID) if 'StudyInstanceUID' in dataset else None
        sop_instance_uid = str(dataset.SOPInstanceUID) if 'SOPInstanceUID' in dataset else None

        sql = """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 = StoragePresentationContexts

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)

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

今日推荐

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

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1514407.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

代码随想录算法训练营第七天| 454.四数相加II、383.赎金信、15.三数之和、18.四数之和

系列文章目录 目录 系列文章目录454.四数相加II使用HashMap法 383.赎金信哈希解法(数组) 15.三数之和双指针法 18.四数之和双指针法 454.四数相加II 题解:该题和1.两数之和的方法是一样的,这个题的难点在于key和value分别是什么。…

发那科数控机床FanucCNC(NCGuide)仿真模拟器配置和数据采集测试

开发日记3.12 此篇用于记录发那科数控机床(Fanuc CNC)采集程序开发中,用虚拟机做测试时,虚拟机的配置和使用以支持采集软件开发和测试。 配置虚拟机使用仿真软件 下载VMware15 「链接:https://pan.xunlei.com/s/VNsl9Gmb14ANBiiNlsT7vA2LA…

Day15 面向对象进阶——接Day14

Day15 面向对象进阶——接Day14 文章目录 Day15 面向对象进阶——接Day14一、访问修饰符二、Object三、深入String的equals()方法四、final 一、访问修饰符 1、含义:修饰类、方法、属性,定义使用的范围 2、经验: 2.1.属性一般使用private修…

springboot3 打包报错32-bit architecture x86 unsupported或者 returned non-zero result

springboot3 打包异常情况处理记录 在测试springboot3 native打包时候遇到的异常,百度和谷歌上方法都无法解决我的问题,最后记录一下我最后的原因和解决方案。 前置要求:自己处理好vs的相关内容后 报错一: [1/7] Initializing…

Vue3 前端生成随机id( 生成 UUID )

效果展示 封装工具(代码展示) 重新创建一个文件**/utils/someTools.js**,并在里面写入如下代码。 function Tools() {}Tools.prototype.guid function () {return xxxxxxxx-xxxx-4xxx-yxxx-xxxxxxxxxxxx.replace(/[xy]/g, function (c) {v…

Android SDK 开发 云托管

开发SDK,拒绝重复造轮子。 本文陈述两种SDK开发方式,第一种AAR方式;第二种远程依赖方式。 具体步骤分为: 一、如何开发SDK? 二、如何打包AAR? 三、如何打包AAR? 四、如何进行SDK远程托管? 五、如…

工业物联网平台在水务环保、暖通制冷、电力能源等行业的应用

随着科技的不断发展,工业物联网平台作为连接物理世界与数字世界的桥梁,正逐渐成为推动各行业智能化转型的关键力量。在水务环保、暖通制冷、电力能源等行业,工业物联网平台的应用尤为广泛,对于提升运营效率、降低能耗、优化管理等…

toJSONString空值被忽略解决办法 toJSONString过程中时间格式丢失问题解决办法

toJSONString空值被忽略解决办法 原代码: GcGarbageBuildingDO data gcGarbageBuildingService.get(id); result.putAll(JSON.parseObject(JSON.toJSONString(data), Map.class));查询出来的data中部分字段值为null,在toJSONString的过程中会把null值…

基于STM32的感应开关垃圾桶

1.定时器介绍 1.1 工作原理 使用精准的时基,通过硬件的方式,实现定时功能。 1.2 定时器分类 基本定时器(TIM6~TIM7)通用定时器(TIM2~TIM5)高级定时器(TIM1和TIM8) 1.3 通用定时器…

【C语言】Windows下的C语言线程编程详解

文章目录 1. 头文件1.1 windows.h1.2 process.h 2. 创建线程3. 线程同步3.1 线程同步方式3.1 互斥量(Mutex)3.2 事件(Event) 4. 线程的结束与资源管理5.线程池(简要) 在Windows平台下,C语言提供…

隧道技术和代理技术(三)

隧道技术 知识点 -隧道技术:解决不出网协议上线的问题(利用出网协议进行封装出网) -代理技术:解决网络通讯不通的问题(利用跳板机建立节点后续操作) 内环境示意图,方便理解 思路:…

【网络原理】UDP协议的详细解析

文章目录 🎄UDP协议概念🌲UDP协议端格式🎄UDP的特点🌸无连接🌸不可靠传输🌸面向数据报🌸缓冲区🌸全双工🌸大小受限 🍀基于UDP的应用层协议🎍扩展问…

(BAT向)Java岗常问高频面试汇总:MyBatis 微服务 Spring 分布式 MySQL等

看面试题可以是为了面试,也可以是对自己学到的东西的一种查漏补缺,更加深刻的去了解一些核心知识点 Spring面试高频问题 问题一:谈 需要zi料 绿色徽【vip1024b】 谈你对spring IOC 和 DI 的理解,它们有什么区别? **问题…

全球首个AI程序员Devin诞生,真的不再需要程序员了吗?

前言 今年到底是怎么回事?前有Sora(首个文生视频大模型)独领风骚,后有Mistral(号称世界第二语言大模型)横刀立马,甚至Claude3的热度都还没有过,今天一则更炸裂的消息就爆出,世界上第一个AI程序…

代码贴--链表--数据机构

本博客将记录链表代码(单链表)&#xff0c;后续其他链表和其他数据结构内容请看我的其他博客 头文件(SList.h) #pragma once #include<iostream> #include<bits/stdc.h> using namespace std;typedef int SLTDataType;struct SListNode {int data;struct SListNo…

windows系统玩游戏找不到d3dx9_43.dll缺失,无法启动此程序的解决方法

今日&#xff0c;我们要深入讨论d3dx9_43.dll文件的重要性及其缺失问题。最近&#xff0c;我也遇到了这个文件丢失的困扰&#xff0c;因此想借此机会与大家分享如何解决d3dx9_43.dll缺失的问题。 一.电脑d3dx9_43.dll丢失会提示什么&#xff1f; 关于电脑提示d3dx9_43.dll丢失…

【人力资源开发】某主题公园人力资源开发管理咨询项目纪实

虽然很多企业将“人事部”改为“人力资源部”&#xff0c;但是企业的人力资源管理水平却仍停留在“人事管理”的阶段。该主题公园也是如此。随着公园的不断发展&#xff0c;其人力资源管理问题逐渐显露&#xff0c;而管理者也不清楚问题的根源在哪里&#xff0c;只能采取“头疼…

AWS入门实践-AWS CLI工具的使用介绍

AWS CLI&#xff08;Amazon Web Services Command Line Interface&#xff09;是一个强大的工具&#xff0c;它允许您直接从命令行与AWS服务进行交互。这不仅可以加快许多任务的处理速度&#xff0c;而且还可以通过脚本自动化。 一、AWS CLI工具的安装 1、Windows 安装下载…

java-双列集合

什么是双列集合&#xff1f; 集合中每次存的数据是成对存入的 以及它的特点是什么&#xff1f; 特别注意&#xff1a;键不可重复&#xff0c;值可以 Map是双列集合的顶层接口 Map 它有哪些方法呢&#xff1f; Map的常用API 添加 添加操作的代码如下 我们要明白一些细节&…

【项目笔记】java微服务:黑马头条(day02)

文章目录 app端文章查看&#xff0c;静态化freemarker,分布式文件系统minIO1)文章列表加载1.1)需求分析1.2)表结构分析1.3)导入文章数据库1.3.1)导入数据库1.3.2)导入对应的实体类 1.4)实现思路1.5)接口定义1.6)功能实现1.6.1)&#xff1a;导入heima-leadnews-article微服务&am…