Metoda szacowania naprężeń kory kości udowej zwłok podczas badania złamań z wykorzystaniem cyfrowej korelacji obrazu

Podsumowanie

W tym protokole, naprężenia powierzchni kości udowej są szacowane podczas testów złamań przy użyciu techniki cyfrowej korelacji obrazu. Nowatorstwo metody polega na zastosowaniu stochastycznego wzoru plamek o wysokim kontraście na powierzchni kości udowej, starannie dobranym oświetleniu, szybkim przechwytywaniu wideo i cyfrowej analizie korelacji obrazu do obliczeń odkształceń.

Streszczenie

Ten protokół opisuje metodę wykorzystującą cyfrową korelację obrazu do oszacowania obciążenia kory mózgowej na podstawie szybkich obrazów wideo powierzchni kości udowej zwłok uzyskanych z testów mechanicznych. Ta metoda optyczna wymaga tekstury wielu kontrastujących ze sobą znaków powierniczych na jednolitym białym tle w celu dokładnego śledzenia deformacji powierzchni w miarę obciążania próbki. Bezpośrednio przed badaniem powierzchnię będącą przedmiotem zainteresowania w widoku z kamery maluje się białym podkładem na bazie wody i pozostawia do wyschnięcia na kilka minut. Następnie na białe tło ostrożnie nakłada się czarną farbę, zwracając szczególną uwagę na równomierny rozmiar i kształt kropel. Oświetlenie jest starannie zaprojektowane i ustawione tak, aby zapewnić optymalny kontrast tych znaków przy jednoczesnej minimalizacji odbić dzięki zastosowaniu filtrów. Obrazy uzyskano dzięki szybkiemu przechwytywaniu wideo z prędkością do 12 000 klatek/s. Kluczowe obrazy przed i w momencie pęknięcia są wyodrębniane, a deformacje są szacowane między kolejnymi klatkami w starannie dobranych oknach przesłuchań w określonym obszarze zainteresowania. Odkształcenia te są następnie wykorzystywane do obliczania czasowego odkształcenia powierzchniowego podczas testu pękania. Dane dotyczące odkształcenia są bardzo przydatne do identyfikacji inicjacji złamania w obrębie kości udowej oraz do ostatecznej walidacji modeli wytrzymałości na złamania bliższej części kości udowej pochodzących z ilościowej analizy elementów skończonych opartej na tomografii komputerowej (QCT/FEA).

Wprowadzenie

Cyfrowa korelacja obrazów (DIC) to metoda przetwarzania obrazu, która jest używana w obecnym protokole do szacowania pełnego odkształcenia powierzchni pola próbek do badań kości udowej zwłok na podstawie obrazów sekwencji czasowej uzyskanych podczas testów mechanicznego złamania. Technika ta została po raz pierwszy opracowana i zastosowana w eksperymentalnej analizie naprężeń w latach 80-tych XX wieku i doświadczyła gwałtownego wzrostu użycia w ostatnich latach^1,2,3. Ma kilka kluczowych zalet w porównaniu z bardziej tradycyjnymi podejściami do montażu czujników tensometrycznych na konstrukcji, w tym zwiększony rozkład przestrzenny pola odkształceniowego, drobniejsze długości mierników dzięki zwiększonej rozdzielczości kamery oraz unikanie problemów z przyczepnością lub zgodnością kleju tensometrycznego. Główną zaletą DIC dla tkanek biologicznych, takich jak kości, jest to, że można go stosować do nieregularnych geometrii składających się z wysoce niejednorodnych właściwości materiału^4,5. Jego główną wadą w porównaniu z tradycyjnymi metodami akwizycji odkształceń jest to, że wymaga drogich, szybkich kamer wideo o wystarczającej rozdzielczości do pomiaru obszaru zainteresowania, aby uzyskać wystarczające próbkowanie przestrzenne i czasowe w celu dokładnego oszacowania pól odkształceń.

