第(1/3)頁(yè)

    林葉很快就進(jìn)入了狀態(tài)。

    吃飽了飯，又喝了一杯咖啡提神，此刻精神十分亢奮。

    “離散模型中使用代數(shù)迭代法（art）對(duì)網(wǎng)格化的離散數(shù)據(jù)直接建模，用矩陣迭代的方法對(duì)像素網(wǎng)格進(jìn)行處理，實(shí)現(xiàn)圖像的重構(gòu)。”

    這是林葉的思路。

    至于數(shù)學(xué)模型具體如何，林葉還在思考。

    這個(gè)離散模型，最重要的還是迭代問(wèn)題。

    fbp（濾波反投影）算法具有的缺陷與限制，只能應(yīng)用代數(shù)迭代才能完美解決。

    “這里得先引入投影矩陣的概念...。”

    林葉嘀咕說(shuō)道。

    時(shí)間一點(diǎn)一點(diǎn)的過(guò)，林葉空白的a4紙上現(xiàn)在已經(jīng)寫(xiě)滿(mǎn)了數(shù)學(xué)公式。

    r_ij={1，i號(hào)射線通過(guò)j號(hào)像素中心//0，其他。

    記x=[x1，x2，...，xj]^t為j維圖像矢量，x_j表示圖像上第j個(gè)像素點(diǎn)的吸收強(qiáng)度，

    p=[p1，p2...pi]^t為i維投影矢量，p_i表示第i條射線所經(jīng)過(guò)的所有像素的總射線投影;

    rx=p;

    ...

    林葉看著草稿紙上最為核心的迭代模型，

    x^(k+1)=x^(k)+λ^(k)(q_ikr_ikx^(k)/||r_ik||^2)

    林葉嘀咕道：

    “這里還差一點(diǎn)東西，那就是降噪算法。不論何種重建算法，不可避免會(huì)帶入偽影點(diǎn)。

    組委會(huì)給的數(shù)據(jù)含有噪聲，成像質(zhì)量可能并不會(huì)達(dá)到最好的要求。這個(gè)其余學(xué)校的建模選手也應(yīng)該會(huì)認(rèn)識(shí)到這個(gè)問(wèn)題才對(duì)。”

    “這里必須得用算法進(jìn)行降噪。我得查一查這方面的文獻(xiàn)。”

    林葉用筆壓著筆記本電腦前的草稿紙，然后查閱起了文獻(xiàn)。

    “這種問(wèn)題應(yīng)該在國(guó)外醫(yī)學(xué)sci期刊上比較好尋找合適的算法。”

    看了一個(gè)多小時(shí)，林葉總算找到了一種合適的算法。

    “這種算法研究一下能不能改進(jìn)。”

    林葉現(xiàn)在數(shù)學(xué)水平已經(jīng)很高了，但是想要在一些頂級(jí)sci期刊上動(dòng)刀別人的算法模型，

    可能性不大。

    “看來(lái)這里只能引用了，這里用這種算法進(jìn)行降噪處理已經(jīng)算是完美答卷了，強(qiáng)行改只會(huì)畫(huà)蛇添足。”

    林葉有些喪氣，

    “早知道這段時(shí)間不忙著做第六個(gè)任務(wù)了。”
    第(1/3)頁(yè)