Nothing Special   »   [go: up one dir, main page]

TWI768575B - 三維影像動態矯正評估與矯具輔助設計方法及其系統 - Google Patents

三維影像動態矯正評估與矯具輔助設計方法及其系統 Download PDF

Info

Publication number
TWI768575B
TWI768575B TW109142659A TW109142659A TWI768575B TW I768575 B TWI768575 B TW I768575B TW 109142659 A TW109142659 A TW 109142659A TW 109142659 A TW109142659 A TW 109142659A TW I768575 B TWI768575 B TW I768575B
Authority
TW
Taiwan
Prior art keywords
dimensional
human body
spine
displacement
curve
Prior art date
Application number
TW109142659A
Other languages
English (en)
Other versions
TW202223919A (zh
Inventor
顏伯甫
黃鐘儀
林尚一
于天立
吳宗明
蔡宗汶
Original Assignee
財團法人工業技術研究院
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 財團法人工業技術研究院 filed Critical 財團法人工業技術研究院
Priority to TW109142659A priority Critical patent/TWI768575B/zh
Priority to CN202011489944.1A priority patent/CN114595522A/zh
Priority to US17/308,413 priority patent/US11475633B2/en
Publication of TW202223919A publication Critical patent/TW202223919A/zh
Application granted granted Critical
Publication of TWI768575B publication Critical patent/TWI768575B/zh

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • G06T17/205Re-meshing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00Computing arrangements based on specific mathematical models
    • G06N7/01Probabilistic graphical models, e.g. probabilistic networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • G06T19/20Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/41Medical
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/20Indexing scheme for editing of 3D models
    • G06T2219/2004Aligning objects, relative positioning of parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/20Indexing scheme for editing of 3D models
    • G06T2219/2016Rotation, translation, scaling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computer Graphics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Public Health (AREA)
  • Data Mining & Analysis (AREA)
  • Medical Informatics (AREA)
  • Computational Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Architecture (AREA)
  • Pathology (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Biomedical Technology (AREA)
  • Databases & Information Systems (AREA)
  • Probability & Statistics with Applications (AREA)
  • Mathematical Physics (AREA)
  • Artificial Intelligence (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Footwear And Its Accessory, Manufacturing Method And Apparatuses (AREA)
  • Testing Of Optical Devices Or Fibers (AREA)
  • Orthopedics, Nursing, And Contraception (AREA)

Abstract

一種三維影像動態矯正評估與矯具輔助設計方法,包括下列步驟。取得人體之三維掃描資訊取得該人體之複數個二維影像進行辨識,計算該些二維影像像素數值,並合成原始三維脊椎曲線。合成該些人體之二維影像與該三維掃描資訊。以影像變形預測矯正體型方法產生該人體之變形後體型。以脊椎材料特性與力學模型預估方法根據變形後體型預測矯具對該人體之施力位置、施力方向及施力大小的參數值。

Description

三維影像動態矯正評估與矯具輔助設計方法及其 系統
本發明是有關於一種三維影像動態矯正評估與矯具輔助設計方法及其系統。
目前脊椎側彎矯具製作仍普遍採用傳統石膏取模方式。然而,繃帶的纏繞方式與鬆緊不一會造成取模誤差,且取模時間較長,患者不易保持固定姿勢造成取模誤差。因此,石膏取模耗時費力又難保證取得正確的軀幹前後左右對稱性,脊椎側彎矯具也因此無法獲得正確的製作依據,後續患者的穿著矯正功能與舒適性也受到影響,有待進一步改善。
本發明係有關於一種三維影像動態矯正評估與矯具輔助設計方法及其系統,用以改善現行療程患者矯治過程的舒適性及矯治效果。
根據本發明之一方面,提出一種三維影像動態矯正評估與矯具輔助設計方法,包括下列步驟。取得人體之三維掃描資訊。取得該人體之複數個二維影像進行辨識,計算該些二維影像像素數值,並合成原始三維脊椎曲線。合成人體之該些二維影像與該三維掃描資訊。以影像變形預測矯正體型方法產生該人體之變形後體型。以脊椎材料特性與力學模型預估方法根據變形後體型預測矯具對該人體之施力位置、施力方向及施力大小的參數值。
根據本發明之一方面,提出一種三維影像動態矯正評估與矯具輔助設計系統,包括三維掃描裝置、運算單元、影像變形預測矯正體型單元、以及脊椎材料特性與力學模型預估單元。三維掃描裝置用以取得人體之三維掃描資訊。運算單元用以取得該人體之複數個二維影像進行辨識,計算該些二維影像像素數值,並合成原始三維脊椎曲線。運算單元用以合成人體之該些二維影像與該三維掃描資訊。影像變形預測矯正體型單元用以產生該人體之變形後體型。脊椎材料特性與力學模型預估單元根據該變形後體型預測矯具對該人體之施力位置、施力方向及施力大小的參數值。
為了對本發明之上述及其他方面有更佳的瞭解,下文特舉實施例,並配合所附圖式詳細說明如下:
100:三維影像動態矯正評估與矯具輔助設計系統
101:三維掃描資訊
102:深度影像
103:二維影像
104:人體之外型影像
105:原始三維脊椎曲線
106:矯正後三維脊椎曲線
107:彈性位能
108:位移量
109:彈性位能最小值
109A:變形後體型
110:三維掃描裝置
111:施力位置
112:施力方向
113:施力大小
114:三角網格
115:矯具
116:有限元素模型
130:運算單元
140:影像變形預測矯正體型單元
150:脊椎材料特性與力學模型預估單元
C:冠狀面
S:矢狀面
T:椎體旋轉
P:壓力
V:椎體
d,e,f,g,h,i,j,k,l,d’,e’,f’,g’,h’,i':點
g1,g2,i2,i3,g1’,g2’,h1’,i1’,i2’,i3’:中繼點
第1A圖繪示依照本發明一實施例的三維影像動態矯正評估與矯具輔助設計系統的示意圖; 第1B圖繪示依照本發明一實施例的三維影像動態矯正評估與矯具輔助設計方法的示意圖;第2圖繪示依照本發明一實施例的將原始三維脊椎曲線影像以及矯正後三維脊椎曲線影像疊合人體之外型影像的示意圖;第3圖繪示辨識人體之三維脊椎類型的示意圖;第4圖繪示矯具對人體之施力位置、施力方向及施力大小的示意圖;第5A至5D圖繪示利用影像變形預測矯正體型方法計算三角網格上的端點的彈性位能、位移量、彈性位能最小值及變形後體型的示意圖;第6圖繪示依照本發明一實施例的影像變形預測矯正體型方法的示意圖;第7圖繪示依照本發明一實施例的脊椎材料特性與力學模型預估方法的示意圖;第8圖繪示各椎體與椎間盤之楊氏模數(E值)的示意圖;及第9A至9F圖繪示另一實施例中利用影像變形預測矯正體型方法計算三角網格上的端點的彈性位能、位移量、彈性位能最小值及變形後體型的示意圖。
以下係提出實施例進行詳細說明,實施例僅用以作為範例說明,並非用以限縮本發明欲保護之範圍。以下是以相同/類似的符號表示相同/類似的元件做說明。
請參照第1A及1B圖,其中第1A圖繪示依照本發明一實施例的三維影像動態矯正評估與矯具輔助設計系統100的示意圖,第1B圖繪示依照本發明一實施例的三維影像動態矯正評估與矯具輔助設計方法的示意圖。
依照本發明之一實施例,此系統100包括一三維掃描裝置110、一運算單元130、一影像變形預測矯正體型單元140以及一脊椎材料特性與力學模型預估單元150。
三維掃描裝置110用以產生複數個深度影像102,對深度影像102進行建模,以產生人體之三維掃描資訊101,以建構人體之外型影像104。如第1B圖所示,在步驟S110中,對深度影像102進行建模,以取得人體之三維掃描資訊101。三維掃描裝置110例如由多台深度攝影機或光達模組所形成,除了藉由機構設計及光機電整合,設置完整包覆待掃描肢體軀幹的掃描空間以掃描人體之外型外,還配置運算及建模軟體,可以將掃描所取得的深度影像102,計算深度資訊,並且將影像疊合建模,使得三維掃描裝置110可以在很短的時間內建立人體三維模型。
運算單元130取得X光片的二維影像103,以產生原始三維脊椎曲線105及矯正後三維脊椎曲線106。原始三維脊椎曲線105例如是使用患者當下所拍攝的X光片疊合生成的影像,而矯正後三維脊椎曲線106例如是利用患者所拍攝的X光片,預估可能的矯正曲線疊合生成的影像。也就是說,運算單元130進行X光片的二維影像103處理,例如,運算單元130使用文字偵測(Text detection)實現X光片尺標數字及 尺標傾斜自動校正,並使用光學文字辨識(Optical Character Recognition,OCR)技術自動辨識二維影像103中之尺標數字並計算像素實際數值。如第1B圖所示,在步驟S120中,取得人體之複數個二維影像103進行辨識,計算此些二維影像103像素數值,並合成原始三維脊椎曲線105。
合成後的原始三維脊椎曲線105可與資料庫中的資料比對,以辨識原始三維脊椎曲線105為何種的脊椎類型。若判斷為需矯正之脊椎類型,判斷此脊椎為何種彎曲型態以規劃矯正計畫,並且合成出預測矯正後三維脊椎曲線106,將此原始三維脊椎曲線105、矯正後三維脊椎曲線106與三維掃描裝置110掃描之人體之外型影像104進行疊合。如第1B圖所示,在步驟S130中,合成人體之這些二維影像103與三維掃描資訊101。
另外,影像變形預測矯正體型單元140用以產生人體之變形後體型109A。當運算單元130判斷原始三維脊椎曲線105為需矯正之脊椎類型,影像變形預測矯正體型單元140可進行三維脊椎曲線與身體外型的變形,顯示原始三維脊椎曲線105、矯正後三維脊椎曲線106及人體之外型影像104的變化過程及結果,並可藉由大數據累積來修正此變形模型。如第1B圖所示,在步驟S140中,以影像變形預測矯正體型方法產生人體之變形後體型109A。
此外,脊椎材料特性與力學模型預估單元150用以預測矯具115對人體之施力位置111、施力方向112及施力大小113的參數值。當分析人員取得需矯正之人體脊椎之變形模型之後,脊椎材料特 性與力學模型預估單元150可根據生物力學模型及前後身體變化結果評估矯具施力位置111、施力方向112及施力大小113的參數值,以獲得最適當的施力參數,作為後續矯治管理與輔護具設計之依據。如第1B圖所示,在步驟S150中,以脊椎材料特性與力學模型預估方法根據變形後體型109A預測矯具對人體之施力位置111、施力方向112及施力大小113的參數值。
請參照第2圖及第3圖,其中第2圖繪示依照本發明一實施例的將原始三維脊椎曲線105影像以及矯正後三維脊椎曲線106影像疊合人體之外型影像104的示意圖,第3圖繪示辨識人體之三維脊椎曲線類型的示意圖。透過上述的三維影像動態矯正評估與矯具輔助設計系統100及其方法,可幫助骨科及復健科醫師預測矯正後體型並進行動態矯正規劃評估。在一實施例中,動態矯正規劃可包括預估矯正期間及動態輔助矯正設計。矯正進展與預期效果評估可協助判斷患者所需矯正期間,也可藉由矯正前後三維體型之變形模型,合理計算每次調整可改善之角度及所需時間。此外,動態矯正設計可依據矯正前後三維體型之變形模型,計算各階段的變形後體型109A及矯具115對人體之施力位置111、施力方向112及施力大小113,並歸納分類及分析矢狀面S的75個偏差樣態動態矯正、冠狀面C及椎體旋轉T等三面向偏差模型,即可估算各椎體所需矯正範圍及矯正所需強度。
請參照第4圖,其繪示矯具115對人體之施力位置111、施力方向112及施力大小113的示意圖。後續進行矯具輔助設計時,可搭配有限元素法分析人體脊椎的生物力學模型,對於三維脊椎側彎矯 正及矯具設計提供更大的幫助。對於脊椎側彎矯具設計的最佳化量化資訊,例如計算矯具115施以多大的壓力P使能產生預定中的側彎角度減少,施力點的方向,墊片的面積大小,各椎體的受力應力分布狀況等等,均可藉由生物力學模型評估以智慧化三維影像動態矯正所獲得的矯治規劃方案,經數值計算後可獲得「側彎角度施力大小」、「椎間盤應力」、「薦椎應力」等結果。以供醫師及矯具設計工程師將數值分析結果,進行綜合評量以決定個別患者最佳三維側彎矯治角度及用以進行達成矯治方案的脊椎側彎參數設計。
請參照第4圖,以特定患者脊椎側彎為基礎建立有限元素模型116,過程中由患者的電腦斷層掃描醫學影像為數據進行幾何建置。將患者的脊椎幾何轉換為有限元素模型116進行相關生物力學分析,但是生物力學模型中若採用且高品質的三維六面體元素,將會嚴重影響計算效率。為了解決這些問題,本實施例利用具有人體之冠狀面、矢狀面的X光片(二維影像)將高品質的三維模型降階為兩個冠狀面C與矢狀面S二維等效有限元素模型。本實施例利用此模型進行脊椎側彎矯正過程的受力情況分析,以此分析結果回饋至矯具設計。目的以模擬和迭代地改善脊椎側彎矯具設計,可在冠狀面C、矢狀面S兩個平面上計算矯治變形量所需的施力位置111、施力大小113、施力方向112與各椎體V所受之應力狀況,以作為矯具設計時之量化參數值。
請參照第5A至5D圖及第6圖,其中第5A至5D圖繪示利用影像變形預測矯正體型方法計算三角網格上的端點的彈性位能 107、位移量108、彈性位能最小值109及變形後體型109A的示意圖,第6圖繪示依照本發明一實施例的影像變形預測矯正體型單元140的各個部分執行的步驟的示意圖。
請參照第5A圖,將身體外型每一個三角網格114的每個端點(d,e,f)各自連接到原始三維脊椎曲線105上最近的點(g,h,i)或節點,假設每條線(f-e,f-d,f-i)由彈簧牽引並遵守虎克定律(F(x)=kx),其中k為彈性係數,x為位移量108。如第6圖所示,以影像變形預測矯正體型單元140的硬體部分及/或軟體部分執行各步驟為例,在步驟S210中,將身體外型的複數三角網格114的端點各自連結到原始三維脊椎曲線105上最近的點,各自計算此些三角網格114上的端點的彈性位能107。
請參照第5B圖,當脊椎矯正過程中,身體外型由原始三維脊椎曲線105逼近矯正後三維脊椎曲線106時,使得點(g,h,i)逼近點(j,k,l),各自計算三角網格114的每個端點(d,e,f)的彈性位能(Ep(x)=1/2kx2)。請參照第5C圖,當i逼近l時,以f為端點的三條線(f-e,f-d,f-i)會移動到以f’為端點的三條線(f’-e,f’-d,f’-l)而得到三角網格114的每個端點的最小彈性位能,以產生f與f’之位移量。如第6圖所示,在步驟S220中,當原始三維脊椎曲線105逼近矯正後三維脊椎曲線106時,各別計算此些三角網格114上的端點的彈性位能107使得彈性位能107達彈性位能最小值109,以產生這些端點的位移量108。
請參照第5D圖,依據相對應之身體外型參數設定,以迭代運算計算所有三角網格114的每個端點的彈性位能107,使得整體彈性位能達最小值,得到新三角網格114’(d’,e’,f’)並根據每個端點的位移量108產生 變形後體型。如第6圖所示,在步驟S230中,依據身體外型參數,迭代運算計算此些三角網格114上的端點的彈性位能107使整體位能達最小值,得到新三角網格114’並根據各自端點的位移量108產生變形後體型109A。在一實施例中,身體頸椎、胸椎及腰椎給予相同的力會有不同的變形量,上述的身體外型參數係指依據身體的力學模型計算變形的難易程度而給予不同的變形參數(例如楊氏係數)。
如上所述,透過上述第5A至5D圖及第6圖中的影像變形預測矯正體型方法可準確地預測欲矯正人體之各自端點的位移量108,據以產生矯正前後三維體型之變形模型,做為變形後體型109A。上述的矯正不限定一次完成脊椎曲線的矯正,可依照醫師及矯具設計工程師的判斷多階段完成脊椎曲線的矯正。也就是說,矯正後三維脊椎曲線106為複數張不同階段的矯正後三維脊椎曲線所形成。
請參照第7及8圖,第7圖繪示依照本發明一實施例的脊椎材料特性與力學模型預估單元150的各個部分執行的步驟的示意圖,第8圖繪示各椎體與椎間盤之楊氏模數(E值)的示意圖。
椎體與椎間盤材料特性與生物力學模型預估方法例如為透過生物力學模擬配合實驗迴歸分析進行迭代求得。在一實施例中,實驗部分例如為患者以吊單槓受自重力作用下,對脊椎進行冠狀面的X光成像,各脊體與椎間盤受到自重力作用後所產生的位移(
Figure 109142659-A0305-02-0011-6
)作為目標函數。在生物力學模型中以患者吊單槓受自重力作用之條件進行建模,配合相關的椎體V與椎間盤幾何形狀來進行最佳化回歸分析以評估椎體V與椎間盤之剛度。如第7圖所述,以脊椎材料特性與力 學模型預估單元150的硬體部分及/或軟體部分執行各步驟為例,步驟S310中,以人體吊單槓受自重力作用下,對人體脊椎進行X光成像,人體之各椎體與椎間盤受到自重力作用後所產生的位移作為目標函數,進行最佳化回歸分析以評估椎體V與椎間盤之剛度。
最佳化回歸分析如下:藉由楊氏係數(E值)的迭代計算可得到每節椎體V與椎間盤元素節點的位移(δl),計算各節點之位移(δl)與目標函數的位移(
Figure 109142659-A0305-02-0012-7
)進行誤差分析,楊氏係數(E值)迭代收斂確認條件為此些節點的位移均方根值達到最小值(誤差約為0.1至0.3),計算後所得知各椎體V與椎間盤之楊氏係數(E值)如第8圖所示。如第7圖所述,步驟S320中,以脊椎拉直時施力狀態的迭代計算可得到人體每節椎體V與椎間盤元素節點的位移,對這些節點的位移與目標函數的位移進行誤差分析,使這些位移的均方根值達到最小值,獲得最適當的施力位置111、施力方向112及施力大小113的參數值。
請參照第9A至9F圖,其繪示另一實施例中利用影像變形預測矯正體型方法計算三角網格114上的端點的彈性位能107、位移量108、彈性位能最小值109及變形後體型109A的示意圖。
請參照第9A圖,依照楊氏係數代入的中繼點密度,在原始三維脊椎曲線105上產生例如相對應密度的點。此外,請參照第9B圖,依照點雲密度在身體內部產生若干中繼點,中繼點例如可依人體構造產生或均勻產生。請參照第9C圖,將身體外型每一個三角網格114的每個端點(d,e,f)各自連接到最近的中繼點(g1,h1,i1),並尋找可能的中繼點,例如:中繼點(g1) 再以中繼點(g2)連接到原始脊椎曲線105上最近的點(g),中繼點(i1)再以中繼點(i2、i3)連接到原始脊椎曲線105上最近的點(i)。請參照第9D圖,假設連接的每條線由彈簧牽引並遵守虎克定律(F(x)=kx),目前的長度為最低位能態。請參照第9E圖,當脊椎矯正過程中,依使用者需求調整或旋轉脊椎,身體外型由原始三維脊椎曲線105逼近矯正後三維脊椎曲線106時,節點(g,h,i)逼近節點(g’,h’,i'),使得各中繼點(g1,g2,h1,i1,i2,i3)移動至新的錨點(g1’,g2’,h1’,i1’,i2’,i3’),各自計算三角網格114的每個端點及中繼點的彈性位能107而使彈性位能107達到彈性位能最小值109,以產生這些端點的位移量108及旋轉量。請參照第9E及9F圖,依據身體外型參數,迭代運算計算此些三角網格114上的端點及中繼點的彈性位能107使整體彈性位能達最小值,得到新三角網格114’(d’,e’,f’),以產生變形後體型109A。
如上所述,透過上述第9A至9D圖中的影像變形預測矯正體型方法可準確地預測欲矯正人體之各自端點的位移量108及旋轉量,據以產生矯正前後三維體型之變形模型,做為變形後體型。
本發明上述實施例的三維影像動態矯正評估與矯具115輔助設計方法,可透過三維掃描取得三維掃描資訊,並合成人體之二維影像與三維掃描資訊,藉以辨識人體之三維脊椎曲線類型。除了可提供設計者準確的矯正後三維脊椎曲線,更可累積大數據給分析人員提供所需要的動態矯正評估與矯具115輔助設計,克服傳統石膏取模方式製作矯具115的缺點,以改善現行療程患者矯治過程的舒適性及矯治效果。
綜上所述,雖然本發明已以實施例揭露如上,然其並非 用以限定本發明。本發明所屬技術領域中具有通常知識者,在不脫離本發明之精神和範圍內,當可作各種之更動與潤飾。因此,本發明之保護範圍當視後附之申請專利範圍所界定者為準。
S110-S150:步驟

Claims (14)

  1. 一種三維影像動態矯正評估與矯具輔助設計方法,包括:取得人體之三維掃描資訊;取得該人體之複數個二維影像進行辨識,計算該些二維影像像素數值,並合成原始三維脊椎曲線;合成該人體之該些二維影像與該三維掃描資訊;以影像變形預測矯正體型方法產生該人體之變形後體型;以及以脊椎材料特性與力學模型預估方法根據該變形後體型預測一矯具對該人體之施力位置、施力方向及施力大小的參數值,其中該影像變形預測矯正體型方法包括:將該人體之該些二維影像疊合生成該原始三維脊椎曲線以及矯正後三維脊椎曲線,並將該原始三維脊椎曲線以及該矯正後三維脊椎曲線疊合該人體之外型影像;在該人體內部建立複數個中繼點,將身體外型的複數三角網格的端點各自連接到最近的中繼點,並將該些中繼點各自連接到該原始三維脊椎曲線上最近的點;及各自計算該些三角網格上的端點及該些中繼點的彈性位能。
  2. 如請求項1所述之三維影像動態矯正評估與矯具輔助設計方法,其中該影像變形預測矯正體型方法包括: 當該原始三維脊椎曲線逼近該矯正後三維脊椎曲線時,計算該些三角網格上的端點及該些中繼點的彈性位能,使該彈性位能達最小值,以產生該些端點的位移量及旋轉量;以及依據身體外型參數,迭代運算計算該些三角網格上的端點及該些中繼點的彈性位能使整體彈性位能達彈性位能最小值,得到新三角網格並根據各自端點的位移量產生該變形後體型。
  3. 如請求項1所述之三維影像動態矯正評估與矯具輔助設計方法,其中該脊椎材料特性與力學模型預估方法包括:以該人體吊單槓受自重力作用下,取得該人體之該些二維影像,該人體之各脊體與椎間盤受到自重力作用後所產生的位移作為目標函數,進行最佳化回歸分析以評估脊體與椎間盤之剛度;以脊椎拉直時施力狀態的迭代計算可得到該人體每節椎體與椎間盤元素節點的位移,對該些節點的位移與該目標函數的位移進行誤差分析,使該些節點的位移均方根值達到最小值,獲得該施力位置、該施力方向及該施力大小的參數值。
  4. 如請求項3所述之三維影像動態矯正評估與矯具輔助設計方法,其中該最佳化回歸分析包括藉由楊氏係數的迭代計算得到每節椎體與椎間盤元素節點的位移,該楊氏係數迭代收斂確認條件為該些節點的位移均方根值達到最小值。
  5. 如請求項1所述之三維影像動態矯正評估與矯具輔助設計方法,其中取得該人體之該些二維影像進行辨識包括使用文字偵測或光學文字辨識自動辨識該些二維影像。
  6. 如請求項3所述之三維影像動態矯正評估與矯具輔助設計方法,其中該脊椎材料特性與力學模型預估方法包括以有限元素法分析該人體脊椎的生物力學模型,利用具有該人體之冠狀面、矢狀面的該些二維影像將三維模型降階為兩個冠狀面與矢狀面二維等效有限元素模型。
  7. 如請求項1所述之三維影像動態矯正評估與矯具輔助設計方法,其中該動態矯正評估包括依據矯正前後三維體型之變形模型,計算各階段的變形後體型、該施力位置、該施力方向及該施力大小。
  8. 如請求項1所述之三維影像動態矯正評估與矯具輔助設計方法,其中該矯具輔助設計包括依據矯正前後三維體型之變形模型,計算該矯具施加的壓力使能產生預定中的側彎角度減少、施力點的方向、墊片的面積大小及各椎體的受力應力分布。
  9. 如請求項2所述之三維影像動態矯正評估與矯具輔助設計方法,其中該矯正後三維脊椎曲線為複數張不同階段的矯正後三維脊椎曲線。
  10. 一種三維影像動態矯正評估與矯具輔助設計系統,包括:三維掃描裝置,用以取得人體之三維掃描資訊;運算單元,用以取得該人體之複數個二維影像進行辨識,計算該些二維影像像素數值,並合成原始三維脊椎曲線,該運算單元用以合成該些人體之二維影像與該三維掃描資訊; 影像變形預測矯正體型單元,用以產生該人體之變形後體型;以及脊椎材料特性與力學模型預估單元,根據該變形後體型預測矯具對該人體之施力位置、施力方向及施力大小的參數值,其中該影像變形預測矯正體型單元包括:第一部分,用以將該人體之該些二維影像疊合生成該原始三維脊椎曲線以及矯正後三維脊椎曲線,並將該原始三維脊椎曲線以及該矯正後三維脊椎曲線疊合該人體之外型影像;第二部分,用以在該人體內部建立複數個中繼點,將身體外型的複數三角網格的端點各自連接到最近的中繼點,並將該些中繼點各自連接到該原始三維脊椎曲線上最近的點;及第三部分,用以各自計算該些三角網格上的端點及該些中繼點的彈性位能。
  11. 如請求項10所述之三維影像動態矯正評估與矯具輔助設計系統,其中該影像變形預測矯正體型單元包括:當該原始三維脊椎曲線逼近該矯正後三維脊椎曲線時,該第三部分計算該些三角網格上的端點及該些中繼點的彈性位能,使該彈性位能達最小值,以產生該些端點的位移量及旋轉量;以及第四部分,依據身體外型參數,迭代運算計算該些三角網格上的端點及該些中繼點的彈性位能使整體彈性位能達最小值,得到新三角網格並根據各自端點的位移量產生該變形後體型。
  12. 如請求項10所述之三維影像動態矯正評估與矯具輔助設計系統,其中該脊椎材料特性與力學模型預估單元包括:第一部分,以該人體吊單槓受自重力作用下,取得該人體之該些二維影像,該人體之各脊體與椎間盤受到自重力作用後所產生的位移作為目標函數,進行最佳化回歸分析以評估脊體與椎間盤之剛度;以及第二部分,以脊椎拉直時施力狀態的迭代計算可得到該人體每節椎體與椎間盤元素節點的位移,對該些節點的位移與該目標函數的位移進行誤差分析,使該些節點的位移均方根值達到最小值,獲得該施力位置、該施力方向及該施力大小的參數值。
  13. 如請求項12所述之三維影像動態矯正評估與矯具輔助設計系統,其中該最佳化回歸分析包括藉由楊氏係數的迭代計算得到每節椎體與椎間盤元素節點的位移,該楊氏係數迭代收斂確認條件為該些節點的位移均方根值達到最小值。
  14. 如請求項12所述之三維影像動態矯正評估與矯具輔助設計系統,其中該脊椎材料特性與力學模型預估單元以有限元素法分析該人體脊椎的生物力學模型,利用具有該人體之冠狀面、矢狀面的該些二維影像將三維模型降階為兩個冠狀面與矢狀面二維等效有限元素模型。
TW109142659A 2020-12-03 2020-12-03 三維影像動態矯正評估與矯具輔助設計方法及其系統 TWI768575B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
TW109142659A TWI768575B (zh) 2020-12-03 2020-12-03 三維影像動態矯正評估與矯具輔助設計方法及其系統
CN202011489944.1A CN114595522A (zh) 2020-12-03 2020-12-16 三维影像动态矫正评估与矫具辅助设计方法及其系统
US17/308,413 US11475633B2 (en) 2020-12-03 2021-05-05 Method and system for 3D image dynamic correction evaluation and orthotics auxiliary design

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW109142659A TWI768575B (zh) 2020-12-03 2020-12-03 三維影像動態矯正評估與矯具輔助設計方法及其系統

Publications (2)

Publication Number Publication Date
TW202223919A TW202223919A (zh) 2022-06-16
TWI768575B true TWI768575B (zh) 2022-06-21

Family

ID=81802865

Family Applications (1)

Application Number Title Priority Date Filing Date
TW109142659A TWI768575B (zh) 2020-12-03 2020-12-03 三維影像動態矯正評估與矯具輔助設計方法及其系統

Country Status (3)

Country Link
US (1) US11475633B2 (zh)
CN (1) CN114595522A (zh)
TW (1) TWI768575B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140169647A1 (en) * 2011-08-08 2014-06-19 Instytut Biologii Doswiadczalnej Im. M. Nenckiego Pan Method and a system for processing an image comprising dendritic spines
CN106214302A (zh) * 2016-07-13 2016-12-14 广东康沃森医疗科技有限责任公司 一种脊柱侧弯矫形支具的制备方法
US20170143426A1 (en) * 2012-03-12 2017-05-25 Nuvasive, Inc. Systems and methods for performing spinal surgery
TW201926253A (zh) * 2017-11-24 2019-07-01 財團法人工業技術研究院 三維模型建構方法及其系統
CN110634554A (zh) * 2018-06-22 2019-12-31 王靖维 脊椎影像注册方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE602007013502D1 (de) 2006-09-18 2011-05-12 Mediguide Ltd Verfahren und System zur Navigation durch ein verschlossenes röhrenförmiges Organ
JP6122864B2 (ja) 2011-11-18 2017-04-26 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. ライブ画像と解剖学的構造表示とのペアリング
TWM443156U (en) 2011-11-23 2012-12-11 Chien Hui Chuan System for three-dimensional shape and size measurement
DE102012111385B4 (de) 2012-11-23 2018-05-03 Diers Engineering Gmbh Bestimmen der räumlichen Lage und Orientierung der Wirbelkörper der Wirbelsäule
US11311402B2 (en) 2013-05-06 2022-04-26 Aspen Medical Partners, Llc Scoliosis brace
TWI605805B (zh) 2016-11-22 2017-11-21 wen-hua Xu Scoliosis correction clothes
EP3651637A4 (en) 2017-07-12 2021-04-28 K2M, Inc. SYSTEMS AND METHODS FOR MODELING SPINE AND TREATMENT OF SPINE BASED ON SPINE MODELS
CN108986891A (zh) 2018-07-24 2018-12-11 北京市商汤科技开发有限公司 医疗影像处理方法及装置、电子设备及存储介质
CN110327146A (zh) 2019-07-02 2019-10-15 武汉联影医疗科技有限公司 一种矫形器设计方法、装置和服务器
CN111311742B (zh) 2020-03-27 2023-05-05 阿波罗智能技术(北京)有限公司 三维重建方法、三维重建装置和电子设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140169647A1 (en) * 2011-08-08 2014-06-19 Instytut Biologii Doswiadczalnej Im. M. Nenckiego Pan Method and a system for processing an image comprising dendritic spines
US20170143426A1 (en) * 2012-03-12 2017-05-25 Nuvasive, Inc. Systems and methods for performing spinal surgery
CN106214302A (zh) * 2016-07-13 2016-12-14 广东康沃森医疗科技有限责任公司 一种脊柱侧弯矫形支具的制备方法
TW201926253A (zh) * 2017-11-24 2019-07-01 財團法人工業技術研究院 三維模型建構方法及其系統
CN110634554A (zh) * 2018-06-22 2019-12-31 王靖维 脊椎影像注册方法

Also Published As

Publication number Publication date
CN114595522A (zh) 2022-06-07
US20220180606A1 (en) 2022-06-09
TW202223919A (zh) 2022-06-16
US11475633B2 (en) 2022-10-18

Similar Documents

Publication Publication Date Title
US11443494B2 (en) Method for estimating at least one of shape, position and orientation of a dental restoration
CN110599528B (zh) 一种基于神经网络的无监督三维医学图像配准方法及系统
Bucki et al. A fast and robust patient specific finite element mesh registration technique: application to 60 clinical cases
CN112638312A (zh) 使用深度学习的自动化正畸治疗计划
JP6633353B2 (ja) 脊柱変形矯正固定術支援装置、脊柱変形矯正固定術支援方法、プログラム、及び脊柱変形矯正固定術に使用するロッド製造方法
CN102651145B (zh) 股骨三维模型可视化方法
WO2015200114A1 (en) Additive smoothing of sharp concave edges on designed 3d printable polygonal mesh models
RU2700114C1 (ru) Устройство, система визуализации и способ коррекции медицинского изображения молочной железы
CN108198235A (zh) 一种三维超声重建方法、装置、设备及存储介质
CN115272582A (zh) 身体建模的系统和方法
CN115116586A (zh) 一种基于联合配准的可变形统计图谱构建方法
EP3921809A1 (en) Method for generating dental models based on an objective function
TWI768575B (zh) 三維影像動態矯正評估與矯具輔助設計方法及其系統
Gutiérrez-Gil et al. Bone remodelling-based numerical evaluation of personalized masticatory forces from CT-scans
US9142042B2 (en) Methods and systems to produce continuous trajectories from discrete anatomical shapes
JP5943267B2 (ja) 個人モデルデータの生成方法、生成プログラム、および生成システム
Hansen et al. Adaptive parametrization of multivariate B-splines for image registration
Makki et al. A new geodesic-based feature for characterization of 3D shapes: application to soft tissue organ temporal deformations
KR102183079B1 (ko) 안와 임플란트 모델링 및 제작을 위한 방법 및 장치
CN116264995A (zh) 脊柱侧弯矫形支具的个性化3d建模方法、计算机设备
Ballit et al. Visual sensor fusion with error compensation strategy toward a rapid and low-cost 3D scanning system for the lower residual limb
Kumar et al. Improving visibility of stereo-radiographic spine reconstruction with geometric inferences
Coevoet et al. Introducing interactive inverse FEM simulation and its application for adaptive radiotherapy
Baka et al. Correspondence free 3D statistical shape model fitting to sparse X-ray projections
KR102714028B1 (ko) 교정치료 정보 제공 시스템 및 이의 동작 방법