Ocena kory kostnej za pomocą ultradźwiękowych fal kierowanych: badanie odtwarzalności w zdrowej populacji

Podsumowanie

Tutaj szczegółowo prezentujemy protokół pomiarowy urządzenia ultradźwiękowego z dwukierunkową transmisją osiową (BDAT) i testujemy go w badaniu odtwarzalności, biorąc pod uwagę 14 zdrowych uczestników i 3 operatorów. Wiarygodność, mierzona za pomocą wewnątrzklasowych współczynników korelacji (ICC), była dobra lub doskonała dla czterech interesujących parametrów.

Streszczenie

Kruche złamania są nadal ogólnoświatowym obciążeniem zdrowotnym w kontekście starzenia się populacji. W szczególności oczekuje się, że w latach 2020-2050 globalna liczba złamań szyjki kości udowej podwoi się. Dlatego tak ważne jest wykrywanie pacjentów zagrożonych złamaniami kruchości w skali populacji. Obecnie złotym standardem jest podwójna absorpcjometria promieniowania rentgenowskiego (DXA), zapewniająca powierzchniową gęstość mineralną kości (aBMD). Urządzenia ultradźwiękowe, zwykle bardziej przenośne i tańsze niż urządzenia rentgenowskie, stanowią interesującą alternatywę DXA jako narzędzia przesiewowe. Jednak zależność od operatora jest zwykle uznawana za ich główną wadę. W niniejszej pracy szczegółowo przedstawiono protokół pomiarowy urządzenia ultradźwiękowego z dwukierunkową transmisją osiową (BDAT). Dedykowana sonda ultradźwiękowa jest umieszczana w jednej trzeciej dystalnego promienia niedominującego przedramienia za pomocą konwencjonalnego żelu sprzęgającego. Interfejs prowadzony dostarcza w czasie quasi-rzeczywistym (około 2 Hz) cztery interesujące parametry: prędkości pierwszego przychodzącego sygnału (vFAS) i tryb A0 (vA0), grubość kory mózgowej (Ct.Th) i porowatość (Ct.Po), a także cztery parametry jakościowe. Operator powoli przesuwa sondę w miejscu pomiaru, uważnie obserwując informacje zwrotne dostarczane przez interfejs, aż do znalezienia stabilnej pozycji i rozpoczęcia serii 10 akwizycji. Po uzyskaniu co najmniej czterech spójnych serii pomiar kończy się i generowany jest automatyczny raport. Pomiar trwa zwykle około 5 minut. Aby określić solidność tego protokołu, przeprowadzono badanie odtwarzalności wśród 3 operatorów (jednego eksperta i dwóch nowicjuszy) i 14 zdrowych uczestników (6 kobiet, 8 mężczyzn, 21-53 lata). Współczynniki korelacji wewnątrzklasowej (ICC) uznano za dobre dla vA0 (0,76), Ct.Po (0,80) lub doskonałe dla Ct.Th (0,87) i vFAS (0,91). Stwierdzono, że odchylenia standardowe wynoszą mniej niż 10% całkowitych zakresów w praktyce klinicznej.

Wprowadzenie

Osteoporoza i związane z nią złamania nadal stanowią poważny problem zdrowia publicznego¹. W szczególności oczekuje się, że do 2050 r. liczba złamań szyjki kości udowej na świecie podwoi się². Łamliwość kości jest spowodowana powolnym i cichym procesem demineralizacji i utraty masy kostnej bez większych znaków ostrzegawczych przed złamaniem łamliwości. Obecnie złotym standardem w wykrywaniu pacjentów zagrożonych złamaniami łamliwości jest podwójna absorpcjometria rentgenowska (DXA), zapewniająca obraz rentgenowski 2D o niskiej rozdzielczości ze skalibrowanym szarym pikselem³. Na podstawie tego obrazu możliwe jest wyodrębnienie powierzchniowej gęstości mineralnej kości (aBMD w ^g.cm-2) w różnych obszarach zainteresowania związanych z głównymi miejscami złamań kruchości: kręgosłupa, nadgarstka i biodra. Wartość aBMD zmniejsza się wraz ze wzrostem szybkości pękania kruchości³. Co więcej, normalizacja T-score, w odniesieniu do normalnej zdrowej populacji, pozwala na porównanie pacjentów mierzonych za pomocą urządzeń proponowanych przez różnych producentów. Światowa Organizacja Zdrowia zaproponowała DXA T-score w celu zdefiniowania diagnostyki osteoporozy w trzech etapach: normalnym (T-score < -1), osteopenicznym (-1 < T-score < -2,5) i osteoporotycznym (T-score < -2,5)⁴.

DXA ma kilka ograniczeń: jego wielkość, stosunkowo wysoki koszt, potrzeba dedykowanego pomieszczenia i zdolność do rozróżniania między pękniętymi i niezłamanymi, a także dostępność w wielu krajach, takich jak Ameryka Łacińska, są umiarkowane⁵. W związku z tym istnieje zapotrzebowanie na alternatywy DXA jako narzędzia przesiewowe do szacowania ryzyka pęknięć kruchości⁶. Jednak niektóre alternatywy DXA, takie jak ilościowa tomografia komputerowa i jej pochodne⁷, obrazowanie metodą rezonansu magnetycznego (MRI)⁸, są również nieporęczne i nie są powszechnie dostępne. Ultrasonografia ilościowa (QUS) stwarza potencjał dla przenośnych, wytrzymałych i łatwych w użyciu urządzeń do badań przesiewowych. Opracowano różne urządzenia do oceny kości korowej, związane z różnymi częstotliwościami w zakresie od kilku kHz do kilku MHz i różnym położeniem przetworników w transmisji, dyfuzją retro ⁹, echem impulsowym¹⁰ oraz transmisją osiową, w której przetworniki są ustawione w jednej linii z osią kości długiej, takiej jak promień i kość piszczelowa. Niektóre urządzenia zapewniają substytuty BMD¹¹, podczas gdy inne zapewniają "klasyczne" parametry ultradźwiękowe, takie jak prędkości¹² lub współczynniki tłumienia⁹, a nawet parametry geometryczne i materiałowe, na przykład grubość kory mózgowej, porowatość lub rozkład wielkości porów⁹. Jednak do dnia dzisiejszego QUS nie udało się jeszcze znaleźć szerokiego zastosowania w praktyce klinicznej do oceny kości, częściowo ze względu na brak homogenizacji między urządzeniami i zależność od operatora¹³.

