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;
 
}