Jak napisać pracę magisterską z fizyki medycznej i dozymetrii

Fizyka medyczna to jeden z nielicznych kierunków, gdzie praca magisterska musi spełniać wymagania jednocześnie naukowe i kliniczne. Nie wystarczy poprawna metodologia i przegląd literatury – w większości tematów trzeba przeprowadzić pomiary na sprzęcie produkcyjnym (akceleratory liniowe, skanery CT/MRI, detektory dozymetryczne) lub zaimplementować algorytm obliczeniowy i porównać go z wynikami systemu planowania leczenia (TPS). Oznacza to, że promotor, instytucja przyjmująca i sam student muszą wcześniej ustalić dostęp do aparatury, zakres pomiarów i obowiązujące protokoły. Ten artykuł pokazuje, jak zaplanować i napisać taką pracę krok po kroku.

Specyfika kierunku – co wchodzi w zakres fizyki medycznej

Fizyka medyczna (ang. medical physics) obejmuje zastosowanie metod fizycznych w diagnostyce i terapii. W praktyce dyplomowej najczęściej pojawiają się cztery obszary.

Radioterapia i planowanie leczenia. Tutaj mieszczą się prace o weryfikacji planów leczenia (Patient-Specific QA), optymalizacji dawek w technikach IMRT (Intensity-Modulated Radiation Therapy), VMAT (Volumetric Modulated Arc Therapy) i SBRT (Stereotactic Body Radiation Therapy), a także charakterystyce wiązek promieniowania produkowanych przez akceleratory liniowe (linac). Promotorzy z tego obszaru pracują zazwyczaj na uczelniach medycznych lub bezpośrednio w zakładach fizyki medycznej szpitali onkologicznych.

Dozymetria promieniowania jonizującego. Dozymetria zajmuje się pomiarami dawki pochłoniętej, mocy dawki i rozkładów dawki w fantomach i tkankach. Prace mogą dotyczyć kalibracji detektorów, porównania różnych typów detektorów, dozymetrii in vivo z użyciem diod, filmów dozymetrycznych (Gafchromic EBT3, EBT-XD) lub termoluminescencyjnych dozymetrów (TLD, LiF:Mg,Ti).

Obrazowanie medyczne: CT, MRI, PET, SPECT, USG. Prace z tego obszaru dotyczą oceny jakości obrazu (NPS, MTF, CNR, SNR), optymalizacji protokołów dawkowania, rekonstrukcji obrazu, analizy artefaktów lub fuzji modalności (np. PET/CT, MRI/PET). Są to tematy interdyscyplinarne – wymagają znajomości fizyki detektorów i podstaw przetwarzania sygnałów.

Ochrona radiologiczna. Praca może dotyczyć szacowania dawek dla personelu i pacjentów, badania skuteczności osłon, analizy ekspozycji przy zabiegach radiologii interwencyjnej lub oceny narażenia zawodowego.

Instytucje i środowisko naukowe

Najważniejsze ośrodki realizujące tematy dyplomowe z fizyki medycznej w Polsce:

• Centrum Onkologii im. Marii Skłodowskiej-Curie (CO-I) z oddziałami w Warszawie, Krakowie, Gliwicach i Gdańsku – największy pracodawca fizyków medycznych w Polsce, regularnie przyjmuje studentów na staże i pomiary do prac magisterskich
• Centrum Medyczne Kształcenia Podyplomowego (CMKP) – organizuje specjalizację z fizyki medycznej, wydaje dokumenty potwierdzające staż
• Zakłady fizyki medycznej przy szpitalach onkologicznych i radioterapeutycznych (m.in. Szpital Uniwersytecki w Krakowie, Klinikum Bochum, WCO Poznań)
• Polskie Towarzystwo Fizyki Medycznej (PTFM) – organizacja zrzeszająca fizyków medycznych, publikuje standardy kliniczne i organizuje konferencje

