El bileği ve distal radyoulnar eklemleri incelemek için bidüzlemal videoradyografi

Özet

Bidüzlemar videoradyografi (BVR), iskelet kemiklerinin ve implantların üç boyutlu hareketini anlamak için gelişmiş bir görüntüleme tekniğidir. Distal üst ekstremitenin yoğunluğa dayalı görüntü hacimlerini ve videoradyografilerini birleştiren BVR, el bileği ve distal radioulnar eklemin in vivo hareketini ve ayrıca eklem artroplastilerini incelemek için kullanılır.

Özet

İskelet kinematiğinin in vivo olarak doğru ölçümü, normal eklem fonksiyonunu, patolojinin etkisini, hastalığın ilerlemesini ve tedavilerin etkilerini anlamak için çok önemlidir. İskelet hareketini anlamak için deri yüzeyi belirteçlerini kullanan ölçüm sistemleri, normal ve patolojik kinematik hakkında önemli bilgiler sağlamıştır, ancak özellikle dinamik aktiviteler sırasında bu sistemler kullanılarak doğru artrokinematik elde edilemez. Son yirmi yılda, iki düzlemli videoradyografi (BVR) sistemleri, birçok araştırmacının günlük yaşam aktiviteleri sırasında eklemlerin iskelet kinematiğini doğrudan incelemesini sağlamıştır. Distal üst ekstremite için BVR sistemlerini uygulamak için, bir denek belirlenmiş bir görevi yerine getirirken, distal radiusun ve elin videoradyografileri kalibre edilmiş iki X-ışını kaynağından elde edilir. Üç boyutlu (3D) rijit gövde pozisyonları, her bir BVR görünümüne en uygun 3D model projeksiyon kayıtları aracılığıyla videoradyograflardan hesaplanır. 3D modeller, bağımsız olarak elde edilen bilgisayarlı tomografi verilerinden elde edilen spesifik kemiğin yoğunluğa dayalı görüntü hacimleridir. Grafik işlemci birimlerini ve yüksek performanslı bilgi işlem sistemlerini kullanan bu model tabanlı izleme yaklaşımının, el bileği ve distal radyoulnar eklem biyomekaniğini değerlendirmede hızlı ve doğru olduğu gösterilmiştir. Bu çalışmada ilk olarak, el bileği ve distal radyoulnar eklem kinematiğinin değerlendirilmesinde in vitro optik hareket yakalama sistemi ile BVR'nin milimetre altı ve milimetre altı uyumunu belirleyen önceki çalışmaları özetledik. Ayrıca, el bileği ekleminin dönme merkezi davranışını hesaplamak, implant bileşenlerinin birbiri üzerindeki artikülasyon paternini değerlendirmek ve önkolun pronosupinasyonu sırasında ulnar varyansın dinamik değişimini değerlendirmek için BVR'yi kullandık. Gelecekte, karpal kemikler, düz panel X-ışını dedektörleri, daha fazla X-ışını kaynağı (yani, çok düzlemli videoradyografi) veya gelişmiş bilgisayarlı görme algoritmalarının eklenmesiyle daha ayrıntılı olarak yakalanabilir.

Giriş

İskelet kinematiğinin in vivo olarak doğru ölçümü, sağlıklı ve yerine konmuş eklem fonksiyonunu, patolojinin etkisini, hastalığın ilerlemesini ve tedavilerin etkilerini anlamak için gereklidir. İskelet kinematiğinin eklem yüzeyinde noninvaziv olarak ölçülmesi (artrokinematik), eklem patolojilerini ve osteoartrit gibi hastalıkları anlamak için çok önemlidir, ancak teknik olarak zordur. Daha önce, iskelet hareketini anlamak için deri yüzeyi belirteçlerini kullanan teknikler, sağlıklı ve patolojik kinematik hakkında önemli bilgiler sağlamıştır. Bununla birlikte, özellikle günlük yaşam aktiviteleri gibi dinamik aktiviteler sırasında bu teknikler kullanılarak doğru artrokinematik elde edilemez. Bu optik sistemler, insan hareketi analizindeki ana hata kaynağı olan altta yatan kemiklere göre cilt hareketi nedeniyle doğal olarak doğruluk açısından sınırlıdır ^1,2.

Üç boyutlu (3D) iskelet kinematiğini ölçmek için mevcut son teknoloji yöntemler, görüntü tabanlı izleme, yani çift düzlemli videoradyografi (BVR)³ ve seri bilgisayarlı tomografi (BT) cilt⁴ ve manyetik rezonans görüntüleme (MRI)^5'tir. Normal 3D BT ve MRI tabanlı teknolojiler, dünyadaki birçok hastanede son derece doğru ve erişilebilir olmasına rağmen, eklemlerin dinamik hareketini ölçemezler. Bu eksikliği gidermek için son yıllarda 4D CT tarama⁶ ve dinamik MRI⁷ gibi görüntüleme teknikleri geliştirilmiştir; Bununla birlikte, bu yöntemler hastaları ya yüksek radyasyon dozuna maruz bırakır ya da düşük zamansal çözünürlükten muzdariptir.