Podstawowym zastosowaniem pól odkształceń czasowych uzyskanych z analizy DIC złamań kości jest walidacja szacunków odkształceń w modelach QCT/MES wytrzymałości kości udowej⁵. Taka walidacja jest przedmiotem zainteresowania wielu ortopedycznych grup badawczych, które głównie wykorzystują zdalne pomiary siły i przemieszczenia z ogniw obciążnikowych i przetworników przemieszczenia^6,7,8. Ponadto, analiza obrazu wzoru złamania po pęknięciu została połączona z tymi zdalnymi pomiarami jako dalsze środki walidacji modelu⁹. Niedawno metoda DIC została zastosowana do walidacji modelu MES pęknięć i propagacji pęknięć w bliższej części kości udowej¹⁰. Wykorzystując korelację szczepów między modelami a eksperymentami, uzyska się jeszcze większą pewność co do trafności modeli obliczeniowych proksymalnej kości udowej i jeszcze bardziej zbliży metodę diagnostyczną QCT/MES do zastosowania klinicznego.

Ta praca wyjaśnia szczegółowy protokół uwzględniający niezbędne kroki do analizy DIC w testach złamań bliższej części kości udowej. Procedura obejmowała etapy przygotowania kości polegające na rozpyleniu białej farby na powierzchni kości, a następnie nakrapianie czarnych plam na wysuszonej białej powierzchni kości, metody uzyskiwania obrazów o wystarczającej rozdzielczości przestrzennej i czasowej przy użyciu szybkich kamer wideo oraz proces i narzędzia, których użyliśmy do obliczenia pól odkształceń z tych obrazów. Wyjaśniliśmy również kilka zastrzeżeń, które mogą mieć wpływ na jakość pomiarów.

Protokół

Wszystkie eksperymenty zostały przeprowadzone za zgodą Institutional Review Board. Próbki uzyskano we współpracy z laboratoriami badań anatomicznych.

1. Przygotowanie próbek do testów

1. Rozmrażać kość udową w temperaturze pokojowej przez 24 godziny.
2. Gdy kość udowa czeka w kolejce do testów, usuń okład, który został nałożony przed zamrożeniem i wytrzyj kość udową suchym ręcznikiem, aby usunąć resztki wilgoci, złogi tłuszczu lub tkanki miękkie. Wbij krętarz większy do prefabrykowanego aluminiowego kubka z cementem kostnym.
3. Używając pudełka, aby zatrzymać cząstki w jak największym stopniu, spryskaj kość białym plastikowym podkładem, aby uzyskać cienką, jednolitą powłokę. Należy pamiętać, aby pokryć kość jedną jednolitą warstwą farby, aby uzyskać optymalny kontrast i silne przyleganie do powierzchni kości udowej.
   UWAGA: Grubość nie została zmierzona.
4. Pozostaw farbę do wyschnięcia na co najmniej 5 minut. Jest to ważne, aby uniknąć niezamierzonego zmieszania z nakrapianymi kropelkami w sekcji 2.
5. Owiń kość wilgotną szmatką, aby uniknąć wysuszenia tkanki.

2. Proces nakrapiania

1. W przybliżeniu dodaj 1 część wody do 2 części farby, aby uzyskać najlepsze plamki na kości. Stopniowo dodawaj wodę (dla lepszej mieszanki) do farby akrylowej, aby uzyskać kolor.
2. Zanurz czystą szczoteczkę do zębów w palecie czarnej farby, aby wchłonąć kolor i przesuń szczoteczką, aby uzyskać czarne plamki na białej powłoce.
3. Pozostaw farbę do wyschnięcia na 5 minut przed kontynuowaniem.

3. Akwizycja obrazu

