X=dicomread('abdomen_002.dcm'); X=X-1024; X=double(X); WW=400; WL=40; for i=1:512 for j=1:512 if (X(i,j)>WL+WW/2) Y(i,j)=255; else if(X(i,j)<WL-WW/2) Y(i,j)=0; else if(X(i,j)>=WL-WW/2 & X(i,j)<=WL+WW/2) Y(i,j)=(X(i,j)+WW/2-WL)*255/WW; end end end end end Y=uint8 (Y); imshow(Y);
123
原创
©著作权归作者所有:来自51CTO博客作者东Frank的原创作品,请联系作者获取转载授权,否则将追究法律责任
下一篇:我的友情链接
提问和评论都可以,用心的回复会被更多人看到
评论
发布评论
相关文章
-
123
1 2 3 4 5 6 7 8 9
职场 休闲