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