Wśród technologii QUS proponowanych jako alternatywy dla DXA, transmisja osiowa (AT) ma tę zaletę, że pomiar może być wykonany na przedramieniu, miejscu (i) łatwo dostępnym i (ii) blisko jednego z głównych miejsc złamań kruchości, tj. nadgarstka. Pierwszy proponowany parametr AT zależy od prędkości rozchodzenia się ultradźwięków w warstwie korowej, oznaczanej prędkością dźwięku (SOS) lub prędkością sygnału pierwszego nadejścia (vFAS), w zależności od przetwarzania sygnału i urządzeń, niektóre z nich są komercyjne^12,14 i inne prototypy laboratoryjne^15,16. Ten parametr był w stanie rozróżnić grupy pacjentów ze złamaniami kruchości lub bez nich o wydajności podobnej do BMD w kilku badaniach klinicznych od końca lat 90. XX wieku^14,15. Został również z powodzeniem zastosowany w wieloośrodkowych badaniach podłużnych, wykazując swoje zastosowanie kliniczne i solidność¹². Precyzja vFAS została ulepszona poprzez połączenie dwóch przeciwstawnych kierunków propagacji w celu zmniejszenia odchylenia wynikającego z kąta między sondą a powierzchnią kości^16,17. Ten punkt widzenia został oznaczony jako dwukierunkowy AT (BDAT).

Nawet jeśli vFAS okazał kliniczne zainteresowanie, jego główną wadą, podobnie jak BMD, jest to, że łączy w sobie różne kluczowe cechy kości korowej, takie jak właściwości geometryczne i materiałowe, co sprawia, że jego kliniczna interpretacja nie jest prosta. Dlatego zaproponowano punkt widzenia fali prowadzonej, biorąc pod uwagę jej potencjał wynikający z doskonałej wrażliwości fal kierowanych na właściwości falowodu. Podejście to powinno łączyć przetwarzanie sygnału, modelowanie falowodu i problemy odwrotne i jest szeroko stosowane w badaniach nieniszczących, biorąc pod uwagę na przykład falowody metalowe, takie jak płytki lub rurki¹⁸. W związku z tym od 2010 roku krok po kroku opracowywane jest urządzenie BDAT drugiej generacji, od fantomów naśladujących kości¹⁹ do walidacji ex vivo²⁰ i pomiarów in vivo²¹. Urządzenie zostało pomyślnie przetestowane w badaniach klinicznych we Francji²², Niemczech²³, Wielkiej Brytanii²⁴ i Chile²⁵, i wykazało poprawę wyników pod względem wskaźnika sukcesu i dyskryminacji pacjentów.

To badanie ma na celu zbadanie odtwarzalności obecnego urządzenia ultradźwiękowego BDAT. Po pierwsze, urządzenie i protokół pomiarowy zostaną szczegółowo opisane. Wyniki uzyskane z udziałem 14 uczestników i 3 operatorów zostaną przedstawione i omówione w kontekście populacyjnych badań przesiewowych w celu wykrycia pacjentów zagrożonych złamaniami kruchości.

Zasada pomiaru: przetwarzanie sygnału, interesujące parametry i parametry jakościowe
Dwukierunkowe urządzenie do transmisji osiowej (BDAT) składa się z różnych części, z których główną jest sonda ultradźwiękowa, moduł elektroniczny i komputer. Pełna lista jest szczegółowo opisana w Tabeli Materiałów i zilustrowana w Rysunek 1. Poniżej opisano interesujące nas parametry, parametry jakości pomiaru oraz protokół pomiaru.

vFAS
Po odebraniu próbkowanych sygnałów przez komputer, są one przetwarzane w różnych krokach. Pierwszy krok polega na przetwarzaniu sygnału w dziedzinie czasu, wykrywaniu FAS za pomocą protokołu opisanego wcześniej^16,17. Po uzyskaniu czasu przybycia dla każdego odbiornika możliwe jest wyznaczenie prędkości FAS, później oznaczonej vFAS, która jest średnią harmoniczną prędkości uzyskanych w obu kierunkach rozchodzenia się. Łącząc informacje z obu kierunków propagacji, możliwe jest uzyskanie kąta wartości między kierunkiem sondy a kierunkiem powierzchni kości oraz wyprowadzenie bezstronnej wartości vFAS¹⁶. Ten dwukierunkowy kąt jest później oznaczany jako alfa i jest używany jako parametr jakości pomiaru. To przetwarzanie czasowe pozwala również na oszacowanie grubości tkanki miękkiej między powierzchnią kości a sondą, oznaczoną jako ST.Th²⁶.

Obraz widma fal kierowanych
Drugi etap polega na przetwarzaniu sygnału w domenie Fouriera, z uwzględnieniem częstotliwości czasowych i przestrzennych, oznaczanych f i k. Podejście to jest metodą opartą na SVD, pozwalającą na transformację sygnałów czasoprzestrzennych w funkcję Norm, oznaczoną również jako obraz widma fal kierowanych (GWSI), jak pokazano w Rysunek 2 dla in vivo forearm¹⁹. Metoda łączy w sobie dwie transformaty Fouriera (czas i przestrzeń) oraz dekompozycję wartości osobliwych (SVD), co pozwala na wizualizację szybkości obecności w odbieranych sygnałach (w skali 0-1) modów kierowanych przez warstwę kości korowej. GWSI może być interpretowany jako ulepszenie czasoprzestrzennej transformaty Fouriera, w której każdy piksel jest powiązany z niezależną płaszczyzną częstotliwości f i liczby falowej k. Należy zauważyć, że podejście zostało ulepszone w celu uwzględnienia wpływu tłumienia materiału²⁷ i liniowej zmienności grubości²⁸.

Szczególna uwaga zostanie zwrócona na górną część spektrum, związaną z trybem A0, a także na najniższą część, związaną z najwyższymi wartościami prędkości fazowej, tj. większą niż 4 mm·^μs-1. Ta część odpowiada obszarowi zainteresowania 3 (ROI 3)²⁹. Średnia wartość ROI 3, później oznaczona jako lowk, jest również używana jako parametr jakości. Duża wartość odpowiada zwykłemu falowodom, co pozwala na wyraźne odbicia fal na granicy kości. Jeśli wartość spadnie, może to być spowodowane nieregularnym falowodem lub źle umieszczoną sondą.

Model falowodu
Dyspersja fali prowadzonej, czyli zmiana prędkości fazowej każdego trybu prowadzonego w odniesieniu do częstotliwości, zależy zarówno od właściwości materiału, jak i właściwości geometrycznych falowodu. W związku z tym potencjalnie możliwe jest odzyskanie tych właściwości za pomocą dedykowanego przetwarzania sygnału, modelowania falowodu i odwrotnych schematów problemów. W przypadku BAT model falowodu odpowiada poprzecznej izotropowej płycie swobodnej 2D, w zależności od materiału falowodu i jednego parametru geometrycznego, thickness³⁰. Materiał kostny kory mózgowej jest homogenizowany z uwzględnieniem ustalonych parametrów macierzy kostnej i zmiennej porowatości³¹. Zatem problem odwrotny zależy od dwóch parametrów, oznaczanych grubością kory (Ct.Th) i porowatością kory mózgowej (Ct.Po). Wpływ absorpcji materiału, krzywizny falowodu i otaczających tkanek miękkich nie jest uwzględniany w modelu, nawet jeśli ma wpływ na pomiar. Jednak ich waga w odwrotnym wyniku problemu nie została uznana za determinującą, co oznacza, że mody w dwóch głównych obszarach zainteresowania (A0 i najniższa część) nie są znacząco zmieniane przez krzywiznę i tkanki miękkie³².