1. Zamontuj przygotowaną nakrapianą kość udową w mechanicznej maszynie testowej, wkładając zatkany dystalny koniec do uchwytu i dokręć dwie, aby zabezpieczyć próbkę.
2. Wyreguluj dwa odbłyśniki światła wyładowczego o dużej intensywności tak, aby powierzchnia kości udowej zapewniała jak najwyższe oświetlenie, unikając jednocześnie odbić na obrazie z kamery. Postępuj szybko, wykonując następujące czynności opisane w sekcji 3, aby uniknąć niepożądanego nagrzania próbki przez promieniowanie świetlne przed badaniem.
3. Zmniejsz przysłonę obiektywów szybkiej kamery wideo do oglądania z przodu i z tyłu, tak aby cały obszar zainteresowania kości udowej w polu view był ostry.
4. Wyreguluj odbłyśniki światła, aby jeszcze bardziej poprawić oświetlenie przy jednoczesnym zmniejszeniu odblasków.
5. Ustaw oprogramowanie do akwizycji obrazu tak, aby rejestrowało obraz z prędkością 6000 klatek/s w rozdzielczości 1024 x 512 pikseli. Ustawia łączną liczbę ramek, które mają zostać pobrane po odebraniu sygnału wyzwalającego na 12288. Uzbrój szybkie oprogramowanie wideo do akwizycji obrazu, gdy sygnał wyzwalający jest odbierany z systemu testowego. Po zakończeniu testowania obraz wideo znajduje się w pamięci bufora kamery.
6. Korzystając z oprogramowania do akwizycji obrazów, zapisz wideo na dysku, określając żądaną ścieżkę i nazwę pliku, a następnie klikając "Zapisz". Przygotuj się na czekanie od 5 do 40 minut na zakończenie tego procesu, w zależności od liczby ramek do zapisania.

4. Przygotowanie obrazu

1. Utwórz osobne katalogi robocze dla widoków z przodu i z tyłu kości udowej.
2. Użyj oprogramowania do analizy wideo, aby otworzyć odpowiednie nagranie wideo z dużą prędkością i zanotuj numery referencyjne klatki kluczowej na 1) początku ruchu siłownika ramy obciążeniowej i 2) klatce bezpośrednio po zdarzeniu pęknięcia.
3. Aby próbkować w dół nieskompresowaną sekwencję obrazu TIFF z szybkiego wideo, otwórz i uruchom skrypt "mov_frames.m" w katalogu roboczym dla odpowiedniej strony kości udowej.
   1. W wyświetlonym oknie dialogowym wprowadź numer klatki końcowej określony w kroku 4.2 o rozmiarze kroku od 25 do 40. Kliknij "Wyodrębnij ramki" i sprawdź katalog roboczy, aby upewnić się, że pliki *.tiff zostały poprawnie wyodrębnione.

5. Tworzenie siatki elementów skończonych

1. Użyj zewnętrznego programu do tworzenia siatki elementów skończonych, aby utworzyć siatkę elementów skończonych. Obliczanie odkształceń 2D na podstawie różnicowych wektorów przemieszczenia za pomocą metody elementów skończonych. Zaimportuj początkowy wyodrębniony obraz *.tiff do preprocesora oprogramowania elementów skończonych jako szablon do tworzenia splajnu.
2. Znajdź dwa łatwe do zidentyfikowania punkty odniesienia na obrazie, które znajdują się w przeciwległych rogach kadru i zapisz ich współrzędne X i Y (zostaną one ostatecznie użyte w kroku 6.1). Współrzędne te są dowolne w oparciu o konwencję używaną przez oprogramowanie MES do importowania obrazu *.tiff. Współrzędne tych punktów zostaną użyte do zarejestrowania węzłów siatki elementów skończonych z odpowiadającymi im pikselami obrazów wideo w kroku 6.2.
3. W oprogramowaniu do edycji obrazów otwórz ten sam obraz, który został zaimportowany do preprocesora oprogramowania elementów skończonych, i zapisz wartości kierunków X i Y pikseli skojarzonych z punktami określonymi w kroku 5.2. Zostaną one ostatecznie użyte w kroku 6.1.
4. W module "Szkic" programu do tworzenia siatki elementów skończonych użyj narzędzia splajn, aby obrysować zamknięty przekrój reprezentujący obszar zainteresowania. Sprawdź, czy obszar nie jest zbyt duży, aby powierzchnia kości przesunęła się poza obszar przed złamaniem z powodu rotacji.
5. Przygotuj zamknięty przekrój utworzony w kroku 5.4 do tworzenia siatki, obsiewając krawędzie globalnym rozmiarem siatki 1 mm w menu "Instancja części źródłowej".
6. W sekcji "Przypisz elementy sterujące siatką" ustaw kształt elementu na czworokątny.
7. Ułóż siatkę na zamkniętym odcinku.
8. Eksportuj siatkę do pliku ASCII bazy danych siatki składającego się ze współrzędnych węzłowych i definicji elementów.
9. Po otwarciu wynikowego pliku wejściowego elementów skończonych w edytorze tekstu skopiuj blok węzła zawierający numery węzłów i współrzędne do nowego pliku tekstowego i zapisz go jako "nodes.txt". Powtórz te czynności dla bloku elementu i zapisz nowy plik tekstowy jako "elements.txt".