Certyfikacja zawodowa fizyka medycznego w Polsce przebiega przez ścieżkę: tytuł magistra + co najmniej 4 lata specjalizacji klinicznej + egzamin państwowy, zgodnie z rozporządzeniem Ministra Zdrowia w sprawie specjalizacji w dziedzinie fizyki medycznej. Tytuł PY (Physicien Medical) lub MPhys uznawany jest przez EFOMP (European Federation of Organisations for Medical Physics).

Tematy prac – konkretne przykłady

Poniżej sześć tematów, które faktycznie pojawiają się w polskich i europejskich uczelniach. Każdy z nich wymaga dostępu do aparatury klinicznej lub oprogramowania TPS.

1. Weryfikacja planów radioterapii techniką IMRT przy użyciu systemu ArcCHECK i oprogramowania SNC Patient

Student mierzy dawkę pochłoniętą za pomocą cylindrycznego detektora matrycowego ArcCHECK (Sun Nuclear Corporation) i porównuje wyniki z dawką obliczoną w TPS metodą Gamma Analysis. Kryterium akceptacji: 3%/3 mm lub 2%/2 mm (zależnie od protokołu ośrodka). Praca wymaga dostępu do akceleratora i systemu Eclipse lub Monaco.

2. Analiza dawki pochłoniętej w tkankach krytycznych podczas radioterapii raka szyjki macicy metodą VMAT

Celem jest ocena ryzyka powikłań popromiennych: dawka na pęcherz moczowy (V40Gy, V70Gy) i odbytnę (V50Gy, D2cc) według zaleceń GEC-ESTRO. Praca analizuje DVH (Dose-Volume Histogram) dla serii planów i porównuje je z danymi klinicznymi lub literaturą.

3. Ocena jakości obrazu CT do planowania radioterapii

Student mierzy parametry jakości obrazu (NPS – Noise Power Spectrum, MTF – Modulation Transfer Function, CNR – Contrast-to-Noise Ratio) dla różnych protokołów dawkowania CT na skanerze stosowanym do planowania. Cel: znalezienie kompromisu między dawką promieniowania a jakością obrazu wystarczającą do delineacji struktur.

4. Charakterystyka kliniczna akceleratora liniowego Elekta Versa HD

Praca polega na pomiarach Output Factors (OFs), PDD (Percentage Depth Dose) i profili dawek (Profiles) w fantomie wodnym (np. IBA Blue Phantom2) dla energii fotonów 6 MV i 10 MV FFF (Flattening Filter Free) oraz elektronów 6-18 MeV. Wyniki wchodzą do bazy danych TPS lub służą kontroli jakości aparatu.

5. Dozymetria in vivo z użyciem diod półprzewodnikowych podczas teleterapii

Celem jest weryfikacja, czy dawka podawana pacjentowi zgadza się z dawką planowaną. Student kalibruje diody krzemowe, mierzy dawkę wejściową i wyjściową na grupie pacjentów i analizuje odchylenia. Prace tego typu wprost wpisują się w wymogi QA (Quality Assurance) ośrodka.

6. Porównanie sekwencji MRI do delineacji tarczówki i narządów krytycznych w raku prostaty

Student porównuje, jak różne sekwencje MRI (T1w – T1-weighted, T2w – T2-weighted, DWI – Diffusion-Weighted Imaging, DCE – Dynamic Contrast Enhanced) wpływają na precyzję wyznaczania GTV (Gross Tumor Volume), CTV (Clinical Target Volume) i OAR (Organs At Risk). Analiza obejmuje inter-observer variability (zgodność między lekarze delineującymi tę samą strukturę na różnych sekwencjach).

Metodologia – sprzęt, detektory, oprogramowanie

Wybór metodologii zależy od tematu, ale kilka elementów pojawia się w prawie każdej pracy dozymetrycznej.

Fantomy dozymetryczne

