QT DICOM医疗影像转换

已于 2023-07-16 07:42:27 修改
阅读量705 收藏 9
点赞数 1
文章标签： qt
于 2023-05-31 22:16:54 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u012156872/article/details/130978168
版权

时间紧任务重，直接上代码吧：

1、引言：借助开源库DCMKT进行转换

2、库文件引用：

#include "dcmtk/config/osconfig.h"
#include "dcmtk/dcmdata/dctk.h"         //DcmFileFormat
#include "dcmtk/dcmimgle/dcmimage.h"    //DicomImage
#include "dcmtk/dcmimgle/diregbas.h"   //DiRegisterBase
#include "dcmtk/dcmimage/diregist.h"   //DiRegisterBase
#include "dcmtk/dcmdata/libi2d/i2d.h"   //I2DImgSource
#include "dcmtk/dcmdata/libi2d/i2dbmps.h"   //I2DBmpSource
#include "dcmtk/dcmdata/libi2d/i2djpgs.h"   //I2DJpegSource
#include "dcmtk/dcmdata/dcpxitem.h"   //DcmPixelItem
#include "dcmtk/dcmdata/dcpath.h"   //DcmPixelItem
#include "dcmtk/dcmdata/dcrledrg.h"   //DcmRLEDecoderRegistration
#include "dcmtk/dcmjpeg/djdecode.h"   //DJDecoderRegistration
#include "dcmtk/dcmjpeg/dipijpeg.h"   //DiJPEGPlugin

3、DCOM影像的显示：