6. Zarejestruj siatkę FE z szybkimi obrazami wideo i przeprowadź cyfrową analizę korelacji obrazu

1. W nowej sesji utwórz 2-elementowe wektory wierszy o nazwach ab1 i ab2 z wartościami określonymi w kroku 5.2 oraz px1 i px2 z wartościami określonymi w krokach 5.3, wpisując te nazwy wektorów w wierszu polecenia. Zapisz obszar roboczy jako "points.mat".
2. Uruchom skrypt "convert_imagesize.m", aby zarejestrować punkty z siatki elementów skończonych z wyodrębnionym obrazem wideo o dużej szybkości.
3. Uruchom skrypt "rrImageTrackGui.m". Załaduj pierwszy obraz ("p01.tif") i wprowadź numer ostatniego pliku *.tiff, który został wyodrębniony, jako całkowitą liczbę obrazów do przetworzenia.
4. Załaduj siatkę utworzoną w kroku 5.7, upewniając się, że opcja siatki jest ustawiona na "Czytaj z pliku" i kliknij "Akceptuj". Siatka elementów skończonych powinna pojawić się nad obrazem kości.
5. Określ wartości śledzenia w oparciu o poniższe wskazówki dotyczące parametrów śledzenia i kliknij "Kontynuuj" (pamiętając, że wartości parametrów to rozmiar obrazu, tekstura i stopień zachodzących deformacji, które należy dokładnie przetestować w każdym przypadku z osobna).
   1. Użyj początkowego rozmiaru jądra 21. Rozmiar jądra, n, to rozmiar okna n x n (gdzie n jest liczbą nieparzystą) pikseli, który jest używany do korelacji krzyżowej i określania wektora deformacji dla tego obszaru, który będzie używany do obliczeń odkształceń.
   2. Użyj początkowego rozmiaru subpiksela wynoszącego 4. Rozmiar subpiksela, m, to rozmiar podokna (2m+1) x (2m+1), na którym oblicza się deformację subpikseli, zakładając jednorodne odkształcenie w tym podoknie.
   3. Użyj początkowego współczynnika gładkości równego 2. Współczynnik gładkości to stopień wygładzenia zastosowany do pola przemieszczenia w śledzonych lokalizacjach przed obliczeniem odkształceń.
   4. Użyj początkowego współczynnika maxMove równego 10. Współczynnik maxMove to maksymalna liczba pikseli, o jaką dowolny węzeł może znajdować się w odległości od swojej trajektorii w stosunku do trajektorii sąsiada. Pomaga to uniknąć złego śledzenia do deformacji.
   5. Użyj początkowego współczynnika smoothGrid wynoszącego 15. Współczynnik smoothGrid to rozmiar siatki (nieco grubszy niż siatka śledzonych węzłów), który jest używany do wygładzania.
