刚刚写出来的c# 的医学图像系统。PS：目前只支持隐式（大部分的CT，DR，CR）已经补上DICM图片

hzr0071 · 发表于 2011-4-23 12:26:17

不知道论坛里有没有医学出身的。不知道放哪里合适。
优化了算法。调窗宽窗位还是挺快的。
暂时还没有其他功能

(原文件名:a.JPG)
dimcreader的exe文件
点击此处下载 ourdev_632940DTDCII.rar(文件大小:7K) (原文件名:dicmreaderexe.rar)
c#源码 PS：数组转图片部分是copy的，懒得写了
点击此处下载 ourdev_632941YA4EPY.rar(文件大小:45K) (原文件名:dicmreader.rar)
点击此处下载 ourdev_632942T5HIDM.pdf(文件大小:794K) (原文件名:DICOM格式白皮书06(dicm数据字典).pdf)
点击此处下载 ourdev_632943NJ42P4.pdf(文件大小:2.23M) (原文件名:VB实现DICOM到BMP的转换.pdf)
点击此处下载 ourdev_632986TWTX7W.rar(文件大小:8.40M) (原文件名:DICM图像.rar) 这是一张飞利浦DR的胸片

有关DICOM的介绍在8楼。

durgy · 发表于 2011-4-23 13:09:22

good

xingliu · 发表于 2011-4-23 13:42:13

mark

hzr0071 · 发表于 2011-4-23 13:46:25

终于有人回帖了。。。。。虽然只是mark。。。。。。还是很欣慰的。，。。

semonpic · 发表于 2011-4-23 13:50:41

我是医工出身，呵呵，你这个不错，去研究研究什么重建算法。应该不错

huzcmail · 发表于 2011-4-23 14:30:25

好帅，希望楼主来点科普知识

linghu2 · 发表于 2011-4-23 15:02:05

LZ,准备开始讲课吧,我等学习

skynet · 发表于 2011-4-23 15:36:18

看不懂,凑热闹

hzr0071 · 发表于 2011-4-23 17:04:37

DICOM（Digitalimaging and Communications in Medicine）医学数字图像与通讯。
分为医学数字图像和信息通讯两个方面。PS：号称dicom通讯卡单买10w一个，不知道有什么好的。普通tcp/ip卡才十几块一个。
x线医学数字图像为高灰阶图像。

(原文件名:x线管.jpg)
上图为x线管。
结构和电子管差不多。电子打在阳极上产生x线。
x线具有穿透效应。经过人体时衰减。由于人体的密度不同，对x线的吸收也不同。最终由探测器测出不同的x线强度形成影像。
ct值：Ct值的含义是： CT值代表X线穿过组织被吸收后的衰减值.
规定 CT值=1000×（u—u水）/ u水
水的ct值是0，上限1000（骨头），下限-1000（空气）。
这并不代表灰度值只有2000，实际图像的灰度值在100000以上（取决于探测器和显示器）。

因为普通显示器（在windows系统下）只有256个灰阶（医用的显示器现在有1000以上的），
图像灰度值远远大于显示器所能显示的范围。所以要进行开窗显示。

窗宽（Window Width）：表示所显示信号强度值的范围。
窗位（Window Level）又称窗水平。是图像显示过程中代表图像灰阶的中心位置。

(原文件名:开窗显示.JPG)
上图为开窗显示的办法。

就写到这里吧。
PS：到医院还是要数字的x线片好。
PPs：拿到片子还是让放射科医生看。。。。绝对比临床医生看的好。

yanrz · 发表于 2011-4-23 17:12:53

厉害，学习了一点，呵呵，楼主继续上课！！！！

hzr0071 · 发表于 2011-4-23 17:20:36

回复【9楼】yanrz yanrz
-----------------------------------------------------------------------

晚上再写了。先去吃饭了。

hasea · 发表于 2011-4-23 17:28:29

这个话题不多见啊，速速学习！

hzr0071 · 发表于 2011-4-23 17:37:18

开窗显示对诊断至关重要。
比如

(原文件名:骨头.jpg)
在高窗位，小窗宽的时候可以清晰地显示骨头的结构，甚至可以看清骨小梁。