Yeni bilgisayarlı görme algoritmalarını ve geleneksel x-ray sistemlerini birleştiren BVR'nin hayvanlarda ve insanlarda birden fazla eklem için doğru olduğu gösterilmiştir; işaretleyici tabanlı veya model tabanlı izleme algoritmaları ile çözülür. İşaretleyici tabanlı yaklaşımlar, kemiklere veya yumuşak dokuya yerleştirilen tantal boncukları izler ve hayvan ve in vitro testler için idealdir. Bununla birlikte, in vivo insan araştırmaları için engelleyici bir şekilde istilacıdırlar. Neyse ki, model tabanlı izleme algoritmalarındaki iyileştirmeler uygun bir alternatif sunuyor. İnsanlarda model tabanlı BVR izleme yaklaşımları, CT veya MRI ile elde edilen hacimsel görüntü setlerinin statik bir duruşta hazırlanmasını ve iki X-ışınının görüş alanındaki ilgi hareketlerinin yakalanmasını içerir. Model tabanlı izleme uygulamalarının çoğu daha sonra statik BT veya MR görüntülerinden kemik veya implantın dijital olarak yeniden yapılandırılmış radyografilerini (DRR) oluşturur ve DRR'ler ile videoradyograflar arasındaki benzerliği gösteren ölçümleri kullanarak bunları özelliği geliştirilmiş videoradyograflarla eşleştirir⁸. Bu işleme kemik veya implantın "izlenmesi" denir.

Kemiklerin veya implantların izlenmesinin birincil çıktı değişkenleri, eklem kinematiği, bağ uzamaları ^9,10, kıkırdak kalınlığı ¹¹ için bir vekil olarak eklem aralığı, eklem teması ^12,13 ve diğer biyobelirteçlerin hesaplanabildiği rijit cisim kinematiğidir. Son zamanlarda, el bileği biyomekaniği, total el bileği artroplastisi (TWA) ve distal radioulnar eklem (DRUJ) hesaplamasında model tabanlı izleme BVR'nin doğruluğunu ^{belgeledik14,15}. Aşağıdaki bölümde, çeşitli görevler sırasında iskelet bileğinin, total el bileği artroplastisinin ve distal radioulnar eklemin hareketini incelemek için bu onaylanmış yöntemin ayrıntılı bir protokolü sunulmaktadır. Kemiklerin ve implantların yoğunluğa dayalı görüntü hacimlerini BT görüntü hacimlerinden segmentlere ayırıyor, bu kısmi görüntü hacimlerini videoradyografiler içinde takip ediyor ve bu yöntemin güçlü yönlerini ve sınırlılıklarını göstermek için rotasyon merkezi, temas paterni ve ulnar varyans gibi sonuçları belirliyoruz.

Protokol

Bu çalışma, AAHRPP tarafından akredite edilmiş bir IRB olan Lifespan - Rhode Island Hastanesi Kurumsal İnceleme Kurulu (IRB) tarafından onaylanmıştır. Toplam 16 hasta kurumsal kılavuzlara göre imzalı bilgilendirilmiş onam verdi.

1. Veri toplama

1. Bilgisayarlı Tomografi (BT)
   1. Numuneleri veya denekleri BT için hazırlayın.
      NOT: Doğruluk değerlendirmesi için total el bileği artroplastisi yapılmış dört sağlam kadavra örneğinden^14,15, 6 sağlam ön kol (70.5 ± 12.3 yıl, 4 sağ, 2 kadın) ve 1 önkol (49 yıl, sağ, kadın) kullanıldı. İn vivo çalışma için¹⁶, sağlıklı bilekleri olan 10 gönüllü (57.0 ± 5.2 yıl, sekiz kadın, dokuz sağ el dominant) ve total el bileği artroplastisi uygulanan 6 hasta (74.7 ± 5.6 yıl, iki kadın, iki sağ, dört baskın bilek) dahil edildi.
   2. Her bir deneği⁴ görüntülemek için, 25 cm'lik bir veri toplama çapı ile distal önkolun BT taramasını alın ve 16 kesitli bir klinik BT tarayıcı kullanarak Bone Plus evrişim algoritmasını kullanarak 20 cm'lik bir görüş alanı ile yeniden yapılandırın. 80 kVp ve 80 mA tüp ayarlarını kullanın. Önkol maruziyetini distal radius yüzeyinin yaklaşık 7 cm proksimalinden distal falankslara kadar sınırlayın.
      NOT: Bu çalışmada görüntünün voksel boyutları önkol enine düzleminde 0.39 × 0.39^mm2 ve önkol uzun ekseni boyunca 0.625 mm olarak bulunmuştur.
   3. Elde edilen BT görüntülerini Digital Imaging and Communications in Medicine (DICOM) formatında HIPPA onaylı bir protokol kullanarak laboratuvarın bilgisayarına aktarın.