QString strFilePath = QFileDialog::getOpenFileName(this, "打开对话框", QDir::currentPath(),"DICOM图像文件(*.dcm)");
    if (!strFilePath.isEmpty()){
        ui->lineEdit_dcm->setText(strFilePath);
        DicomImage *dcmImage = new DicomImage(ui->lineEdit_dcm->text().toLocal8Bit().toStdString().c_str());
        if (dcmImage != NULL)
        {
            if (dcmImage->getStatus() == EIS_Normal)
            {
                EP_Interpretation PhotometricInterpretation = dcmImage->getPhotometricInterpretation();
                if ((PhotometricInterpretation == EPI_Monochrome1) || (PhotometricInterpretation == EPI_Monochrome2))
                {
                    dcmImage->setMinMaxWindow();
                    ulong size = dcmImage->getWidth() * dcmImage->getHeight();
                    uchar *data = new uchar[size];
                    dcmImage->setMinMaxWindow();
                    dcmImage->getOutputData(data, size, 8);
                    QImage image = QImage(data, dcmImage->getWidth(), dcmImage->getHeight(), QImage::Format_Indexed8);
                    QImage image2 = image.copy();
                    image2 = image2.scaled(ui->label_dcm->width(), ui->label_dcm->height(), Qt::IgnoreAspectRatio, Qt::SmoothTransformation);
                    ui->label_dcm->setPixmap(QPixmap::fromImage(image2));

                    delete[] data;
                    data = nullptr;
                }
                else if (PhotometricInterpretation == EPI_RGB){
                    DcmFileFormat dicomfile;//将dcm格式文件读取出来解析相关类
                    OFCondition oc = dicomfile.loadFile(ui->lineEdit_dcm->text().toLocal8Bit().data());
                    if(oc.bad()){
                        QMessageBox::information(this, "提示",  "转换失败！");
                        delete dcmImage;
                        dcmImage = nullptr;
                        return;
                    }

                    DcmDataset* pDataset = dicomfile.getDataset();//得到dcm文件的数据

                    int nWidth, nHeight, nBandCount, nWindowWidth, nWindowCenter,nRescaleSlope,nRescaleIntercept,nBitsAllocated,nHighBit,
                        nNumberOfFrames;
                    OFString strInstanceNumber, strWidth, strHeight,
                        strWindowWidth, strWindowCenter, strSamplesPerPixel,
                        strRescaleSlope,strRescaleIntercept,strBitsAllocated,strHighBit,strNumberOfFrames;

                    //读取dcm格式文件里相关元素值
                    pDataset->findAndGetOFString(DCM_InstanceNumber, strInstanceNumber);
                    pDataset->findAndGetOFString(DCM_Columns ,strWidth);
                    pDataset->findAndGetOFString(DCM_Rows ,strHeight);
                    pDataset->findAndGetOFString(DCM_WindowWidth ,strWindowWidth);
                    pDataset->findAndGetOFString(DCM_WindowCenter ,strWindowCenter);
                    pDataset->findAndGetOFString(DCM_SamplesPerPixel ,strSamplesPerPixel);
                    pDataset->findAndGetOFString(DCM_RescaleSlope ,strRescaleSlope);
                    pDataset->findAndGetOFString(DCM_RescaleIntercept,strRescaleIntercept);
                    pDataset->findAndGetOFString(DCM_BitsAllocated,strBitsAllocated);
                    pDataset->findAndGetOFString(DCM_HighBit,strHighBit);
                    pDataset->findAndGetOFString(DCM_NumberOfFrames,strNumberOfFrames);

                    nWidth = atoi(strWidth.c_str());
                    nHeight = atoi(strHeight.c_str());
                    if(strSamplesPerPixel.empty())
                    {
                        strSamplesPerPixel = "1";//若没有说明哪种类型图像，默认为灰度图像
                    }

                    if(strRescaleSlope.empty())
                    {
                        strRescaleSlope = "1";//像素转换的比例系数，没有则默认存储为1
                    }

                    if(strRescaleIntercept.empty())
                    {
                        strRescaleIntercept = "0";//像素转换的比例参数，没有则默认存储为0
                    }

                    nBandCount = atoi(strSamplesPerPixel.c_str());
                    nWindowWidth = atoi(strWindowWidth.c_str());
                    nWindowCenter = atoi(strWindowCenter.c_str());
                    nRescaleSlope = atoi(strRescaleSlope.c_str());
                    nRescaleIntercept = atoi(strRescaleIntercept.c_str());
                    nBitsAllocated = atoi(strBitsAllocated.c_str());
                    nHighBit = atoi(strHighBit.c_str());
                    nNumberOfFrames = atoi(strNumberOfFrames.c_str());

                    E_TransferSyntax opt_oxfer = pDataset->getOriginalXfer();
                    //图像存储的压缩格式
                    if(EXS_RLELossless == opt_oxfer)
                    {
                        E_TransferSyntax opt_oxfer = EXS_RLELossless;
                        OFBool opt_uidcreation = OFFalse;
                        OFBool opt_reversebyteorder = OFFalse;
                        DcmRLEDecoderRegistration::registerCodecs(opt_uidcreation, opt_reversebyteorder);
                        DcmXfer original_xfer(opt_oxfer);
                        if (!original_xfer.isEncapsulated()){
                            QMessageBox::information(this, "提示",  "转换失败！");
                            delete dcmImage;
                            dcmImage = nullptr;
                            return;
                        }

                        oc = pDataset->chooseRepresentation(EXS_LittleEndianExplicit, NULL);
                        if (oc != EC_Normal){
                            QMessageBox::information(this, "提示",  "转换失败！");
                            delete dcmImage;
                            dcmImage = nullptr;
                            return;
                        }
                    }

                    else if(opt_oxfer >= EXS_JPEGProcess1 && opt_oxfer <= EXS_JPEGLSLossy)
                    {
                        E_DecompressionColorSpaceConversion opt_decompCSconversion = EDC_guessLossyOnly;
                        E_UIDCreation    opt_uidcreation = EUC_default;
                        DJDecoderRegistration::registerCodecs(opt_decompCSconversion, opt_uidcreation);
                        DcmXfer original_xfer(opt_oxfer);
                        if (!original_xfer.isEncapsulated()){
                            QMessageBox::information(this, "提示",  "转换失败！");
                            delete dcmImage;
                            dcmImage = nullptr;
                            return;
                        }

                        oc = pDataset->chooseRepresentation(EXS_LittleEndianExplicit, NULL);
                        if(oc != EC_Normal){
                            QMessageBox::information(this, "提示",  "转换失败！");
                            delete dcmImage;
                            dcmImage = nullptr;
                            return;
                        }
                    }

                    DcmElement* element = NULL;
                    OFCondition result = pDataset->findAndGetElement(DCM_PixelData, element);//获取文件中存储的像素元素
                    long nCount = nWidth * nHeight;//unsigned char
                    if(nBandCount == 3)
                        nCount *= 3;//若是RGB彩色图像则有R,G,B三种色彩值进行存储

                    /*获取的像素值需要进行转化
                     * 实际像素值 = RescaleSlope * 获取像素值 + RescaleIntercept，
                     * 若BitsAllocated超过8会超出颜色值的表示范围0~255需要对其进行转化.
                     */
                    unsigned char *pRgbData = new unsigned char[nCount];
                    if(8 == nBitsAllocated)
                    {
                        Uint8* pixData = NULL;
                        result = element->getUint8Array(pixData);
                        if(result.good() && pixData != NULL)
                        {
                            if(pRgbData != NULL)
                            {
                                for(long i = 0; i < nCount; i++)
                                {
                                    while(nBitsAllocated - nHighBit > 0)
                                    {
                                        pixData[i] &= ~(1 << nHighBit);
                                        nHighBit++;
                                    }

                                    pRgbData[i] = nRescaleSlope*pixData[i] + nRescaleIntercept;
                                }
                            }

                            QImage image = QImage((uchar*)pRgbData, nWidth, nHeight, nBandCount==1 ? QImage::Format_Indexed8 : QImage::Format_RGB888);
                            QImage image2 = image.copy();
                            image2 = image2.scaled(ui->label_dcm->width(), ui->label_dcm->height(), Qt::IgnoreAspectRatio, Qt::SmoothTransformation);
                            ui->label_dcm->setPixmap(QPixmap::fromImage(image2));
                        }
                    }
                    else if(16 == nBitsAllocated)//若颜色值存储为16位，则有可能超出颜色的表示范围(0~255)此时需要对其进行转化
                    {
                        Uint16* pixData = NULL;
                        result = element->getUint16Array(pixData);
                        if(result.good() && pixData != NULL)
                        {
                            if(pRgbData != NULL)
                            {
                                //short *pshort = (short *)pixData;
                                for(long i = 0; i < nCount; i++)
                                {
                                    while(nBitsAllocated - nHighBit > 0)
                                    {
                                        pixData[i] &= ~(1 << nHighBit);
                                        nHighBit++;
                                    }

                                    int data= nRescaleSlope*pixData[i] + nRescaleIntercept;
                                    //float fdata = float (nRescaleSlope*pixData[i] + nRescaleIntercept);
                                    if(data < nWindowCenter - nWindowWidth/2)
                                    {
                                        pRgbData[i] = 0;
                                    }
                                    else if(data > nWindowCenter + nWindowWidth/2)
                                    {
                                        pRgbData[i] = 0xFF;
                                    }
                                    else
                                    {
                                        double dValue = (data + nWindowWidth/2.0f - nWindowCenter)*255.0f/nWindowWidth;

                                        int nGray = 0;
                                        if(dValue < 0)
                                        {
                                            nGray = 0;
                                        }
                                        else if(dValue > 255)
                                        {
                                            nGray = 255;
                                        }
                                        else
                                        {
                                            nGray = (int)dValue;
                                        }

                                        pRgbData[i] = nGray;
                                    }
                                }

                                QImage image = QImage((uchar*)pRgbData, nWidth, nHeight, nBandCount==1 ? QImage::Format_Indexed8 : QImage::Format_RGB888);
                                QImage image2 = image.copy();
                                image2 = image2.scaled(ui->label_dcm->width(), ui->label_dcm->height(), Qt::IgnoreAspectRatio, Qt::SmoothTransformation);
                                ui->label_dcm->setPixmap(QPixmap::fromImage(image2));
                            }
                        }
                    }

                    if(NULL != pRgbData)
                    {
                        delete []pRgbData;
                        pRgbData = NULL;
                    }

                    if(EXS_RLELossless == opt_oxfer)
                    {
                        DcmRLEDecoderRegistration::cleanup();
                    }

                    else if(opt_oxfer >= EXS_JPEGProcess1 && opt_oxfer <= EXS_JPEGLSLossy)
                    {
                        DJDecoderRegistration::cleanup(); // deregister JPEG codecs
                    }
                }
            }
            else {
                QMessageBox::information(this, "提示",  QString("Error: cannot load DICOM image (") + QString(DicomImage::getString(dcmImage->getStatus())) + QString(")"));
            }

            delete dcmImage;
            dcmImage = nullptr;
        }
    }

