カラービットマップ

pixelDataからカラービットマップを取得
BGR > RGB の変換が必要

{
if(PhotometricInterpretation != "RGB")
    retern;
//dataSet: DicomDataSet
Bitmap bitmap;
int columns = dataSet.Columns;
int rows = dataSet.Rows;
byte[] pixelData = dataSet.PixelData;
bitmap = new Bitmap(columns, rows, PixelFormat.Format24bppRgb);
BitmapData bitmapData = bitmap.LockBits(new Rectangle(0, 0, columns, rows), 
ImageLockMode.WriteOnly, PixelFormat.Format24bppRgb);
int padding = bitmapData.Stride - columns;
byte* pBitmap = (byte*)bitmapData.Scan0;
int j = 0, i = 0;
for (int y = 0; y < rows; y++)
{
    j = y * columns;
    for (int x = 0; x < columns; x++)
    {
        i = j + x;
        pBitmap[i * 3]     = pixelData[i * 3 + 2];   //G
        pBitmap[i * 3 + 1] = pixelData[i * 3 + 1];   //B
        pBitmap[i * 3 + 2] = pixelData[i * 3];       //R
    }
    i += padding;
}
bitmap.UnlockBits(bitmapData);
return bitmap;
}

.net 4.0から Parallel.For が使えるようになったので一部訂正。2015-08-12

fixed (byte* pData = pixelData)
{
    byte* rgb = pData;
    int stride = bmpData.Stride;
    Parallel.For(0, rows, y =>
    {
        for (int x = 0; x < columns; x++)
        {
            int pos = (x * 3) + y * stride;
            pBitmap[pos]     = rgb[pos + 2];   //G  
            pBitmap[pos + 1] = rgb[pos + 1];   //B
            pBitmap[pos + 2] = rgb[pos];       //R
        }
    });
}

マルチフレーム

NumberOfFrames タグ (0028,0008)があるとマルチフレーム画像。
PixelData(7FD0,0010)の長さをNumberOfFramesで割った値が一枚当たりの画像の PixelData 。

List<byte[]> pixelDataList = new List<byte[]>();

int numberOfFrames = dicomData.NumberOfFrames;

int length = dicomData.PixelData.Length / numberOfFrames;

for(int k=0; k < numberOfFrames;k++)

{

    byte[] buf;

    Array.Copy(dicomData.PixelData, numberOfFrames * length, buf, 0, length);

    pixelDataList.Add(buf);

}

RescaleIntercept タグ と RescaleSlopeタグ

他院からのCDで Rescale Interceptタグ と　Rescale Slopeタグがある DICOM Fileを散見。
シーメンスのCTなど。
画像データに傾斜と切片を付加して8bitグレースケールを得る。

fixed (byte* pData = pixelData)

{

    int pixel;

    short* ptr = (short*)pData;             //short型ポインタ

    byte* pBitmap = (byte*)bitmapData.Scan0;

    int j = 0, i = 0;

    int difference = max - min;

    if (rescaleIntercept == 0.0d )

    {

        for (int x = 0; x < _columns; x++)

        {

            pixel = (int)(((*ptr + rescaleIntercept - windowCenter) / windowWidth + 0.5d) * 255.0d);

            if (pixel > 255) pixel = 255;

            else if (pixel < 0) pixel = 0;

            *pBitmap = (byte)pixel;             ptr++;             pBitmap++;

        }

        pBitmap += padding;

    }

    else

    {

        for (int y = 0; y < _rows; y++)

        {

            j = y * _columns;

            for (int x = 0; x < _columns; x++)

            {

                pixel = (int)((((int)ptr[j + x] + rescaleIntercept - windowCenter) / windowWidth + 0.5d) * 255.0d);

                if (pixel > 255) pixel = 255;

                if (pixel < 0) pixel = 0;

                pixel = (int)(rescaleSlope * pixel + rescaleIntercept);

                pBitmap[i] = (byte)pixel;

                i++;

           }

           i += padding;

       }

    }

}

16 bit グレースケール

CT, MR, CR, DR などの医用画像は殆どが16bit グレースケール画像。
PCでは8bitしか表示できないので、これを変換。
また Bitmap.SetPixel　がとても遅いため　LockBits　を使用。

using System;
using System.Drawing;
using System.Drawing.Imaging;

DICOMオブジェクト

int columns = dicomData.Columns;
int rows = dicomData.Rows;
int windowCenter = dicomData.WindowCenter;
int windowWidth = dicomData.WindowWidth;

ビットマップ
 1byteが1pixelのグレースケール用ColorPaletteを作成する。
 3byte=1pixelのグレースケール（R=G=B）より容量減か。

byte[] pixelData = dicomData.PixelData;Bitmap bitmap;
ColorPalette colorPalette;
for (int i = 0; i < 256; i++)
    colorPalette.Entries[i] = Color.FromArgb(i, i, i);
Bitmap bitmap = new Bitmap(columns, rows, PixelFormat.Format8bppIndexed);

データの挿入 ()

public unsafe Bitmap Grayscale16To8()

{

int padding = bitmapData.Stride - columns;

int min = windowCenter - windowWidth / 2;

int max = windowCenter + windowWidth / 2;

fixed (byte* pData = pixelData)

{

    int pixel;

    short* ptr = (short*)pData;

    byte* pBitmap = (byte*)bitmapData.Scan0;

    int j = 0, i = 0;

    for (int y = 0; y < _rows; y++)

    {

        j = y * _columns;

        for (int x = 0; x < _columns; x++)

        {

            pixel = (int)ptr[j + x];   

            if (pixel <= min)

                pixel = 0;

            else if (pixel >= max)

                pixel = 255;

            else

                pixel = (int)((pixel - min) * 255 / (max - min));

            pBitmap[i] = (byte)pixel;

            i++;

        }

        i += padding;

    }

}

bitmap.UnlockBits(bitmapData);

return bitmap;

}