TWI803848B

TWI803848B - 一種帕金森氏症上肢震顫分析裝置與方法

Info

Publication number: TWI803848B
Application number: TW110111425A
Authority: TW
Inventors: 巫建興; 林家宏
Original assignee: 國立勤益科技大學
Priority date: 2021-03-29
Filing date: 2021-03-29
Publication date: 2023-06-01
Also published as: TW202237029A

Abstract

本發明提出一種帕金森氏症上肢震顫分析裝置與方法，特別是，本發明係一種陣列式雷達微波束之時域特徵參數於帕金森氏症的上肢震顫量化分析裝置。

Description

一種帕金森氏症上肢震顫分析裝置與方法

本發明提出一種帕金森氏症上肢震顫分析裝置與方法，特別是關於一種陣列式雷達微波束之時域特徵參數於帕金森氏症的上肢震顫量化分析裝置。

在醫學領域中，「帕金森氏症(Parkinson's disease,PD)」是一種會嚴重影響中樞神經系統的慢性神經退化疾病，而帕金森氏症多好發於中老年人的身上，且有逐年普及上升的趨勢。典型的帕金森氏症主要靠症狀診斷，而帕金森氏症的主要表現包括了運動症候群，其中運動症候群係數種症狀類似的運動障礙的合稱。

在帕金森氏症的數種運動障礙的主要表現症狀中，其中「原發性震顫(Essential Tremor,ET)」係為主要的表現方式。原發性震顫常會伴隨著語調和輕微步態異常，病患通常先由上肢開始(主要影響上肢)，逐漸影響到頭，下肢，軀幹，發聲和面部肌肉，主要表現為「姿態性震顫(Postural Tremor,PT)」。該「姿態性震顫」亦同時含有運動性(Kinetic Tremor)，意向性(Action Tremor)，或靜止性震顫(Resting Tremor)，而震顫的頻率約為4Hz至8Hz。

而前述「姿態性震顫」的震顫是一種反射性的，通常會迅速出現及消失(僅持續數秒)，而隨著震顫幅度的增加，會變得難以控制震顫，甚至會影響到正常工作。「姿態性震顫」的震顫幅度，震顫頻率在不同的動作，不同的姿態下，常會有所變動，但仍可靠自我抑制震顫。

在臨床醫學上，醫師可以由震顫的頻率，震顫的時間，以及震顫的振幅強度，藉以評估震顫的嚴重程度。通常帕金森氏症發病時，其頻率約為8Hz至12Hz，而隨著病態發展和年齡的增加，震顫頻率會逐漸降低，而震顫幅度則會逐漸增加。

於傳統的臨床醫學領域，對於帕金森氏症(PD)已有早期診斷的方法與相關技術，包括了「抽血」，「腦部電腦斷層掃描(CT)」技術，「磁振造影(MRI)」技術，「功能性磁振造影(Functional MRI,FMRI)」技術，亦或是「腦部正子掃描」(PET)。

如前述，其中「抽血」係藉由偵測腦部不同區域的血氧濃度相關參數，進而評估大腦內執行各項功能所對應的腦區。

而如前述，「磁振造影」技術和「功能性磁振造影」技術，係藉由強磁場造影，得到腦部內不同組織之對比影像，再利用電腦的計算取得影像，藉由臨床醫師判讀相關病灶，並做出最適當的臨床診斷。

但於臨床實際運用上，藉由所取得的影像並無法完全瞭解腦部的實際運作狀況，主要由於中腦部位的黑質細胞數量並不多，即使黑質細胞全部死亡，該部位的神經結構仍無明顯變化，故而，即使使用了「磁振造影」技術，「功能性磁振造影」技術，或是「腦部電腦斷層掃描」技術等影像掃描技術，仍舊無法精確分辨帕金森氏症疾病與正常腦部的差異，對於疾病早期診斷並無太大的幫助。

再者，由於「磁振造影」技術，「功能性磁振造影」技術，或是「腦部電腦斷層掃描」技術等影像掃描技術的設備相當昂貴，且設備體積占地甚為鉅大，而於進行檢查時，更需要花費相當長久的時間，加上使用的成本過高，且無法進行早期篩檢，且前述的數種磁振造影設備，僅能適用於醫療環境下，並無法適用於帕金森氏症患者的居家照護使用。

除了前述儀器設備外，另有一種「腦部AV133正子掃描」(PET)，其係一種以放射性追蹤劑，根據放射性訊號的減弱，精準評估帕金森氏症的嚴重程度，得以進行區隔正常老化或帕金森氏症，而可作為正確評估藥效，以及預後或是疾病嚴重程度的基準，然而，該種「腦部AV133正子掃描」技術亦有前述設備的使用成本問題與使用空間問題，仍舊無法適用於帕金森氏症患者的居家照護使用。

請參考圖1，以下係說明於傳統技術領域，一種都卜勒(Doppler)雷達微波的感測基本原理。因該感測基本原理可使用於本發明，故於此處提出基本原理說明，其中該都卜勒(Doppler)雷達微波的感測基本原理，係經由天線發射入射電波到目標物，碰到(待測)物體產生反射波，反射波的頻率會隨著(待測)物體移動的狀態而改變。

仍請參考圖1，倘使用於本發明之一種帕金森氏症上肢震顫分析裝置於感測時，若受測者上肢發生規則或不規則震顫時，入射波碰到震顫的上肢，則波長與相位會改變，隨著顫抖頻率及顫抖振幅大小，將會有明顯的差異。

續請參考圖1，使用於本發明一種帕金森氏症上肢震顫分析裝置中，具有雷達微波之入射電波，係由振盪器(Oscillator)101提供載波中頻訊號(IF Frequency)，透過諧振電路(Resonant Circuit)102決定振盪頻率，而具有負電阻特性的主動元件(Negative Resistance Circuit)103，會將直流功率轉換成射頻功率，再由升頻元件104升頻至射頻頻率，而透過天線TX(Transmit Antenna)發射出入射電波信號。

如圖1所示，藉由壓控振盪器105以改變振盪頻率，在特定的頻率下，使迴授增益值為｜H(S)｜=1，以及相角∠H(S)=0°(其中當S=jω _IF，ω為頻率)時，即X為輸入訊號，H為放大器，Y為輸出，而產生轉移函數，則可形成振盪物理現象。

如圖1所示，在接收天線RX(Receive Antenna)接受反射電波信號後，先由降頻元件降至中頻，振盪器101會提供本地訊號(Local Signal)，並由混頻器(Mixer)將反射電波信號與本地信號混合，如(01)所示：V(t)=A _IF cos(ω _IF t)+A _LO cos(ω _LO t) (01)

其中(01)式的符號表示為：V為反射波訊號；t為時間；A _IF為中頻訊號的振幅；A _LO 為本地訊號的振幅；ω_IF=2πf _IF為中頻訊號的頻率；ω_LO=2πf _LO為本地訊號(Local Signal,LoS)的頻率(Local Frequency,Lo)。

仍如圖1所述，之後，在混頻器輸出後，可產生頻率為f _LO±f _IF的輸出訊號，再送入後端的數位訊號處理器(DSP)，以進行訊號前置處理，特徵參數萃取，以及分析等應用。

續如圖1所示，當將前述傳統技術的都卜勒雷達微波，應用於非接觸生理訊號的感測時，假設有單一都卜勒雷達微波的入射電波信號T(t)，則可表示為：T(t)=A _T cos(ω _T t) (02)

而經過空氣傳輸一段d距離後，碰到震顫的上肢後反射，其反射信號R(t)，可表示為：

其中λ為雷達微波發射頻率在空氣中之波長。

而上肢震顫造成的位移x(t)，可表示為A _x sin(ω _x t)，將信號T(t)及信號R(t)混頻後之基頻訊號B(t)為：

透過頻率分析後，可求得混頻後之頻率成份為ωx，此頻率即為上肢震顫所造成的頻率。

故而如前述，當透過都卜勒雷達微波感測原理的運用，可使得雷達微波運用於非接觸感測生理狀態的偵測，且雷達微波除可用於動態物體的偵測，目前亦已應用於無人機/載具避撞控制，智慧交通監控，以及入侵偵測等相關領域的應用。

有鑑於解決前述問題，可達到早期發現早期治療的效果，本發明提出一種帕金森氏症上肢震顫分析裝置與方法。

本發明之一種帕金森氏症上肢震顫分析裝置，包含了以下組件：壓控振盪器，功率放大器，兩個濾波器(Filter)包含第一濾波器以及第二濾波器，天線發射端，混波器，接收天線，類比數位轉換器，以及數位信號處理器。

本發明之一種帕金森氏症上肢震顫分析方法，包含了以下的步驟：首先，擷取反射中頻訊號，進行訊號數位濾波處理，進行訊號放大，執行時域特徵參數萃取，以及，進行訊號時域分析。接著，執行頻域特徵參數萃取，以及，進行訊號頻域分析。繼續，進行訊號映射轉換，進行單點灰階或是彩色特徵參數轉換，以及，進行二維特徵圖樣分析，藉以評估上肢震顫的嚴重程度。

本發明係一種陣列式雷達微波束之時域特徵參數於帕金森氏症的上肢震顫量化分析裝置與方法，可藉由非侵入的感測方式，發展適用於帕金森氏症病患，於居家自我評估及照護使用的智能偵測裝置。

本發明一種帕金森氏症上肢震顫分析裝置與方法，係以個人化為中心之智能感測，智慧篩檢，行動顯示，以及個人化決策支援等功能之照護使用裝置，以達到預防與早期診斷的目的。

本發明一種帕金森氏症上肢震顫分析裝置與方法，其具有上肢動態震顫量測，震顫程度量化，以及自動數值資料分析等預期性的功能。

本發明一種帕金森氏症上肢震顫分析裝置與方法，係應用都卜勒雷達微波感測基本原理，進行帕金森氏症(PD)上肢震顫的偵測，而於感測時，正常受測者的上肢是固定不動的，且反射電波的頻率會和發射電波的頻率非常接近。

本發明一種帕金森氏症上肢震顫分析裝置與方法，係應用雷達微波感測裝置，以雷達微波偵測動態物體，而可偵測受測者之上肢動態的細微變化。

本發明一種帕金森氏症上肢震顫分析裝置與方法，除可偵測移動物體之外，雷達微波亦可以測量物體的範圍和徑向速度。

本發明之一種帕金森氏症上肢震顫分析裝置與方法，其優點係可不受環境因素(白晝夜、溫度、濕度)的影響，並藉由非接觸的感測方式，能夠檢測目標物體的動態等相關參數。

為讓上述目的，技術特徵，以及優點能更明顯易懂，下文係以較佳實施例配合所附圖式進行詳細說明。

101:振盪器

102:諧振電路

103:主動元件

104:升頻元件

105:壓控振盪器

201:壓控振盪器

202:功率放大器

203:第一濾波器

204:發送天線

205:混波器

206:接收天線

207:第二濾波器

208:低雜訊放大器

210:類比數位轉換器

211:數位信號處理器

310:平面式雷達微波陣列感測裝置

311:雷達微波感測器元件

31:背後瓣波

32:主要瓣波

33:兩側瓣波

a:間距

有關本發明之前述及其它許多優點，於以下配合參考圖示之一個較佳實施例的詳細說明中，將可更清楚呈現，其中如下：圖1係說明於傳統技術領域，一種都卜勒雷達微波的感測基本原理。

圖2係本發明之一種「帕金森氏症上肢震顫分析」裝置。

圖3A係本發明之平面式雷達微波陣列感測裝置之作動說明。

圖3B係本發明之平面式雷達微波陣列感測裝置的基本組成。

圖3C係本發明之平面式雷達微波陣列感測裝置的每個單一的雷達微波感測器元件。

圖4係本發明之一種「帕金森氏症上肢震顫分析」方法。

圖5A係本發明一種「帕金森氏症上肢震顫分析」的初始狀態顯示圖。

圖5B係本發明一種「帕金森氏症上肢震顫分析」的結果顯示圖。

以下請參照所附圖式說明與敘述，以對本發明之實施形態據以示意圖進行描述。且於示意圖式中，相同之元件符號表示相同之元件，而為求清楚說明，元件之大小或厚度可能誇大顯示。

如圖2所示，係本發明之一種「帕金森氏症上肢震顫分析」裝置，包括壓控振盪器(Power Source and Oscillator)201，功率放大器(Amplifier)202，兩個濾波器(Filter)包含第一濾波器203以及第二濾波器207，發送天線(TX)204，混波器(Mixer)205，以及接收天線(RX)206。

另外，仍如圖2所示，尚包括類比數位轉換器(Analog-to-Digital Converter,ADC)210，以及數位信號處理器(Digital Signal Processor,DSP)211。

如圖2所示，本發明係以壓控振盪器201之輸出訊號，透過功率放大器202的放大後，經由第一濾波器203，而分離成兩個同相位且振幅大小相同的信號，分別傳送至發送天線204，以及混波器205，而進入前述混波器205的信號，即為本地信號(LoS)。

續如圖2所示，當發射電波信號偵測到物體後反射之訊號，經由接收天線206接受後，經過第二濾波器207，再透過低雜訊放大器208進行放大，再饋入至混波器205，當以混波器205降頻後，會產生基頻訊號(IF)。在時域上，可作為觀察到生物體之呼吸或心跳的週期性相關應用，進而計算出呼吸及心跳之頻率。

仍如圖2所示，本發明之雷達微波感測裝置可分為兩部分，其一部分為前述的無線電頻率發射單元，係透過射頻技術產生雷達微波，再透過發送天線204加以發射電波，並以接收天線206接收反射波，進行差頻分析。

猶如圖2所示，另外一部分為嵌入式處理單元，係將混頻器輸出則產生的中頻訊號，送至類比數位轉換器210，以及數位信號處理器211，藉以進行數位化，以及時域/頻域參數分析，接者，再進行計算應用分析，以及上肢震顫量化分析，藉以完成整個偵測過程。

仍如圖2所示，本發明藉由時域反射電波信號，透過類比數位轉換器(Analog-to-Digital Converter,ADC)210，轉換成數位訊號 Bn後，其中n=1,2,3,…,N，N為取樣點數，可藉由零交越點(Zero-Crossing,ZC)，以及振幅差((Amplitude Difference,AD)的總合，Willison振幅(Willison Amplitude,WANP)，以及振盪訊號的平均能量(Mean Energy,ME)等參數，直接以時域訊號取得時域參數(Time-domain Parameter)，觀察上肢震顫的嚴重程度，如(05)式至(10)式所示如下：

零交越點：

振幅差總合：

Willison振幅：

振盪訊號的平均能量：

當零交越點(ZC)，振幅差(AD)總合，Willison振幅(WAMP)，以及，振盪訊號的平均能量(ME)等物理量，均具有明顯的變化時，則可顯示出上肢有明顯的震顫的症狀。

本發明以前述四個參數所呈現的量化刻度，得以進行正常，原發性震顫(ET)，以及帕金森氏症(PD)等相關震顫的篩選。值得注意的是，本發明係應用都卜勒雷達微波感測基本原理，進行帕金森氏症(PD)上肢震顫的偵測，而於感測時，正常受測者的上肢是固定不動的，且反射電波的頻率會和發射電波的頻率非常接近。而本發明具有一使用上的優勢，即通常雷達微波感測具有「非侵入」及「非接觸」的特點，故而當使用雷達微波感測受測者/使用者時，不一定需要穿戴在身上。

如圖3A所示，係本發明之平面式雷達微波陣列感測裝置310系統之作動說明，亦即為圖2所示之接收天線(RX)206之作動說明。而仍如圖3A所示，本發明之平面式雷達微波陣列感測裝置310系統，包含一平面式雷達微波陣列與一可動結構元件(圖3A中未式)，而該可動結構元件係以馬達帶動，可改變雷達微波陣列的掃描角度，而可達到作動的效果。

仍如圖3A所示，本發明係以「平面式雷達微波陣列感測裝置」310系統，發射出「陣列式雷達波束」，以集中波束進行上肢手部的震顫偵測，再以具有接收反射波功能的「接收器單元」(圖3A中未示)以接收反射波。

如圖3B係本發明之平面式雷達微波陣列感測裝置310系統的雷達微波感測器元件311，如在前述裝置310平面的平台上，可裝置M×M個(亦及複數個)的雷達陣列，最少為一個2×2的雷達陣列(若當M=2；則有4個雷達微波感測器元件)。而該平台可提供一個方位角上M×M個波束的陣列，在同一時間可將偵測範圍，從點擴展到區域波束所能覆蓋的範圍，藉由嵌入式系統快速切換M×M個驅動波束的控制(近似同步)，此種切換驅動方式可有效率控制及減少後功率放大器的損耗，使用陣列式掃描可增加掃描的覆蓋範圍，及降低方向所造成影響。

仍如圖3B，當M×M個雷達微波感測元件，裝設在前述裝置310系統平面的平台上，每個雷達微波感測器元件311，具有包含發射電波與接收反射電波的功能，以快速切換驅動發射電波，應用M×M個驅動波束，以偵測上肢震顫的狀況。通常雷達微波束的指向與等相位面垂直，在均勻口徑的幅射情況下，3dB波束寬度的寬度因子k=0.886，則每個波束寬度：

其中：λ為發射波的波長；M為陣元個數(M=4；2×2個雷達微波感測器元件)； a為陣元間距；θ為掃描角度，亦即波束指向方向。

如圖3B，以及(11)式所示，本發明之平面式雷達微波陣列感測裝置310系統可在有限的範圍內增加陣元的間距a，以減少陣元數量來達到相同的口徑降低成本的目標，並集中雷達微波束進行手部震顫的偵測，再以馬達控制集中雷達波束指向的範圍。

而如圖3C，係本發明之平面式雷達微波陣列感測裝置310系統的每個單一的雷達微波感測器元件311，其可發射單一雷達波束，包括有背後瓣波(Back lobe)31，主要瓣波(Main lobe)32，和兩側瓣波(Side lobe)33。

猶如圖3A，圖3B，以及圖3C所示，本發明係以單點雷達微波擴充成陣列式雷達微波束，並藉由可動機構改變陣列波束的掃描角度，以微控制器切換波束，進而增加其空間中的掃描範圍，對於移動物體更能夠精準偵測目標物的動態等相關參數。

如圖4所示，本發明之一種「帕金森氏症上肢震顫分析」方法，亦係本發明之一種時域/頻域特徵參數萃取與訊號處理方法，其步驟包括了以下：首先根據混頻器輸出所產生的中頻訊號，由原始時域的反射中頻訊號，進行類比數位訊號轉換，亦及，首先擷取反射中頻訊號401，再進行訊號數位濾波處理402，接著進行訊號放大403，再執行時域特徵參數萃取404，繼續，當進行訊號時域分析405步驟時，係直接從時域的反射訊號，萃取零交越點(ZC)，振幅差(AD)總合，Willison振幅(WANP)，以及振盪訊號的平均能量(ME)等相關時域特徵參數，以進行上肢震顫量化與嚴重程度的評估。而隨著上肢震顫的嚴重程度增加，四個時域參數的物理變化數值，會嚴重程度增加而隨著增加，藉此量化上肢震顫的等級，分成正常(Normal)，原發性震顫(ET)，姿態性震顫(PT)，以及帕金森氏症(PD)所造成的震顫。

仍如圖4所示，或是執行頻域特徵參數萃取406，再進行訊號頻域分析407，即將時域的反射訊號，進行時域至頻域的轉換，以取得頻域的參數(能量幅度與頻率)，藉由高頻與低頻成份的能量幅度變化與頻率分佈範圍，藉此量化上肢震顫的等級，得以分成前述的正常，原發性震顫，姿態性震顫，及帕金森氏症所造成的震顫。而頻域參數會隨著上肢震顫的嚴重程度增加，會逐漸出現高頻成份，即能量幅度在高頻範圍會逐漸增加。

猶如圖4所示，將每個陣元所取得將時域的反射訊號，進行訊號映射轉換408，再進行單點灰階或是彩色的特徵參數轉換409，接續進行二維特徵圖樣分析(重組二維特徵圖樣)410，而訊號的轉換式如下：

其中於(12)式中，NB參數可改變轉換訊號的解析度(NB=8,10,12)；而於(12)式中，abs(．)為計算絕對值運算子；仍於(12)式中，max(．)為取最大值運算子；再於(12)式中，B_max為絕對值之最大值；猶於(12)式中，映射轉換訊號Φ的刻度介於0至2^NB-1之間，再將每個陣元所取得映射轉換訊號Φ，轉換成單點的灰階或彩色的特徵參數，亦即二維特徵圖樣分析。

圖5A係本發明一種「帕金森氏症上肢震顫分析」的初始狀態顯示圖，可看出一開始僅能顯示所擷取的波型而已。

圖5B係本發明一種「帕金森氏症上肢震顫分析」的結果顯示圖，其係藉由最後重組成M×M的二維特徵圖樣，即藉由二維特徵圖樣的灰階或彩色圖樣的變化，以顯示上肢震顫的量化分析，藉以如圖4步驟411的評估上肢震顫的嚴重程度411。

依照前述，本發明係一種陣列式雷達微波束之時域特徵參數於帕金森氏症的上肢震顫量化分析裝置，可藉由非侵入的感測方式，發展適用於帕金森氏症病患，於居家自我評估及照護使用的智能偵測裝置。且本發明一種帕金森氏症上肢震顫分析裝置，係以個人化為中心之智能感測，智慧篩檢，行動顯示，以及個人化決策支援等功能之照護使用裝置，以達到預防與早期診斷的目的。

此外，本發明一種帕金森氏症上肢震顫分析裝置具有上肢動態震顫量測，震顫程度量化，以及自動數值資料分析等預期性的功能。甚且，本發明一種帕金森氏症上肢震顫分析裝置係應用都卜勒雷達微波感測基本原理，進行帕金森氏症上肢震顫的偵測，而於感測時，正常受測者的上肢是固定不動的，且反射電波的頻率會和發射電波的頻率非常接近。

本發明係應用雷達微波感測裝置，以雷達微波偵測動態物體，而可偵測受測者之上肢動態的細微變化。本發明除可偵測移動物體之外，雷達微波亦可以測量物體的範圍和徑向速度。而本發明之優點，可不受環境因素(白晝夜、溫度、濕度)的影響，並藉由非接觸的感測方式，能夠檢測目標物體的動態等相關參數。

以上所述僅為本發明之較佳實施例而已，並非用以限定本發明之申請專利範圍；凡其它未脫離本發明所揭示之精神下所完成之等效改變或修飾，均應包含在下述之申請專利範圍內。