$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Osteoporoz ve buna bağlı frajilite kırıkları halen önemli bir halk sağlığı sorunu oluşturmaktadır1. Özellikle, dünya çapında kalça kırığı sayısının 2050 yılına kadar iki katına çıkması beklenmektedir2. Kemik kırılganlığı, kırılganlık kırığı olayından önce önemli uyarı işaretleri olmaksızın yavaş ve sessiz bir demineralizasyon ve kemik kaybı sürecinden kaynaklanır. Kırılganlık kırığı riski taşıyan hastaları tespit etmek için mevcut altın standart, kalibre edilmiş gri piksel3 ile 2D, düşük çözünürlüklü bir X-ışını görüntüsü sağlayan çift X-ışını absorpsiyometrisidir (DXA). Bu görüntüden, ana kırılganlık kırığı bölgeleriyle ilişkili farklı ilgi bölgelerinde alansal kemik mineral yoğunluğunu (g.cm-2'de aBMD) çıkarmak mümkündür: omurga, bilek ve kalça. Kırılganlık kırılma oranı arttıkça aBMD değeri azalır3. Ayrıca, normal sağlıklı bir popülasyona göre T-skoru normalizasyonu, farklı üreticiler tarafından önerilen cihazlarla ölçülen hastaların karşılaştırılmasına izin verir. DXA T-skoru, Dünya Sağlık Örgütü tarafından osteoporoz tanısını üç aşamada tanımlamak için önerilmiştir: normal (T-skoru <-1), osteopenik (-1 < T-skoru < -2.5) ve osteoporotik (T-skoru < -2.5)4.
DXA çeşitli sınırlamalar sunar: boyutu, nispeten yüksek maliyeti, özel bir odaya olan ihtiyacı ve kırıklı ve kırıksız arasında ayrım yapma yeteneğinin yanı sıra Latin Amerika gibi çok sayıda ülkede bulunabilirliği orta düzeydedir5. Bu nedenle, frajilite kırılma riski tahmini için tarama araçları olarak DXA alternatiflerine ihtiyaç vardır6. Bununla birlikte, kantitatif bilgisayarlı tomografi ve türevleri7, manyetik rezonans görüntüleme (MRI)8 gibi bazı DXA alternatifleri de hacimlidir ve yaygın olarak bulunmaz. Kantitatif ultrason (QUS), taşınabilir, sağlam, kullanımı kolay tarama cihazları için potansiyel sunar. Kortikal kemik değerlendirmesi için, birkaç kHz'den birkaç MHz'e kadar değişen farklı frekanslarla ilişkili farklı cihazlar ve transdüserlerin transdüser ve tibia gibi uzun bir kemiğin ekseni ile hizalandığı transdüser, retro difüzyon9, darbe eko10 ve eksenel transmisyonda farklı transdüser konumlandırması ile ilişkili farklı cihazlar geliştirilmiştir. Bazı cihazlar aBMD vekilleri11 sağlarken, diğerleri hızlar12 veya zayıflama katsayıları9 gibi "klasik" ultrasonik parametreler ve hatta kortikal kalınlık, gözeneklilik veya gözenek boyutu dağılımı9 gibi geometrik ve malzeme parametreleri sağlar. Bununla birlikte, bugüne kadar, QUS, kısmen cihazlar arasındaki homojenizasyon eksikliği ve operatör bağımlılığı nedeniyle, kemik değerlendirmesi için klinik uygulamada yaygın olarak kullanılmayı henüz başaramamıştır13.
DXA alternatifleri olarak önerilen QUS teknolojileri arasında, eksenel transmisyon (AT), ölçümün (i) kolayca erişilebilir ve (ii) kırılganlık kırıklarının ana bölgelerinden birine, yani el bileğine yakın bir bölge olan önkolda yapılabilmesi avantajına sahiptir. Önerilen ilk AT parametresi, kortikal tabakadaki ultrasonik yayılma hızına, ses hızının (SOS) veya ilk varış sinyalinin (vFAS) hızına, sinyal işlemeye ve cihazlara bağlı olarak, bazıları ticari12,14 olanlar ve diğer laboratuvar prototipleri15,16. Bu parametre, 1990'ların sonlarından bu yana çeşitli klinik çalışmalarda KMY ile benzer performanslara sahip frajilite kırığı olan veya olmayan hasta grupları arasında ayrım yapabilmiştir14,15. Ayrıca, klinik uygulamasını ve sağlamlığını gösteren çok merkezli boylamsal çalışmalar için başarıyla uygulanmıştır12. vFAS hassasiyeti, prob ile kemik yüzeyiarasındaki açı 16,17 arasındaki açıyı azaltmak için iki zıt yayılma yönünün birleştirilmesiyle geliştirilmiştir. Bu bakış açısı iki yönlü AT (BDAT) olarak ifade edilmiştir.
vFAS klinik ilgi göstermiş olsa bile, BMD'ye benzer şekilde ana dezavantajı, geometrik ve malzeme özellikleri gibi farklı temel kortikal kemik özelliklerini bir araya getirmesi ve klinik yorumunu basit hale getirmesidir. Bu nedenle, kılavuzlu dalgaların dalga kılavuzu özelliklerine olan ince duyarlılığı nedeniyle potansiyeli göz önünde bulundurularak, kılavuzlu dalga bakış açısı önerilmiştir. Bu yaklaşım, sinyal işleme, dalga kılavuzu modelleme ve ters problemleri birleştirmelidir ve büyük ölçüde, örneğin plakalar veya tüpler gibi metalik dalga kılavuzları dikkate alınarak tahribatsız muayenede kullanılır18. Bu nedenle, 2010 yılından bu yana, kemik taklit eden fantomlardan19'dan ex vivo doğrulama20'ye ve in vivo ölçümlere21 kadar ikinci nesil bir BDAT cihazı adım adım geliştirilmiştir. Cihaz, Fransa22, Almanya23, Birleşik Krallık24 ve Şili25'teki klinik çalışmalarda başarıyla test edilmiş ve başarı oranı ve hasta ayrımcılığı açısından iyileştirici sonuçlar göstermiştir.
Bu çalışma, mevcut BDAT ultrasonik cihazının tekrarlanabilirliğini araştırmayı amaçlamaktadır. İlk olarak, cihaz ve ölçüm protokolü detaylandırılacaktır. 14 katılımcı ve 3 operatör ile elde edilen sonuçlar, frajilite kırığı riski taşıyan hastaların tespiti için popülasyon taraması açısından sunulacak ve tartışılacaktır.
Ölçüm prensibi: sinyal işleme, ilgilenilen parametreler ve kalite parametreleri
İki yönlü eksenel iletim (BDAT) cihazı, ana ultrasonik prob, elektronik modül ve bilgisayar olmak üzere farklı parçalardan oluşur. Tam liste Malzeme Tablosunda detaylandırılmıştır ve Şekil 1'de gösterilmiştir. Aşağıda, ilgilenilen parametreler, ölçüm kalitesi parametreleri ve ölçüm protokolü açıklanmaktadır.
vFAS (İngilizce)
Örneklenen sinyaller bilgisayar tarafından alındıktan sonra, farklı adımlar izlenerek işlenir. İlk adım, daha önce açıklanan protokolü kullanarak FAS'ı tespit eden zaman alanındaki sinyal işlemeden oluşur16,17. Her alıcı için varış süresi elde edildikten sonra, her iki yayılma yönünde elde edilen hızların harmonik ortalaması olan ve daha sonra vFAS olarak gösterilen FAS hızını belirlemek mümkündür. Her iki yayılma yönünden gelen bilgileri birleştirerek, prob ve kemik yüzey yönleri arasındaki değer açısını elde etmek ve tarafsız bir vFAS değeri16 elde etmek mümkündür. Bu çift yönlü açı daha sonra alfa olarak gösterilir ve ölçüm kalitesinin bir parametresi olarak kullanılır. Bu zamansal işlem aynı zamanda kemik yüzeyi ile prob arasındaki yumuşak doku kalınlığının tahmin edilmesine de izin verir, ST.Th26 olarak gösterilir.
Güdümlü dalga spektrumu görüntüsü
İkinci adım, f ve k ile gösterilen zamansal ve uzamsal frekansları göz önünde bulundurarak Fourier alanında sinyal işlemeden oluşur. Yaklaşım, uzaysal-zamansal sinyallerin Norm fonksiyonuna dönüştürülmesine izin veren SVD tabanlı bir yöntemdir, ayrıca Şekil 2'de bir in vivo önkol19 için gösterildiği gibi kılavuzlu dalga spektrumu görüntüsünü (GWSI) de gösterir. Yöntem, iki Fourier dönüşümünü (zaman ve uzay) ve bir tekil değer ayrıştırmasını (SVD) birleştirerek, kortikal kemik tabakası tarafından yönlendirilen modların alınan sinyallerdeki (0-1 ölçeğinde) varlık oranını görselleştirmeye izin verir. GWSI, uzaysal-zamansal Fourier Dönüşümünün bir geliştirmesi olarak yorumlanabilir, her piksel bağımsız bir f frekansı ve k dalga sayısı düzlemi ile ilişkilendirilir. Yaklaşımın, malzeme zayıflamasının27 ve doğrusal kalınlık değişiminin28 etkisini hesaba katmak için geliştirildiğini unutmayın.
A0 modu ile ilişkili spektrumun üst kısmına ve ayrıca en yüksek faz hızı değerleri, yani 4 mm·μs-1'den büyük olan en düşük kısma özellikle dikkat edilecektir. Bu kısım ilgilenilen bölgeye karşılık gelir 3 (ROI 3)29. Daha sonra lowk olarak gösterilen ROI 3'ün ortalama değeri de bir kalite parametresi olarak kullanılır. Büyük bir değer, kemik arayüzlerinde net dalga yansımalarına izin veren düzenli bir dalga kılavuzuna karşılık gelir. Değer düşerse, bunun nedeni düzensiz bir dalga kılavuzu veya yanlış yerleştirilmiş bir prob olabilir.
Dalga kılavuzu modeli
Kılavuzlu dalga dağılımı veya her bir kılavuzlu modun faz hızının frekansa göre değişimi, dalga kılavuzunun hem malzeme hem de geometrik özelliklerine bağlıdır. Bu nedenle, özel sinyal işleme, dalga kılavuzu modelleme ve ters problem şemaları kullanarak bu özellikleri elde etmek potansiyel olarak mümkündür. BDAT durumunda, dalga kılavuzu modeli, dalga kılavuzu malzemesine ve bir geometrik parametreye, kalınlığa30 bağlı olarak 2D enine izotropik serbest bir plakaya karşılık gelir. Kortikal kemik materyali, kemik matriksi için sabit parametreler ve değişken gözeneklilik31 dikkate alınarak homojenize edilir. Bu nedenle, ters problem, kortikal kalınlık (Ct.Th) ve kortikal gözeneklilik (Ct.Po) olarak gösterilen iki parametreye bağlıdır. Malzeme absorpsiyonu, dalga kılavuzu eğriliği ve çevredeki yumuşak dokuların etkileri, ölçümü etkilese bile modelde dikkate alınmaz. Bununla birlikte, ters problem sonucu üzerindeki ağırlıkları belirleyici bulunmamıştır, bu da iki ana ilgi alanındaki (A0 ve en alt kısım) modların eğrilik ve yumuşak dokular tarafından önemli ölçüde değişmediği anlamına gelir32.
Ters problem
Başlangıçta, ters problem iki adıma ayrıldı: birincisi, deneysel kılavuzlu dalga dağılımını çıkarın ve ikincisi, dalga kılavuzu modeliyle karşılaştırın. Bu bakış açısı gürültü ve mod etiketlemesi ile sınırlıydı 30,32. Bu nedenle, Norm fonksiyonu bakış açısının bir uzantısı olarak bu sınırlamaların üstesinden gelmek için özel bir yaklaşım önerilmiştir. Her bir düzlem dalgasını bağımsız olarak düşünmek yerine, yalnızca dalga kılavuzu modeli tarafından sağlanan olası kılavuzlu dalgalar dikkate alınır20. Bu, model parametre alanında, yani Ct.th - Ct.Po düzleminde ifade edilen ters problem görüntüsüne yol açar (Şekil 2 sağ altta). En uygun modele maksimum konum verilirken, nihai ikincil tepe noktaları (gri noktalı ters problem görüntüleri ile gösterilir), açık gri çizgili deneysel modlarla fk karşılaştırmasında belirtilen belirsiz çözümlere karşılık gelir. Daha önce olduğu gibi, piksel değeri yapı ile normalleştirilir ve bu durumda, alınan sinyallerde belirli bir dalga kılavuzu modelinin varlığını yansıtır. Maksimum değer (maksimum ile gösterilir) ve ikinci maksimum ile fark (fark ile gösterilir) de kalite parametreleri olarak kullanılır.
Ters problem başlangıçta çevrimdışı hesaplama için, yani sinyaller elde edildikten sonra, model dalga sayılarının kesin değerleri kullanılarak önerilmiştir. Bu yaklaşım, ex vivo 20,33 ve in vivo 21,34,35 çalışmaları göz önüne alındığında hem radius hem de tibia bölgeleri için doğrulanmıştır. Bu hesaplamaları insan-makine arayüzüne (HMI) dahil etmek için, seyrek bir matris bakış açısı kullanılarak gerçek zamanlı uygulama ile uyumlu yaklaşık bir versiyon önerilmiştir36.
vA0
GWSI'dan, serbest plakanın veya Lamb modelinin33,35 ilk antisimetrik modu A0 ile ilişkili en yavaş kılavuzlu modun hızını çıkarmak da mümkündür. Yönlendirilmiş dalga spektrumunun üst kısmı, vA0 hızının değerini sağlayan eğim ile doğrusal olarak yaklaştırılabilir (Şekil 2, sol alt).
Parametre özeti:
Son olarak, ilgilenilen dört parametre ölçülür: (i) vFAS: İlk Gelen Sinyalin hızı (m·s-1); (ii) vA0: en yavaş kılavuzlu modun hızı (m·s-1); (iii) Ct.Th: kortikal kalınlık (mm); ve (iv) Ct.Po: kortikal gözeneklilik (%).
Dört kalite parametresi dikkate alınır: (i) alfa: çift yönlü açı (°); (ii) lowk: GWSI'nin en düşük kısmının ortalama değeri (0 ile 1 arasında normalleştirilmiş değer); (iii) maks: ters problem fonksiyonunun maksimumu (0 ile 1 arasında normalleştirilmiş değer); ve (iv) fark: ters problem fonksiyonunun birinci ve ikinci maksimumları arasındaki fark (0 ile 100 arasında normalleştirilmiş değer).
Tüm bu parametrelerin yanı sıra iki kılavuzlu dalga spektrumu görüntüsü (bir eşit yayılma yönü) ve ters problem görüntüsü, HMI tarafından yaklaşık 2 Hz'lik bir kare hızıyla "gerçek zamanlı" olarak görüntülenir. Tipik bir örnek Şekil 3'te gösterilmiştir. Aşağıdaki bölümde, bu parametreleri kullanma yöntemi ayrıntılı olarak açıklanmaktadır. Ana fikir, operatörün probu ölçüm alanında yavaşça hareket ettirmesi, sabit bir konum bulana ve bir dizi 10 alıma başlayana kadar arayüzün farklı parçaları tarafından sağlanan geri bildirimi dikkatlice gözlemlemesidir. En az dört tutarlı seri elde edildiğinde ölçüm sona erer ve otomatik bir rapor oluşturulur.