Re: PET中,蒙特卡罗怎么产生那些随机数啊

bigben
bigben 2004-12-07 字数 667

所谓projector就是根据图像值计算投影值的一个过程(或者说算子,等等)

里面最重要的因素就是图像上的某一个像素(体素)对投影中的某一个像素的贡献值

怎么确定的问题。area-weighted按我的理解就是投影上的某像素,其投影线覆盖

图像中某像素的面积是多少,就按这个面积来计算贡献值。我一会贴一个图,可能

便于你理解。

在Monte Carlo模拟中,改变总共跟踪的photon history的数目,可以控制投影中

获得的总计数。这两个基本上是成正比的(考虑到统计涨落),比例因子就是系统的探测效率。

回忆一下概率里面方差是怎么计算的?

D(x)=E(x-E(x))^2

现在就等于你有8000个图像的样本,来求其方差,呵呵。

我的理解,蒙卡是用来产生多组投影数据样本的,你对这多组数据用MLEM迭代算法得到

多组重建图像。所以两个不是一回事,呵呵。

【 在 yjhhgp (yjhhgp) 的大作中提到: 】

: 标  题: Re: PET中,蒙特卡罗怎么产生那些随机数啊

: 发信站: BBS 水木清华站 (Tue Dec  7 11:19:19 2004), 站内

: 谢谢

: 还有就是在PET里面,phantom是怎样产生的?

: 以及如何用蒙卡计算图像的方差啊?

: 我看的一篇文献里面一段是这样:

: 为了衡量ML-EM算法的噪音特性(如方差,协方差),需要和蒙卡计算的噪音特性进行比较

: 无噪音投影数据是由simple area-weighted projector(不知道这个是什么东东)产生,噪音数据是加入泊松噪音。

: we simulated several projection count levels ranging from 20000 to 1500000 total projection counts.这个counts怎么控制啊?

: 有没有哪位用过这种projector?

: 蒙卡计算噪音特性的样本大小为8000副图像,迭代次数为50,总得投影数为500000

: 这句话的意思是不是:蒙卡计算方差的方法是不是对这8000副图像的方差进行平均啊?

: 还有在蒙卡里也用的上迭代吗?

: 说得没有条理,不好意思

: 希望哪位大哥能帮帮!

: 【 在 xct (XCT) 的大作中提到: 】

: : 语言里面有radom函数啊

: : 不过这算是伪随机

: : 所以现在有的cpu有了硬件随机数产生,对于蒙卡的准确度有意义

: --

MedImaging 医学影像
2 个回复
bigben
bigben 2004-12-07

版权属于:http://www.eecs.umich.edu/~fessler

【 在 bigben (zeze) 的大作中提到: 】

: 所谓projector就是根据图像值计算投影值的一个过程(或者说算子,等等)

: 里面最重要的因素就是图像上的某一个像素(体素)对投影中的某一个像素的贡献值

: 怎么确定的问题。area-weighted按我的理解就是投影上的某像素,其投影线覆盖

: ...................

bigben
bigben 2004-12-07

原来是这篇文献啊。Dr. Tsui和Dr. Barrett都是大牛啊,呵呵。Dr. Tsui

来我们组做报告的时候讲过一点和这个有关的内容。

关于Hoffman phantom,我前面贴过一个,嗯,可能是用xia0maer的ID贴的。

projector不是软件,你可以理解为为一种由图像计算投影的方法。我已经解释过了。如果还不清楚,可以看文献当中列的两个reference: (Huesman er a1 1977)和(Frey and Tsui

1991)。具体的计算过程需要你自己去实现。

图像中各个像素的值就是你的phantom啊,呵呵。

在这篇文献里,作者是把原始的phantom经过projector运算首先得到一个noise free

projection,然后用Monte Carlo方法作多次Possion抽样得到那8000组noisy projection。

然后用MLEM方法重建(好像是迭代50次吧),然后利用重建的结果来做统计分析。

如果你对Monte Carlo抽样的方法不是很清楚,可以先看看有关的书籍,呵呵。

【 在 yjhhgp (yjhhgp) 的大作中提到: 】

: 标  题: Re: PET中,蒙特卡罗怎么产生那些随机数啊

: 发信站: BBS 水木清华站 (Tue Dec  7 15:20:37 2004), 站内

: 还有就是在PET里面,phantom是怎样产生的?

: : 以及如何用蒙卡计算图像的方差啊?

: :

: : 我看的一篇文献里面一段是这样:

: : 为了衡量ML-EM算法的噪音特性(如方差,协方差),需要和蒙卡计算的噪音特性进行比较

: : 无噪音投影数据是由simple area-weighted projector(不知道这个是什么东东)产生,噪音数据是加入泊松噪音。

: 所谓projector就是根据图像值计算投影值的一个过程(或者说算子,等等)

: 里面最重要的因素就是图像上的某一个像素(体素)对投影中的某一个像素的贡献值

: 怎么确定的问题。area-weighted按我的理解就是投影上的某像素,其投影线覆盖

: 图像中某像素的面积是多少,就按这个面积来计算贡献值。我一会贴一个图,可能

: 便于你理解。

: *******这个projector是软件还是什么?如是软件的话什么地方可以找到。

: 另外,之前是不知道各个象素的值因此也没法按照你所说的得到投影啊.

: simple area-weighted projector这个好像可以产生均匀的Phantom(这个什么意思?:))

: : we simulated several projection count levels ranging from 20000 to 1500000 total projection counts.这个counts怎么控制啊?

: 在Monte Carlo模拟中,改变总共跟踪的photon history的数目,可以控制投影中

: 获得的总计数。这两个基本上是成正比的(考虑到统计涨落),比例因子就是系统的探测效率。

: ***********要作这种模拟需要什么软件或者其他。我以前没有用过蒙卡。

: : 有没有哪位用过这种projector?

: :

: : 蒙卡计算噪音特性的样本大小为8000副图像,迭代次数为50,总得投影数为500000

: : 这句话的意思是不是:蒙卡计算方差的方法是不是对这8000副图像的方差进行平均啊?

: : 还有在蒙卡里也用的上迭代吗?

: 回忆一下概率里面方差是怎么计算的?

: D(x)=E(x-E(x))^2

: 现在就等于你有8000个图像的样本,来求其方差,呵呵。

: 我的理解,蒙卡是用来产生多组投影数据样本的,你对这多组数据用MLEM迭代算法得到

: 多组重建图像。所以两个不是一回事,呵呵。

: :

: **********文献的意思好像这8000个图像不是用MLEM法重建的

: 如您有空就麻烦您看下这篇文章,我们图书馆没有这个的电子版

: 谢谢

: Wilson D W,Tsui B M and Barret  1994 Noise properties of EM algorithm :Monte Carlo simulations   Phys.Med.Biol. 39  847-72.

: --