2. Bidüzlemli videoradyografi (BVR)
   1. Deney düzeneği (Şekil 1). Distal yarıçapı, distal ulnayı, karpal kemikleri ve metakarpalları yakalayabilen bir görüş alanı çapının yanı sıra çeşitli bilek duruşlarında kemik tıkanıklıklarını en aza indirecek bir perspektif sağlamak için, görüntü yoğunlaştırıcıları (II'ler) ön kolların altına yerleştirin ve aralarındaki açıyı 110°'ye ayarlayın (X-ışını kaynakları, yere göre 55° oryantasyonda karşılık gelen II'lere doğru aşağıyı hedefler).
      1. Kemiklerin/implantın görünürlüğüne bağlı olarak 65 ila 75 kV arasında bir pozlama ayarı ve 80 mA ile sürekli modda X-ışınları oluşturun. Her iki X-Ray-II çiftinin kaynak-görüntü mesafelerini yaklaşık 130 cm'de ve kaynak-el mesafelerini 90 cm'de ayarlayın.
         NOT: Bu çalışmada, denekler laboratuvarda bir x-ışını koruyucu kurşun önlük, bir tiroid yaka kalkanı boyun örtüsü ve bir dijital dozimetre giydiler. Tüm araştırmacılar deney sırasında kurşun koruyucu bir camın arkasında duruyorlardı.
      2. II'ler arasına bir sehpa veya benzeri bir yapı yerleştirin ve bunu denekler için öğelerin (örn. kapı kolu, çekiç, sürahi) yerleştirilmesi için kullanın.

Şekil 1. Deney düzeneği. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

2. Bozulma ızgarası (Şekil 2A). Dış elektromanyetik alanlar, X-ışını demetlerinin (yani elektron yollarının) davranışını etkilediğinden, radyografik görüntüler doğrusal olmayan iğne yastığı ve S-şekilli bozulmalardan muzdariptir¹⁷. Bu bozulmaları gidermek için, düz bir hassas delikli sac levha parçasının görüntüsünü kullanın (yani, bozulmayan ızgara)^18,19. II'lerin üzerine bozulma olmayan ızgaralar yerleştirin ve videolarını videoradyografi verilerinden önce ve sonra yakalayın. Bu görüntüleri her iki X-ışını kaynağı için 8 bit TIFF görüntüleri olarak kaydedin.
3. Kalibrasyon (Şekil 2B). Nesnelerin (kemikler veya implantlar) 3D koordinatlarını hesaplamak için X-ışını kaynaklarının kalibrasyonu ^gereklidir20. X-ışını kaynağının 3B uzaydaki konumu ve yönü ile görüntü alanı (yani radyograflar) arasında bir ilişki oluşturmak için, düzenli aralıklarla yerleştirilmiş radyo-opak kürelere sahip bilinen bir geometriye sahip bir 3B nesnenin görüntülerini kullanın. Akrilik levhalardan ve çelik küreli ara parçalardan yapılmış kübik bir nesne (yani kalibrasyon küpü), BVR ayarları¹⁹ için etkili ve doğrudur.
   1. Kalibrasyon küpünün¹⁸ radyografilerini çeşitli yönlerde, deneyi gerçekleştirmeden önce 4 kez ve deneyden sonra 4 kez yakalayın. Deney sonrası kalibrasyonlar, deney düzeneğinin deney boyunca değişmeden kaldığını doğrulamaya hizmet eder. Her iki X-ışını kaynağı için tüm görüntüleri 8 bit TIFF görüntüleri olarak kaydedin.

Şekil 2. A) Bozulma ızgarası. B) Kalibrasyon küpü ve referans öğeleri. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

4. Görevlerin açıklaması. Fleksiyon-ekstansiyon, radyal-ulnar deviasyon ve dolaşım sırasında bilek hareketini yakalamak için çeşitli görevler tanımlayın¹⁶. DRUJ hareket görevlerini, bir T-kollu kapı kolunun saat yönünde ve saat yönünün tersine dönmesi olarak tanımlanan bir tutma pozisyonundan tam pronasyona veya tam supinasyona ön kol dönüşü olarak tanımlayın²¹. Bileğin nötr pozisyonunda (yani nötr statik görev) videoradyografilerini alın. Nötr pozu, elin arkası ön kolun arkası ile aynı hizada olacak şekilde yakalayın.
   NOT: Tüm görevler için başlangıç duruşunda, omuz addüksiyondaydı, dirsek eklem seviyesinde desteklendi ve önkol kısıtlanmamış ve nötr rotasyonundaydı. Fleksiyon-ekstansiyon ve radyal-ulnar deviasyon görevleri, avuç içinin sırasıyla elin volar-dorsal ve radial-ulnar taraflarına göre dönmesi olarak tanımlandı. Circumduction, her yönde maksimum aktif hareket açıklığı elde etmeyi amaçlayan bilek hareketi olarak tanımlandı.
