[發明專利]圖像配準方法、裝置及設備和存儲介質

有效

• 申請號: 201810556548.2
• 申請日: 2018-05-31
• 公開（公告）號: CN108830886B
• 公開（公告）日: 2022-04-01
• 發明（設計）人: 于海晶;李微微
• 申請（專利權）人: 東軟醫療系統股份有限公司
• 主分類號: G06T7/32
• 分類號: G06T7/32;G06T7/33
• 代理公司: 北京博思佳知識產權代理有限公司 11415
• 代理人: 林祥

1.一種圖像配準方法，其特征在于，包括：

獲取待配圖像及預定的標準圖像；所述標準圖像為超聲圖像庫中反映被檢體的預定標準切面的圖像；所述待配圖像為超聲設備生成的實時圖像；

分別對所述待配圖像與所述標準圖像進行小波變換，對應得到第一近似分量圖像和第二近似分量圖像；

對所述第一近似分量圖像進行至少一次幾何變換，并計算每次變換后的近似分量圖像與所述第二近似分量圖像的互信息量；

基于互信息量計算結果，獲取最大互信息量為配準結果；

如果所述最大互信息量不在預定的匹配范圍內，根據所述最大互信息量對應的幾何變換參數，確定超聲探頭的位置調整信息。

2.根據權利要求1所述的方法，其特征在于，分別對所述待配圖像與所述標準圖像進行小波變換，對應得到第一近似分量圖像和第二近似分量圖像，包括：

采用二維小波分解算法，分別對所述待配圖像與所述標準圖像進行至少兩層分解，得到分別包含至少兩層近似分量的第一近似分量圖像和第二近似分量圖像。

3.根據權利要求2所述的方法，其特征在于，對所述第一近似分量圖像進行至少一次幾何變換，包括：

分別對所述第一近似分量圖像的各層近似分量進行多次幾何變換；

計算每次變換后的近似分量圖像與所述第二近似分量圖像的互信息量，包括：

計算每次變換后的第N層近似分量，與所述第二近似分量圖像的第N層近似分量的互信息量，為第N層近似分量的互信息量；N為1至M之間的整數，M為分解層數。

4.根據權利要求3所述的方法，其特征在于，所述幾何變換參數包括旋轉角度和/或平移距離；

分別對所述第一近似分量圖像的各層近似分量進行多次幾何變換，包括：

按照尺寸從小到大的順序，對所述第一近似分量圖像的各層近似分量進行排序；

按照排序，順次對所述第一近似分量圖像的每層近似分量進行多次幾何變換；

其中，尺寸最小的第M層近似分量的幾何變換參數，在預定的參數范圍內；

第K-1層近似分量的幾何變換參數，至少由前一層近似分量的最大互信息量對應的幾何變換參數決定；K為2至M之間的整數。

5.根據權利要求4所述的方法，其特征在于，所述幾何變換參數包括旋轉角度和平移距離；

第M層近似分量在水平方向的尺寸為MA，在垂直方向的尺寸為MB；

第M層近似分量的旋轉角度在-10度與10度之間，水平方向的平移距離在-MA/2與-MA/2之間，垂直方向的平移距離在-MB/2與-MB/2之間；

第K層近似分量的最大互信息量對應的旋轉角度為θK，對應的水平方向的平移距離為SK，對應的垂直方向的平移距離為HK；

第K層近似分量的下一層近似分量的旋轉角度，在θK-2度與θK+2度之間；

第K層近似分量的下一層近似分量的水平方向的平移距離，在2SK-10與2SK+10之間；

第K層近似分量的下一層近似分量的垂直方向的平移距離，在2HK-10與2HK+10之間。

6.根據權利要求5所述的方法，其特征在于，針對同一層近似分量進行多次幾何變換時，相鄰兩次幾何變換的幾何變換參數的差值相等。

7.根據權利要求1所述的方法，其特征在于，根據所獲取的幾何變換參數，確定超聲探頭的位置調整信息，包括：

根據超聲探頭所在坐標系與超聲圖像坐標系之間的映射關系，將所獲取的幾何變換參數映射到超聲探頭所在坐標系，得到所述位置調整信息。

8.根據權利要求1所述的方法，其特征在于，所述方法還包括：

如果所述最大互信息量在預定的匹配范圍內，則基于所述實時圖像進行解剖結構的識別。