Odwrotny problem
Początkowo problem odwrotny został podzielony na dwa kroki: po pierwsze, wyodrębnić eksperymentalną dyspersję fal kierowanych, a po drugie, porównać z modelem falowodu. Ten punkt widzenia był ograniczony przez szum i etykietowanie trybów^30,32. W związku z tym zaproponowano specjalne podejście mające na celu przezwyciężenie tych ograniczeń jako rozszerzenie punktu widzenia funkcji normy. Zamiast rozpatrywać każdą falę płaską niezależnie, brane są pod uwagę tylko możliwe fale kierowane dostarczone przez model falowodu²⁰. Prowadzi to do odwrotnego obrazu problemu, wyrażonego w dziedzinie parametrów modelu, tj. płaszczyzny Ct.th - Ct.Po (Rysunek 2 w prawym dolnym rogu). Najlepiej dopasowanemu modelowi podaje się maksymalną pozycję, podczas gdy ostateczne wtórne piki (oznaczone odwróconymi obrazami problemu z szarą kropką) odpowiadają niejednoznacznym rozwiązaniom, wskazanym w porównaniu f-k z modami eksperymentalnymi z jasnoszarymi liniami. Tak jak poprzednio, wartość piksela jest znormalizowana konstrukcyjnie i odzwierciedla, w tym przypadku, obecność jednego konkretnego modelu falowodu w odbieranych sygnałach. Wartość maksymalna (oznaczona jako maksimum) i różnica w stosunku do drugiego maksimum (oznaczona jako różnica) są również używane jako parametry jakościowe.

Odwrotny problem został pierwotnie zaproponowany do obliczeń offline, tj. po pobraniu sygnałów, przy użyciu dokładnych wartości numerów fal modelu. To podejście zostało zweryfikowane zarówno dla stron radius jak i tibia, biorąc pod uwagę ex vivo^20,33 oraz in vivo^21,34,35 badań. Aby uwzględnić te obliczenia w interfejsie człowiek-maszyna (HMI), zaproponowano przybliżoną wersję, kompatybilną z aplikacją czasu rzeczywistego, wykorzystującą punkt widzenia rzadkiej macierzy³⁶.

vA0
Z GWSI możliwe jest również wyodrębnienie prędkości najwolniejszego trybu prowadzonego, powiązanego z pierwszym trybem antysymetrycznym A0 wolnej płyty lub modelu Lamba^33,35. Górna część widma fal kierowanych może być aproksymowana liniowo, przy czym nachylenie zapewnia wartość prędkości vA0 (rysunek 2 na dole po lewej).

Podsumowanie parametrów:
Na koniec mierzone są cztery interesujące parametry: (i) vFAS: prędkość pierwszego przychodzącego sygnału (m·^s-1); (ii) vA0: prędkość najwolniejszego trybu kierowania (m·^s-1); (iii) Ct.Th: grubość kory mózgowej (mm); oraz (iv) Ct.Po: porowatość kory mózgowej (%).

Pod uwagę brane są cztery parametry jakości: (i) alfa: kąt dwukierunkowy (°); (ii) lowk: średnia wartość najniższej części GWSI (wartość znormalizowana między 0 a 1); (iii) max: maksimum odwrotnej funkcji problemu (wartość znormalizowana między 0 a 1); oraz (iv) diff: różnica między pierwszym a drugim maksimum odwrotnej funkcji problemu (wartość znormalizowana z zakresu od 0 do 100).

Wszystkie te parametry, jak również dwa obrazy widma fal kierowanych (jeden par kierunek rozchodzenia się) i obraz odwrotnego problemu, są wyświetlane w "czasie rzeczywistym" przez HMI, z częstotliwością odświeżania około 2 Hz. Typowy przykład przedstawiono na rysunku 3. W dalszej sekcji szczegółowo opisano metodę korzystania z tych parametrów. Główną ideą jest to, że operator powoli przesuwa sondę w miejscu pomiaru, uważnie obserwując informacje zwrotne dostarczane przez różne części interfejsu, aż do znalezienia stabilnej pozycji i rozpoczęcia serii 10 akwizycji. Po uzyskaniu co najmniej czterech spójnych serii pomiar kończy się i generowany jest automatyczny raport.

Protokół

Badanie zostało zatwierdzone przez Komitet Etyczny Universidad de Valparaiso, Chile, pod numerem protokołu CEC213-20. Uczestnikom udzielono pisemnej świadomej zgody. Przeprowadzono wywiad telefoniczny w celu ustalenia włączenia/wykluczenia uczestników. Badanie zostało zarejestrowane pod następującym numerem referencyjnym: NCT05424536.

1. Konfiguracja urządzenia

1. Umieść główne części urządzenia na dużym stole.
   1. Na dużym stole ułóż obok siebie następujące części: transformator izolacji elektrycznej, moduł elektroniczny i laptop. Upewnij się, że przed tymi częściami jest wystarczająco dużo miejsca, aby później z łatwością ułożyć przedramię uczestnika.
   2. Ostatecznie umieść laptopa bezpośrednio na module elektronicznym w przypadku ograniczonej przestrzeni, uważając, aby nie zasłaniać okien wentylacyjnych modułu, które są wyraźnie wskazane.
2. Podłączyć transformator izolacji elektrycznej.
   1. Podłącz transformator izolacji elektrycznej do domowego źródła zasilania w pomieszczeniu za pomocą dedykowanego.
3. Podłączyć moduł elektroniczny.
   1. Podłącz moduł elektroniczny do transformatora izolacji elektrycznej za pomocą dedykowanego zasilającego.
   2. Naciśnij przycisk ON-OFF na transformatorze, aby zasilić moduł.
4. Podłącz laptopa.
   1. Podłącz laptopa do modułu za pomocą dedykowanego uniwersalnej magistrali szeregowej (USB), aby przesłać zdigitalizowane odebrane sygnały do komputera w celu dalszego przetwarzania.
   2. Jeśli laptop wymaga zasilania, podłącz jego zasilający do elektronicznego transformatora izolacyjnego.