5. Videoradyograf alımı. Kamera deklanşör hızı 500 μs olarak ayarlandığında 200 Hz hızında video radyograflar elde edin. 2 tekrar elde etmek için her görev için iki saniyelik görüntüleme kaydedin. Her X-ışını kaynağından alınan videoradyografları 8 bit formatında CINE video olarak saklayın.
   NOT: Bu çalışmada, radyografiler görüntü boyutunda 1760×1760 piksel idi ve piksel başına yaklaşık 0.22mm×0.22mm çözünürlüğe sahipti.
   NOT: Mevcut X-ışını kaynaklarının görüş alanını büyütme (FOV) seçeneği vardır. Ancak bu çalışmada donanımın sunduğu büyütülmüş mod kullanılmamıştır. Görev tanımına bağlı olarak, bu özellik daha yüksek çözünürlüklü çıktılar elde etmenizi sağlayabilir. Bu özellik kullanılırsa, Bozulmayı Giderme Izgarası ve Kalibrasyon adımları bu modda yakalanmalıdır.

3. Optik hareket yakalama (OMC) - doğruluk çalışması
   1. Önceki çalışmalarda açıklanan optik hareket yakalama (OMC) veri toplama protokolünü kullanın ^14,15. Her biri dört adet 9,5 mm çapında retro-yansıtıcı işaretleyici küreye sahip iki işaretleyici kümesini üçüncü metakarpal ve yarıçaplı kemiklere doğrudan takın. 8 optik hareket yakalama kamerasını kullanarak işaret hareketini yakalayın. OMC veri alımının başlangıcını aktif bir düşük harici tetikleyici ile BVR ile senkronize edin.

2. Veri İşleme

1. Bilgisayarlı tomografi (BT)
   1. DICOM anonimleştirme. DICOM'ları anonimleştirin ve 3D modelleme için görüntü işleme yazılımına aktarın (örneğin, Mimics Materialise). Mimics Materialise yazılımında, Dosya araç çubuğundan Anonimleştir penceresini açın. Yeni pencere açıldığında sol sekmeden Tüm Görselleri Anonimleştir'i seçin ve sağ altta bulunan Anonymize butonuna basın.
   2. Segmentasyon (Şekil 3). Daha önce açıklanan bir metodolojiyi izleyerek ilgilenilen kemikleri veya implantları^{segmentlere ayırın} 16,22,23.
      1. Kısaca, Segment araç çubuğundan Threshold (Eşik ) seçeneğini belirleyin ve pencere açıldığında Bone için önceden tanımlanmış eşik kümelerini seçin. Eşikli maske oluşturulduktan sonra, Bölünmüş Maske seçeneğini kullanın ve Kemiklerin veya implant bileşenlerinin Bölgelerini manuel olarak boyayın.
      2. Tüm bölgeler renklendiğinde, Tamam'a tıklayın ve yazılımın her kemik ve implant bileşeni için maskeleri otomatik olarak oluşturmasını bekleyin. Amaçlanan bölgeleri tamamen kapladıklarından emin olmak için tüm maskeleri kontrol edin. Değilse, bölgeleri manuel olarak ve tek tek değiştirmek için Maskeleri Düzenle seçeneğini kullanın.
         NOT: El bileği eklemini incelemek için 2. ve 3^{. metakarpallar} ve radius kemikleri bölümlere ayrılmalı ve diğer kemiklerden ayrılmalıdır. İsteğe bağlı olarak, kapitat ve diğer kemikler işlem sonrası aşamalar için bölümlere ayrılabilir.
         NOT: DRUJ'u incelemek için, yarıçap ve ulna bölümlere ayrılmalı ve diğer kemiklerden ayrılmalıdır.
         NOT: İmplantlar için, radyal bileşen ve vidaları ile karpal bileşen segmentlere ayrılmalıdır.

Şekil 3. El bileğinin bilgisayarlı tomografi görüntüsü ve radius, üçüncü metakarpal ve ulna'nın yeniden yapılandırılmış modelleri. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

3. Segmentasyon - yalnızca radyal bileşen. CT görüntülerinde radyal bileşenin yüzeyinin yoğunluğa dayalı çizgi artefaktları nedeniyle, radyal bileşen modeli görüntülerden bölümlere ayrılamaz. Bu nedenle, ekilen radyal bileşenin 3B modelini oluşturmak için bir CAD modeli veya 3B yüzey tarayıcı kullanın. Bu modeli, Mimics Materialise default'un yöntemini kullanarak radyal bileşenin kaba segmentasyonuna kaydedin.
   NOT: Bu çalışmada 0.1 mm çözünürlüğe sahip bir yüzey tarayıcısı kullanılmıştır. İmplant bir masaya yerleştirildi ve 3D tarayıcı, tüm dış yüzeyleri yakalamak için implantın etrafında manuel olarak hareket ettirildi.