(原文件名:软组织.jpg)
在低窗位，大窗宽的时候可以清晰地显示软组织的影像。比如有化脓性骨髓炎的时候可以看到有网状软组织。
我这没有ct的DICM片子，全是转成的bmp的，所以没法给大家演示，原理是一样的。过几天弄几张。
PS：这个人的手腕好像是骨折了。考赖氏骨折。

skynet · 发表于 2011-4-23 17:42:29

楼主专业人士阿,启动详解

zxq6 · 发表于 2011-4-23 18:05:27

楼主，顶你，以前做过放射科工作站，感觉你技术很NB的说。

hzr0071 · 发表于 2011-4-23 18:07:53

回复【14楼】zxq6 小马哈
-----------------------------------------------------------------------

我是学的影像诊断，做医生的。实在的说。。

fanwt · 发表于 2011-4-23 18:30:34

感觉这个实用性最强

wcm_e · 发表于 2011-4-23 18:34:41

flyforyou85 · 发表于 2011-4-23 18:46:44

像灰度值远远大于显示器所能显示的范围。所以要进行开窗显示。
==============================================================

不错，期待LZ详细的往下讲，上面的那句不是很清楚，说详细点吧！

srygg · 发表于 2011-4-23 18:58:15

楼主的意思是不是有了这个软件能看清更多的细节？

118139 · 发表于 2011-4-23 19:00:15

有点灰度增强的效果。。

118139 · 发表于 2011-4-23 19:04:12

我看过
脑CT拍出来的片子

一根根血管都是立体的真是漂亮。
只显示血管，其它没有，是不是经过处理。。

hzr0071 · 发表于 2011-4-23 19:05:15

回复【18楼】flyforyou85
-----------------------------------------------------------------------

医学图像的灰阶很多（10000以上），而普通显示器的灰阶只有256（每个像素由RGB构成，32bit的还有alpha值。每个颜色由8bit构成，最大0xff，当RGB相同时显示灰阶，所以只有256个灰阶）。
开窗显示：比如窗位为30000，窗宽为20000，那么显示的灰度范围为(30000-10000)至(30000+10000),不足20000的位0，超过40000为255，
中间的20000个灰阶分为256份在显示器上显示。

(原文件名:开窗公式.jpg)

由于人的眼睛不能分辨灰度相近的两个灰阶，所以可以减小窗宽以增加对比度。
如果要显示不同密度的组织，就可以降低窗宽，以显示以窗位为中心的大范围组织。

hzr0071 · 发表于 2011-4-23 19:07:46

回复【21楼】118139
-----------------------------------------------------------------------

那个叫血管重建。很高的科技含量，包括边缘检查，最大密度投影，和容积表面成像。
有人专门做那个，好像一次100块RMB，不过不同的人使用同一软件做出的效果也不一样。影响诊断结果。

hzr0071 · 发表于 2011-4-23 19:20:53

既然有人关注，那就多说一点。
以前照x片的时候，医生拿一个小板子放在身体下面，其实那是增感屏，里面有胶片。

(原文件名:增感屏.JPG)

增感屏里面的物质把x线转变为荧光，使胶片感光。然后胶片再拿去洗（洗下来的全是银哦）。
增感物质最好的是稀土的。（稀土真是TMD好东西，一会再说其他用途）。

(原文件名:非晶硒探测器.JPG)

这是现在的数字x线机（DR）的探测器。和现在的相机的传感器很像。这是非晶硒的，不过也有CCD的。
非晶硒的探测器大概100wRMB一张（30cmx40cm左右）。
对温度敏感，必须在25℃左右使用和存放。夏天不开空调一会就会坏掉。

zx5111030 · 发表于 2011-4-23 19:27:18

楼主去万东吧。

hzr0071 · 发表于 2011-4-23 19:28:19

回复【25楼】zx5111030
-----------------------------------------------------------------------

想去来着。。。。不过后来打听一下好像待遇不好。。。
而且投简历去，人家不理我
我还是去当医生吧

STM32_Study · 发表于 2011-4-23 19:40:16

有意思，很感兴趣啊～～～