5. Podłącz sondę ultradźwiękową.
   1. Podłącz sondę ultradźwiękową do modułu za pomocą dedykowanego otworu na znajdującego się z przodu modułu. Istnieją dwie różne sondy do pomiaru miejsc pomiaru przedramienia i nogi. W tym badaniu brana jest pod uwagę tylko promień (sonda przedramienia).
6. Podłącz przełącznik pedału.
   1. Umieść przełącznik pedału na podłodze blisko stóp, biorąc pod uwagę pozycję operatora podczas pomiaru uczestnika. Podłącz przełącznik pedału do komputera za pomocą USB. Użyj pedału, aby rozpocząć serię akwizycji.

2. Instalacja uczestnika

1. Ustaw uczestnika.
   1. Poproś uczestnika, aby usiadł przed operatorem z gołym przedramieniem opartym na stole przed wcześniej zainstalowanym urządzeniem (patrz Rysunek 3).
      UWAGA: Mierzona jest strona kontralateralna (tj. lewa strona dla uczestnika praworęcznego).
2. Zaznacz miejsce pomiaru (jedna trzecia dystalnego promienia).
   1. Zmierz długość promienia za pomocą linijki od promieniowego rylcowatego rylca (koniec kości blisko nadgarstka) do łokcia.
   2. Podziel tę długość przez trzy.
   3. Zaznacz miejsce pomiaru, tj. jedną trzecią dystalnego promienia, za pomocą pisaka, mierząc jedną długość promienia na trzecią część nadgarstka.
3. Uruchom oprogramowanie HMI.
   1. Uruchom oprogramowanie HMI, klikając odpowiednią ikonę na pulpicie laptopa.
4. Dodaj dane uczestnika.
   1. Dodaj dane uczestnika (zanonimizowany identyfikator, lateralność, zmierzone miejsce, identyfikator operatora, płeć itp.) za pomocą wyskakującego okienka, które otwiera się automatycznie po uruchomieniu oprogramowania.
5. Dodaj żel ekologiczny.
   1. Dodaj żel ekologiczny z przodu sondy i w miejscu pomiaru, oznaczonym na przedramieniu uczestnika, aby zapewnić propagację fali ultradźwiękowej.
6. Umieść sondę w kontakcie z przedramieniem.
   1. Umieść sondę w kontakcie z przedramieniem, tak aby środek sondy znajdował się na znaku, który wykonano wcześniej w kroku 2.2.

3. Szukam stabilnej pozycji

UWAGA: HMI wyświetla cztery interesujące parametry: dwie prędkości, vFAS i vA0, oraz dwie odwrotne wartości problemu, grubość kory (Ct.Th) i porowatość kory (Ct.Po). HMI wyświetla również cztery parametry jakościowe oznaczone jako alpha, lowk, max i diff. Parametry te zostały szczegółowo opisane we wprowadzeniu.

1. Rozpocznij wizualizację w czasie rzeczywistym.
   1. Rozpocznij wizualizację w czasie rzeczywistym, klikając przycisk Start w prawym dolnym rogu interfejsu oprogramowania. Czas między dwoma kolejnymi wyświetleniami wartości wynosi około 0,5 s.
2. Znajdź stabilną wartość vFAS.
   1. Powoli wyreguluj pozycję sondy, obserwując wartość parametru vFAS wyświetlaną w konkretnym przypadku interfejsu. Normalne wartości wahają się od około 3800 m∙^s-1 do około 4200 m∙^s-1.
   2. Jeśli zostanie znaleziona stabilna pozycja, upewnij się, że odchylenie vFAS jest mniejsze niż około 40 m∙^s-1 między dwoma kolejnymi obliczeniami.
3. Dostosuj kąt dwukierunkowy.
   1. Powoli wyreguluj pozycję sondy, obserwując wartość dwukierunkową (parametr jakości alfa) wyświetlaną w konkretnym przypadku interfejsu.
   2. Dostosuj pozycję sondy, delikatnie zwiększając nacisk po jednej stronie sondy, aż wartość bezwzględna kąta będzie mniejsza niż 2°, aby poprawić równoległość między sondą a powierzchnią kości.
4. Znajdź stabilną wartość vA0.
   1. Powoli wyreguluj pozycję sondy, obserwując wartość parametru vA0 wyświetlaną w określonym przypadku interfejsu. Normalne wartości wahają się od około 1500 m∙^s-1 do około 1900 m∙^s-1.
   2. Jeśli zostanie znaleziona stabilna pozycja, należy upewnić się, że odchylenie vA0 jest mniejsze niż około 40 m∙^s-1 między dwoma kolejnymi obliczeniami.
   3. W przypadku trudności obserwuj widma obrazu fali kierowanej wyświetlane w prawej kolumnie interfejsu. Upewnij się, że górna część widma jest widoczna jako linia ciągła, której nachylenie zapewnia wartość vA0.
5. Obserwuj odwrotny obraz problemu.
   1. Obserwuj odwrotny obraz problemu, który pojawia się automatycznie po ustabilizowaniu dwóch prędkości (vFAS i vA0) oraz wartości kąta.
   2. Upewnij się, że obraz pokazuje co najmniej jedno maksimum, oznaczone wyraźnym pikselem, a ostatecznie jedno lub kilka drugorzędnych maksimów, oznaczonych innym kolorem. Trzy brakujące parametry jakościowe (max. diff, lowk) są automatycznie obliczane w czasie rzeczywistym.
6. Popraw obraz problemu odwrotnego.
   1. Powoli dostosuj pozycję sondy, obserwując odwrotne maksima obrazu.
   2. Znajdź najwyższe możliwe pierwsze maksimum i najniższe możliwe maksimum pomocnicze, patrząc na odpowiednie przypadki interfejsu (wartości max i diff).
   3. W przypadku trudności obserwuj obraz widma fal prowadzonych wyświetlany w prawej kolumnie interfejsu. Upewnij się, że dolna część widma jest pokryta kilkoma ciągłymi liniami, tak długimi, jak to możliwe, związanymi z modami wysokiej prędkości fazowej i niską jakością parametru, tak wysokim, jak to możliwe.
7. Znajdź stabilną pozycję.
   1. Po znalezieniu akceptowalnego obrazu problemu odwrotnego ustabilizuj pozycję sondy. Upewnij się, że między dwoma kolejnymi obliczeniami nie są widoczne żadne znaczące zmiany obrazu odwrotnego problemu.

4. Akwizycja danych

1. Rozpocznij serię 10 przejęć.
   1. Po znalezieniu stabilnej pozycji rozpocznij serię 10 akwizycji, naciskając stopą przełącznik pedału.
   2. Zachowaj jak największą stabilność podczas 10 przejęć, trwających około 5 sekund.