4. Kısmi görüntü hacimleri. İzleme yazılımı, yoğunluğa dayalı bir görüntü hacmi olarak girilmek üzere bölümlere ayrılmış bir ilgi alanı gerektirir. Bu kısmi hacmi oluşturmak için, Segment araç çubuğundan Threshold (Eşik) seçeneğini belirleyerek tüm gri değerleri içeren görüntü hacminin tam maskesini oluşturun. Eşik değerlerini minimumdan maksimum gri değerlere ayarlayın.
   1. Ardından, Boole Operatörleri seçeneğini kullanın ve oluşturulan tam maskeden ilgilenilen kemiğin veya implantın segmentli hacmini çıkarın . Kayıt yazılımı yalnızca gri tonlamalı görüntü hacimlerini kabul ettiğinden, bu kısmi birimin maskesinin rengini siyah (RGB: 0,0/0,0/0,0) olarak değiştirin. DICOM'u Dışa Aktar'ı seçerek maskeyi DICOM olarak dışa aktarın... yazılımda. Çıktı, izole edilen kemiğin veya implantın bir görüntü hacmidir.
5. Üç boyutlu modeller. En Uygun Çözünürlük ile 3B Hesapla seçeneğini belirleyerek segmentli maskeleri üçgen örgülü 3B yüzey modellerine dönüştürün. Bu modelleri, köşelerin (noktaların) 3B konumunu ve bağlantı sıralarını içeren VRML formatında dışa aktarın. Modeller, koordinat sistemi oluşturma ve daha fazla görselleştirme için oluşturulur.

2. Videoradyografi Ön işleme
   1. XMALab¹⁸ İçe Aktarma. Her konu için XMALab yazılımında yeni bir deneme sürümü oluşturun. Klasör yollarını, Yeni Deneme Sürümü iletişim kutusundaki Bozulmayı Giderme Izgarası görüntülerine ve 8 kalibrasyon görüntüsü setine atayın.
      NOT: X-ışını kaynakları çekim boyunca ters çevrilirse, XMALab yazılımında çevirme seçeneği işaretlenmelidir.
   2. Bozulma. Radyografların deformasyonunu gidermek için, deneme sürümü oluşturulduktan ve yazılım ızgara görüntülerini Bozulmayı Giderme çalışma alanında açtıktan sonra, Bozulmayı Kaldır düğmesine tıklayın. XMALab, ızgaraların bozulmuş görüntülerini otomatik olarak belirler ve görüntüler. Tüm deliklerin merkezleri yazılım tarafından algılandıktan sonra, sol üst menüdeki çalışma alanını Kalibrasyon çalışma alanı olarak değiştirin. Daha fazla bilgi ve olası sorunları giderme konusunda yardım için başvuru^24'e bakın.
   3. Kalibrasyon. Bu aşamada, XMALab yazılımı kalibrasyon matrislerini belirler. XMALab'de X-ışını kaynakları, odak uzaklığı, kaynağın merkezi ve kaynaktan görüntüye mesafeler gibi içsel ve dışsal parametrelere sahip iğne deliği kameraları olarak modellenir ve kalibrasyon aşaması bu metrikleri belirleyen bir süreçtir¹⁸. Kalibre etmek için, tüm görüntü setleri ve her görünüm için kalibrasyon küpüne yerleştirilen 4 referans öğesinin her birini seçin. Kurulumda bu referanslar Piramit, Jack, Dome ve Helix idi (Şekil 2B). Yazılım noktaları otomatik olarak bulmazsa, Hesaplama Kalibrasyonu'na basın. Tüm kalibrasyon görüntüleri kalibre edildikten sonra Kamera Kurulumunu Optimize Et (Kamera Kurulumunu Optimize Et) öğesini seçin.
   4. Kalibrasyon Matrislerinin Dışa Aktarılması. Optimizasyondan sonra, 8 kalibrasyon setinin tümünü (yani, deneyden önce yakalanan 4 set ve deneyden sonra yakalanan 4 set) sıralayın ve bozulmamış en küçük hataya sahip görüntü setini seçin. Her kalibrasyon penceresinin sağ üst köşesindeki bilgi seçeneğine basarak bu değeri bulun. En düşük hataya sahip görüntüde, Dosya-Dışa Aktar menüsünden MayaCams 2.0 seçeneğini seçerek görüntüleri dışa aktarın. Daha fazla bilgi ve olası sorunları gidermek için başvuru^20'ye bakın.
   5. Videoradyografi Bozulma. Çalışma alanını İşaretçi izleme olarak değiştirin ve + simgesine basarak CINE biçiminde yakalanan videoradyografları yükleyin. Yükleme işlemi tamamlandıktan sonra, Dosya Dışa Aktar menüsünden Bozulmamış Deneme Görüntüleri'ni seçerek bozulmamış deneme görüntülerini dışa aktarın.