6. Zaznacz punkt pomocniczy, który ma wokół siebie znaczny kontrast, unikając jednocześnie obszarów z odblaskami lub rozmyciem. Sprawdź ten punkt, klikając "Sprawdź przewodnik" i sprawdź, czy pik korelacji jest silny (co najmniej dwukrotnie większą amplitudę) w porównaniu do swoich sąsiadów. Kliknij "Akceptuj" i "Wykonaj śledzenie", gdy jesteś zadowolony. Może to być długotrwały proces obliczeniowy, w którym przemieszczenie różnicowe jest obliczane dla czasowej sekwencji obrazu.
7. Po wykonaniu kroku 6.6 kliknij "Animuj". Po zakończeniu animacji kliknij opcję "Write strains (post processing software)", wpisz *.exe, a następnie wybierz opcję writeStrainRR_simple.exe¹¹. To obliczy szczepy. Zamknij graficzny interfejs użytkownika.

7. Przetwarzanie końcowe danych o przemieszczeniach i odkształceniach

1. Aby uzyskać odkształcenie w stosunku do numeru ramki, uruchom polecenie "analyzeFailurePrecursor.m" z wiersza polecenia z argumentem wejściowym o rozmiarze kroku (wybierz 20–30). Szczyty będą oznaczać uszkodzenie kości, a największy szczyt będzie odpowiadał ramie w pobliżu globalnej niewydolności kości.
2. Aby utworzyć pliki filmowe szczepów, uruchom "makeMovies.m" w oknie poleceń z argumentami (numVars, endstep, flag).
   UWAGA: Argument liczbaZmienne jest zdefiniowany jako 1 - 3 to przemieszczenia, 4 - 6 to składowe odkształcenia _xx, _yy i _xy, 7 i 8 to dwie główne odkształcenia i odkształcenia von Misesa, a 9 to energia odkształcenia. Koniec argumentu jest ostatnią klatką, która ma zostać uwzględniona w filmie.
   1. Ustaw opcjonalną flagę argumentu na 1, aby tworzyć filmy tylko dla jednostki określonej dla argumentu numVars i na 0, aby tworzyć filmy wszystkich zmiennych.

8. Dostrajanie i udoskonalanie wyników

1. Jeśli śledzenie DIC dało słabe wyniki, takie jak nieciągłe pole odkształcenia, które wykracza poza założenia mechaniki kontinuum, określ, co się dzieje i dlaczego śledzenie kończy się niepowodzeniem. Powtórz sekcję 6, zwracając szczególną uwagę na dostosowanie parametrów śledzenia. Drugorzędną opcją może być powrót do oprogramowania elementów skończonych i stworzenie bardziej jednolitej i możliwie drobniejszej siatki.
2. Jeśli śledzenie DIC dało rozsądne wyniki, utwórz dokładniejszą serię obrazów dla DIC. Korzystając z numerów referencyjnych klatek kluczowych z kroku 4.2 i odpowiadającej im liczby klatek na sekundę wideo, zidentyfikuj odstępy między klatkami dla trzech różnych reżimów zainteresowania w teście pękania, mając na uwadze, że punkty nie powinny przesuwać się między klatkami o więcej niż 6 pikseli.
   UWAGA: W początkowym segmencie badania, gdy naprężenia w kości udowej narastają powoli, odstępy między ramami będą stosunkowo duże (na przykład dla szybkości przemieszczenia 100 mm/s odstęp między ramkami dla tej części wynosi 3333 μs). W przypadku pośredniej części badania bliżej ramy pękającej odkształcenie rośnie szybciej i potrzebne są mniejsze odstępy między ramami (1667 μs dla szybkości przemieszczenia 100 mm/s). W przypadku ostatniej części tuż przed pęknięciem odstęp między ramami jest najmniejszy (16,7 μs przy prędkości przemieszczenia 100 mm/s).
