matlab——开窗显示法读取dicom图像
直接读取dicom图像,对比度很低,利用开窗显示法,读取的图片更为清晰,更利于医学分析
ct值和灰度值之间的变换公式:ct= pixel * RescaleSlope +RescaleIntercept
- 将ct值转换为灰度值的开窗显示
info=dicominfo('img');%读取dicom图像的元数据
wc=info.WindowCenter;%读取窗位
ww=info.WindowWidth;%读取窗宽
rs=info.RescaleSlope;
ri=info.RescaleIntercept;
lc=wc-ww/2;
hc=wc+ww/2;
lc_pixel=(lc-ri)/rs;%将CT值转换为灰度值
hc_pixel=(hc-ri)/rs;
subplot(1,2,1),imshow(I),title('原始图像');
subplot(1,2,2),imshow(I,[lc_pixel,hc_pixel]),title('开窗显示');
结果显示
- 将灰度值转换为CT值的开窗显示
I=dicomread('img');
I_ct=