3. Model Tabanlı İzleme (Şekil 4)
   1. Radyografilerin/DRR'lerin Geliştirilmiş Özellikleri (Şekil 4A-B). Filtreleri kullanarak görüntü özelliklerini geliştirmek, doğru bir eşleşme elde etmek için gerekli bir adımdır. Yazılımda kontrast (yoğunluk algılama), Sobel (kenar algılama), Gauss (bulanıklaştırma/yumuşatma) ve Keskinleştirme (kenarları kalınlaştırma) dahil olmak üzere dört filtrenin parametrelerini tanımlayın. Bu filtreler, sol üst widget'ta gösterilen Rad Renderer veya DRR Renderer nesnelerine sağ tıklanarak seçilebilir.
      NOT: Bu çalışmada, radyografilerdeki kemik görüntülerinin kenarları, ölçek faktörü 3 ve karışım değeri 0.4 olan bir Sobel filtresi kullanılarak geliştirilmiştir. Ek olarak, radyografların kalitesini daha da iyileştirmek için alfa (görüntü kontrastı için) 2.5 ve beta (görüntü parlaklığı için) 0.9 olan bir kontrast filtresi kullanıldı. DRR'leri radyografla eşleştirmek için, 0.1 karışım değeri ve 1.7 ölçek faktörüne sahip bir Sobel filtresi ile 0.35'lik bir ışın yoğunluğu değeri kullanıldı.
   2. Manuel başlatma. DRR ve radyografiler arasında mümkün olan en iyi eşleşmeyi elde etmek için Pivot noktasını kullanarak her kemiğin ve bileşenin DRR'sini hareket ettirin. Pivot noktası , kullanıcıların seçilen DRR'yi sol alt widget'tan döndürmesine ve çevirmesine yardımcı olur. Bu adım gereklidir, çünkü optimizasyonu gerçekleştirmeden önce kemiğin veya implantın konumu en az bir kez ve en azından ilk kare için başlatılmalıdır. DRR'yi çeviri modunda hareket ettirmek için W tuşuna basın ve DRR'yi döndürmek için E düğmesine basın.
      NOT: Pivot noktası konumunu değiştirmek genellikle gerekli değildir, ancak araç çubuğunun Pivotu Taşı seçeneği kullanılarak mümkündür.
   3. Optimizasyon. DRR ve radyograflar arasında en iyi eşleşmeyi bulmak için C'ye basın veya araç çubuğundan Akımı İzle düğmesine tıklayın. Yazılımın varsayılan prosedürü, normalleştirilmiş çapraz korelasyon (NCC) maliyet fonksiyonunu ve hem parçacık sürüsü optimizasyon yöntemini hem de yokuş aşağı simpleks prosedürlerini^{kullanır 26,27}. Bu işlem aynı zamanda eşleştirme veya izleme adımı olarak da adlandırılır. Optimizasyon prosedürü, İzleme İletişim Kutusu aracılığıyla da yürütülebilir. Bunu yapmak için, araç çubuğundan İzleme İletişim Kutusu seçeneğine basın. Yeni pencere açıldığında, Çerçeveden ve Çerçeveye'yi 0 olarak değiştirin ve Tamam'a basın.
      NOT: NCC, radyografiler ve DRR'ler arasındaki benzerliği ölçen bir fonksiyondur ve 0'a yakın bir NCC değeri daha iyi bir eşleşme göstermiştir. Ancak bu değere ulaşmak tamamen görüntü filtrelerine ve yakalama özelliklerine bağlıdır. NCC'nin değeri "Terminal Penceresi" nde görülebilir.
   4. Eşleştirme işlemi - ilk geçiş. İlk kare otomatik olarak optimize edildikten sonra, İzleme İletişim Kutusunu açın, Kareden 0 olarak değiştirin ve Çerçeveye Kadar'ı son kare numarasıyla değiştirin (bu sağ altta bulunabilir). Yazılımın varsayılan eşleştirme işlemi, eşleştirme işlemi için DRR'nin başlangıç konumu olarak Geçerli kareyi kullanır. Optimizasyon işlemi, başlangıç konumunun 3 mm ve 3°'lik bir bölgesinde parçacık sürüsü optimizasyonu (PSO) yöntemini kullanır.
   5. Eşleştirme Süreci - Son Geçiş (Şekil 4C). Otomatik eşleştirme işlemi, zayıf bir başlatma veya filtreler nedeniyle başarısız olabileceğinden, optimizasyonun doğru şekilde yürütüldüğünden emin olmak için tüm kareleri manuel olarak kontrol edin. DRR ve radyografiler doğru şekilde eşleşmezse, bölüm 2.3.2 ve 2.3.3'te açıklanan adımları tekrarlayın.
      NOT: Bileğin izlenmesi ilk olarak, kemik örtüşmesinin etkilerini azaltmak ve eşleştirme doğruluğunu artırmak için ikinci ve üçüncü metakarpalları (MC2-MC3) içeren bir modelle gerçekleştirildi. Radyografilerde birleşik DRR pozisyonu ve rotasyonu bulunduktan sonra, eşleşen ve optimize edilmiş MC2-MC3 konumları MC3 kısmi görüntü hacmine dönüştürüldü ve MC3'ün başlangıç pozisyonları tohumlandı. Bu tohumlanmış pozisyonlar kullanılarak, MC3 izlendi. MC3'e benzer şekilde, tüm görevler sırasında videoradyografilerdeki yarıçap, ulna, radyal bileşen ve karpal bileşen de izlendi.

