Kalibrace diagnostických displejů v radiologii (3)

Při pravidelném testování diagnostických displejů je potřeba proměřit alespoň několik bodů GSDF křivky, aby se dalo zjistit, zda funguje displej z hlediska jasu podle požadavků GSDF křivky. American Association of Physicists in Medicine (AAPM) ve svém starším doporučení TG18 doporučovala jako dostatečné pro QC proměření alespoň 18 bodů. Nicméně toto doporučení bylo primárně vytvořeno pro CRT displeje, nikoliv LCD displeje, které v dnešní praxi převažují. A právě pro LCD displeje není těchto 18 proměřených bodů dostatečných k tomu, aby se prokázal souhlas s GSDF křivkou. Ukázka nesouhlasu je uvedena na obr. 1.

Obr. 1: Výřez měření GSDF křivky pro prvních 20 hodnot [8]

V případě, že byl proměřen jas pouze prvního (nejtmavšího) a posledního (nejsvětlejšího) bodu, displej splnil požadavky GSDF křivky v rámci tolerance. Ale v případě, že byly změřeny i body mezi těmito okrajovými hodnotami, ukázalo se, že displej nesplňuje požadavky GSDF křivku v rámci tolerance, hlavně pro body 3-13. V novějším doporučení AAPM TG270 se proto již doporučuje proměřit alespoň 52 bodů.

Proveďme si nyní proměření displeje alespoň pro těch 18 bodů dle TG18.

Máme DICOM obrazy s hloubkou 12 bitů, tedy 2^12 = 4096 stupňů šedi, tedy 4096 P-hodnot. Obrazů máme pouze 18. Abychom těmito obrazy rovnoměrně pokryli interval P-hodnot, které jsou v rozsahu 0-4095, mají jednotlivé obrazy následující P-hodnoty: 0, 240, 480, 720, 960, 1200, 1440, 1680, 1920, 2160, 2400, 2640, 2880, 3120, 3360, 3600, 3840 a 4080. Ukázka prvního obrazu a pak každého druhého je uvedena na obr. 2.

Obr. 2: Obrazy dle TG18 pro kontrolu zobrazení dle GSDF křivky

Nejprve se změří maximální L_max a minimální L_min kalibrovaný jas. Nechť je L_min = 0,98 cd/m^2 a L_max = 535,2 cd/m^2 (reálné hodnoty odpovídající jednomu z našich displejů). Tomu dle vztahu (1) z předešlého příspěvku odpovídají hodnoty JND indexů JND_min = 70,67 a JND_max = 716,10. Pak JND_tot = 716,10 – 70,67 = 645,43. Nyní si stanovíme vztah mezi JND indexy a P-hodnotami pomocí JND_ave:

Pomocí vztahu (3) opět z předešlého příspěvku, který ale pro úplnost uvádím i níže, dopočítáme hodnoty JND indexů pro všech zbývajících 16 P-hodnot (první a poslední již máme).

Vztah (3) z předešlého příspěvku:

Pro jednotlivé P-hodnoty dostáváme následující hodnoty JND indexů:
Pro P = 0: JND_P(0) = 70,67 + 0*0,158 = 70,67
Pro P = 240: JND_P(240) = 70,67 + 240*0,158 = 108,59
Pro P = 480: JND_P(480) = 70,67 + 480*0,158 = 550,83
Pro P = 720: JND_P(720) = 70,67 + 720*0,158 = 184,43

Pro P = 4080: JND_P(4080) = 70,67 + 4080*0,158 = 715,31

Máme tedy hodnoty všech 18 JND indexů, ke kterým ze vztahu (2) opět z předešlého příspěvku dopočítáme jas L z GSDF křivky. Hodnoty jsou následující (tučně zvýrazněné hodnoty jasu jsou změřené hodnoty jasů):

Dále změříme jas pro každý z 18 obrazů fotometrickou sondou. Změřené hodnoty jasu jsou následující:

Porovnáme-li hodnoty změřeného jasu a dopočítaného jasu z GSDF křivky ze dvou předešlých tabulek, musí pro splnění požadavků platit, že hodnoty pro tentýž bod křivky se od sebe liší pouze v rámci tolerance, tj. pro primární diagnostické displeje se mohou od sebe lišit maximálně o ±10 %. Liší-li se více, pak je potřeba displej překalibrovat, protože nevyhovuje GSDF křivce. Ukázka GSDF křivky s dopočítanými hodnotami jasu včetně křivek tolerance ±10 %, ale také ±20 % jsou uvedeny na obr. 3. Jsou zde uvedeny také hodnoty změřeného jasu.

