Badanie to zostało zatwierdzone przez Komitet Etyki Medycznej, Wydział Stomatologii, Universiti Malaya, Kuala Lumpur (DF OS1921/0082(P)). Ponieważ badania obejmowały dane z retrospektywnej i anonimizowanej tomografii komputerowej (CT) bez użycia tkanek biologicznych lub płynów, nie wymagała się zgody pacjenta. Metodologia tego badania opierała się na procedurach określonych przez Hariri i in.9 , z kilkoma modyfikacjami. Używany sprzęt i oprogramowanie są wymienione w Tabeli Materiałów.
1. Wybór obiektu i pobieranie danych z tomografii komputerowej
Retrospektywne tomografie komputerowej czaszki pacjentów pediatrycznych z diagnozą syniosynostozy czaszki (SC) zostały pobrane z cyfrowego archiwum obrazowego Kliniki Czaszkowo-Twarzowej w University Malaya Medical Centre (UMMC), obejmujące okres od listopada 2015 do grudnia 2022 roku. Dodatkowo uwzględniono prospektywne tomografie komputerowe uzyskane między majem 2023 a styczniem 2024 roku.
Potencjalni uczestnicy byli identyfikowani na podstawie dokumentacji Kliniki Czaszkowo-Twarzowej i ręcznie przesiewani przez badacza zgodnie z wcześniej ustalonymi kryteriami kwalifikacyjnymi. Kryteria włączenia obejmowały: (1) potwierdzoną diagnozę syndromicznej kraniosynostozy przez chirurgów czaszkowo-twarzowych, (2) dostępność pełnej dokumentacji medycznej oraz (3) wiek <12 lat. Osoby kontrolne również miały mniej niż 12 lat i wymagały wykonania kompletnych tomografii czaszki i twarzy bez wypadłych anomalii czaszkowo-twarzowych. Kryteria wykluczające stosowane dla wszystkich kohort obejmowały wiek ≥12 lat, niesyndromową lub izolowaną kraniosynostozę, niepełną dokumentację kliniczną lub obrazową, historię operacji czaszkowo-twarzowych lub hipoplazję śródtwarzy związaną z innymi zespołami. Po badaniu ostateczna lista kwalifikujących się uczestników została przekazana Jednostce Badawczej Obrazowania Biomedycznego w celu uzyskania zestawów danych z niekontrastowej tomografii czaszki. Dane obrazowe były udostępniane na płytach kompaktowych (CD) w formacie Digital Imaging and Communications in Medicine (DICOM), z protokołami akwizycji standaryzowanymi zgodnie z wytycznymi działu (grubość osi ≤1,00 mm, pełne pokrycie czaszkowo-twarzowe oraz spójna rozdzielczość wokseli). Jednostka obrazowania biomedycznego zweryfikowała kompletność zbioru danych, brak istotnych artefaktów oraz zgodność ze standardami DICOM przed wydaniem. Wszystkie zweryfikowane zbiory danych zostały bezpiecznie zarchiwizowane i następnie przygotowane do importu do oprogramowania do 3D medycznej obróbki obrazów w celu trójwymiarowej rekonstrukcji i analizy. Wybrano łącznie trzydzieści uczestników, podzielonych równomiernie na trzy kohorty: nieoperowaną grupę SC (SCNO), grupę operacyjną SC (SCO) oraz grupę kontrolną (n = 10 na grupę).
2. 3D rekonstrukcja i identyfikacja zabytków
Wszystkie importowane zbiory danych były wyświetlane w płaszczyźnie osiowej, strzałkowej i koronalnej, wraz z trójwymiarowym podglądem. Metadane w plikach DICOM umożliwiały automatyczne wyodrębnianie parametrów akwizycji (grubość fragmentów, rozdzielczość wokseli i liczba fragmentów). Aby wygenerować trójwymiarowy model objętościowy, najpierw stworzono maskę kością, a następnie ją zrekonstruowano, wybierając Segmentację → Obliczanie 3D w Mimikach. Segmentacja została dodatkowo udoskonalona za pomocą funkcji thresholding, dostępnej przez Segmentation → Thresholding na głównym pasku narzędzi. Zmienność definicji kości lub obecność nadmiernego szumu można zaobserwować, gdy podczas segmentacji zastosowano nieodpowiednie zakresy Jednostki Hounsfield (HU). Problem ten można rozwiązać poprzez zastosowanie standardowego progu "Kość (CT)" dostępnego w oprogramowaniu obrazowania medycznego 3D oraz dopracowanie masek segmentacyjnych za pomocą funkcji Edit Mask , aby wykluczyć artefakty i poprawić ciągłość strukturalną. Proces ten zapewniał, że w masce podświetlano tylko struktury kostne, wykluczając sąsiednie tkanki miękkie i zmniejszając szumy obrazowe lub artefakty.
Anatomiczne punkty orientacyjne były identyfikowane i oznaczane za pomocą narzędzia do tworzenia punktów, znajdującego się w sekcji Pomiary → Utworzenie punktu. Dla każdego punktu orientacyjnego badacz obracał i powiększał rekonstrukcję 3D, aby zoptymalizować widoczność, a następnie umieszczał punkt bezpośrednio na stanowisku anatomicznym jednym kliknięciem myszy. Zarejestrowano następujące punkty orientacyjne: sella (S) w środku sella turcica, Nasion (N) na styku szwu czołowo-nosowego, Basion (Ba), definiowany jako przedni środek kości potylicznej w synhondrosis klinowo-potylicznej, oraz prawy i lewy szew jarzmako-szczękowy (ZMR, ZML) na krawędzi podoczodołu. Drobna zmienność w pozycjonowaniu punktów orientacyjnych może wystąpić, gdy punkty są definiowane z jednej perspektywy obserwacyjnej. Aby zwiększyć dokładność, każdy punkt orientacyjny był potwierdzany pod wieloma kątami widzenia poprzez obracanie trójwymiarowej rekonstrukcji i porównanie z płaszczyznami osiowymi, strzałkowymi i koronalnymi przed ostatecznym umieszczeniem. Każdy punkt orientacyjny był automatycznie przechowywany w drzewie projektu pod listą "Pomiary". Pełny zestaw współrzędnych i etykiet punktów orientacyjnych został wyeksportowany za pomocą funkcji Eksport pomiarów , tworząc zbiory danych w formacie .csv lub .xls. Te punkty orientacyjne zostały następnie wykorzystane jako punkty odniesienia do liniowych pomiarów podstawy czaszki i szkieletu śródtwarzy.
3. Liniowe pomiary zmiennych czaszkowych i śródtwarzowych
Liniowe pomiary czaszkowo-twarzowe wykonywano w oprogramowaniu do przetwarzania obrazów medycznych 3D, wykorzystując funkcję pomiaru odległości. Po trójwymiarowej rekonstrukcji i umieszczeniu punktów orientacyjnych narzędzie pomiarowe było dostępne przez menu Analiza → Pomiary → Odległość między punktami . W tej funkcji podczas rekonstrukcji 3D ręcznie wybrano dwa punkty orientacyjne anatomiczne, a oprogramowanie automatycznie obliczało euklidesową odległość między wybranymi współrzędnymi, a wyniki były wyświetlane w oknie Wyniki pomiarów.
Uzyskano następujące pomiary liniowe: SN (przednia długość nasady czaszki), zdefiniowana jako odległość między sellą (S) a nosem (N); SBa (tylna długość podstawy czaszki), definiowana jako odległość między sellą (S) a basionem (Ba); NBa (całkowita długość podstawy czaszki), definiowana jako odległość między nasionem (N) a basionem (Ba); oraz ZMR - ZML (szerokość szczęki), definiowana jako odległość międzyzygomatyczna między prawym a lewym szwem jarzmatematyczno-szczękowym. Każdy pomiar wykonywano bezpośrednio na modelu 3D, stosując narzędzia do zoomu i rotacji w razie potrzeby, aby zapewnić dokładne rozmieszczenie punktów.
Po zakończeniu lista pomiarów była automatycznie aktualizowana w drzewie projektu oprogramowania. Cały zbiór danych został wyeksportowany za pomocą funkcji Eksport Pomiarów i zapisany w .xls formacie, uwzględniając wszystkie etykiety punktów orientacyjnych z odpowiadającymi im wartościami liniowymi. Te eksportowane wartości zostały następnie wykorzystane do analizy statystycznej.
4. Zastosowanie modelu regresji
Eksportowane pomiary zostały zaimportowane do programu arkusza kalkulacyjnego do dalszej analizy. Stosując wzór regresji Haririego-Ros-Nor, obliczono przewidywane wartości zarówno dla całkowitej długości podstawy czaszki, jak i szerokości szczęki. Zastosowano następujące równania:

Przewidywane wartości uzyskane z tych wzorów porównywano następnie z odpowiadającymi zmierzonymi wartościami dla każdego badanego, aby ocenić spójność i odchylenia morfologii czaszki poprzez obliczenie odchylenia standardowego zbioru danych. Następnie zbiór danych jest gotowy do analizy statystycznej.
5. Analiza statystyczna
Wszystkie zmierzone i przewidywane wartości zostały zebrane w jeden zbiór danych i przeanalizowane za pomocą oprogramowania do analizy statystycznej. Najpierw obliczono statystyki opisowe, aby podsumować parametry podstawy i środkowej części twarzy czaszki we wszystkich grupach. Różnice między grupami oceniano za pomocą nieparametrycznego testu Mann–Whitney U. Dodatkowo przeprowadzono analizę korelacji Pearsona, aby ocenić zależność między wiekiem a pomiarami czaszki w każdej grupie. Opracowano wizualne reprezentacje danych, w tym wykresy pudełkowe do porównań grup oraz wykresy liniowe dla trendów korelacji, aby wspierać i wyjaśniać interpretacje statystyczne.