楼主继续～～～

hzr0071 · 发表于 2011-4-23 20:28:52

在这个论坛里电子技术都是高人，电子方面我就不班门弄斧了，那就说一些图像重建的方法。
DSA（数字减影血管造影）：

(原文件名:脑DSA.jpg)

这是一张椎-基动脉系的DSA造影，包括了椎动脉，基底动脉，大脑后动脉，小脑下前动脉和小脑下后动脉。
DSA以对血管清晰地影像呗医生们称为诊断的金标准。
DSA有可以应用在介入上，以帮助医生以微创（甚至无创）进行心脏的冠脉放支架，心脏心房、心室壁缺损的修补等手术。

(原文件名:DSA机器.JPG)

上次有一朋友去给医院修DSA，人家问他减影是怎么减出来的，他是搞维修的，不清楚这些，使得人家不打信任他。

那DSA的成像机制是什么呢？

(原文件名:DSA成像原理.JPG)

对比剂：打进人体的，密度很大（胶体）或者很小（气体）的物质，以增加组织与周围的对比度。
在打对比剂之前先照一张片子。然后打入对比剂，然后再照一张，用后一张逐个像素的减去前一图片同位置像素的灰度值。
这样，在打入对比剂的地方，灰度值不一样（一般用密度很高的对比剂），相减出来的图片就是DSA图片。

哈哈，很简单吧。不过很有用哦。。。。。

jackielau · 发表于 2011-4-23 20:42:03

厉害，继续开讲！！搬个小板凳

hzr0071 · 发表于 2011-4-23 20:58:31

(原文件名:CT机器.jpg)

(原文件名:ct.jpg)

(原文件名:ct探测器.JPG)

有些人很纳闷（我让也很纳闷），为什么CT转一圈就能照出人体的断层图像呢。
其实原理很简单。让我们先看一下CT图像重建的反投影法。
接：32楼

wukong · 发表于 2011-4-23 20:58:33

发一些人体造影的特写看看

hzr0071 · 发表于 2011-4-23 21:06:00

(原文件名:反投影法示意图.JPG)
反投影法：先在各个角度对射，获得一维的数据，然后叠加，减去基数，就获得了断层图像。
简单吧。哈哈哈、、、、、、

(原文件名:反投影法.JPG)
方法简单，但是问题也来了，叠加然后减去基数后，会发现图像周围有毛刺，
这就是为什么老的CT扫描大脑颅后窝的时候会有很多星形
伪影。
为了解决这一问题，人们开始想办法。办法之一就是转化为频率空间后滤波，然后再转换回来。

(原文件名:反投影法矫正.JPG)

hzr0071 · 发表于 2011-4-23 21:07:53

回复【31楼】wukong
-----------------------------------------------------------------------

我找好看。心脏和脑的比较多

jamesyu · 发表于 2011-4-23 21:08:55

Cool

118139 · 发表于 2011-4-23 21:12:05

刚才说错了，是核磁共振拍出来的图像很立体，摄像机就像在脑里拍的一样
CT的图像太平面很烂。。不在一个档次。。

有没有核磁共振拍的图，看一次都要赞叹一次。。经典图像

hzr0071 · 发表于 2011-4-23 21:37:42

回复【35楼】118139
-----------------------------------------------------------------------

磁共振的图像也是平面的。因为磁共振的特殊性，对大脑、腰椎和膝盖的成像有不可替代的地位。
不过普通磁共振不能扫心脏、肺等组织，而且由于大部分医生不清楚磁共振的成像机制和磁共振多序列成像的复杂性，
磁共振不是所有医生都能看的懂得，能看出病来也是考经验。
CT或者磁共振只是诊断工具，有不同的用处。不能够比较谁好谁差。

立体图像比如

(原文件名:心脏冠脉重建.jpg)

这是64排CT扫描出的图像，然后三维重建（容积再现VR）出来的图像

oldtom · 发表于 2011-4-23 21:57:18

这个坛里面，N人N多啊!

这里有没有搞卫星的？

afei8856 · 发表于 2011-4-23 22:04:51

看着真清楚啊呵呵 mark