2. Kontroluj jakość serii.
   1. Przyjrzyj się średnim i odchyleniom standardowym interesujących Cię parametrów, które są automatycznie obliczane i wyświetlane w wyskakujących okienkach pojawiających się po zakończeniu serii.
   2. Jeżeli odchylenia standardowe są niższe niż ustalone progi, należy wziąć pod uwagę szeregi. Wręcz przeciwnie, odrzuć serię.
   3. Odpowiedź na pytanie zadane w drugim wyskakującym okienku, z pytaniem, czy operator chciałby przerwać, czy kontynuować serię akwizycji tego samego uczestnika.
3. Zmień położenie sondy.
   1. Rozpocznij ponownie poprzednie kroki (od kroków 2.1 do 3.2), aby znaleźć bardziej stabilne pozycje i zdobyć więcej serii 10 przejęć. Na koniec, jeśli to konieczne, pozwól uczestnikowi odpocząć między dwoma zmianami pozycji.
   2. Tak jak poprzednio, średnie i odchylenia standardowe parametrów będących przedmiotem zainteresowania są obliczane automatycznie dla każdej serii.
   3. Spójrz na wyskakujące okienko wyników, aby sprawdzić, czy ostatnia nabyta seria została zachowana, czy odrzucona. Pomiar uczestnika kończy się, gdy zarejestrowane zostaną co najmniej cztery spójne serie. Serie outsiderów są automatycznie odrzucane.
4. Uzyskaj wartości końcowe.
   UWAGA: Końcowe wartości interesujących nas parametrów są uzyskiwane automatycznie, biorąc pod uwagę średnią średnich uzyskanych za pomocą spójnego szeregu.
5. Sprawdź plik PDF z raportem automatycznym.
   1. Sprawdź, czy wartości końcowe są zgłaszane w pliku PDF raportu automatycznie i natychmiast generowane po wybraniu opcji zatrzymania w wyskakującym okienku. Przykład jest pokazany w Rysunek 4.
      UWAGA: Plik PDF znajduje się w tym samym folderze, co nieprzetworzone dane, które mogą być później ponownie analizowane w trybie offline.
6. Zweryfikuj drugi dokładny raport.
   1. Zweryfikuj drugi dokładny raport wygenerowany przy użyciu dokładnych wartości modelu falowodu do obliczenia problemu odwrotnego, a nie wartości przybliżonych, jak w przypadku pierwszego raportu automatycznego. Drugie wygenerowanie raportu zajmuje mniej niż 5 minut. Przykłady są pokazane w Rysunek 5 i Rysunek 6.
   2. Sprawdź, czy raport automatyczny zgadza się z dokładnym raportem. Usuń serie, które nie zostały automatycznie wyeliminowane, aby zachować spójność serii.

Wyniki

Przeprowadzono badanie odtwarzalności z udziałem 3 operatorów (jeden ekspert, dwóch nowicjuszy) i 14 zdrowych uczestników (6 kobiet, 8 mężczyzn, 21-53 lata). Początkujący operatorzy byli szkoleni przez około 3 godziny, aby zrozumieć i przećwiczyć protokół akwizycji. Następnie uczestnicy byli mierzeni przez 2 tygodnie w sierpniu 2023 roku. Każdy pomiar był wykonywany niezależnie. Wszyscy operatorzy byli zaślepieni, tj. jeden operator nie znał wyników uzyskanych przez dwóch pozostałych.

Powtarzalność wewnątrzoperatorska
Rysunek 5 ilustruje powtarzalność wewnątrzoperatorską dla spójnego przypadku: 4 serie po 10 przejęć nabytych na jednym uczestniku i jednym operatorze. Pierwsza kolumna odpowiada odwrotnym funkcjom problemu, podczas gdy druga kolumna pokazuje eksperymentalną dyspersję modów kierowanych w porównaniu z najlepiej dopasowanym modelem. Każda linia figury odpowiada udanej serii. Liczba zachowanych nabytków podana jest w tytule: 7 na 10 dla pierwszej serii i 10 na 10 dla kolejnych. Średnie i odchylenia standardowe oraz cztery parametry (vFAS, vA0, Ct.Th i Ct.Po) są podane dla każdej serii. Dodatkowo pokazane są również parametry jakościowe: kąt dwukierunkowy (alfa), maksimum funkcji odwrotnego problemu (max), różnica bezwzględna z drugim maksimum (diff) oraz średnia wartość najniższej części GWSI (lowk).

Odchylenia standardowe wewnątrz serii są niskie, około 0,02 mm dla grubości kory, mniej niż 0,5% dla porowatości kory i mniej niż 20 m∙^s-1 dla dwóch prędkości, co wskazuje, że pozycje sondy zasięgu były stabilne. Następnie można zauważyć, że średnie wartości uzyskane dla każdej serii są bardzo zbliżone, szczególnie dla wartości grubości w zakresie od 3,4 mm do 3,5 mm oraz wartości vFAS w zakresie od 4040 m∙^s-1 do 4070 m∙^s-1. Należy zauważyć, że różnica 40 m∙^s-1 odpowiada różnicy 1% dla średniej wartości 4000 m∙^s-1. Większe wahania obserwuje się dla porowatości kory mózgowej, wynoszącej od 8% do 12%, oraz prędkości vA0, wynoszącej od 1700 m∙^s-1 do 1740 m∙^s-1. W tym konsekwentnym przypadku prawie wszystkie akwizycje są spójne, czyli blisko siebie. Nie ma prawie żadnych niejasności w końcowym wyniku czterech interesujących nas parametrów.

Drugi przypadek jest zilustrowany w Rysunek 6 dla mniej regularnego przypadku. W tym przypadku trzy parametry, vFAS, vA0 i Ct.Po, są bardzo stabilne, z wartościami odpowiednio około 4120 m∙^s-1, 1700 m∙^{s-1 i} 10%. Przypadek grubości kory mózgowej jest trudniejszy, ponieważ obserwuje się dwa niejednoznaczne rozwiązania, 2,6 mm i 3,5 mm, zgodne z niewielką wartością różnicy (różnica mniejsza niż 0,5%) między dwoma pierwszymi maksimami odwrotnej funkcji problemu. W poprzednim regularnym przypadku różnica ta wynosiła od 1 % do 3 %. Niejednoznaczność jest usuwana przez analizę ekspercką, w tym przypadku przyglądającą się zgodności między eksperymentalnymi i teoretycznymi trybami kierowanymi (prawa kolumna). W przypadku najmniejszej grubości zgodność jest lepsza w najniższej części widma (dwie pierwsze serie). Dla dwóch ostatnich serii istnieje tryb teoretyczny z bardzo małą liczbą punktów eksperymentalnych, około 0,5 MHz, co wskazuje na gorszą zgodność w porównaniu z poprzednią serią. Co więcej, parametr diff (0,1%) jest mniejszy niż wartości dwóch pierwszych serii (0,4% i 0,2%). W takim przypadku wybór zachowanego rozwiązania (2,6 mm) nie jest jeszcze zautomatyzowany i nadal potrzebny jest ekspert. Trzej operatorzy borykali się jednak z podobnymi problemami i zdecydowali się na podobne rozwiązania, bliskie 2,6 mm.