3. Opcjonalnie tylko do celów dokumentacyjnych: Korzystając z informacji z kroku 8.2, utwórz sformatowane wpisy danych w pliku ASCII o nazwie "steps.txt" zawierającym wiersz danych dla każdego odstępu między klatkami. Format każdego wiersza będzie początkową klatką tego reżimu oddzieloną liczbą klatek do pominięcia (na podstawie kroku 8.2), oddzieloną dwukropkiem, a następnie końcową klatką tego reżimu (tj. format "1:20:200" poinstruuje oprogramowanie do ekstrakcji, aby wyodrębnić klatkę 1 do klatki 200 w krokach co 20).
   1. Bezpośrednio po tym oznaczeniu wstaw zakładkę i oznacz zakres numerów wyodrębniania obrazu (w przykładzie "1:20:200" pełne oznaczenie wiersza to "1:20:200 1:11" bez cudzysłowów). Powtórz te czynności dla pozostałych dwóch reżimów testowych, aby w pliku "steps.txt" znajdowały się trzy wiersze informacji. Ten plik służy jako zapis tego, w jaki sposób obrazy zostały wyodrębnione z oryginalnego wideo o dużej szybkości.
4. Uruchom ponownie kod mov_Frames.m, tym razem określając wiele reżimów odstępów między ramkami w oknie dialogowym. Wprowadź numery klatek i rozmiary kroków określone w kroku 8.2, aby określić parametry początkowe, końcowe i pomijania klatek oczekiwane w oknie dialogowym narzędzia. Upewnij się, że robisz to w nowym katalogu, w przeciwnym razie oryginalne obrazy zostaną nadpisane.
5. Powtórz sekcje 6 i 7 i przejrzyj wyniki w celu poprawy. Każda kość udowa może wymagać różnych dodatkowych iteracji w zależności od charakteru złamania, wzoru plamek i oświetlenia. Powtarzając krok 6.5, zachowaj te same ustawienia, z wyjątkiem zmniejszenia maxMove do 6 (z 10).

Wyniki

Przed procesem nakładania wzoru kropkowego kość udową oczyszcza się z nadmiaru tłuszczu i tkanek miękkich, a krętarz większy osadza się w aluminiowym kubku. Podczas twardnienia polimetakrylanu metylu (PMMA) kość owija się tkaniną nasączoną solą fizjologiczną, aby zapobiec wysuszeniu tkanki. Po utwardzeniu PMMA kość ponownie oczyszcza się bezpośrednio przed natryskiwaniem (Rysunek 1). Następnie powierzchnię kości natryskuje się lub maluje pędz...

Dyskusja

Wprowadziliśmy protokół konsekwentnego przygotowywania próbek kości udowej do obrazowania o wysokim kontraście podczas testów złamań, które następnie wykorzystano do oszacowania rozkładu odkształceń w pełnym polu za pomocą DIC. Protokół ten zapewniał odpowiedni kontrast tekstury czarnych plamek śledzących na jednolitym białym tle na powierzchni kości. Postępując zgodnie z tym protokołem, z powodzeniem powtórzyliśmy oszacowanie szczepów za pomocą analizy DIC dla osiemdziesięciu dziewięciu kości udowych.

Materiały

Lista materiałów użytych w tym artykule