pjdu · 发表于 2011-4-23 22:12:34

mark

flor · 发表于 2011-4-23 22:17:59

不懂

eiglxl · 发表于 2011-4-23 22:32:41

fy024 · 发表于 2011-4-23 22:46:50

高深学科~~~~

lcgforward · 发表于 2011-4-23 22:52:00

长见识了,谢谢!

STM32_Study · 发表于 2011-4-24 01:02:24

64排CT在国内还不算多，要比较好的医院才有呢

gzhuli · 发表于 2011-4-24 01:31:54

回复【37楼】oldtom
这个坛里面，n人n多啊!
这里有没有搞卫星的？
-----------------------------------------------------------------------

说不准还有打卫星的呢……

eworker · 发表于 2011-4-24 07:30:33

专业

1ongquan · 发表于 2011-4-24 08:34:27

还是想不通CT

cody · 发表于 2011-4-24 08:53:25

强人,关注一下,虽然看不懂

wjhhh · 发表于 2011-4-24 10:08:28

胡诌一下， 12楼“手腕好像是骨折”是这里吗？

(原文件名:未命名.JPG)

hzr0071 · 发表于 2011-4-24 11:49:19

回复【49楼】wjhhh
-----------------------------------------------------------------------

就是这里，桡骨远端近关节面，好像尺骨远端也有点。这种桡骨远端近关节面2.5CM以内的骨折称为克莱氏骨折

回复【47楼】1ongquan
-----------------------------------------------------------------------
那一会我自己画个图，好让大部分人能看懂。
就是光照和影子的原理，原始的方法不加傅里叶变换，应该挺好理解的

ndust · 发表于 2011-4-24 12:09:40

newID · 发表于 2011-4-24 14:15:18

搬凳子听课。如果涉及到专业的术语或缩写，希望楼主解释一下。

hzr0071 · 发表于 2011-4-24 14:16:38

假设一物体横截面具有4种密度，分别为1,2,3，4，用x线在0°、45°、90°、135°进行投照，在探测其上就会获得该物体的“影子”--一维数组。如图：

(原文件名:x线投影0°.jpg)

(原文件名:x线投影45°.jpg)

(原文件名:x线投影90°.jpg)

(原文件名:x线投影135°.jpg)
然后利用各个一维数据，在对应的角度进行反投影，得到二维数据。如图：

(原文件名:0°反投影.jpg)

(原文件名:45°反投影.jpg)

(原文件名:90°反投影.jpg)

(原文件名:135°反投影.jpg)
然后把获得的各个角度的反投影进行叠加。
得到一个叠加的二维数据，然后减去基数10（0°或者90°时获得的一行或一列的和），然后除去最大公约数，就得到了
原始物体的断面数据。

(原文件名:反投影叠加.jpg)
反投影法矫正在楼上。
PS：CT不要是不是就做，做一次全身CT剂量相当于一次放疗剂量。会掉光头发的哦、、、、、

stdioh · 发表于 2011-4-24 16:41:32

太牛了，过来慢慢听课。

jackielau · 发表于 2011-4-24 17:15:49

明白了，原理很简单，以前很神秘！

继续小板凳

qdsywh · 发表于 2011-4-24 19:25:53

请楼主更新，继续听课

CadProE · 发表于 2011-4-24 22:59:48

原理讲的挺通俗易懂的

lbn2006 · 发表于 2011-4-24 23:32:49

mark

icetea · 发表于 2011-4-25 03:07:18

好贴留名
话说透视骨骼的效果核磁共振效果应该是最好吧，
我老爸04年做腰椎间盘突出手术室就是靠核磁共振最终明确手术部位的。
那是做CT的价格是180左右，核磁共振是400多，估计是现在应该过800了。
话说楼主那张MRI是彩色的，我老爸那会还是黑色的不过清晰度和CT比就像1000万像素点和300万像素点DC的差别一样。
看样子楼主应该是三甲以上的医院了，我记得一台这样的设备都是相当昂贵的，我这边04年那会就两家医院有这种设备。

kneken · 发表于 2011-4-25 07:55:53

mark

leoyang · 发表于 2011-4-25 08:15:43