Niezawodność poza operatorem
Wszystkie wyniki dla 4 interesujących parametrów uzyskanych przez 3 operatorów z 14 uczestnikami są pokazane w Rysunek 7. Współczynniki korelacji wewnątrzklasowej (ICC) zostały obliczone zgodnie ze wzorem i wcześniej opublikowanym kodem Matlaba^37,38,39. ICC jest powszechnie stosowany do oceny wiarygodności skal pomiarowych, w szczególności do zastosowań biomedycznych. Wartości ICC w zakresie od 0,75 do 0,9 są zwykle związane z dobrą niezawodnością, podczas gdy wartości ICC powyżej 0,9 są uważane za doskonałą niezawodność. Najniższą wartość dla ICC (0,76) uzyskano dla parametru vA0. Odchylenie standardowe było równe 17 m∙^s-1, co stanowi około 7% zakresu pomiarowego rzędu 250 m∙^s-1. Podobne wartości zaobserwowano dla Ct.Po z ICC równym 0,80 i odchyleniem standardowym 1,1%, około 10% zakresu. Doskonałą rzetelność (ICC około 0,9) uzyskano dla dwóch pozostałych parametrów, Ct.Th i Ct.Po, przy odchyleniu standardowym mniejszym niż 10% zakresu.

Rysunek 1: różne części urządzenia ultradźwiękowego Bi-Directional Axial Transmission (BDAT). Prototyp składa się z izolacji elektronicznej (1), przełącznika nożnego (2), dwóch sond (3.1 i 3.2), modułu elektronicznego (4), komputera (5), oraz linijki (6). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Typowa akwizycja na przedramieniu in vivo. Dwie funkcje Norm (po jednej na kierunek propagacji) są pokazane w górnym wierszu obrazu. Są one również oznaczone jako obraz widma fal prowadzonych (GWSI). Z maksimów tych obrazów możliwe jest wyodrębnienie eksperymentalnych trybów kierowanych (niebieskie i czerwone kropki) w porównaniu z najlepiej dopasowanym modelem (lewy dolny obraz podrzędny). Najlepiej dopasowany model jest parametryzowany przez dwie wartości, grubość kory mózgowej (Ct.Th) i porowatość (Ct.Po), odpowiadające maksymalnej pozycji funkcji problemu odwrotnego (prawy dolny obraz podrzędny). Ich wartości są pokazane w tytule każdego panelu. Łącznik vA0 jest pokazany jako linia przerywana (po lewej). Obrazy widma fal kierowanych i obraz problemu odwrotnego są znormalizowane (tj. wartość piksela mieści się w zakresie od 0 do 1) przez construction^19,20. Wartości dwóch zmierzonych prędkości, vFAS i vA0, są wskazane w tytule prawego górnego obrazu podrzędnego. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Interfejs człowiek-maszyna (HMI). HMI pokazuje w rzeczywistości dwa GWSI (po jednym na kierunek propagacji), odwrotny obraz problemu, interesujące parametry i parametry jakościowe. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Przykład automatycznego raportu. W raporcie wskazano dane uczestnika i operatora oraz końcowe wartości interesujących nas parametrów i parametrów jakościowych. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5: Przykład drugiego raportu dla spójnego przypadku, taki sam jak raport automatyczny pokazany w Rysunek 4. Rysunek przedstawia 4 serie po 1 uczestnika i 1 operatora: odwrotne obrazy problemów (lewa kolumna) i eksperymentalne liczby falowe w porównaniu z najlepiej dopasowanym modelem (prawa kolumna). Przedstawione wartości odpowiadają średniej i odchyleniu standardowemu utrzymywanych przejęć w ciągu 10 serii przejęć. Liczba przechowywanych nabytków jest wskazana w tytule prawej kolumny, na przykład (7/10) dla pierwszej serii. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Przykład drugiego raportu dla niejednoznacznego przypadku grubości kory mózgowej. Rysunek przedstawia 4 serie po 1 uczestnika i 1 operatora: odwrotne obrazy problemów (lewa kolumna) i eksperymentalne liczby falowe w porównaniu z najlepiej dopasowanym modelem (prawa kolumna). Przedstawione wartości odpowiadają średniej i odchyleniu standardowemu w stosunku do utrzymywanych przejęć w ciągu 10 serii akwizycji. Liczba przechowywanych nabytków jest wskazana w tytule prawej kolumny, na przykład (10/10) dla pierwszej serii. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 7: Niezawodność między operatorami. Wyniki przedstawiono dla 4 interesujących parametrów uzyskanych przez 3 operatorów z 14 uczestnikami. Wartości uzyskane przez 3 operatory (oś y) są porównywane ze średnią wartością 3 operatorów (oś x). Współczynniki korelacji wewnątrzklasowej (ICC) są wskazane w tytułach. Wskazane jest również σ odchylenia standardowego i współczynnik zmienności CV. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Dyskusja

Punktem krytycznym dla powodzenia pomiaru jest prawidłowe ustawienie sondy w stosunku do kości. Położenie względem powierzchni kości udało się rozwiązać poprzez zintegrowanie z interfejsem prowadzącym dwukierunkowej korekcji^16,17. Wykazano, że bez tej poprawki można osiągnąć kilkuprocentowy błąd dla parametrów vFAS¹⁶. Poprawka ta okazała się kluczowa dla dokładnego rozróżnienia pacjentów, ponieważ cały zakres pomiarowy (około 3800-4200 m∙^s-1) wynosi około 10% średniej wartości¹⁷, około 4000 m∙^s-1. Odtwarzalność stwierdzona w tym badaniu (odchylenie standardowe 24 m∙^s-1) była podobna do poprzedniego badania¹⁷, wskazując na odchylenie standardowe 20 m∙^s-1 i współczynnik zmienności (CV) około 0,5%. Podobny współczynnik zmienności (0,5 %) stwierdzono przy niższej częstotliwości, tj. 0,3 MHz⁴⁰. Należy zauważyć, że doskonała wiarygodność parametru vFAS (ICC = 0,91) okazała się wyższa niż ostatnio uzyskana za pomocą innego urządzenia AT w starszej populacji (ICC = 0,77)⁴¹. Należy zauważyć, że populacja w tym przypadku składała się nie tylko ze osób starszych, ale także samych kobiet i niektórych osób z innymi objawami klinicznymi zwiększonej porowatości kości⁴¹.