Şekil 4. A) Kemiklerin dijital olarak yeniden yapılandırılmış radyografileri (DRR'ler) ile bir X-ışını kaynağının yakalanan radyografisi. B) Gelişmiş (filtrelenmiş) radyografi ve DRR'ler. C) Optimizasyon işleminden sonra eşleşen DRR'ler. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

6. Veri azaltma. Kare hızına ve radyografik görüntü kalitesine bağlı olarak, hareketli ortalama yöntemi (5 karelik bir açıklık ile) gibi filtreleme yaklaşımları kullanılabilir. Bunu yapmak için, araç çubuğu menüsünden İzleme-Gelişmiş İletişim Kutusu'nu seçin ve Düzgün düğmesine tıklayın. Bu yaklaşım, bu çalışmadaki tüm yakalamalar için kullanılmıştır.
7. Kemik/implant pozisyonlarının dışa aktarılması. Her kemik radyografilerle eşleştirildikten ve tüm çerçeveler optimize edildikten sonra, İzlemeyi Kaydet düğmesine tıklayarak ve matris formatı çıktısını seçerek X-ışını dünyasındaki DRR'lerin 4x4 dönüşüm matrislerini dışa aktarın. Bu dönüşüm matrisi, kalibre edilmiş X-ışını dünyası koordinat sisteminde kemik/implantın kısmi görüntü hacminin dönüşünü ve ötelenmesini tanımlar ve ek işleme ve hareket hesaplaması için kullanılabilir. Kısmi görüntü hacminin koordinat sistemi, anlamlı biyomekanik analiz için anatomik koordinat sistemine aktarılmalıdır.

3. Veri Analizi

1. Koordinat Sistemleri (Şekil 5). Kinematiği raporlamak için, anatomik ve geometrik işaretlerine dayalı olarak kemikler ve implant bileşenleri için koordinat sistemleri (CS) oluşturun. Bu CS'ler oluşturulduktan sonra, yeni CS'deki modelleri, yinelemeli en yakın nokta algoritması kullanılarak Bölüm 2.1.4'te oluşturulan yüzey modellerine kaydedin. Bu bire bir kayıt, kinematik hesaplamaları için kullanılacak yeni bir 4x4 dönüştürme matrisini hesaplar.
   NOT: Bu çalışmada, yarıçap, ulna, üçüncü metakarpal, karpal bileşen, radyal bileşen, polietilen başlık ve kapitat için CS, daha önce açıklanan yöntemler kullanılarak^{oluşturulmuştur} 16,22,23.

Şekil 5. Kemiklerin ve implant bileşenlerinin koordinat sistemleri. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

2. Kinematik. Bu çalışmada, bilek hareketini MC3'ün yarıçapa göre hareketi olarak tanımlayın ve kinematik veriler nötr bilek pozisyonuna göre rapor edilmiştir.
   1. DRUJ analizi için, kinematiği ulnar CS'deki CS yarıçapının nötr poza göre göreceli hareketi olarak tanımlamak için sarmal hareket ekseni (HAM) parametrelerini²⁸ hesaplayın.
   2. Fleksiyon-ekstansiyon ve radyal-ulnar sapmada el bileğinin dönme merkezini (COR) belirlemek için HAM analizini kullanın.
   3. COR'un tek boyutlu proksimal-distal kaymasını tanımlamak için, sağlıklı bilekler için öngörülen COR'u, vida dönme ekseninden en kısa mesafeye sahip olan kapitat kemiğin merkezi ekseni üzerinde bir nokta olarak tanımlayın¹⁶.
   4. HAM vida eksenini kullanarak ulnar koordinat sistemindeki HAM rotasyonunu ayrıştırın ve önkolun pronasyonunu (pozitif) ve supinasyonunu (negatif) hesaplayın. Nötr poz, pronasyon-supinasyon (PS) sıfır olduğunda tanımlandı.
   5. Son olarak, ulnar varyansı (UV) ulnar fovea ile yarıçap üzerindeki sigmoid çentik arasındaki mesafe olarak tanımlayın.
