[發明專利]一種調制傳遞函數的計算方法、裝置及電子設備

有效

• 申請號: 202111489689.5
• 申請日: 2021-12-08
• 公開（公告）號: CN114202521B
• 公開（公告）日: 2023-03-24
• 發明（設計）人: 贠曉帆;魏國亮;丁佳;呂晨翀
• 申請（專利權）人: 北京醫準智能科技有限公司
• 主分類號: G06T7/00
• 分類號: G06T7/00;G06T5/30;G06F17/14
• 代理公司: 北京樂知新創知識產權代理事務所(普通合伙) 11734
• 代理人: 秦然

1.一種調制傳遞函數的計算方法，其特征在于，包括：

獲取至少一張主動脈CT圖像；

確定所述主動脈CT圖像中主動脈的中心；

計算以所述中心為起點的射線對應的CT值，所述射線的終點在與所述主動脈相鄰的組織內；

根據所述CT值生成徑向基函數f(r)，f(r)＝C，C為該CT值，r為主動脈的半徑；

根據所述徑向基函數、所述射線在線性成像系統中的理論成像計算調制傳遞函數；可以用階梯函數f(x，y)來表示以中心為起點的射線在線性成像系統中的理論成像，即當(x，y)∈S時，f(x，y)＝1，否則等于0，S的數學表達式為所述根據所述徑向基函數、所述射線在線性成像系統中的理論成像計算調制傳遞函數，包括：對所述徑向基函數進行傅里葉變換；對所述射線在線性成像系統中的理論成像進行傅里葉變換；確定所述徑向基函數的傅里葉變換值與所述理論成像的傅里葉變換值的比值，為調制傳遞函數。

2.根據權利要求1所述的調制傳遞函數的計算方法，其特征在于，所述確定所述主動脈CT圖像中主動脈的中心，包括：

確定所述主動脈CT圖像中的主動脈；

對所述主動脈進行形態學腐蝕處理，得到所述主動脈的中心。

3.根據權利要求2所述的調制傳遞函數的計算方法，其特征在于，所述確定所述主動脈CT圖像中主動脈，包括：

對所述主動脈CT圖像進行圖像分割，得到所述主動脈。

4.根據權利要求1所述的調制傳遞函數的計算方法，其特征在于，所述計算以所述中心為起點的射線對應的CT值，包括：

以水平方向的正方向為采樣起點，按照預設角度間隔對以所述中心為起點的射線進行均勻采樣，得到多條射線；

對于每一條射線：按照預設距離間隔對該射線進行均勻采樣，得到多個采樣點；確定每個采樣點對應的CT值；將各個采樣點對應的CT值進行加和，得到該射線的CT值；

將各條射線的CT值進行平均計算，得到以所述中心為起點的射線對應的CT值。

5.根據權利要求4所述的調制傳遞函數的計算方法，其特征在于，所述確定每個采樣點對應的CT值，包括：

對于每一個采樣點：

查找到與每個采樣點相鄰的多個像素點；

根據多個所述像素點的CT值確定所述采樣點的CT值。

6.一種調制傳遞函數的計算裝置，其特征在于，包括：

獲取單元，用于獲取至少一張主動脈CT圖像；

確定單元，用于確定所述主動脈CT圖像中主動脈的中心；

第一計算單元，用于計算以所述中心為起點的射線對應的CT值，所述射線的終點在與所述主動脈相鄰的組織內；

生成單元，用于根據所述CT值生成徑向基函數f(r)，f(r)＝C，C為該CT值，r為主動脈的半徑；

第二計算單元，用于根據所述徑向基函數、所述射線在線性成像系統中的理論成像計算調制傳遞函數；可以用階梯函數f(x，y)來表示以中心為起點的射線在線性成像系統中的理論成像，即當(x，y)∈S時，f(x，y)＝1，否則等于0，S的數學表達式為所述根據所述徑向基函數、所述射線在線性成像系統中的理論成像計算調制傳遞函數，包括：對所述徑向基函數進行傅里葉變換；對所述射線在線性成像系統中的理論成像進行傅里葉變換；確定所述徑向基函數的傅里葉變換值與所述理論成像的傅里葉變換值的比值，為調制傳遞函數。

7.根據權利要求6所述的調制傳遞函數的計算裝置，其特征在于，所述確定單元包括：

確定子單元，用于確定所述主動脈CT圖像中的主動脈；

處理子單元，用于對所述主動脈進行形態學腐蝕，得到所述主動脈的中心。

8.一種電子設備，其特征在于，包括：

至少一個處理器；以及與所述至少一個處理器通信連接的存儲器；其中，所述存儲器存儲有可被所述至少一個處理器執行的指令，所述指令被所述至少一個處理器執行，以使所述至少一個處理器執行如權利要求1-5任意一項所述的調制傳遞函數的計算方法。

9.一種計算機可讀存儲介質，其特征在于，所述計算機可讀存儲介質存儲有計算機指令，所述計算機指令用于使計算機執行如權利要求1-5任意一項所述的調制傳遞函數的計算方法。