Fantomat wodny (water phantom) to podstawowe narzędzie do charakterystyki wiązek promieniowania. W pracy magisterskiej najczęściej używa się:

• IBA Blue Phantom2 – trójwymiarowy skaner wodny do pomiarów PDD i profili dawek
• CIRS (Computerized Imaging Reference Systems) – fantomy tkankopodobne (tissue-equivalent) do pomiarów heterogenicznych
• Gammex RMI – fantomy do CT i MRI, wzorce gęstości elektronowej

Do pomiarów in vivo lub na fantomach stałych stosuje się detektory matrycowe: MatriXX Evolution (IBA), ArcCHECK (Sun Nuclear), Delta4 (ScandiDos).

Detektory jonizujące i inne

Pomiary dawki absolutnej (absolute dosimetry) wymagają kalibrowanej komory jonizacyjnej:

• Farmer chamber (np. FC65-G, IBA) – do pomiarów referencyjnych w wiązkach fotonów i elektronów
• Semiflex 3D (PTW Freiburg) – do skanowania profili w fantomie wodnym
• PinPoint (PTW) – do pomiarów małych pól (small fields, < 3×3 cm)

Oprócz komór jonizacyjnych stosuje się:

• Diody krzemowe (np. IBA EDP-30) – do dozymetrii in vivo i small fields
• Filmy dozymetryczne Gafchromic EBT3 i EBT-XD – wysoka rozdzielczość przestrzenna, czytane za pomocą skanera flatbed (Epson Expression 10000XL + FilmQA Pro)
• TLD (termoluminescencyjne dozymetry, LiF:Mg,Ti) – do pomiarów dawki na personel i pacjentów, odczytywane czytnikiem TLD (np. Harshaw 5500)
• MOSFET (Metal-Oxide-Semiconductor Field-Effect Transistor) – miniaturowe detektory in vivo

Systemy planowania leczenia (TPS)

Wszystkie obliczenia dawki przed leczeniem wykonuje TPS. W Polsce dominują trzy platformy:

• Eclipse (Varian Medical Systems) – z algorytmami AAA (Analytic Anisotropic Algorithm) i Acuros XB
• Monaco (Elekta) – algorytm Monte Carlo (XVMC)
• RayStation (RaySearch Laboratories) – obsługuje wiele algorytmów, w tym Monte Carlo

Do prac stricte obliczeniowych lub walidacyjnych używa się niezależnych kodów Monte Carlo: EGSnrc (National Research Council Canada) lub Geant4 (CERN). Wymagają one modelu głowicy akceleratora (phase space file) i znaczącej mocy obliczeniowej.

Analiza statystyczna i miary jakości

W pracy dozymetrycznej pojawiają się następujące miary i analizy:

• Gamma Analysis – metoda porównania dwóch rozkładów dawki, kryterium: Delta D (procentowa różnica dawki) i Delta r (odległość od zgodności, Distance to Agreement). Standardowe kryterium: 3%/3 mm; bardziej rygorystyczne: 2%/2 mm lub 1%/1 mm dla SBRT
• DVH (Dose-Volume Histogram) – zbiorcze zestawienie dawki w danym narządzie lub tarczówce. DVH kumulatywny (C-DVH) i różniczkowy (D-DVH) służą do różnych analiz
• Parametry DVH: V95% (objętość tarczówki pokryta co najmniej 95% dawki przepisanej), D98% (dawka pokrywająca co najmniej 98% objętości), Dmean (dawka średnia), D2cc (dawka na 2 cm3 danego narządu).

Niepewność pomiarową należy podawać zgodnie z GUM (ang. Guide to the Expression of Uncertainty in Measurement, ISO/IEC Guide 98-3). Praca powinna zawierać budżet niepewności: niepewność standardowa kalibracji (typ B), niepewność statystyczna pomiarów (typ A) i łączna niepewność rozszerzona U (k=2, poziom ufności 95%).