LZ作为一个医生的话，能画出PCB来，确实比较NB了，呵呵，厉害
我搞过两年的图像处理，基本的图像处理原理还是看得懂，不过应用在医学方面的，这次算是相对集中的见识了，谢谢LZ科普

alengend · 发表于 2011-4-25 09:27:14

感谢楼主的无私奉献，谢谢！

guaipizi · 发表于 2011-4-25 09:48:36

楼主真牛逼

lee345 · 发表于 2011-4-25 10:02:49

知识，学习

abao_wjx · 发表于 2011-4-25 10:23:31

好厉害的牛人啊！！收下真想研究研究这个呢！！

bangbangji · 发表于 2011-4-25 10:32:15

这个相当要顶啊

fanwt · 发表于 2011-4-25 10:35:45

顶楼主，科普呀

wuxi_stl · 发表于 2011-4-25 10:40:57

学习，

skynet · 发表于 2011-4-25 10:54:34

楼主,有空去给这贴看看
http://www.ourdev.cn/bbs///bbs_content.jsp?bbs_sn=4694974&bbs_id=9999
貌似论坛就你1个医生噢

hzr0071 · 发表于 2011-4-25 10:56:18

可能好多人对CT图像和MRI(磁共振)图像有误会。
认为MRI图像是立体的或者是彩色的。
这是不正确的。下面有两张MRI的图像，大家可以看一下。MRI有不同的序列，脑脊液是低信号（黑色）的是T1加权像（A），脑脊液是高信号（白色）的是T2加权图像（B）。
MRI有很多的序列，可以形成不同的图像，为诊断疾病提供依据。

(原文件名:MIR图像.jpg)

(原文件名:MRI图像.png)