4、DCOM影像的转换：

QString strFilePath = QFileDialog::getSaveFileName(this, "保存对话框", QDir::currentPath(),"DICOM图像文件(*.dcm)");
    if (!strFilePath.isEmpty()){
        QString strDcmFilePath = strFilePath;//ui->lineEdit_raw->text().replace(".raw",".dcm");
        if (!strDcmFilePath.endsWith(".dcm", Qt::CaseInsensitive))
            strDcmFilePath += ".dcm";

        QFile file(ui->lineEdit_raw->text()/*"template.raw"*/);
        if (file.open(QIODevice::ReadOnly)){
            int width = 0;
            int height = 0;
            file.read((char*)&width, sizeof(int));
            file.read((char*)&height, sizeof(int));
            QByteArray pixelData = file.readAll();
            file.close();

            Uint16 samplesPerPixel = 3;
            Uint16 rows = height;
            Uint16 cols = width;
            Uint16 bitsAlloc = 8;
            Uint16 bitsStored = 8;
            Uint16 highBit = 7;
            OFString photoMetrInt = "RGB";

            DcmFileFormat fileFormat;
            OFCondition oc = fileFormat.loadFile("template.dcm");//DCOM影像有些标签需要填写，这里为了省事，直接从其他影像文件获取吧
            if (oc.bad()){
                QMessageBox::information(this, "提示",  "转换失败！");
                return;
            }

            DcmDataset *dataSet = fileFormat.getDataset();
            dataSet->putAndInsertUint16(DCM_SamplesPerPixel, samplesPerPixel);//3
            dataSet->putAndInsertString(DCM_NumberOfFrames, "1");   //1
            dataSet->putAndInsertUint16(DCM_Rows, rows);    //3200
            dataSet->putAndInsertUint16(DCM_Columns, cols); // 4000
            dataSet->putAndInsertUint16(DCM_BitsAllocated, bitsAlloc);  //8
            dataSet->putAndInsertUint16(DCM_BitsStored, bitsStored);    //8
            dataSet->putAndInsertUint16(DCM_HighBit, highBit);  //7
            dataSet->putAndInsertOFStringArray(DCM_PhotometricInterpretation, "RGB");//RGB

            dataSet->putAndInsertUint8Array(DCM_PixelData, OFreinterpret_cast(Uint8*, pixelData.data()), pixelData.size());

            DcmFileFormat dicomfile(dataSet);
            oc = dicomfile.saveFile(strDcmFilePath.toLocal8Bit().toStdString().c_str());
            if (oc.bad()){
                QMessageBox::information(this, "提示",  "转换失败！");
                return;
            }
            else{
                QMessageBox::information(this, "提示",  "转换成功！");
            }
        }
    }