Drugim wyzwaniem jest również prawidłowe ustawienie sondy w stosunku do osi kostnej i dwóch interfejsów w celu uzyskania modów kierowanych, w szczególności tych o dużej prędkości fazowej lub niskiej liczbie falowej. Mody te są bliskie rezonansowi i związane z wieloma ścieżkami odbicia. Jeśli wyrównanie nie jest prawidłowe, fale zostaną rozproszone poza przestrzenny zakres pomiarowy sondy. Wręcz przeciwnie, jeśli sonda jest prawidłowo ustawiona, te tryby prędkości wysokofazowej pojawią się jako ciągłe gałęzie w ROI3. W przypadku pierwszej wersji obecnego interfejsu HMI to wyrównanie zostało rozwiązane poprzez przyjrzenie się GWSI²⁰ w czasie rzeczywistym. Podejście to okazało się jednak niewystarczające do osiągnięcia solidnego pomiaru: w pilotażowym badaniu klinicznym stwierdzono, że odsetek niepowodzeń wynosił około²⁰%. Włączenie parametrów jakościowych w czasie rzeczywistym pozwala na zmniejszenie liczby parametrów do około 10% w drugim badaniu klinicznym²³. Następnie, włączenie odwrotnego obrazu problemu do HMI "w czasie rzeczywistym" pozwala na dalszą poprawę, przy obecnym wskaźniku awaryjności dla pomiaru promienia wynoszącym około 5%²⁵. Należy zauważyć, że początkowa awaria z pierwszym urządzeniem mierzącym tylko vFAS wyniosła około 15%¹⁷.

Stwierdzono, że wiarygodność grubości kory mózgowej jest podobna do poprzedniego badania²². W poprzedniej wersji tego samego prototypu stwierdzono również, że ICC jest bliskie 0,9, z odchyleniem standardowym około 0,1 mm i CV około 3%. Obserwuje się jednak wyraźną poprawę porowatości kory mózgowej: ICC wzrosło z 0,6²² do 0,8, a odchylenie standardowe zmniejszyło się z 1,5% do 1%. Najtrudniejszym parametrem jest vA0 ze względu na bliskość jego prędkości do prędkości tkanek miękkich, wynoszącej około 1500 m∙^s-1. Jeśli sonda jest prawidłowo ustawiona, tryb A0 wydaje się unikalny i ciągły. Wręcz przeciwnie, wydaje się nieciągły i/lub wielokrotny z powodu sprzężenia z otaczającymi tkankami miękkimi. Efekt ten jest szczególnie silny u osób o małej (mniej niż 4 mm) lub dużej (ponad 10 mm) grubości tkanek miękkich (oznaczanej ST.Th).

Jak wspomniano wcześniej, kluczową kwestią jest znalezienie właściwej pozycji sondy. Pozycja powinna być jednak nie tylko poprawna, ale także stabilna, aby można było wyznaczyć średnią i odchylenie standardowe dla serii 10 przejęć. W większości przypadków znalezienie stajni zajmuje mniej niż 1 minutę, a pełny pomiar pacjenta trwałby około 5 minut. Nawet jeśli standardowy protokół opisany w tym badaniu jest dobrze dopasowany do większości pacjentów, niektóre osoby są trudniejsze do zmierzenia; Znalezienie prawidłowych pozycji jest trudne, ale możliwe, ale prawie niemożliwe jest znalezienie stabilnej. W takim przypadku operator może zdecydować się na nagrywanie dłuższych serii, do 200 akwizycji. Najlepsze przejęcia są później określane offline na podstawie parametrów jakościowych. To filtrowanie powinno być stosowane w czasie rzeczywistym w przyszłości. W praktyce operator rejestruje więcej serii i/lub akwizycji niż 4 serie po 10 akwizycji protokołu stanu idealnego. Jednak zwykły czas pomiaru wynosi około 5 minut dla jednego miejsca, przedramienia lub nogi. Jeśli po kilku minutach HMI nie wykryje żadnej prawidłowej pozycji, pomiar zostaje zatrzymany i jest uważany za niepowodzenie. W przypadku obecnego urządzenia i protokołu stwierdzono, że wskaźnik awaryjności wynosi mniej niż 5%²⁵.

Obecne urządzenie ma różne ograniczenia:
orazWymiary i waga: Obecne urządzenie BDAT jest przenośne: mieści się w standardowym bagażu i waży około 25 kg. Waga ta jest jednak duża w stosunku do najnowszych urządzeń ultradźwiękowych. Można rozważyć nową konstrukcję elektroniczną, w której sonda i przetwarzanie sygnału pozostaną takie same. Możliwe jest jednak przenoszenie obecnego urządzenia, szczególnie w przypadku pacjentów o ograniczonej sprawności ruchowej leżących w łóżku lub w domu.

orazSzybkość akwizycji: Obecna częstotliwość odświeżania wynosi około 2-4 Hz, co oznacza, że eksploracja pomiaru jest powolna w porównaniu z rzeczywistą akwizycją w czasie rzeczywistym, tj. wyższa niż 25 Hz. W przyszłości można to poprawić, biorąc pod uwagę szybsze komputery, szybszą analizę danych i transmisję między elektroniką a komputerem. Zwiększenie szybkości akwizycji poprawiłoby łatwość obsługi pomiaru, w szczególności poszukiwanie prawidłowego położenia sondy.

orazGrubość tkanki miękkiej: Obecne podejście jest ograniczone przez dużą warstwę tkanek miękkich, zwykle przekraczającą 10 mm. W tym przypadku pierwszy przychodzący sygnał jest powiązany ze ścieżką tkanek miękkich, a nie z kością korową. W związku z tym vFAS i związany z nim kąt dwukierunkowy nie mogą być używane. Podobnie, vA0 jest bardzo trudny do zmierzenia w przypadku dużych warstw tkanek miękkich. Bez tych dwóch prędkości nie można wykonać problemu odwrotnego. W przyszłości można by zastosować inne metody korekcji dwukierunkowej, na przykład za pomocą technik obrazowania. Pacjenci z dużymi warstwami tkanek miękkich są zwykle związani z otyłością i wskaźnikiem masy ciała (BMI) większym niż 30 ^kg.m-2.

orazRegularność falowodu: Podejście do problemu odwrotnego zakłada regularny falowód z wieloma ścieżkami propagacji. U pacjentów z osteoporozą wewnętrzny interfejs korowy może być nieregularny, co oznacza słabe obrazy widma fal kierowanych, szczególnie w dolnej części. Pacjenci ci są zwykle kojarzeni z dużą niejednoznacznością roztworu. Jeśli tkanka miękka lub złe ustawienie nie mogą być uważane za źródło słabego obrazu widma, a wartość parametru lowk jest niska, falowód jest przypuszczalnie nieregularny i brane jest pod uwagę rozwiązanie o najmniejszej grubości. Można również stosować podejścia oparte na uczeniu maszynowym, które nie wymaga modelowania fizycznego²⁹.