Obr. 3: Změřené luminance a GSDF křivka s tolerancemi

Pro testování se však nepoužívá přímo porovnání jednotlivých měřených bodů s GSDF křivkou (tak jak je na obr. 3), protože splnění či nesplnění tolerance je hůře viditelné, ale pro zjednodušení se stanoví relativní změna d = dL/L jasu v daném intervalu (vždy pro dva sousední body) z následujícího vztahu (1):

kde L_k je jas pro P-hodnotu k, L_k+q je jas pro P-hodnotu k+q. Ke křivce GSDF jsou dopočítány tolerance ±10 %, případně ±20 %, a z graficky vynesených bodů je pak ihned zřejmé, je-li hodnota jasu bodu v toleranci či nikoliv. Hodnoty dL/L stanovené z GSDF křivky včetně tolerance ±10 % a ±20 % a hodnoty dL/L ze změřeného jasu jsou uvedeny na obr. 4.

Obr. 4: Hodnoty dL/L z GSDF křivky včetně tolerancí a hodnoty ze změřené luminance

Ukázka displeje nevyhovujícího GSDF křivce ani v rámci tolerance je uvedena na obr. 5. Jednalo se o běžný kancelářský displej.

Obr. 5: Ukázka hodnot displeje nevyhovujícího GSDF křivce

Většina displejů, i běžných kancelářských, splní požadavky GSDF křivky v nějaké její části, většinou prostřední, avšak nesplní požadavky v krajních částech, tj. v oblasti velmi tmavých a velmi světlých stupňů šedi. To je také případ displeje na obr. 5. A to je právě to, v čem se odlišují diagnostické displeje od ostatních displejů.

Subjektivní kontrast vnímaný lidským okem v tmavé části GSDF křivky je velmi ovlivněn okolním jasem a okolním osvětlením. Na obr. 6 jsou uvedeny křivky změřené při různém okolním jasu. Okolní jas je ovlivňen mnoho faktory, samozřejmě okolním osvětlením, ale např. také výmalbou místnosti (popisovny) (bílá vs. černá), zda má lékař bílý plášť. V USA je na některých pracovištích dokonce zakázáno, aby lékař vstupoval do popisovny v bílém plášti a stěny popisoven jsou natřeny černou barvou. Report AAPM TG270 doporučuje, aby při kalibraci displeje bylo zohledněno okolní osvětlení, přičemž by mělo platit, že okolní jas, resp.okolní osvětlení, nebude vyšší než 1/4 minimálního kalibrovaného jasu L_min. S vyšším okolním jasem se samozřejmě předpokládá, že bude vyšší minimální kalibrovaný jas L_min. V případě, že nebude při kalibraci brán v potaz okolní jas a okolní osvětlení, se může stát, že vizuální sytém člověka nerozpozná kontrast na úrovni tmavých odstínů šedi a může tak dojít k přehlédnutí patologie.

Obr. 6: Vliv okolního jasu na subjektivní kontrast vnímaný okem [8]

Použitá literatura
[1] Pianykh OS. Digital Imaging and Communications in Medicine (DICOM). Springer, 2012.
[2] Fetterly KA, Blume HR, Flynn MJ, et al. Introduction to grayscale calibration and related aspects of medical imaging grade liquid crystal displays. Journal of Digital Imaging, 2008; 21(2): 193-207.
[3] http://dicom.nema.org/medical/dicom/current/output/pdf/part14.pdf
[4] https://siim.org/page/displays_chapter3
[5] http://www.otpedia.com/entryDetails.cfm?id=226
[6] http://fyzika.jreichl.com/main.article/view/535-fotometricke-veliciny
[7] https://en.wikipedia.org/wiki/Luminance
[8] https://www.aapm.org/pubs/reports/RPT_270.pdf
[9] Silosky MS, Marsh RM. Performance characteristics and quality assurance considerations for displays used in interventional radiology and cardiac catheterization facilities. J Appl Clin Med Phys 2018; 19(5): 708-717

Napsat komentář

Vaše e-mailová adresa nebude zveřejněna. Vyžadované informace jsou označeny *