老狼主
关注
  • 1
    点赞
  • 踩
  • 9
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
博客
TightVNC H264编解码（一）
12-31 1548
经过尽二个月的不断摸索，TightVNC终于支持H264编解码了，前期真正编写H264编解码器只花了一周左右时间，但是测试发现效果并不是太理想（帧率太低，延迟大），这大概是所有远程控制软件要处理的通病吧。作为一位经验丰富的程序员，造成这种问题的原因，有以下几点：1、抓屏？2、H264编码？3、TCP网络发送？4、H264解码？5、视频渲染？为了排除到底哪个环节出了问题，看来不摸透底层事务流程是行不通的，不得已摸透代码画出了下面这个简陋的流程图。Server端事务处理流程图.
博客
QT HttpServer
07-31 855
资源从网上收集，包含httpserver、sslserver、http-parser三个模块，解决了编译过程中由于路径引用错误问题。下载链接：https://download.csdn.net/download/u012156872/20664799，下载编译即可应用到自己项目中，简单快速，不用造车。...
博客
QT 界面设计篇（水波纹进度条QProgressBarWater）
11-30 1666
效果图：头文件声明（QProgressBarWater.h）：#ifndef QPROGRESSBARWATER_H#define QPROGRESSBARWATER_H#include <QWidget>#include <QPainter>#include <QPaintEvent>class QProgressBarWater : public QWidget{ Q_OBJECTpublic: explic.
博客
EasyPR编译指南
08-02 3221
1、下载源码https://github.com/liuruoze/EasyPRDownloadZIP或者git（https://github.com/liuruoze/EasyPR.git）将代码解压到硬盘目录下。例如F:\车牌识别\EasyPR（车牌识别）\EasyPR2、双击解决方案文件EasyPR.sln，用VSIDE（这里使用VS2019）打开，提示点击确定。3、选择libeasypr项目，编译出现如下错误这是因为没有配置OpenCV环境所致。4、...
博客
C++ linux下FFmpeg转发摄像头视频流(续)
06-19 149
本篇文章介绍使用v4l2接口捕获摄像头图像数据进行ffmpeg转发。
博客
C++ linux下FFmpeg转发摄像头视频流
06-19 205
在liunx下，捕获摄像头视频流有2种方案，一种是使用ffmpeg捕获流，另一种是使用v4l2接口捕获视频流。前者使用方法简单，但是有后遗症，并不是所有型号摄像头都能够获取成功，这就是我写这篇文章的目的。老规矩，直接粘贴代码(这里虽然是使用Qt代码格式，稍加变换，即可支持其他平台)。源文件申明（PushVideoThread.cpp）
博客
QT MySQL驱动
05-14 127
在某些场合QT自带的QSQLDatabase驱动有效，但是open的时候，程序直接崩溃了，原因未知。不得已将上次使用过的。根据QT接口进行改动，进行无缝衔接。
博客
QT GDAL加载显示tif图像
05-11 242
2、加载tiff图像文件。在此之前记得引用头文件。3、GDAL库资源释放。1、GDAL库初始化。
博客
QT vlc音视频开发踩坑
05-11 192
为此我检查错误耗时将近1整天时间，甚至于都一度怀疑人生（之前在VC++下开发一切顺利）。SDK初始化一切顺利，最终在调用如下代码的时候始终失败，返回NULL值。问题解决，原来如此简单啊！
博客
QT 在线安装失败
05-07 231
QT在线安装使用的国外服务器，所以导致安装极其缓慢或者失败，唯一能想到的解决方法就是切换到国内的镜像源（除非翻墙）。
博客
QT 无标题栏窗口的拖动和拉伸
12-14 513
如下2种方案，2种方案都不完美（拉伸窗口都会抖动），但是如果一定要使用，强烈推荐第1种方案（可减小拉伸时窗口抖动）
博客
C++ libVNC开源库
12-11 608
VNC具体用途，不用过多介绍了。如果需要实时更新服务器端画面，while循环就可以了。
博客
QT QCefView自编译
12-09 574
6、回到QCefView目录，执行批处理generate-win-x86_64.bat即可，出现如下信息表示生成成功。5、设置QT版本，修改配置文件QCefView/cmake/QtConfig.cmake，如下。4、请根据对应的版本修改编译文件CefViewCore/CMakeLists.txt，如下。2、下载QCefViewCore，放到QCefViewCore目录下，地址。请注意最新版本兼容性可能存在问题，我下载的版本是。1、下载QCefView，地址。3、下载Cef库文件，地址。
博客
Ubuntu 配置共享文件夹
10-26 273
安装虚拟机经常会遇到已经共享了文件夹，但是就是找不到目录文件，搜集网上资料，汇总一下指令，留作以后备用吧：vmware-hgfsclient在文件最后一行添加如下内容：
博客
C++ OleDB数据库驱动（MySql驱动姊妹篇）
10-11 240
这里只介绍OelDB调用数据的方法。OleDB因为兼容ODBC和ADO，所以受大众追捧。大家也都非常熟悉了，这里就直接上代码吧。
博客
Visual Assist X插件安装失败总结
09-27 1077
2、安装成功以后，却显示插件无法启用，或者是插件启动了，但是实则没有启动，看不到插件选项。1、如果是运行VA_X_SETUP.exe无反应，无法启动，请右键属性，修改兼容模式。3、此时重启应该就正常了。
博客
常用时间日期格式符号定义
09-27 2324
编程过程中经常会用到时间日期格式化，这里是一些常见的日期时间格式字符定义的列表：年份：yyyyyy月份：MMMMMM天数：ddd小时：hhh分钟：mmm秒数：sss毫秒：zzzz星期几：dddddddAP时区：t。
博客
C++ linux下使用X11实现屏幕截图（续）
09-18 291
使用libpng库将XGetImage获取到的内容保存为PNG格式的文件。
博客
C++ MFC中如何根据不同文件编码格式读取和显示文件内容
09-17 804
这里假设VC++的字符集编码格式是Unicode，要显示的文本内容是strResult。上一篇，介绍了写文件，这一篇我要介绍的如何读取文件内容并显示出来。3、读取unicode-le编码格式的文本文件。4、读取unicode-be编码格式的文本文件。2、读取utf-8编码格式的文本文件。1、读取ANSI编码格式的文本文件。
博客
C++ MFC中如何根据不同文件编码格式创建和写文本内容
09-17 572
这里假设VC++的字符集编码格式是Unicode，要写入的文本内容是m_strTEXT。3、保存unicode-le编码格式的文本文件。4、保存unicode-be编码格式的文本文件。2、保存utf-8编码格式的文本文件。1、保存ANSI编码格式的文本文件。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值