Normy i protokoły – co cytować w pracy

Dobra praca z dozymetrii musi powołać się na właściwe dokumenty normalizacyjne. Poniżej lista z opisem zastosowania.

Normy IAEA:

• IAEA TRS-398 (Technical Reports Series No. 398, „Absorbed Dose Determination in External Beam Radiotherapy”) – protokół dozymetrii absolutnej wiązek fotonów i elektronów z akceleratorów liniowych. Podstawa każdego pomiaru dawki bezwzględnej w fantomie wodnym. Konieczny przy obliczaniu Output Factors i PDD
• ICRU Report 83 („Prescribing, Recording, and Reporting Photon-Beam Intensity-Modulated Radiation Therapy (IMRT)”) – definiuje jak opisywać dawkę i objętości leczenia w IMRT; każda praca o planach IMRT musi się na niego powołać

Protokoły AAPM (American Association of Physicists in Medicine):

• TG-51 (Task Group 51, „Protocol for Clinical Reference Dosimetry of High-Energy Photon and Electron Beams”) – dozymetria absolutna, alternatywa dla TRS-398 stosowana w Ameryce Północnej, cytowana także w pracach porównawczych
• TG-142 („Quality Assurance of Medical Accelerators”) – wymagania dotyczące kontroli jakości akceleratorów liniowych: dzienne, tygodniowe i miesięczne testy dozymetryczne i mechaniczne

Prawo polskie i dyrektywy EU:

• Rozporządzenie Ministra Zdrowia z dnia 27 marca 2008 r. w sprawie minimalnych wymagań dla jednostek ochrony zdrowia ubiegających się o zawarcie umów na udzielanie świadczeń zdrowotnych w zakresie radioterapii onkologicznej (Dz.U. 2008 nr 59 poz. 365) – reguluje wymagania sprzętowe i kadrowe zakładów radioterapii
• Ustawa z dnia 29 listopada 2000 r. Prawo atomowe (Dz.U. 2000 nr 116 poz. 1216 ze zm.) – podstawa prawna nadzoru nad źródłami promieniowania jonizującego w Polsce. Reguluje obowiązki jednostki organizacyjnej wykonującej działalność z promieniowaniem
• Dyrektywa EURATOM 2013/59/Euratom z dnia 5 grudnia 2013 r. ustanawiająca podstawowe normy bezpieczeństwa w dziedzinie ochrony przed promieniowaniem jonizującym – implementowana do prawa polskiego, reguluje dawki graniczne dla personelu i pacjentów

Normy PN-EN:

• PN-EN 62083:2009 (IEC 62083) – wymagania dotyczące bezpieczeństwa systemów planowania leczenia (TPS) stosowanych w radioterapii. Cytuj przy opisie walidacji oprogramowania TPS
• PN-EN IEC 60601-2-1 – wymagania szczegółowe dla bezpieczeństwa akceleratorów cząstek stosowanych w radioterapii

Dokumenty PTFM:

PTFM (Polskie Towarzystwo Fizyki Medycznej) wydaje wytyczne kliniczne adaptujące normy IAEA i AAPM do warunków polskich. Warto sprawdzić najnowsze zalecenia przed złożeniem pracy.

Struktura pracy – jak podzielić rozdziały

Praca magisterska z fizyki medycznej ma zazwyczaj 60-90 stron (bez załączników). Promotorzy w CO-I i na wydziałach fizyki oczekują następującego układu:

Wstęp (5-10% objętości). Uzasadnienie tematu: dlaczego ten problem jest ważny klinicznie lub metodologicznie. Cytuj dane epidemiologiczne (np. liczba pacjentów leczonych radioterapią rocznie w Polsce, źródło: ROIS), statystyki sprzętowe (IAEA DIRAC) i lukę badawczą (gap), którą praca wypełnia. Unikaj pisania wstępu jako ogólnego kursu z fizyki promieniowania – to błąd, który recenzenci określają jako „encyklopedyzm”.

