[發(fā)明專利]基于OCTA視網(wǎng)膜圖像的黃斑中心無(wú)血管區(qū)體積計(jì)算方法及系統(tǒng)

有效

• 申請(qǐng)?zhí)枺?/td>: 202010094886.6
• 申請(qǐng)日: 2020-02-17
• 公開(kāi)（公告）號(hào): CN111292311B
• 公開(kāi)（公告）日: 2022-08-09
• 發(fā)明（設(shè)計(jì)）人: 陳強(qiáng);徐秋卓
• 申請(qǐng)（專利權(quán)）人: 南京理工大學(xué)
• 主分類號(hào): G06T7/00
• 分類號(hào): G06T7/00;G06T7/62;A61B3/10;A61B3/14
• 代理公司: 南京理工大學(xué)專利中心 32203
• 代理人: 馬魯晉

1.一種基于OCTA視網(wǎng)膜圖像的黃斑中心無(wú)血管區(qū)體積計(jì)算方法，其特征在于，所述方法包括以下步驟：

步驟1，采集OCTA視網(wǎng)膜圖像；

步驟2，設(shè)置包含黃斑中心無(wú)血管區(qū)域即FAZ區(qū)域的限制區(qū)域；具體過(guò)程包括：

步驟2-1，分割OCTA視網(wǎng)膜圖像的內(nèi)界膜ILM和外叢狀層OPL下邊界；

步驟2-2，將ILM作為限制區(qū)域的上邊界，將OPL下邊界下方lμm處作為限制區(qū)域的下邊界；

步驟3，拉平所述限制區(qū)域；具體過(guò)程包括：

步驟3-1，獲取OPL下邊界三維全局最低點(diǎn)，并將該最低點(diǎn)下方lμm處的行坐標(biāo)作為下邊界基準(zhǔn)；

步驟3-2，將限制區(qū)域內(nèi)的每一列移動(dòng)至下邊界基準(zhǔn)，實(shí)現(xiàn)每一列下邊界對(duì)齊，獲得初步拉平限制區(qū)域；

步驟3-3，獲取所述初步拉平限制區(qū)域的上邊界三維全局最高點(diǎn)，并將該最高點(diǎn)作為上邊界基準(zhǔn)；

步驟3-4，將所述初步拉平限制區(qū)域內(nèi)的每一列變換為固定長(zhǎng)度，該固定長(zhǎng)度為所述下邊界基準(zhǔn)、上邊界基準(zhǔn)之間的距離長(zhǎng)度；

步驟4，在拉平后的限制區(qū)域內(nèi)在軸向方向依次生成若干幀一個(gè)像素厚度的投影圖像，并在每一幀投影圖像上標(biāo)注FAZ區(qū)域，記為標(biāo)簽圖像；

步驟5，將所述標(biāo)簽圖像還原至未拉平限制區(qū)域前的原始層位置；具體過(guò)程包括：

反向執(zhí)行步驟3的反操作，將標(biāo)簽圖像還原為符合OCTA視網(wǎng)膜圖像原始層結(jié)構(gòu)分布的標(biāo)簽圖像；

步驟6，統(tǒng)計(jì)還原后的標(biāo)簽圖像中標(biāo)簽像素的個(gè)數(shù)，由此求取FAZ區(qū)域的體積。

2.根據(jù)權(quán)利要求1所述的基于OCTA視網(wǎng)膜圖像的黃斑中心無(wú)血管區(qū)體積計(jì)算方法，其特征在于，所述l的取值為OCTA圖像采集設(shè)備劃定的深層血管復(fù)合體DVC下邊界值。

3.根據(jù)權(quán)利要求1所述的基于OCTA視網(wǎng)膜圖像的黃斑中心無(wú)血管區(qū)體積計(jì)算方法，其特征在于，步驟3-4中將所述初步拉平限制區(qū)域內(nèi)的每一列變換為固定長(zhǎng)度，具體采用線性插值法。

4.根據(jù)權(quán)利要求1所述的基于OCTA視網(wǎng)膜圖像的黃斑中心無(wú)血管區(qū)體積計(jì)算方法，其特征在于，步驟6中所述求取FAZ區(qū)域的體積，所用公式為：

V＝n×p

式中，V表示FAZ區(qū)域的體積，n表示標(biāo)簽像素的個(gè)數(shù)，p表示圖像的分辨率參數(shù)。

5.一種基于OCTA視網(wǎng)膜圖像的黃斑中心無(wú)血管區(qū)體積計(jì)算系統(tǒng)，其特征在于，所述系統(tǒng)包括：

采集模塊，用于采集OCTA視網(wǎng)膜圖像；

區(qū)域設(shè)置模塊，用于設(shè)置包含黃斑中心無(wú)血管區(qū)域即FAZ區(qū)域的限制區(qū)域；具體過(guò)程包括：

步驟2-1，分割OCTA視網(wǎng)膜圖像的內(nèi)界膜ILM和外叢狀層OPL下邊界；

步驟2-2，將ILM作為限制區(qū)域的上邊界，將OPL下邊界下方lμm處作為限制區(qū)域的下邊界；

區(qū)域變換模塊，用于拉平所述限制區(qū)域；具體過(guò)程包括：

步驟3-1，獲取OPL下邊界三維全局最低點(diǎn)，并將該最低點(diǎn)下方lμm處的行坐標(biāo)作為下邊界基準(zhǔn)；

步驟3-2，將限制區(qū)域內(nèi)的每一列移動(dòng)至下邊界基準(zhǔn)，實(shí)現(xiàn)每一列下邊界對(duì)齊，獲得初步拉平限制區(qū)域；

步驟3-3，獲取所述初步拉平限制區(qū)域的上邊界三維全局最高點(diǎn)，并將該最高點(diǎn)作為上邊界基準(zhǔn)；

步驟3-4，將所述初步拉平限制區(qū)域內(nèi)的每一列變換為固定長(zhǎng)度，該固定長(zhǎng)度為所述下邊界基準(zhǔn)、上邊界基準(zhǔn)之間的距離長(zhǎng)度；

標(biāo)簽圖像生成模塊，用于在拉平后的限制區(qū)域內(nèi)在軸向方向依次生成若干幀一個(gè)像素厚度的投影圖像，并在每一幀投影圖像上標(biāo)注FAZ區(qū)域，記為標(biāo)簽圖像；

圖像還原模塊，用于將所述標(biāo)簽圖像還原至未拉平限制區(qū)域前的原始層位置；具體過(guò)程包括：

反向執(zhí)行區(qū)域變換模塊的反操作，將標(biāo)簽圖像還原為符合OCTA視網(wǎng)膜圖像原始層結(jié)構(gòu)分布的標(biāo)簽圖像；

體積計(jì)算模塊，用于統(tǒng)計(jì)還原后的標(biāo)簽圖像中標(biāo)簽像素的個(gè)數(shù)，由此求取FAZ區(qū)域的體積。