使用dcmtk获取dicom影像的像素值

#include "widget.h"

#include "ui_widget.h"

#include "dcmtk/config/osconfig.h"    /* make sure OS specific configuration is included first */

#include "dcmtk/dcmdata/dctk.h"

#include "dcmtk/dcmdata/dcpxitem.h"

#include "dcmtk/dcmimgle/dcmimage.h"

#include 

uint getPixelValue(DicomImage *image,long x, long y, EP_Representation &r)

{

    if (image) {

        const DiPixel* pixel = image->getInterData();

        if (pixel && (x < (long)image->getWidth()) && (x >= 0)

                && (y < (long)image->getHeight()) && (y >= 0)) {

            r = pixel->getRepresentation();

            switch (r) {

            case EPR_Sint8:

            case EPR_Uint8:

                return *((char*)(pixel->getData()) + (y * image->getWidth() + x));

            case EPR_Sint16:

            case EPR_Uint16:

                return *((short*)(pixel->getData()) + (y * image->getWidth() + x));

            case EPR_Sint32:

            case EPR_Uint32:

                return *((int*)(pixel->getData()) + (y * image->getWidth() + x));

            }

        }

    }

    r = (EP_Representation)-1;

    return 0;

}

Widget::Widget(QWidget *parent) :

    QWidget(parent),

    ui(new Ui::Widget)

{

    ui->setupUi(this);

    DicomImage* dcmImage = new DicomImage("D:/dcm/dcm_rc1/57DF1DB1");

    Sint16 *pixelData = (Sint16 *)(dcmImage->getOutputData(16));

    for(int i=0;i<100;i++)

    qDebug()<<(int)pixelData[i];

    double min,max;

    dcmImage->getMinMaxValues(min,max);

    qDebug()< 
  
    const DiPixel* pixel=dcmImage->getInterData();
 
  
    qDebug()<getRepresentation();
 
  
    int x=0;
 
  
    int y=0;
 
  
    uint data=*((short*)(pixel->getData()) + (y * dcmImage->getWidth() + x));
 
  
    qDebug()<<"data is:"<<(int)data;
 
  
    EP_Representation r;
 
  
    uint data2=getPixelValue(dcmImage,x,y,r);
 
  
    qDebug()<<"data2 is:--------------------------"<<(int)data2;
 
  
    for(int i=0;i<20;i++)
 
  
    {
 
  
        for(int j=0;j<20;j++)
 
  
        {
 
  
            //uint data3=getPixelValue(dcmImage,i,j,r);
 
  
            //qDebug()<<(qint16)data3;
 
  
            uint data=*((short*)(pixel->getData()) + (i * dcmImage->getWidth() + j));
 
  
            qDebug()<<"data is:"<<(qint16)data;
 
  
        }
 
  
    }
 
  
}
 
  
 
 
  
Widget::~Widget()
 
  
{
 
  
    delete ui;
 
  
}