Przegląd literatury (15-20% objętości). Omów stan wiedzy, cytując prace opublikowane w ciągu ostatnich 10 lat w czasopismach takich jak Medical Physics (AAPM), Radiotherapy and Oncology (ESTRO), Physics in Medicine and Biology (IOP). Sprawdź impact factor i wybieraj artykuły recenzowane. Bazy: PubMed, Scopus, Web of Science.

Materiał i metody (20-30% objętości). To serce pracy. Opisz sprzęt (z numerami seryjnymi), detektory (z certyfikatami kalibracji), protokoły pomiarowe (z powołaniem na TRS-398, TG-142 lub inny dokument normalizacyjny), oprogramowanie (z numerem wersji) i procedurę analizy danych. Każde zdanie w tej sekcji powinno być powtarzalne przez kogoś, kto ma dostęp do podobnego sprzętu.

Wyniki (20-25% objętości). Tylko dane, bez interpretacji. Tabele z wartościami liczbowymi (z niepewnościami), wykresy PDD, profili dawek, krzywych Gamma Analysis, DVH. Format: każda tabela i rycina z numerem, podpisem i opisem osi.

Dyskusja (15-20% objętości). Porównanie wyników z literaturą i wymaganiami protokołów. Tu tłumaczysz, dlaczego wyniki są takie a nie inne, co mogło wpłynąć na odchylenia od wartości referencyjnych, jakie są ograniczenia pracy.

Wnioski (2-3 strony). Krótkie, ponumerowane, odpowiadające wprost na pytania badawcze ze wstępu. Każdy wniosek powinien wynikać z danych – nie z ogólnej wiedzy o temacie.

Bibliografia. Format zazwyczaj Vancouver (numerowany) lub APA – sprawdź wymagania uczelni. PubMed generuje gotowe cytowania w formacie NLM, z którego można korzystać bezpośrednio.

Zarządzanie danymi pomiarowymi

Dane pomiarowe trzymaj w formacie czytelnym dla promotora i recenzenta. Spis praktycznych narzędzi:

• Surowe dane z detektorów matrycowych (MatriXX, ArcCHECK) eksportuj do CSV lub DICOM-RT Dose i przetwarzaj w Pythonie (pydicom, numpy)
• Wyniki skanowania profili z IBA OmniPro Accept lub BlueScannerX zapisuj z metadanymi (data, energia, rozmiar pola, SSD)
• Plany leczenia archiwizuj w formacie DICOM-RT Plan + RT Dose + RT Structure Set (DICOM-RT kompletu); bez tych plików analiza DVH jest niemożliwa
• Zdjęcia filmów dozymetrycznych skanuj przy rozdzielczości >= 150 dpi (dla EBT3 optymalnie 300-400 dpi), tryb 48-bit RGB, bez korekcji kolorów przez skaner

Rejestracja danych i instytucje zewnętrzne

W pracach zawierających dane kliniczne (parametry planów leczenia, dozymetria in vivo na pacjentach) pojawia się obowiązek powołania się na dane zbiorcze lub porównanie z rejestrem.

• IAEA DIRAC (DIRectory of RAdiotherapy Centres) – publiczny rejestr ośrodków radioterapii na świecie, zawiera dane o sprzęcie i liczbie pacjentów. Użyteczny do sekcji „wprowadzenie” i uzasadnienia tematu
• ROIS (Rejestr Ośrodków i Informacja Statystyczna CO-I) – dane epidemiologiczne dotyczące nowotworów i radioterapii w Polsce, prowadzone przez CO-I w Warszawie
• NCI (National Cancer Institute, USA) – statystyki epidemiologiczne, literatura kliniczna. Baza PubMed jest częścią NCI

