DICOM图像像素值、灰度值与CT值

做医学图像处理会涉及DICOM图像，初学者应该要对图像像素有一定了解，在对像素处理中会一些算法会利用的更好。

图像灰度值的概念是什么？灰度也可以认为是亮度，简单说就是色彩的深浅程度。

    实际上在我们的日常生活中，通过三原色色彩深浅的组合，可以组成各种不同的颜色。产品能够展现的灰度数量越多，也就意味着这款产品的色彩表现力更加丰富，能够实现更强的色彩层次。例如三原色16级灰度，能显示的颜色就是16×16×16=4096色。不过目前的产品256级灰度已经非常地普遍了。

　所谓颜色或灰度级指黑白显示器中显示像素点的亮暗差别，在彩色显示器中表现为颜色的不同，灰度级越多，图像层次越清楚逼真。灰度级取决于每个像素对应的刷新存储单元的位数和显示器本身的性能。如每个象素的颜色用16位二进制数表示，我们就叫它16位图，它可以表达2的16次方即65536种颜色。如每一个象素采用24位二进制数表示，我们就叫它24位图，它可以表达2的24次方即16777216种颜色。

      灰度就是没有色彩，RGB色彩分量全部相等。如果是一个二值灰度图象，它的象素值只能为0或1，我们说它的灰度级为2。用个例子来说明吧：一个256级灰度的图象，RGB(100,100,100)就代表灰度为100,RGB(50,50,50)代表灰度为50。

      灰度是指黑白图像中点的颜色深度，范围一般从0到255，白色为255 ，黑色为0，故黑白图片也称灰度图像，在医学、图像识别领域有很广泛的用途

   　彩色图象的灰度其实在转化为黑白图像后的像素值（是一种广义的提法），转化的方法看应用的领域而定，一般按加权的方法转换，R ， G ，B 的比一般为3：6：1。

       任何颜色都有红、绿、蓝三原色组成，假如原来某点的颜色为RGB(R，G，B)，那么，我们可以通过下面几种方法，将其转换为灰度：

    1.浮点算法：Gray=R*0.3+G*0.59+B*0.11

    2.整数方法：Gray=(R*30+G*59+B*11)/100

    3.移位方法：Gray =(R*28+G*151+B*77)>>8;

    4.平均值法：Gray=（R+G+B）/3;

     5.仅取绿色：Gray=G；

   通过上述任一种方法求得Gray后，将原来的RGB(R,G,B)中的R,G,B统一用Gray替换，形成新的颜色RGB(Gray,Gray,Gray)，用它替换原来的RGB(R,G,B)就是灰度图了。

灰度值与像素值的关系
      记录了明亮（暗与亮，黑与白）度信息的2D图像成为亮度(灰度）图像。如果对于一张本身就是灰度图像（8位灰度图像）来说，他的像素值就是它的灰度值，如果是一张彩色图像，则它的灰度值需要经过函数映射来得到。灰度图像是由纯黑和纯白来过渡得到的，在黑色中加入白色就得到灰色，纯黑和纯白按不同的比例来混合就得到不同的灰度值。

    像素中具有RGB三通道的彩色图像来说怎样才能得到灰度图像呢？只要使R=G=B，三者的值相等就可以得到灰度图像。R=G=B=255为白色，R=G=B=0为黑色，R=G=B=小于255的某个整数时，此时就为某个灰度值。

像素值（灰度值）转换为CT值
    CT值的单位是Hounsfield，简称为Hu，范围是-1024-3071。用于衡量人体组织对X射线的吸收率，设定水的吸收率为0Hu。

　在DICOM图像读取的过程中，我们会发现图像的像素值有可能不是这个范围，通常是0-4096，这是我们常见到的像素值或者灰度值，这就需要我们在图像像素值（灰度值）转换为CT值。

　　首先，需要读取两个DICOM Tag信息，（0028|1052）：rescale intercept和（0028|1053）：rescale slope.

　　然后通过公式：

　　　　Hu = pixel * slope + intercept

　　计算得到CT值。

获得三维值：
 

PathDicom = "./PATIENT_DICOM5/"  # 与python文件同一个目录下的文件夹
lstFilesDCM = []
slices = []
number = 0
name = 0
for dirName, subdirList, fileList in os.walk(PathDicom):
    for filename in fileList:
            lstFilesDCM.append(os.path.join(dirName, filename))
            number += 1
RefDs = pydicom.read_file(lstFilesDCM[0])
ConstPixelSpacing = (float(RefDs.PixelSpacing[0]), float(RefDs.PixelSpacing[1]), float(RefDs.SliceThickness))
print(ConstPixelSpacing)