Jak wspomniano we wprowadzeniu, obecnie złotym standardem w wykrywaniu pacjentów zagrożonych złamaniami kruchymi jest DXA, który ma pewne ograniczenia: duże rozmiary, słabą dostępność w niektórych regionach, stosunkowo wysoki koszt i stosunkowo umiarkowaną skuteczność. Pierwsze ograniczenia można złagodzić za pomocą urządzeń ultradźwiękowych, znanych ze swojej atrakcyjnej przenośności i kosztów. Jednak zdolność do skutecznego wykrywania pacjentów zagrożonych powinna być co najmniej równoważna z DXA. W rzeczywistości czasami oczekuje się, że będzie ona wyższa niż DXA, aby uzasadnić dostosowanie większości odniesień (decyzje medyczne, leczenie, koszty, pokoje itp.) związanych ze złotym standardem. Dlatego niektóre urządzenia ultradźwiękowe proponują substytuty BMD 10,11,42. Jednak jedną z wad parametrów klinicznych, takich jak aBMD, a także vFAS, jest integracja różnych właściwości kości korowej. Dlatego też inne urządzenia ultradźwiękowe, w tym BDAT, proponują uzupełniające się punkty widzenia, proponując parametry, które mogą być łatwiej interpretowane przez lekarza i pacjenta, takie jak porowatość kory mózgowej, grubość czy rozkład wielkości porów⁹. Parametry te odzwierciedlają właściwości geometryczne i materiałowe: kość korowa może być potencjalnie oceniana pod względem niezależnych zmian ilości lub jakości. Ten punkt widzenia może być bardzo pomocny w badaniu różnych możliwych przyczyn łamliwości kości. Na przykład wewnątrz- lub zewnątrztorebkowe złamania szyjki kości udowej, tj. miejsca złamania szyjki kości udowej lub krętarza, mają różne pochodzenie medyczne⁴³. Podobnie, możliwe może być stosowanie różnych leków mających na celu różny wpływ na kość korową, a także pod względem ilości lub jakości³.

Należy zauważyć, że precyzja uzyskana za pomocą BDAT dla grubości kory mózgowej (0,1 mm) jest lepsza niż w przypadku innych metod ultradźwiękowych, zwykle powyżej 0,25 mm⁴⁴. Różnica ta wynika częściowo z faktu, że problem odwrotny BDAT uwzględnia połączone różnice geometryczne i materiałowe. Niektóre inne podejścia, takie jak echo impulsowe, uwzględniają unikalne właściwości materiału kostnego dla wszystkich pacjentów^10,44. Ta wartość precyzji, około 0,1 mm (CV około 3%), jest rzeczywiście kluczowa dla precyzyjnego rozróżnienia pacjentów, ponieważ zakres grubości jest mniejszy niż 2 mm. Precyzja porowatości kory mózgowej (1%, CV około 14%) nie jest jeszcze tak dobra, jak w przypadku grubości. Zaobserwowano już jednak znaczną poprawę w stosunku do poprzedniego badania odtwarzalności²² . Można się spodziewać, że podobne ulepszenia mogą zostać osiągnięte w niedalekiej przyszłości dzięki przyszłym ulepszeniom HMI, szczególnie pod względem liczby klatek na sekundę bliższej czasu rzeczywistego.

BDAT może być stosowany na dużą skalę do badań przesiewowych populacji w regionach, w których DXA nie jest powszechnie dostępne. Co więcej, najnowsze wyniki kliniczne wykazały potencjał, że BDAT może być jeszcze bardziej wydajny niż DXA. Wyniki te powinny być jednak potwierdzone poprzez włączenie większej liczby pacjentów. Kolejnym wyzwaniem powinny być badania wieloośrodkowe i/lub podłużne ^11,12. Urządzenie BDAT jest jednak nadal prototypem nadającym się do współpracy naukowej, ponieważ zostało to już zrobione w Niemczech²³ i Wielkiej Brytanii²⁴ . Potrzebne są wysiłki w kierunku uprzemysłowienia nowej generacji urządzeń BDAT, które z pewnością byłyby szybsze i bardziej przenośne.

Materiały

Lista materiałów użytych w tym artykule

Nazwa	Firma	Numer katalogowy	Komentarze
	Notebook komputerowy HP	mod Zbook (16 Go RAM, Inrel Core i7)	do odbierania próbkowanych sygnałów i stosowania kroków przetwarzania sygnału. Wyniki są wyświetlane w czasie quasi rzeczywistym (do 4 na sekundę) za pośrednictwem dedykowanego interfejsu człowiek-maszyna (HMI, BleuSolid, Pomponne, Francja) umożliwiającego prowadzenie pomiaru;
Transformator elektroizolacyjny	REOMED, Solingen, Niemcy	IEC / EN 60601-1	w celu ochrony urządzenia, pacjenta i innych urządzeń przed wszelkimi zagrożeniami elektrycznymi
Moduł elektroniczny	Althaï s, Tours, Francja	we własnym zakresie	wysyłają sygnały wzbudzenia (połowa okresu ujemnego napięcia 170 V) i dyskretnie odbierają sygnały (1024 punkty czasowe na odbiornik przy częstotliwości próbkowania 20 MHz i dynamice 12 bitów) przed wysłaniem ich do komputera. Opóźnienia i wzmocnienie można dostosować za pomocą praw liniowych w celu optymalizacji próbkowania danych w dostępnych zakresach.
Interfejs człowiek-maszyna	BleuSolid, Pomponne, Francja	Nie dotyczy
Przełącznik nożny	Kosa, Niemcy	Przełącznik nożny USB 2	do rozpoczęcia serii akwizycyjnej
Linijka	Westcott, USA	10417	w celu zlokalizowania miejsca pomiaru
Promień sondy ultradźwiękowej	Vermon, Tours, Francja	we własnym zakresie	Częstotliwość centralna 1 MHz, 24 odbiorniki o rastrze 0,8 mm i dwa bloki po 5 nadajników o rastrze 1 mm.
Sonda ultradźwiękowa piszczel	Vermon, Tours, Francja	własna	częstotliwość centralna 0,5 MHz, 24 odbiorniki o rastrze 1,2 mm i dwa bloki po 5 nadajników o rastrze 1,5 mm
Sondy ultradźwiękowe			zaprojektowane zgodnie z geometrią dwukierunkową: pojedyncza tablica odbiorników otoczona dwoma układami nadajników. Trzy układy są wyrównane, odizolowane mechanicznie i elektrycznie w celu zminimalizowania sygnałów sprzęgających.   Sondy są przystosowane do dwóch różnych miejsc, jednej trzeciej dystalnej kości promieniowej i środkowej kości piszczelowej.