Każde badanie na pacjentach lub danych z dokumentacji medycznej wymaga zgody Komisji Bioetycznej. W przypadku prac retrospektywnych (analiza gotowych planów leczenia bez kontaktu z pacjentem) część komisji zwalnia z wymogu pełnej zgody – sprawdź to z promotorem przed planowaniem metodologii.

Typowe błędy – i jak ich unikać

To jest sekcja, którą najchętniej pomijają studenci w recenzjach. Każdy z poniższych błędów pojawia się w pracach przynajmniej raz na pięć składanych dyplomów.

Błędna normalizacja dawki przy kalibracji komory jonizacyjnej. Protokół TRS-398 wymaga uwzględnienia współczynnika jakości wiązki kQ (beam quality correction factor), który zależy od energii i rodzaju komory. Pominięcie kQ lub użycie wartości z tabeli dla innej komory generuje błąd systematyczny kilku procent – co w dozymetrii jest wartością nieakceptowalną.

Użycie fantomu wodnego o złej geometrii dla małych pól (< 3×3 cm). Przy pomiarach małych pól promieniowania (np. dla technik stereotaktycznych) dochodzi do efektu braku równowagi elektronowej (lateral charge equilibrium, LCE). Fantomat wodny zaprojektowany do standardowych pól 10×10 cm nie nadaje się bez modyfikacji do tych pomiarów – wymagane są detektory o małej objętości (PinPoint) i specjalne procedury.

Brak pomiaru in-air przy obliczaniu Output Factors dla małych pól. Output Factor (OF) mierzy się jako stosunek dawki dla danego pola do dawki referencyjnej (10×10 cm). Przy małych polach pomiar bez korekt związanych z rozmiarem detektora i geometrią wiązki daje niepoprawne wyniki. Prace oparte na IAEA TRS-483 („Dosimetry of Small Static Fields Used in External Beam Radiotherapy”) powinny dokumentować procedurę korekt.

Nieuwzględnienie współczynnika rekombinacji ks przy wysokich mocach dawki. Akceleratory FFF (Flattening Filter Free, np. TrueBeam STx, Versa HD) podają dawkę z bardzo wysoką mocą (do 2400 MU/min vs. 600 MU/min w standardzie). Przy takich mocach jonizacji w komorze wzrasta rekombinacja jonów – współczynnik ks koryguje ten efekt. Jego pominięcie prowadzi do niedoszacowania dawki.

Mylenie DVH kumulatywnego z różniczkowym. DVH kumulatywny (C-DVH) pokazuje, jaka objętość otrzymuje dawkę co najmniej X Gy. DVH różniczkowy (D-DVH) pokazuje rozkład dawki w przedziałach. Parametry kliniczne (V40Gy, D2cc) odnoszą się do C-DVH. Błędna interpretacja prowadzi do błędnych wniosków o pokryciu dawką tarczówki i oszczędzeniu narządów.

Brak analizy niepewności pomiarowej. Prace ograniczające się do podania średniej i odchylenia standardowego bez pełnego budżetu niepewności zgodnie z GUM (ISO/IEC Guide 98-3) są recenzowane jako metodologicznie niekompletne. Budżet powinien zawierać: niepewność kalibracji detektora (certyfikat laboratorium wzorcującego), niepewność pozycjonowania detektora, niepewność statystyczną serii pomiarów oraz połączoną niepewność standardową i rozszerzoną.

Brak walidacji kodu Monte Carlo. Jeśli praca zawiera obliczenia MC, recenzent oczekuje porównania wyników z pomiarem referencyjnym lub danymi z literatury dla tych samych warunków geometrycznych. Samo uruchomienie EGSnrc lub Geant4 bez walidacji nie stanowi kompletnej metodologii.

FAQ

Natalia Witek-Dąbrowska

Redaktorka i korektorka z 20-letnim doświadczeniem

Specjalizuję się w korekcie prac naukowych i dyplomowych. Jeśli potrzebujesz pomocy z tekstem, napisz do mnie.