3. Artikülasyon Modeli
   1. Bilek kinematiği ve bileşene özgü mesafe alanlarını²⁹ kullanarak polietilen kapak ile radyal bileşen arasındaki temas modellerini hesaplayın.
   2. Her 3B yüzey modeli için uzaklık alanlarını, yüzeyden işaretli mesafelerin hacimsel bir dizisi olarak hesaplayın.
   3. Her duruş için polietilen kapağın ve radyal bileşenin yüzeyindeki mesafe alanlarının yakınlık değerlerini hesaplayın.
   4. Bileşenler arasındaki temas yamasını elde etmek için, bu yakınlık değerlerini alım sisteminin çözünürlüğüne göre ayarlayın (bu çalışmada mesafe dışlama eşiği 0,70 mm olarak seçilmiştir).
4. İstatistiksel Analiz
   1. BVR ve OMC^{14, 30} arasındaki yanlılığı (ortalama farklar) ve %95 uyum sınırını (sapma ± 1.96SD) göstererek doğruluğu değerlendirmek için Bland-Altman analizini kullanın.
   2. Sağlıklı bileğin dönme merkezi davranışını tanımlamak ve UV ile PS'deki değişim arasında bir ilişki bulmak için genelleştirilmiş doğrusal modeller (GLM) kullanın (UV = p₁×PS² + p₂×PS), burada p₁ ve p₂ , doğrusal regresyon kullanılarak optimize edilmiş katsayılardır. Modelin doğruluğunu değerlendirmek için kök-ortalama-kare-hata (RMSE) ölçümünü kullanın.

Sonuçlar

Model tabanlı izleme için 2B'den 3B'ye görüntü kayıt yazılımının seçimi, kısmen grafik işlemci birimine (GPU) ve yüksek performanslı bilgi işlem (HPC) sistemlerine erişime bağlıdır. Bu programların farklı boru hatları vardır ve şu an itibariyle programlar arasında ortak bir metodoloji yoktur. Bu çalışmada, Brown Üniversitesi^25'te geliştirilen açık kaynaklı bir 2D'den 3D'ye görüntü kayıt programı olan Autoscoper'ı kullanıyor...

Tartışmalar

Biplanar videoradyografi (BVR), el bileği ve distal radioulnar eklemdeki kemik ve implant hareketini milimetre altı ve alt derece doğrulukla ölçmek için kullanılabilen görüntü tabanlı bir yöntemdir. Burada anlattığımız çalışmalarda, sağlıklı bir bilek için öngörülen COR'un doğru bir modelini ve TWA temas modellerini belirlemek için BVR kullanıldı. Bu tür bulgular, yeni nesil total bilek protezlerinin tasarımını bilgilendirebilir ve modellerin hesaplanmas...

Malzemeler

Name	Company	Catalog Number	Comments
3D Surface Scanner	Artec 3D	Artec Space SpiderTM	Luxembourg
Autoscoper	Brown University	https://simtk.org/projects/autoscoper	https://doi.org/10.1016/j.jbiomech.2019.05.040
CT Scanner	General Electric (GE)	Lightspeed 16	Milwaukee, WI, USA
Geomagic Wrap 3D	3DSystems	Version 2017	Rock Hill, SC, USA
Graphics Processing Unit (GPU)	Nvidia	GeForce GTX 1080	CUDA-enabled GPU
High-speed Video Cameras	Phantom	Version 10	Vision Research, Wayne, NJ, USA
Image Intensifier	Dunlee	40 cm diameter	Aurora, IL, USA
ImageJ	Open-source (Brown University)	https://imagej.net/Fiji	https://doi.org/10.1038/nmeth.2019
Matlab	The MathWorks, Inc.	R2017a to R2020a	Natick, MA, USA
Mimics	Materialise	Version 19.0 to 22.0	Leuven, Belgium
Motion Capture Cameras	Qualisys	Oqus 5+	Gothenburg, Sweden
Pulsed X-ray Generators	EMD Technologies	EPS 45–80	Saint-Eustache, Quebec, QC, Canada
Undistortion Grid	McMaster-Carr	9255T641	Steel Perforated Sheet Staggered Holes, 0.048" Thk, 0.125" Hole Dia, 36" X 40"
Wrist Implant (In-vitro Study)	Integra LifeSciences	Universal 2	Plainsboro, NJ, USA
Wrist Implant (In-vivo Study)	Integra LifeSciences	Freedom	Plainsboro, NJ, USA
WristViz	Open-source (Brown University)	https://github.com/DavidLaidlaw/WristVisualizer/tree/master	Open-source software
X-ray Tubes	Varian Medical Systems	Model G-1086	Palo Alto, CA, USA
XMALab	Open-source (Brown University)	https://www.xromm.org/xmalab/	https://doi.org/10.1242/jeb.145383