Nazwa	Firma	Numer katalogowy	Komentarze
Krylon plastic primer biały	Krylon, Peoria, AZ, USA	N/A	Stosowany jako podkład do uzyskania gładkiego białego wykończenia na powierzchni kości
Farba akrylowa akrylowa na bazie wody w kolorze białym i czarnym	Plaid Enterprises (Ceramcoat), Norcross, GA, USA	N/A	Źródło farby dla kolorów białego i czarnego
Miska do mieszania	Niespecyficzne (ogólne)	Nie dotyczy	Służy do mieszania i przygotowania farby
Szczotka piankowa	Linzer Products, Wyandanch, NY, USA	N/A	Służy do nakładania farby na powierzchnię kości
Szczoteczka do zębów	Colgate-Palmolive, Nowy Jork, NY, USA	Twarde włosie	Służy do nakładania odpowiedniego rozmiaru i rozmieszczenia wzoru nakrapianego
Higieniczna żywica ortodontyczna (PMMA)	Patterson Dental, St Paul, MN, USA	H02252	Substancja kontrolowana i może być zakupiona za odpowiednią zgodą
Kenmore Freezer	Sears Holdings, Hoffman Estates, IL, USA	N/A	Służy do utrzymywania środowiska przechowywania -20oC dla próbek kości
Sól fizjologiczna (0,9% chlorek sodu)	Baxter Healthcare, Deerfield, IL, USA	NDC 0338-0048-04	Służy do utrzymywania nawodnionych próbek
Skalpele i skrobaki	Aspen Surgical (Bard-Parker), Kaledonia, MI, USA	Nie dotyczy	Służy do usuwania tkanek miękkich z próbek kości
Okap	dymowy Hamilton Laboratory Solutions, Manitowoc, WI, USA	70532	Używany do wentylacji podczas przygotowywania PMMA do zalewania próbek
Jednostki oświetleniowe	ARRI, Monachium, Niemcy	Nie dotyczy	Potrzebne do oświetlenia celu w celu przechwytywania obrazu
Szybka kamera wideo	Photron Inc., San Diego, CA, USA	Photron Fastcam APX-RS	Służy do przechwytywania szybkich nagrań wideo ze zdarzeń związanych ze złamaniem
Photron FASTCAM Imager and Viewer	Photron Inc., San Diego, CA, USA	Ver.3392(x64)	Służy do nagrywania i przeglądania szybkich nagrań wideo
Obiektyw kamery	Zeiss, Oberkochen, Niemcy	Obiektyw Zeiss Planar L4/50 ZF Potrzebny	do odpowiedniej rozdzielczości obrazu
ABAQUS CAE	Dassault Systemè s, Waltham, MA, USA	Wersje 6.13-4	Służy do definiowania obszaru zainteresowania i tworzenia siatki elementów skończonych
MATLAB	Mathworks, Natick, MA, USA	Wersja 2015b	Służy do przetwarzania obrazu i analizy DIC
TecPlot	TecPlot Inc., Bellevue, WA		Służy do przetwarzania końcowego pól odkształceń
Kalkulator odkształceń Oprogramowanie	Victor Barocas Research Group, University of Minnesota, Minneapolis, MN, USA	http://license.umn.edu/technologies/20130022_robust-image-correlation-based-strain-calculator-for-tissue-systems	Służy do obliczania pola odkształcenia
mov_frames.m	Skrypt Matlab, Mayo Clinic, Rochester, MN, USA	N/A	Służy do próbkowania w dół nieskompresowanych obrazów z szybkich plików wideo
convert_imagesize.m	Skrypt Matlab, Mayo Clinic, Rochester, MN,USA	N/A	Służy do rejestrowania współrzędnych pikseli obrazu ze współrzędnymi siatki
rrImageTrackGui.m	Skrypt Matlab, Mayo Clinic, Rochester, MN, USA	N/A	Służy do wykonywania korelacji krzyżowej obrazu w celu uzyskania deformacji i uruchomienia kalkulatora odkształceń
analyzeFailurePrecursor.m	Skrypt Matlab, Mayo Clinic, Rochester, MN, USA	N/A	Służy do czasowego śledzenia szczytowych
składników odkształceniaSkrypt makeMovies.m	Matlab, Mayo Clinic, Rochester, MN, USA	N/A	Służy do tworzenia przenośnych filmów *.avi przedstawiających składowe deformacji, składowe odkształceń, główne odkształcenia, odkształcenia von Misesa i energię odkształcenia