不论CT还是MRI，都可以进行3D重建，是对多个图片插值后的叠加，所有医学图片上的颜色都是伪彩色，是人为加上去的。
不要认为MRI比CT好，只是有些病需要MRI来成像。有些病用MRI反而会扫不清楚，造成误诊(有些医生向病人传播MRI比CT好的思想是不负责任的）

(原文件名:3D重建.jpg)

这一张图片必须用CT和MRI扫描，然后叠加才行，软组织MRI的清晰，不过MRI扫不出骨骼和钙化灶。
PS：临床好多医生不会看MRI图像，还得靠我们放射科的医生，我们下的报告才是诊断的金标准，哈哈哈。
PPS：咱们论坛如果有人能做出国产的MRI谱仪(大概用7、8个DSP，然后几个FPGA，然后是各种信号算法)来，那绝对是中国一大国宝。。。。。
PPPS：图像清楚还得靠咱们搞电子的把机器做好。

hzr0071 · 发表于 2011-4-25 11:02:09

回复【70楼】skynet 青青草原
-----------------------------------------------------------------------

我已经给他回了，好像是真菌感染。不过我不是皮肤科的。所以也不敢确定。最好还是去做病理。

albert_w · 发表于 2011-4-25 14:15:17

CT的原理知道,只是感觉真实使用下计算量好大啊!

Meteora · 发表于 2011-4-25 16:01:07

厉害

hzr0071 · 发表于 2011-4-25 16:09:12

回复【74楼】albert_w
-----------------------------------------------------------------------

真正的使用都是基于傅里叶空间变换的，用dsp进行计算，计算量还可以，成像很快。

ggyyll8683 · 发表于 2011-4-25 16:10:28

学习，本科毕业的时候做过图像配准。研究生阶段主要做信号了。我也是生物医学工程的

PaulDE · 发表于 2011-4-25 23:01:57

搬个小板凳来坐着看

zgx8 · 发表于 2011-4-25 23:39:14

象科普讲座，学习了！

zgx8 · 发表于 2011-4-25 23:41:01

我认为学了医学知识可以终身受益，现在看病难呀

au730 · 发表于 2011-4-26 00:36:54

mark

10000years · 发表于 2011-4-26 09:16:52

mark长见识

fsclub · 发表于 2011-4-26 09:18:27

又懂医学又懂电子，这个才是强人啊！

fsclub · 发表于 2011-4-26 09:18:58

又懂医学又懂电子，这个才是强人啊！

Gorgon_Meducer · 发表于 2011-4-26 09:21:37

搬板凳~

radiation · 发表于 2011-4-26 14:48:21

楼主很专业，顶起来。

yusufu · 发表于 2011-4-26 16:26:45

听着真爽~

peak · 发表于 2011-4-26 17:27:27

回复【楼主位】hzr0071
-----------------------------------------------------------------------

楼主，你好，请问你的邮箱多少？我想请教你问题。我的是hf-zhu@qq.com

shaozh · 发表于 2011-4-26 17:59:36

学习了！

hzr0071 · 发表于 2011-4-26 21:32:25

感谢大家这么支持。这是我的第一个回复这么多的帖子。
关于CT的这是最后一个帖子。本来想再讲一下MRI，不过MRI的东西太多，而且关系到磁物理，量子物理，还有旧的傅里叶变换，所以以后整理出资料来再说吧。
先说一下CT的应用在实际的比较快的方法（还多机子是反投影的，效果还是挺好的）。

(原文件名:100_0631.JPG)

(原文件名:100_0632.JPG)

中心切片定理：密度函数f（x，y）在某一方向上的投影函数g θ （R）的一维傅里叶变换函数 G θ （p），是密度函数f（x，y）的二维傅里叶变换F（p，θ）在p，θ平面上沿同一方向过原点直线上的值。

这样的话，我们就可以把各个方向的投影的傅里叶变换算出来，然后放在一起，然后再插值，然后进行二维傅里叶变换。
这些用DSP实现不是问题，速度超级快。不过也有很多缺点。好像是出现函数拖尾什么的，信号学的不好，想不起来了。

上次说医学图像的灰度值有100000以上，这个说法有点错误。这个根据扫描的机器而定。在dicm图片的0x0028,0101处储存了图像的灰度储存的位数，最大有15位的（0x0000000f），也就是最大灰度32768（philips的以一台DR），普通的CT的灰度在12位左右。
新的阅读器已经支持了几款机器的图片，支持图像自动大小，和图像最大灰度测定，不过还没加辅助功能，等完全做完再放出来吧。

(原文件名:QQ截图未命名.jpg)

gwh1128 · 发表于 2011-4-26 22:47:56

不错啊

ma_gaap_7 · 发表于 2011-4-26 23:25:09

水一下:

CT進行3D重建的精度有多高?

可以output什麼fromat?

ndt2000 · 发表于 2011-4-27 10:05:40

够历害

af00 · 发表于 2011-4-27 10:49:25

楼主太强了，顶你

hzr0071 · 发表于 2011-5-11 10:55:44

DICM阅读器最新版本到这里下载。
http://www.ourdev.cn/bbs/bbs_content.jsp?bbs_sn=4728354&bbs_page_no=1&search_mode=3&search_text=hzr0071&bbs_id=9999

linghu2 · 发表于 2011-5-11 11:28:45

回复【76楼】hzr0071

真正的使用都是基于傅里叶空间变换的，用dsp进行计算，计算量还可以，成像很快。
-----------------------------------------------------------------------

速度具体多少,DSP型号什么?

reflecter · 发表于 2011-5-11 11:31:13

回复【15楼】hzr0071
回复【14楼】zxq6 小马哈
-----------------------------------------------------------------------
我是学的影像诊断，做医生的。实在的说。。
-----------------------------------------------------------------------
会写程序的医生？

msq-tc · 发表于 2011-5-16 12:00:58

要学习下!

stone21com · 发表于 2011-5-25 16:45:25

标记下,这个对于我来说是扫盲,呵呵

armok · 发表于 2011-6-29 23:25:41

cool !

diyer_zhou · 发表于 2011-6-30 12:35:10

记号，有空看

刚刚写出来的c# 的医学图像系统。PS：目前只支持隐式（大部分的CT，DR，CR） 已经补上DICM图片

阿莫论坛20周年了！感谢大家的支持与爱护！！

刚刚写出来的c# 的医学图像系统。PS：目前只支持隐式（大部分的CT，DR，CR）已经补上DICM图片