Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
YurkinPhD.pdf
Скачиваний:
67
Добавлен:
28.03.2016
Размер:
4.03 Mб
Скачать

 

800

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

700

 

 

 

 

теория

 

 

 

 

 

 

V = 92 мкм3

 

)

 

 

 

 

 

эксперимент

 

 

 

 

 

D = 6.33 мкм

 

)I(θ

600

 

 

 

 

V = 110 мкм3

 

 

 

 

 

β = 90°

 

 

. w(θ

500

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D = 7.6 мкм

 

 

 

 

 

χ

2

= 31

 

 

интен

 

 

 

 

 

β = 90°

 

 

 

 

 

 

 

 

 

400

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Модиф.

300

 

 

 

 

χ2=28

 

 

 

 

 

 

 

 

 

 

 

200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

700

 

 

 

 

V = 92 мкм3

 

 

 

 

 

V = 110 мкм3

 

(θ )

 

 

 

 

 

D = 6.75 мкм

 

 

 

 

 

D = 6.75 мкм

 

600

 

 

 

 

 

 

 

 

 

β = 80°

 

 

 

 

 

 

β = 90°

 

 

 

 

 

 

 

 

)I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. w(θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

500

 

 

 

 

χ

2

= 33

 

 

 

 

 

 

χ2 = 36

 

 

интен

400

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Модиф.

300

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

700

 

 

 

 

V = 100 мкм3

 

 

 

 

 

V = 100 мкм3

 

)

 

 

 

 

 

D = 8.28 мкм

 

 

 

 

 

D = 6.84 мкм

 

)I(θ

600

 

 

 

 

β = 90°

 

 

 

 

 

 

β = 80°

 

 

. w(θ

500

 

 

 

 

χ2 = 30

 

 

 

 

 

 

χ2 = 30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

интен

400

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Модиф.

300

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

010

15

20

25

30

35

 

40

45

15

20

25

30

35

40

45

50

 

 

 

Угол рассеяния θ, градусы

 

 

 

 

Угол рассеяния θ, градусы

 

 

 

 

 

Рис. 38. Экспериментальные и теоретические модифицированные индикатрисы зрелых

эритроцитов. Также приведены χ 2 расстояния.

 

 

 

 

 

 

 

 

 

 

3.2.4. Характеризация эритроцитов

Практически невозможно использовать МДД для подгонки экспериментальных индикатрис, ввиду большого времени вычислений даже для современной программы и аппаратных средств (см. подраздел 2.4.2). В качестве подхода к решению обратной задачи светорассеяния мы вычислили индикатрисы нескольких двояковогнутых дисков с разными диаметрами и объёмами, представленными в таблице 13, создав тем самым малую базу данных. В ячейках таблицы указаны значения η для соответствующих диаметра и объёма, а прочерк означает, что для этих параметров вычисления не проводились. Для каждой пары диаметра и объёма вычислялись четыре индикатрисы

148

Таблица 13. Параметры эритроцитов для предварительных вычислений.a

D, мкм

 

 

V, мкм3

 

 

86

92

100

105

110

6.08

0.638

6.33

0.565

0.605

6.51

0.556

0.604

0.665

6.75

0.466

0.499

0.542

0.569

0.596

6.84

0.521

0.573

7.01

0.484

0.532

7.60

0.327

0.349

0.380

0.399

0.418

8.28

0.294

0.418

0.323

a Приведены отношения толщины к диаметру, прочерки означают, что вычисления не проводились.

Количество

250

200

150

100

50

0

60

70

80

90

 

Угол ориентации β, градусы

 

Рис. 39. Распределение зрелых эритроцитов по углу ориентации, полученное χ2 методом.

для β = 60°, 70°, 80° и 90° – всего было вычислено 92 индикатрисы. Эти теоретические индикатрисы использовались для характеризации эритроцитов χ2 методом, как описано в подразделе 3.2.2.

Несколько представительных примеров результатов χ2 метода представлены на рис. 38 – видно хорошее согласие между теоретическими и экспериментальными индикатрисами. Средняя интенсивность (по всему диапазону углов) измеренных индикатрис соответствует теоретическим индикатрисам для значений β равных 70° и

90° [рис. 36(б)], и превосходит среднюю интенсивность для меньших β. Используя все экспериментальные индикатрисы, для которых значение χ2 было меньше заданного уровня (примерно 10% всех экспериментальных индикатрис), построено распределение зрелых эритроцитов по углу ориентации, показанное на рис. 39. Можно заключить, что ориентация эритроцитов в капилляре СПЦ близка к перпендикулярной (β = 90°), что согласуется с предыдущими результатами [225]. При этом рис. 39 даёт оценку отклонения от этой предпочтительной ориентации.

149

 

800

 

 

 

 

 

 

 

 

 

700

 

 

 

 

 

 

(а)

 

 

 

 

 

 

 

 

 

 

)

600

 

 

 

 

 

 

 

 

)I(θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(θ

500

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

интен

400

 

 

 

 

 

эритроцит

 

 

 

 

 

 

 

 

Модиф.

300

 

 

 

 

 

 

 

 

 

 

 

диск

 

 

200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

350

 

 

 

 

 

 

(б)

 

 

 

 

 

 

 

 

 

 

)I(θ )

300

 

 

 

 

 

 

 

 

250

 

 

 

 

 

 

 

 

w(θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

интен.

200

 

 

 

 

 

 

 

 

150

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

Модиф

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

50

 

 

 

 

 

 

 

 

 

0

15

20

25

30

35

40

45

50

 

10

 

 

 

 

Угол рассеяния θ, градусы

 

 

 

Рис. 40. Модифицированные индикатрисы двояковогнутого диска и шарового диска с тем же диаметром и объёмом: (а) перпендикулярная и (б) параллельная ориентация.

3.2.5. Приближённые формы

Существует несколько приближений формы нативных эритроцитов. Кроме совсем простых, наиболее популярными являются шаровой диск и сплюснутый сфероид. Их преимуществом по сравнению с реалистичными формами является то, что их можно быстро моделировать с помощью МРГУ. Мы сравнили индикатрисы эритроцита с

D = 7.60 мкм и η = 0.380 (V = 100 мкм3) с приближёнными формами для того же диаметра и объёма, рассматривая две ориентации оси эритроцита: перпендикулярную

(β = 90°) и параллельную (β = 0°) падающему излучению. Результаты для шарового диска показаны на рис. 40, где видно, что он является удовлетворительной моделью только в диапазоне от 10° до 15°.

Результаты сравнения для сплюснутого сфероида (рис. 41) показывают хорошее согласие с реальной формой, но только для перпендикулярной ориентации. Это согласуется с выводами на основе метода граничных элементов [226] и МДИ [228]. Однако обоснованность замены эритроцита сфероидом требует более обширной

150

 

800

 

 

 

 

 

 

 

 

 

700

 

 

 

 

эритроцит

 

(а)

 

 

 

 

 

сфероид

 

 

 

 

 

 

 

 

 

 

)

600

 

 

 

 

 

 

 

 

)I(θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(θ

500

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

интен

400

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Модиф.

300

 

 

 

 

 

 

 

 

200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(б)

 

250

 

 

 

 

 

 

 

 

)

200

 

 

 

 

 

 

 

 

)I(θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. w(θ

150

 

 

 

 

 

 

 

 

интен

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Модиф.

100

 

 

 

 

 

 

 

 

50

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

15

20

25

30

35

40

45

50

 

10

 

 

 

 

Угол рассеяния θ, градусы

 

 

 

Рис. 41. Модифицированные индикатрисы двояковогнутого диска и сплюснутого сфероида с тем же диаметром и объёмом: (а) перпендикулярная и (б) параллельная ориентация.

проверки для различных размеров эритроцита в контексте конкретной задачи светорассеяния.

3.2.6. Выводы

Моделирование светорассеяния зрелыми эритроцитами показало, что индикатриса чувствительна к форме эритроцита, поэтому в общем случае требуется использовать МДД (или другой метод способный работать с реалистичной формой) для изучения влияния характеристик эритроцита на его индикатрису. Однако при ориентации оси эритроцита перпендикулярно падающему излучению он может быть заменён на сплюснутый сфероид, светорассеяние которым можно моделировать с помощью метода расширенных граничных условий (МРГУ). К счастью, гидродинамическая система сканирующего проточного цитометра (СПЦ) доставляет эритроциты в измерительную зону именно в этой ориентации. Это может сильно облегчить решение обратной задачи, например, с помощью параметризации или нейронной сети, так как МРГУ намного быстрее чем МДД. Тем не менее точность подобных алгоритмов необходимо

151

проверять, используя реалистичные индикатрисы, вычисленные, например, с помощью МДД.

Подход на основе базы данных не упрощает форму эритроцита и является перспективным для его точной характеризации. Но малая база данных, построенная в данном разделе, только демонстрирует идею и не может использоваться для построения распределений эритроцитов в измеряемой на СПЦ пробе по диаметру, объёму и другим параметрам. Для того чтобы иметь возможность строить подобные распределения, необходимо составить намного бóльшую базу данных и отрегулировать

χ2 метод, а модель эритроцита должна включать, по крайней мере, вариацию концентрации гемоглобина. Все эти задачи рассматриваются в разделе 3.3.

152

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]