Исследования свойств гексагональных кодирующих коллиматоров для однофотонной эмиссионной томографии

Построение гексагонального кодирующего коллиматора. Полуэмпирическая формула расчета глубинного расширения. Влияние на свойства кодирующего коллиматора циклических перестановок элементов в последовательности. Тексты программ и описание их работы.

Рубрика Коммуникации, связь, цифровые приборы и радиоэлектроника
Вид дипломная работа
Язык русский
Дата добавления 10.06.2012
Размер файла 1,6 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

Реферат

Цель работы: Численно исследовать аппаратные функции кодирующих коллиматоров, построенных на базе псевдослучайных последовательностей, расширенных псевдослучайных последовательностей, троичных последовательностей, расширенных троичных последовательностей. Оптимизировать скорость расчета аппаратных функций гексагональных кодирующих коллиматоров. Исследовать полученные аппаратные функции, выбрать критерии отбора гексагональных кодирующих коллиматоров для различных целей. Оптимизировать итерационный алгоритм, позволяющий проводить восстановление пространственного распределения источников излучения.

Положения, выносимые на защиту: 1) Получена база данных характеристик глубинных аппаратных функций и сфокусированных изображений тестовых пространственных распределений источников излучения, позволяет осуществлять быстрый отбор кодирующих коллиматоров по раличным критериям. 2) Реализовано 2 итерационных алгоритма, позволяющие восстанавливать пространственное распределение источников излучения. 3) Проведено сравнение гексагональных кодирующих коллиматоров, построенных на основе псевдослучайных последовательностей, и построенных на основе троичных последовательностей.

В работе были проведены исследования более 250 кодирующих коллиматоров, работающих по различным схемам измерений. Для них были рассчитаны глубинные аппаратные функции и получены сфокусированные изображения тестовых пространственных распределений источников излучения. Полученные данные были сведены в общую базу данных, позволяющую проводить выбор кодирующих коллиматоров по различным критериям.

Глава 1. Введение в интегрально-кодовые системы измерений

В 1968 году Диком [1] была предложена идея использования кодирующих коллиматоров. Первоначально камера-обскура заменяется на большое количество пинхолов, расположенных случайным образом. Рисунок 1 дает простое представление об этой концепции. Каждый источник излучающего объекта вносит свой вклад в изображение. Последующая обработка позволяет получить восстановленное изображение, которое похоже на первоначальный объект.

Рис. 1 - Концепция использования кодирующих коллиматоров

Первоначальной целью было получить систему, которая поддерживает высокое угловое разрешение одиночного пинхола, но получаемое изображение имеет отношение сигнала к шуму (ОСШ), которое соизмеримо с общей открытой площадью апертуры. Метод, как правило, применяется для рентгеновского изображения, так как большинство источников рентгеновского излучения настолько слабы, что размер одной камеры-обскуры должен быть очень велик для того, чтобы получить разумное ОСШ. Большое отверстие исключает возможность получить хорошее угловое разрешение.

Если в коллиматоре N пинхолов, то изображение состоит из N перекрытых изображений объекта. Использование кодирующего коллиматора (для точечного источника), может улучшить ОСШ примерно в по сравнению с камерой-обскура. Так как N может достигать 100000, то идея улучшить ОСШ выполнима [1]. Кроме того, была поставлена и решена задача оптимизации среднего пропускания кодирующих коллиматоров в зависимости от квантовой статистики полезного сигнала и некодируемого фона [2].

Второй главной целью было проводить томографию, как было предложено Барретом [3]. Точки объекта на разных расстояниях от апертуры будут оставлять тень от апертуры всевозможных размеров на изображении. Можно получить послойные сфокусированные изображения трёхмерного источника излучения. Это свойство кодирующих коллиматоров в частности очень полезно в медицине, но также используется и в промышленности.

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

Интегрально-кодовые системы измерений (ИКСИ) ионизирующих излучений нашли применение в рентгеновской и гамма-астрономии [4], спектрометрии нейтронов по времени пролета [5, 6], радиационной интроскопии [2, 7], радиационной безопасности [8, 9], рентгеновской диагностике и других областях.

Дальнейшее развитие ИКСИ связано с использованием кодирующих коллиматоров для томографической реконструкции трехмерных пространственных распределений радионуклидов без вращения массивной детектирующей системы вокруг объекта исследования. При этом большое значение имеют фокусирующие свойства кодирующих коллиматоров [10].

Среди кодирующих устройств особого внимания заслуживают двумерные многопинхольные кодирующие коллиматоры [11], применение которых в ИКСИ дает возможность получать не только планарное двумерное, но и восстанавливать трехмерное распределение радионуклидов в объектах [12].

Двумерные многопинхольные кодирующие коллиматоры строятся на основе двумерных кодовых таблиц (ДКТ). ДКТ чаще всего строятся на основе одномерных двоичных псевдослучайных последовательностей (ПСП) [13] из нулей и единиц. Также ДКТ могут строиться на основе троичных последовательностей (ТП) [14], состоящих из +1, 0 и -1.

Выполненные исследования показали, что успешное практическое решение таких задач как планирование радиационно-физического эксперимента [15] и планарная эмиссионная томография [12] во многом ограничено очень небольшим набором существующих псевдослучайных последовательностей требуемой размерности и, особенно, с нужным коэффициентом пропускания, необходимым для построения оптимальных кодирующих коллиматоров.

Это обстоятельство стимулировало поиск неизвестных ранее классов последовательностей. В 2007 году был предложен новый класс последовательностей, названных расширенными последовательностями [16, 17], который в несколько раз увеличивает число пригодных для использования в ИКСИ кодирующих коллиматоров. При этом расширенные последовательности можно получать как из ПСП, так и из ТП.

Кодирующие устройства на основе квадратных и прямоугольных таблиц хорошо сочетаются с квадратными и прямоугольными детекторами. Однако при применении круглых детекторов в сочетании с кодирующими устройствами прямоугольной конфигурации существенная часть (36% и более) полезной площади детекторов не используется. По этой причине и вследствие целесообразности дальнейшего расширения класса двумерных кодирующих устройств представляют интерес гексагональные конфигурации, расположение пинхолов в которых определяется как на основе одномерных ПСП [18,19], так и на основе ТП. Аналогично, расширенные последовательности в несколько раз увеличивают число пригодных для использования в ИКСИ гексагональных кодирующих коллиматоров (ГКК).

При существенном увеличении количества ГКК актуальна задача поиска среди них предпочтительных коллиматоров. Существует задача по сравнению между собой ГКК, построенных на основе двоичных последовательностей на основании сравнения аппаратных функций (АФ) кодирующих коллиматоров. Аналогичная задача существует и для ГКК, построенных на основе троичных последовательностей. Кроме того, целесообразно сравнить ГКК, построенные на основе двоичных последовательностей, с ГКК, построенными на основе троичных последовательностей, с целью определения лучших конфигураций.

1.1 Построение ГКК

Псевдослучайная гексагональная конфигурация (ПСГК) строится с помощью ПСП или с помощью расширенной ПСП (РПСП) путем сворачивания ее в гексагональную структуру. Троичная гексагональная конфигурация (ТГК) строится с помощью троичной последовательности (ТП) или с помощью расширенной троичной последовательности (РТП) [13] сворачиванием ее в гексагональную структуру. На длину последовательности накладывается при этом следующее условие: , где - целое число, называемое рангом конфигурации, а - длина последовательности (количество ее элементов). Из последнего условия следует, что делится без остатка на 6.

На практике использование кодирующих устройств на основе ТП и РТП сводится к изготовлению двух коллиматоров, в одном из которых положение открытых ячеек задается позицией +1, а в другом позицией -1 в ТП или РТП; проведению двух измерений и вычитанию результатов второго измерения из результатов первого [13].

Процедуру построения последовательности длиной в базовую и мозаичную ПСГК иллюстрирует рисунке 2. Количество элементов на диагонали базовой конфигурации равно . Большое значение имеет также возможность построения многопинхольного коллиматора на основе двумерной мозаики ее базовой части, что делает матрицу, описывающую процесс кодирования, циклической (матрицей-циркулянтом), как это имеет место и для прямоугольных коллиматоров [17]. Мозаика строится добавлением к базовому шестиугольнику со всех (шести) сторон половин таких же шестиугольников за исключением элементов, расположенных на их диагонали [18, 20]. Поэтому полное число элементов в мозаичной структуре равно .

В работах [21,22] сравнивались многопинхольные коллиматоры со случайным расположением пинхолов с коллиматорами с регулярным и неизбыточным (минимально избыточным) расположением пинхолов. Неизбыточное расположение первоначально было использовано при синтезе антенн [23]. Некоторые двумерные неизбыточные расположения были найдены в работе [24]. Затем неизбыточные расположения были использованы при построении многопинхольных кодирующих коллиматоров.

Рис. 2 - Преобразование ПСП в базовую (заштрихована) и мозаичную (обведена жирными внешними линиями) гексагональные конфигурации с рангом

Примеры мозаичных гексагональных кодирующих коллиматоров разных рангов показаны на рисунке 3.

a

б

Рис. 3 - Мозаичный кодирующий коллиматор построенный на основе ПСП с (а), мозаичный кодирующий коллиматор построенный на основе ТП с (черные ячейки - «0», белые ячейки - «1», серые ячейки - «-1») (б)

Удобно разделить класс ПСП на три подкласса: классические (), ассоциированные () и вырожденные (), где - количество единиц в ПСП. Алгоритм построения РПСП и РТП заключается в добавлении к каждому элементу псевдослучайной последовательности любого одинакового числа нулей , при этом для полученной длины должен выполняться критерий . В табл. 1 указано количество возможных гексагональных конфигураций на основе ПСП и РПСП; ТП И РТП с , свидетельствующее о преобладании последних.

Таблица 1 - Количество возможных гексагональных конфигураций

ПСП 

Классические

Ассоциированные

Вырожденные

Всего

ТП

Нерасширенные

19

19

39

77

8

Расширенные

64

64

77

205

45

Всего

83

83

116

282

53

1.2 Аппаратные функции кодирующих коллиматоров

В связи с увеличением количества возможных кодирующих коллиматоров (КК) является актуальным систематическое исследование их свойств и выбор лучших КК для решения задач вычислительной томографии. При этом для ИКСИ с ГКК также может быть использован метод фокусных плоскостей, исследованный ранее для многопинхольных прямоугольных кодирующих коллиматоров [2]

Удобным инструментом для такого исследования являются глубинные аппаратные функции ИКСИ, характеризующие влияние нефокусных источников на изображение в плоскости фокуса [10], которые определяются по формуле (1)

,(1)

где - координаты источника, - Фурье-образ функции пропускания кодирующего коллиматора, - расстояние между детектором и кодирующим коллиматором.

Аппаратные функции можно получить только численно. Пример относительно хорошей и очень плохой глубинной АФ, характеризующей вклад нефокусных источников излучения в результаты измерений в фокусной плоскости, представлен на рисунке 4. Значение АФ в плоскости фокуса принято за единицу. Коллиматоры, имеющие очень плохие АФ, целесообразно исключить из дальнейшего анализа.

Рис. 4 - Относительно хорошая (сплошная линия) и плохая (штриховая линия) аппаратные функции кодирующих коллиматоров в зависимости от смещения нефокусной плоскости от плоскости фокуса ()

Для характеристики томографических (фокусирующих) свойств ИКСИ могут быть использованы такие параметры, как глубинное разрешение , равное ширине АФ на половине высоты около плоскости фокуса, максимальная амплитуда ложных пиков и интегральные критерии и (рисунок 5).

Рис. 5 - Основные параметры аппаратной функции: глубинное разрешение , максимальная амплитуда ложного пика и интегральные критерии и

1.3 Круговая аппроксимация

Для ПСГК больших размерностей расчет АФ является ресурсоёмкой задачей, поэтому было предложено увеличение скорости расчета АФ за счет замены шестиугольной ячейки круглой. При этом полезная площадь круглого пинхола составляет 91% от площади шестиугольной ячейки. Такое приближение было названо круговой аппроксимацией (КА).

Рис. 6 - Верхние (1), средние(2) и нижние (3) АФ для ПСГК ранга 6 рассчитанные без КА (штриховая линия) и с КА (сплошная линия)

Оказалось, численно полученные аппаратные функции точно и с КА (рисунок 6) примерно одинаковы. Однако аппаратная функция ПСГК, рассчитанная с КА является более гладкой по сравнению с АФ, рассчитанной точно. Кроме того, скорость расчета аппаратной функции с КА для КК больших размерностей существенно выше.

Таким образом, схема измерений преобразуется к виду на рисунке 7.

Рис. 7 - точечный источник; 2 - пинхол; 3 - ячейка детектора; 4 - тень от точечного источника через (2)

Глава 2. Анализ аппаратных функций

2.1 Численное исследование аппаратных функций

Было исследованы 226 ИКСИ построенных на основе многопинхольных гексагональных кодирующих коллиматоров, из которых 53 построены на основе РТП и ТП, и 173 построены на основе РПСП и ПСП. В таблицах 2 и 3 представлены характеристики численно рассчитанных АФ.

Расчеты выполнены при следующих геометрических условиях: размер пинхола , размер ячейки ПЧД , фокусное расстояние . В табл. 2 , - коэффициент расширения ПСП, и - верхнее и среднее значение АФ в мм. соответственно, и - среднее и нижнее значение наибольшего ложного пика в относительных единицах. В табл. 3 - количество открытых пинхолов, - количество +1 в ТП или РТП, - количество --1 в ТП или РТП, - коэффициент расширения ТП, , и - верхнее, среднее и нижнее значение АФ в миллиметрах соответственно, , и - верхнее, среднее и нижнее значение наибольшего ложного пика в относительных единицах.

Таблица 2 - Характеристики численно рассчитанных АФ ГКК, построенных на основе ПСП и РПСП (вырожденные ПСП в таблице не представлены)

Параметры ПСП

 

Параметры АФ

 

 

л

 

 

 

 

 

 

 

7

7

3

1

0

1

-

102,13

114,11

2,14

1,15

0,69

0,01

0,02

7

7

4

2

0

1

-

99,68

112,12

1,54

1,18

0,87

0,01

0,02

19

19

9

4

0

2

125,46

65,89

117,23

1,73

0,97

0,56

0,02

0,04

19

19

10

5

0

2

124,31

56,39

117,61

1,76

1,01

0,59

0,03

0,06

37

37

9

2

0

3

127,11

35,38

113,13

1,71

0,77

0,39

0,07

0,29

37

37

28

21

0

3

126

35,67

115,07

1,69

0,96

0,59

0,07

0,09

91

7

3

1

12

5

-

19,72

40,64

3,56

0,58

0,28

0,11

3,34

91

7

4

2

12

5

59,14

24,59

43,29

2,77

0,67

0,37

0,18

4,10

91

13

4

1

6

5

-

21,2

47,05

2,45

0,57

0,27

0,08

1,89

91

13

9

6

6

5

117

29,2

38,47

2,31

1,2

0,47

0,13

1,31

91

91

10

1

0

5

119

23,41

36,38

1,99

0,65

0,27

0,17

1,51

91

91

81

72

0

5

118,94

23,56

36,37

1,73

0,89

0,57

0,17

0,19

127

127

63

31

0

6

44,41

22,54

31,03

1,18

0,66

0,45

0,28

0,56

127

127

64

32

0

6

44,51

22,71

31,32

1,09

0,66

0,46

0,29

0,57

169

13

4

1

12

7

56,8

15,64

30,93

2,67

0,56

0,33

0,15

6,34

169

13

6

3

12

7

49,37

20,41

27,42

2,47

0,93

0,44

0,18

5,07

217

7

3

1

30

8

47,73

12,2

23,33

3,7

0,55

0,22

0,19

13,74

217

7

4

2

30

8

38,4

14,75

26,15

2,77

0,66

0,23

0,26

14,11

217

31

6

1

6

8

-

13,43

24,73

2,22

0,48

0,29

0,23

8,32

217

31

15

7

6

8

37,47

15,76

22,78

2,35

0,54

0,24

0,27

3,91

217

31

16

8

6

8

34,25

17,01

24,87

2,04

0,61

0,31

0,30

4,07

217

31

25

20

6

8

32,04

17

23,93

2,22

1,31

0,71

0,20

1,74

271

271

135

67

0

9

30,29

15,89

20,83

0,78

0,57

0,45

0,46

0,93

271

271

136

68

0

9

30,25

15,88

20,83

0,78

0,57

0,44

0,46

0,91

331

331

165

82

0

10

25,78

14,28

19,24

0,78

0,55

0,43

0,69

1,38

331

331

166

83

0

10

25,82

14,28

19,24

0,78

0,55

0,41

0,49

0,97

469

7

3

1

66

12

33,34

8,19

15,52

3,72

0,52

0,22

0,24

37,52

469

7

4

2

66

12

26,41

10,07

17,4

2,88

0,61

0,2

0,31

36,35

469

67

33

16

6

12

22,12

11,35

16,38

1,33

0,55

0,29

0,45

6,40

469

67

34

17

6

12

22,92

11,36

16,05

1,4

0,53

0,27

0,45

6,21

547

547

273

136

0

13

19,78

11,22

14,89

0,68

0,57

0,41

0,55

1,10

547

547

274

137

0

13

19,81

11,24

14,89

0,68

0,57

0,42

0,55

1,10

631

631

315

157

0

14

20,44

10,36

13,81

0,79

0,56

0,44

0,57

1,14

631

631

316

158

0

14

20,44

10,36

13,8

0,79

0,56

0,44

0,57

1,14

721

7

3

1

102

15

27,12

6,55

12,55

3,84

0,55

0,21

0,26

62,49

721

7

4

2

102

15

21,36

7,88

14

2,88

0,61

0,18

0,33

59,48

721

103

51

25

6

15

18,14

9,35

12,97

1,43

0,55

0,31

0,52

7,35

721

103

52

26

6

15

17,43

9,26

13,19

1,33

0,57

0,31

0,53

7,35

817

19

9

4

42

16

19,78

7,32

12,41

2,62

0,55

0,18

0,39

35,40

817

19

10

5

42

16

18,77

7,51

12,25

2,32

0,54

0,16

0,41

33,50

817

43

43

43

18

16

18,68

8,58

11,82

1,65

0,5

0,21

0,53

10,07

817

43

43

43

18

16

17,38

9,15

12,78

1,73

0,56

0,21

0,53

10,07

919

919

459

229

0

17

14,75

8,73

11,72

0,65

0,55

0,44

0,63

1,26

919

919

460

230

0

17

14,76

8,73

11,72

0,65

0,55

0,44

0,63

1,26

1027

13

4

1

78

18

22,21

6

11,12

3,17

0,53

0,22

0,31

79,59

1027

13

9

6

78

18

15,44

8

11,05

2,99

0,86

0,34

0,34

38,80

1027

79

39

19

12

18

16,02

7,59

10,73

1,48

0,54

0,34

0,57

15,01

1027

79

40

20

12

18

14,7

7,69

11,12

1,38

0,58

0,33

0,58

14,89

1141

7

3

1

162

19

21,7

5,39

10,2

3,84

0,51

0,21

0,27

102,69

1141

7

4

2

162

19

17,03

6,35

11,11

2,88

0,61

0,18

0,35

99,84

1141

163

81

40

6

19

13,83

7,41

10,45

0,7

0,55

0,33

0,60

8,45

1141

163

82

41

6

19

14,12

7,48

10,4

0,7

0,54

0,35

0,61

8,49

1261

13

4

1

96

20

19,92

5,55

10,29

2,88

0,53

0,21

0,32

100,88

1261

13

9

6

96

20

14,21

7,13

10,09

2,56

0,85

0,37

0,34

47,64

1387

19

9

4

72

21

15,33

5,74

9,68

2,56

0,56

0,16

0,43

66,27

1387

19

10

5

72

21

14,5

5,83

9,52

2,31

0,54

0,15

0,45

62,42

1387

73

9

1

18

21

17,32

5,62

8,95

1,74

0,49

0,2

0,36

55,48

1387

73

64

56

18

21

12,29

6,94

9,65

1,67

1,19

0,63

0,35

7,59

1519

7

3

1

216

22

18,87

4,57

8,65

3,84

0,51

0,19

0,28

141,77

1519

7

4

2

216

22

14,88

5,48

9,8

2,88

0,62

0,16

0,36

136,71

1519

31

6

1

48

22

17,96

5,04

9,29

2,44

0,47

0,18

0,38

96,20

1519

31

15

7

48

22

14,55

5,84

8,75

1,75

0,51

0,18

0,45

45,57

1519

31

16

8

48

22

12,73

6,41

9,53

1,6

0,59

0,16

0,60

56,96

1519

31

25

20

48

22

12,39

6,39

9,17

2,31

0,86

0,52

0,34

20,66

1951

1951

975

487

0

25

10,52

6,04

7,93

0,63

0,55

0,47

0,70

1,40

1951

1951

976

488

0

25

10,52

6,04

7,94

0,63

0,55

0,46

0,70

1,40

2107

7

4

2

300

26

12,64

4,77

8,19

2,88

0,61

0,17

0,37

194,90

2107

43

21

10

48

26

11,71

5,49

7,36

1,44

0,49

0,2

0,64

64,21

2107

43

22

11

48

26

10,86

5,69

7,9

1,36

0,56

0,2

0,65

62,25

2611

7

3

1

372

29

14,45

3,4

7,06

3,84

0,51

0,19

0,29

252,40

2611

7

4

2

372

29

11,38

4,11

7,38

2,88

0,61

0,16

0,37

241,52

2791

2791

1395

697

0

30

8,73

5,2

6,77

0,55

0,55

0,48

0,73

1,46

2791

2791

1396

698

0

30

8,73

5,2

6,78

0,55

0,55

0,48

0,73

1,46

2977

13

4

1

228

31

12,75

3,49

6,99

2,88

0,53

0,2

0,35

260,49

2977

13

9

6

228

31

9,24

4,8

6,67

2,56

0,84

0,33

0,38

125,70

3367

7

3

1

480

33

12,72

3,09

6,51

2,56

0,5

0,19

0,30

336,70

3367

7

4

2

480

33

10,09

3,63

6,74

1,92

0,59

0,14

0,38

319,87

3367

13

4

1

258

33

11,97

3,31

6,72

2,88

0,53

0,2

0,36

303,03

3367

13

9

6

258

33

8,52

4,44

6,38

1,28

0,99

0,18

0,67

250,65

3367

37

9

2

90

33

9,24

3,45

6,35

1,28

0,42

0,14

0,58

216,98

3367

37

28

21

90

33

8,55

4,27

6,35

1,32

0,71

0,37

0,45

53,51

3367

91

10

1

36

33

9,81

3,51

6,26

1,3

0,35

0,15

0,49

164,98

3367

91

81

72

36

33

7,9

4,8

6,32

1,27

0,97

0,56

0,46

19,12

3571

3571

1785

892

0

34

7,9

4,53

6,15

0,61

0,54

0,44

0,76

1,52

3781

19

9

4

198

35

9,33

3,37

6,17

1,71

0,53

0,16

0,65

273,07

3781

19

10

5

198

35

8,86

3,45

6,05

1,54

0,52

0,15

0,49

185,27

3781

199

99

49

18

35

7,93

3,94

5,99

0,7

0,52

0,34

0,73

27,88

3781

199

100

50

18

35

7,58

6,06

4,14

0,7

0,53

0,35

0,73

27,60

3997

7

3

1

570

36

11,72

2,92

6

2,56

0,5

0,19

0,30

399,70

3997

7

4

2

570

36

9,27

3,33

6,33

1,92

0,59

0,14

0,38

379,72

3997

571

285

142

6

36

7,41

4,02

5,89

0,64

0,52

0,4

0,73

10,24

3997

571

286

143

6

36

7,41

4,06

5,89

0,64

0,52

0,41

0,73

10,20

4681

31

6

1

150

39

10,34

2,97

5,71

0,93

0,46

0,17

0,47

366,68

4681

31

15

7

150

39

8,31

3,25

5,36

0,91

0,48

0,18

0,51

159,15

4681

31

16

8

150

39

7,48

3,54

5,71

0,89

0,58

0,16

0,72

210,65

4681

151

75

37

30

39

7,37

3,63

5,43

0,87

0,5

0,33

0,75

46,81

4681

151

76

38

30

39

6,87

3,74

5,58

0,87

0,53

0,35

0,75

46,19

4921

7

3

1

702

40

10,57

2,77

5,23

0,62

0,5

0,19

0,30

492,10

4921

7

4

2

702

40

8,34

3,05

5,81

0,64

0,58

0,16

0,38

467,50

4921

19

9

4

258

40

8,14

3,05

5,47

0,63

0,52

0,16

0,66

360,87

4921

19

10

5

258

40

7,74

3,05

5,38

0,66

0,5

0,15

0,50

246,05

4921

37

9

2

132

40

7,52

2,98

5,35

0,64

0,41

0,14

0,62

339,00

4921

37

28

21

132

40

7,13

3,47

5,41

0,64

0,68

0,37

0,46

80,85

4921

133

12

1

36

40

8,16

3,07

5,39

1,19

0,32

0,12

0,61

250,15

4921

133

33

8

36

40

7,81

3,52

5,3

0,62

0,3

0,16

0,68

101,40

5677

7

3

1

810

43

9,84

2,65

4,62

0,65

0,51

0,19

0,30

567,70

5677

7

4

2

810

43

7,76

2,88

5,4

0,65

0,58

0,16

0,38

539,32

5941

13

4

1

456

44

9,09

2,74

4,86

2,55

0,52

0,2

0,37

549,54

5941

13

9

6

456

44

6,83

3,3

4,97

2,27

0,82

0,33

0,40

264,04

6487

13

4

1

498

46

8,58

2,68

4,51

2,08

0,52

0,2

0,37

600,05

6487

13

9

6

498

46

6,58

3,16

4,68

1,85

0,82

0,32

0,40

288,31

6769

7

3

1

966

47

9,03

2,57

4,04

2,18

0,5

0,19

0,31

699,46

Таблица 3 - Характеристики численно рассчитанных АФ ГКК, построенных на основе ТП и РТП

 

Исходная ТП

Параметры АФ

 

 

 

 

 

 

7

7

4

3

1

0

-

77,9

-

2,1

1,33

0,51

0,01

0,02

91

7

4

3

1

12

58,79

20,43

38,59

2,77

0,57

0,3

0,19

4,32

91

13

9

6

3

6

53,22

27,61

40

2,09

0,62

0,39

0,25

2,53

91

91

36

21

15

0

54,77

25,61

37,53

1,59

0,68

0,42

0,25

0,63

91

91

81

45

36

0

64,79

26,27

37,9

1,19

0,8

0,53

0,26

0,29

127

127

64

36

28

0

46

22,51

32,02

1,29

0,65

0,36

0,27

0,54

169

13

9

6

3

12

39,38

19,17

28,47

2,28

0,6

0,26

0,34

6,38

217

7

4

3

1

30

38,19

12,89

23,07

2,88

0,55

0,24

0,27

14,65

217

31

16

10

6

6

34,45

16,53

24,34

1,87

0,59

0,23

0,3

4,07

217

31

25

15

10

6

31,84

17,01

23,91

1,71

0,47

0,29

0,45

3,91

217

217

64

36

28

0

33,41

16,8

23,83

1,12

0,57

0,33

0,38

1,29

217

217

100

55

45

0

33,08

17,58

23,88

0,93

0,58

0,39

0,45

0,98

469

7

4

3

1

66

26,44

8,54

17,57

2,88

0,54

0,16

0,31

36,35

721

7

4

3

1

102

21,43

6,88

12,6

2,88

0,53

0,2

0,33

59,48

1027

13

9

6

3

78

16,57

7,54

11,37

2,56

0,57

0,18

0,59

67,33

1141

7

4

3

1

162

17,01

5,56

10,23

2,88

0,53

0,16

0,35

99,84

1261

13

9

6

3

96

14,95

6,75

10,38

2,56

0,58

0,17

0,61

85,47

1387

73

64

36

28

18

12,23

6,94

9,65

0,7

0,57

0,42

0,71

15,39

1519

7

4

3

1

216

14,85

4,82

8,73

2,88

0,53

0,15

0,36

136,71

1519

31

16

10

6

48

13,12

6,21

9,3

1,44

0,57

0,15

0,46

43,67

1519

31

25

15

10

48

12,34

6,39

9,16

1,43

0,46

0,27

0,7

42,53

1519

217

64

36

28

6

12,11

6,51

9,19

0,77

0,56

0,32

0,66

15,66

1519

217

100

55

45

6

12,11

6,67

9,1

0,52

0,28

0,2

0,71

10,78

2107

7

4

3

1

300

12,67

3,93

7,44

2,88

0,53

0,14

0,37

194,90

2611

7

4

3

1

372

11,41

3,56

6,98

2,88

0,52

0,16

0,37

241,52

2977

13

9

6

3

228

9,87

4,55

6,87

1,28

0,56

0,16

0,67

221,62

3367

7

4

3

1

480

10,08

3,16

6,42

1,92

0,52

0,15

0,38

319,87

3367

13

9

6

3

258

9,28

4,09

6,53

1,28

0,56

0,16

0,67

250,65

3367

91

36

21

15

36

8,32

4,3

6,3

0,5

0,34

0,18

0,75

70,15

3367

91

81

45

36

36

8,01

4,67

6,31

0,52

0,33

0,22

0,79

32,84

3997

7

4

3

1

570

9,27

2,99

5,96

1,92

0,51

0,15

0,38

379,72

4681

31

16

10

6

150

7,54

3,44

5,59

1,04

0,53

0,14

0,52

152,13

4681

31

25

15

10

150

7,1

3,54

5,5

1

0,43

0,25

0,79

147,92

4681

4681

4096

2080

2016

0

6,81

3,8

5,5

0,7

0,51

0,43

0,79

0,90

4921

7

4

3

1

702

8,35

2,8

5,28

1,67

0,51

0,15

0,38

467,50

4921

133

121

66

55

36

6,69

3,67

5,37

0,4

0,25

0,18

0,82

33,35

5677

7

4

3

1

810

7,76

2,69

4,72

1,66

0,5

0,14

0,38

539,32

5941

13

9

6

3

456

7

3,13

5,07

1,13

0,53

0,16

0,69

455,48

6487

13

9

6

3

498

6,75

3

4,8

0,93

0,53

0,16

0,69

497,34

6769

7

4

3

1

966

7,16

2,58

4,09

1,64

0,49

0,14

0,39

659,98

7651

7

4

3

1

1092

6,76

2,57

3,86

1,58

0,49

0,14

0,39

745,97

7651

1093

729

378

351

6

5,53

3,09

3,98

0,38

0,32

0,27

0,83

8,71

7651

7651

2916

1485

1431

0

 

 

 

 

 

 

 

0,00

7957

73

64

36

28

108

5,36

3,05

3,96

0,6

0,48

0,36

0,84

104,44

8587

31

16

10

6

276

5,78

2,78

3,83

0,93

0,48

0,14

0,76

407,88

8587

31

25

15

10

276

5,37

2,95

3,8

0,65

0,4

0,18

0,82

281,65

8911

7

4

3

1

1272

6,31

2,57

3,71

1,52

0,48

0,14

0,39

868,82

8911

133

121

66

55

66

5,04

3,03

3,75

0,37

0,23

0,17

0,85

62,82

9919

7

4

3

1

1416

6,01

2,56

3,63

1,48

0,48

0,14

0,39

967,10

9919

13

9

6

3

762

5,52

2,81

3,7

0,93

0,5

0,16

0,70

771,48

9919

91

36

21

15

108

5

2,85

3,58

0,46

0,3

0,17

0,81

223,18

9919

91

81

45

36

108

4,82

2,96

3,58

0,46

0,32

0,21

0,85

104,09

2.2 Полуэмпирическая формула расчета глубинного разрешения

Для прямоугольных конфигураций существует полуэмпирическая формула расчета глубинного разрешения как полуширины АФ в области фокуса [10]. Можно предложить аналогичную по структуре формулу (2) и для гексагональных конфигураций, эффективность которой иллюстрирует рис. 8.

, (2)

где - размер пинхола, - размер ячейки детектора, - расстояние от детектора до коллиматора, R - ранг конфигурации.

а) б)

в)

Рис. 8 - Глубинное разрешение усредненных АФ исследованных ГКК: построенных на основе классических и ассоциированных ПСП (а), построенных на основе всех ПСП (б), построенных на основе ТП и РТП (в) : рассчитанное по результатам моделирования (квадратики) и определенное по формуле (2) (сплошная линия) при и

2.3 Зависимость интегрального критерия и нормированного интегрального критерия от среднего пропускания ГКК

Интересно сравнить параметры ГКК построенных на основе двоичных и троичных последовательностей. Основными характеристиками АФ ГКК являются среднее пропускание и интегральный критерий . На рисунках 9 и 10 представлены зависимость интегрального критерия и нормального интегрального критерия соответственно от среднего пропускания.

Рис. 9 - Зависимость интегрального критерия от среднего пропускания ГКК для двоичных последовательностей (звездочки) и для троичных последовательностей (кружки)

Рис. 10 - Зависимость нормированного интегрального критерия от среднего пропускания ГКК для двоичных последовательностей (звездочки) и для троичных последовательностей (кружки)

Из рисунков 9 и 10 видно, что для АФ ГКК построенных на основе двоичных последовательностей максимальное значение интегрального критерия и нормированного интегрального критерия меньше, чем максимальное значение для ГКК построенных на основе троичных. В то же время, большие значения для ТП располагаются в области низкого среднего пропускания, а большие значения для двоичных последовательностей располагаются в районе среднего пропускания равного 0.5.

Кроме того, данный факт важен при непосредственном производстве кодирующего коллиматора - необходимо проделать меньшее число отверстий, что существенно облегчает технологический процесс.

2.4 Влияние на свойства ГКК циклических перестановок элементов в последовательности

Для формирования базовой части КК можно использовать любую циклическую перестановку данной ПСП или ТП. Поскольку при построении мозаичного КК из базового КК ранга не используются пинхолов, среди которых могут быть и открытые, то среднее пропускание мозаичного КК может заметно меняться в зависимости от явного вида исходной ПСП. Среди мозаичных гексагональных конфигураций, построенных на основе ПСП, существуют инвариантные по отношению к повороту на 120° вокруг центрального элемента (см. рис. 3а). Из одной ПСП или ТП можно получить как симметричные конфигурации, так и несимметричные в зависимости от использованной циклической перестановки. Поэтому интересно исследовать влияние среднего пропускания мозаичного КК на параметры АФ, а также сравнить симметричные и несимметричные конфигурации. В таблице 4 приведены основные параметры АФ таких ИКСИ и их вариации. Как видно из таблицы 4 симметричные конфигурации обладают лучшим интегральным критерием , особенно при малых рангах, и практически такими же глубинным разрешением и максимальной амплитудой ложных пиков по сравнению с несимметричными конфигурациями.

Таблица 4 - Основные параметры АФ и их вариации ГКК различных рангов

Ранг

2

3

6

9

Минимум

26

28

227

504

Максимум

31

34

242

519

Количество различных k

6

7

16

16

Минимум

119

122

41,4

28,1

Максимум

130

129

49,8

30,6

Среднее

122

124

44,6

29,4

Для симметричной ГКК

130

129

43,7

30,6

Минимум

54,77

34,95

21,83

15,33

Максимум

65,89

43,75

23,92

16,4

Среднее

58,71

38,11

22,95

15,64

Для симметричной ГКК

59,81

40,08

22,39

16,4

Минимум

1,03

1,01

0,97

0,73

Максимум

1,1

1,02

1,06

0,83

Среднее

1,08

1,01

0,99

0,78

Для симметричной ГКК

1,03

1,01

1,02

0,76

Минимум

1

0,8

0,6

0,6

Максимум

1,3

0,9

0,7

0,6

Среднее

1,1

0,8

0,6

0,6

Для симметричной ГКК

1,3

0,8

0,7

0,6

Минимум

0

0

0,2

0,5

Максимум

0,08

0,09

0,29

0,47

Среднее

0,02

0,06

0,27

0,46

Для симметричной ГКК

0,1

0,1

0,3

0,5

Минимум

0

0,17

0,49

0,91

Максимум

0,16

0,36

0,58

0,95

Среднее

0,03

0,26

0,56

0,92

Для симметричной ГКК

0,16

0,33

0,58

0,95

Инвариантные по отношению к повороту на 120° вокруг центрального элемента (симметричные) ПСГК обладают полезным свойством, заключающемся в том, что при их повороте на 60° меняются местами открытые и закрытые ячейки КК за исключением центральной [18]. Это может быть практически легко использовано для реализации более информативных при решении некоторых задач разностных схем измерений. Судя по табл. 4, симметричная конфигурация практически не имеет других предпочтительных свойств среди всех возможных конфигураций. При этом с увеличением ранга ПСГК диапазон изменения характеристик АФ всех конфигураций уменьшается.

2.5 Исследование локализации ложных пиков АФ

Отрицательной характеристикой аппаратных функции являются наличие ложных пиковв плоскости фокуса от нефокусных источников. Чем больше их число и амплитуда, тем большее влияние они оказывают на восстановленное распределение (изображение) радионуклидов в плоскости фокуса. В целом поведение аппаратной функции является трудно предсказуемым, однако ранее было отмечено, что у большинства прямоугольных КК ложные пики расположены примерно при одних и тех же смещениях [25]. Выполненные исследования показали, что таким же свойством обладают ГКК.

Рис. 11 - Аппаратные функции ГКК разных рангов: (штрих-пунктирная линия), (сплошная линия), (штриховая линия) и предполагаемые положения ложных пиков согласно таблице 4 (пунктирные вертикальные линии)

кодирующий коллиматор расширение программа

Как и для прямоугольных конфигураций, было предположено, что локальные экстремумы в АФ существуют при таких смещениях от фокусной плоскости, при которых тень от целого числа ячеек коллиматора полностью накрывает целое число ячеек детектора . Поскольку для одного и того же смещения может существовать несколько наборов (, ), можно предположить, что при большем количестве таких наборов существование экстремума более вероятно. Часть таких наборов и их количества для предполагаемых положений ложных пиков для ГКК ранга приведена в табл. 5.

Экстремумы в приведенных точках не обязательно существуют, однако все существующие особенности АФ (экстремумы и изломы) появляются только в этих точках.

Таблица 5 - Предполагаемые положения ложных пиков АФ и соответствующие им наборы ячеек

-66,67

1

2

12

-96,55

1

15

2

-63,16

4

13

2

-80,00

1

3

8

-96,30

1

14

2

-60,87

5

16

2

-85,71

1

4

6

-96,00

1

13

2

-58,82

4

12

2

-40,00

2

4

6

-95,65

1

12

2

-53,33

4

11

2

-88,89

1

5

5

-95,24

1

11

2

-50,63

5

14

2

-57,14

5

15

5

-94,74

1

10

2

-47,62

6

16

2

-92,31

1

7

4

-89,66

2

16

2

-47,06

5

13

2

-90,91

1

6

4

-88,00

2

14

2

-46,15

4

10

2

-72,73

2

7

4

-86,96

2

13

2

-42,11

6

15

2

-44,44

3

7

4

-84,21

2

11

2

-35,29

6

14

2

-28,57

3

6

4

-82,35

2

10

2

-31,58

7

16

2

-94,12

1

9

3

-81,48

3

16

2

-30,77

5

11

2

-93,33

1

8

3

-78,26

3

14

2

-26,67

6

13

2

-76,92

2

8

3

-76,19

3

13

2

-23,53

7

15

2

-61,54

3

9

3

-73,68

3

12

2

-18,18

5

10

2

-54,55

3

8

3

-72,00

4

16

2

-15,38

6

12

2

-36,36

4

9

3

-70,59

3

11

2

-13,33

7

14

2

-22,22

4

8

3

-69,57

4

15

2

-11,76

8

16

2

-96,77

1

16

2

 

 

 

 

 

 

 

 

Было получено следующее выражение, описывающее предполагаемое положение особых точек АФ:

.(3)

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

2.6 Семейства гексагональных кодирующих коллиматоров

При анализе положения ложных пиков АФ было выявлено, что ГКК, построенные на основе расширенных ПСП c одинаковой длиной обладают похожими аппаратными функциями.

Легко видеть, что ложные пики не только находятся в одних и тех же местах, но и близки по величине. Можно предположить, что АФ кодирующих коллиматоров ещё бoльших рангов, построенных на основе этих расширенных ПСП, будут обладать близкими к этим характеристиками. Таким образом, можно узнать примерный вид и основные характеристики АФ ГКК больших рангов, без их непосредственного расчета, что существенно упрощает исследования, потому что расчет АФ задача очень трудоемкая.

а)

б)

Рис. 12 - Семейства аппаратных функций, построенных на основе ПСП с (а), (б)

Результаты исследований описанные в первых двух главах были представлены в [19, 30-31].

Глава 3. Реконструкция распределений источников излучения

Основным способом реконструкции трёхмерных распределений источников излучения с применением ИКСИ является метод фокусных плоскостей (МФП), сводящий эту задачу к решению системы линейных алгебраических уравнений (СЛАУ) большого порядка [2, 7, 12, 26]. В МФП объект (трехмерное распределение источников ионизирующего излучения) условно разбивают на плоскостей, параллельных плоскости детектора (рисунок 13). Перемещая детектор с КК перпендикулярно плоскостям объекта, выполняют измерений, в каждом из которых в фокусе находится одна из этих плоскостей. Такая схема измерений позволяет не только оценить распределения источников в фокусных плоскостях (сфокусированные изображения), но и получить необходимые данные для точной реконструкции трёхмерного распределения источников излучения за счет исключения влияния внефокусных плоскостей. Для такого исключения требуется решить систему линейных алгебраических уравнений (СЛАУ), что, вследствие ее большого порядка, требует использования итерационных алгоритмов.

Рис. 13 - Метод фокусных плоскостей: 1 - позиционно-чувствительный детектор (ПЧД), 2 - кодирующий коллиматор, 3 - пространственный источник излучения, 4 - фокусная плоскость

В рентгеновской трансмиссионной томографии для тестирования используют фантом Шеппа-Логана [27], моделирующий сложное пространственное распределение коэффициента поглощения фотонов в среде (в голове пациента). В эмиссионной томографии также можно использовать аналогичный фантом, моделирующий сложное пространственное распределение источников излучения. Визуально фантом представляет собой суперпозицию эллипсоидов, имеющих различную ориентацию (рисунок 14). Внутри каждого отдельного эллипсоида интенсивность излучения имеет некоторое постоянное значение активности. В местах перекрытия эллипсоидов активность источника является суммой активностей всех перекрывающихся эллипсоидов.

а) б)

Рис. 14 - Трехмерный фантом Шеппа-Логана (а), и одно из его сечений (б)

Параметры такого модифицированного фантома Шеппа-Логана, использованного в данной работе, представлены в таблице 6.

Таблица 6 - Параметры модифицированного фантома Шеппа-Логана

Номер эллипсоида

Координаты центра эллипсоида

Полуоси

Угол

(град.)

Активность

1

(0.00; 0.00; 0.00)

(0.69; 0.92; 0.90)

0

1.0

2

(0.00; 0.00; 0.00)

(0.66; 0.87; 0.88)

0

-0.8

3

(-0.22; 0.00; -0.25)

(0.16; 0.41; 0.22)

18

-0.2

4

(0.22; 0.00; -0.25)

(0.11; 0.31; 0.21)

-18

-0.2

5

(0.00; 0.35; -0.25)

(0.21; 0.25; 0.24)

90

0.1

6

(0.00; 0.10; 0.00)

(0.046; 0.046; 0.16)

0

0.1

7

(0.00; -0.10; 0.00)

(0.046; 0.046; 0.046)

0

0.1

8

(-0.08; -0.605; 0.25)

(0.046; 0.023; 0.10)

0

0.1

9

(0.00; -0.606; 0.25)

(0.023; 0.023; 0.10)

0

0.1

10

(0.06; -0.605; 0.25)

(0.023; 0.046; 0.10)

0

0.1

Примечание: координаты расположены в плоскости параллельной плоскости детектора, ось - перпендикулярна плоскости детектора, - угол поворота в плоскости . Отрицательные активности введены для того, чтобы можно было промоделировать снижение активности относительно среднего уровня.

Для начала была решена задача восстановления плоского источника, рисунок 15.

Рис. 15 - Восстановление плоского источника: истинный плоский источник (а), восстановленный плоский источник (б)

Восстановление происходит идеально.

Следующим этапом было восстановление объёмного источника. На рисунке в качестве источника взяты буквы алфавита, расположенные в двух плоскостях, расстояние между которыми 16мм, рисунок 16.

Рис. 16 - Восстановление объёмного источника, истинный объёмный источник (а), восстановленное изображение (б)

Как видно из рисунка 16, восстановление происходит неточно, с артефактами. Это связано с влиянием внефокусных плоскостей. Возникает задача по улучшению восстановленного изображения.

При пространственной модуляции излучения в МФП для решения СЛАУ используются методы простой итерации (МПИ) [28] и скорейшего спуска (МСС) [29], а также модифицированный метод наименьшего направленного расхождения (МНР) [29].

Результаты модельных исследований МСС с использованием фантома Шеппа-Логана приведены на рис. 17.

Рис. 17 - Результаты модельных исследований с использованием фантома Шеппа-Логана: истинный объёмный источник (а), восстановленное изображение после 100 итераций по методу МСС (б)

Помимо визуальной оценки и сравнения полученных изображений, начальное и последующие приближения можно характеризовать средним квадратичным, абсолютным и максимальным отклонениями восстановленного изображения от действительного, соответственно:

(12)

Рис. 18 - Среднеквадратичное отклонение по 100 итерациям МСС

Результаты модельных исследований МНР с использованием фантома Шеппа-Логана приведены на рис. 19.

Рис. 19 - Результаты модельных исследований с использованием фантома Шеппа-Логана: истинный объёмный источник (а), восстановленное изображение после 100 итераций по методу МНР (б)

Зависимость среднеквадратичного отклонения от количества итераций приведена на рисунке 20.

Рис. 20 - Среднеквадратичное отклонение по 100 итерациям МНР

Можно сделать вывод о том, что скорость схождения МНР намного больше, чем скорость схождения МСС.

Глава 4. Тексты программ и описание их работы

Головная программа

%Расчет аппаратных функций ГКК, построенных на основе ПСП

global Params X_d XY_Det Y_d d1;

name_pst='LIST_PST.txt'; %Аппаратные функции, которые необходимо посчитать прописаны в файле

fid = fopen(name_pst,'r');

id = fgetl(fid);

while id~=-1

v = str2num(id(1:4)); %Получаем длину

k = str2num(id(5:8)); %Количество открытых пинхолов

l = str2num(id(9:12)); %Лямбда

n = str2num(id(13:16)); %Коэффициент расширения

maxR=(-3+sqrt(12*v*(n+1)-3))/6; %Из длины получаем Ранг ГКК

naaame='csv\';

Params.A=8; %Ячейка детектора

Params.a=4; %Ячейка КК

Params.L=100;%Расстояние от КК до ПЧД

Params.NumT=2; %Дискретность

if ~(k==v-1)

PSHEX=Make_list_KK_HEX_ex(minR,maxR,Params.a,naaame, n)

%считываем явный вид ПСП из ее параметров v,k,l,n

switch k %выбираем тип (вырожденная, ассоциированная или классическая)

case 1

pp=4;

case v-1

pp=3;

otherwise

if k>v/2

pp=2;

else

pp=1;

end

end

pshex = PSHEX(pp)

Params.R = pshex.R;

Params.MKK = pshex.open; %Записываем центры открытых пинхолов

XY_Det = Make_Centr_Hex(Params.A, pshex.R);

d1 = 0.5*numel(pshex.open);

X_d = repmat(XY_Det(:,1), 1, d1);

Y_d = repmat(XY_Det(:,2), 1, d1);

af=getaf(pshex); %Производится расчет АФ

af_d = min(af); %Находятся верхняя, нижняя, и средняя АФ

af_u = max(af);

af_m = mean(af);

%Отрисовываем аппаратные функции

figure()

cla

hold on

title('Up, Middle, Down')

f=Params.A*Params.L/(Params.A-Params.a);

Z=(-Params.L+Params.NumT:Params.NumT:Params.L);

af=[Z' af_u' af_d' af_m'];

plot(Z,af_d,'b','Linewidth',2 );

plot(Z,af_u,'r','Linewidth',2 );

plot(Z,af_m,'g','Linewidth',2 );

AF=flipud(af); %Приводим АФ к нужному для расчета виду

Анализ верхней АФ

nU = Archive_Array(AF(:,1),AF(:,2));

P = Return_Local_Max_New(nU(:,1),nU(:,2)); %поиск локальных максимумов

if isempty(P)

Glav_Peak = [NaN, NaN];

P = [];

AF_U = [nU(:,1),nU(:,2)];

GU = max(nU(:,2));

else

K = Analys_TU(P);

G = K{1};

Glav_Peak = G;

GU = G(1,2);

AF_U = Cut_AF(nU(:,1),nU(:,2),GU);

P = Return_Local_Max_New(AF_U(:,1),AF_U(:,2));

K = Analys_TU(P);

G = K{1};

Glav_Peak = G;

GU = G(1,2);

P = K{2};

end

if isempty(P)

T(1)=NaN;

TU_ar = [];

Envelope_U = [AF_U(:,1),AF_U(:,2)];

Dol_params = [NaN, NaN, NaN, NaN, NaN, NaN];

else

T(1) = max(P(:,2));

TU_ar = Form_TU_L(P);

Envelope_U = Return_Envelope(AF_U(:,1),AF_U(:,2),TU_ar);

Dol_params = Return_Dols_Params(AF_U(:,1),AF_U(:,2),Glav_Peak);

end

%2) анализ нижней АФ

nU = Archive_Array(AF(:,1),AF(:,3));

P = Return_Local_Max_New(nU(:,1),nU(:,2));

K = Analys_TU(P);

G = K{1};

GD = G(1,2);

P = K{2};

if isempty(P)

T(2)=NaN;

else

T(2) = max(P(:,2));

end

%3) анализ средней АФ

nU = Archive_Array(AF(:,1),AF(:,4));

P = Return_Local_Max_New(nU(:,1),nU(:,2));

K = Analys_TU(P);

G = K{1};

GM = G(1,2);

P = K{2};

if isempty(P)

T(3) = NaN;

else

T(3) = max(P(:,2));

end

Polush = Return_FWHM(AF(:,1),AF(:,2), AF(:,3), AF(:,4), GU, GD, GM);

%Формирование структуры характеристик

S = struct( 'V',str.vkln(1),...

'R',minR,...

'Pinholes', str.vkln(2),...

'AF',AF,...

'Du_1',Polush(1),...

'Dm_1',Polush(3),...

'Dd_1',Polush(2),...

'Tu_1',T(1),...

'Tm_1',T(3),...

'Td_1',T(2),...

'Env_U',Envelope_U,...

'GL_Peak',Glav_Peak,...

'IU_m', Dol_params(1), ...

'TU_array', TU_ar);

plot(S.Env_U(:,1), S.Env_U(:,2), 'k')

name = ['analizNEW\',sv,'_', svn, sk, sl, sn,'_', sA, sa, sL,'_',num2str(Params.ind), '.mat'];

save(name, 'S');

%Сохранение в файл всех характеристик

id = fgetl(fid);

end;

fclose(fid);

getaf.m

%Непосредственный расчет АФ

function massiv_tochkek_AF=getaf(PSHEX)

global XY_Det Params MMB MMBF;

% XY_Det = Make_Centr_Hex(Params.A, PSHEX.R);

FP_A = Params.A*Params.a/(Params.A-Params.a);

Формирование центров ячеек в фокусной плоскости

XY_FP = Make_Centr_Hex(FP_A, PSHEX.R);

%Процедура, необходимая для свертки, при использовании РПСП

MBF=eye(PSHEX.vkln(4)+1);

for i=1:PSHEX.vkln(1)

for j=1:PSHEX.vkln(1)

MMBF((i-1)*(PSHEX.vkln(4)+1)+1:(i-1)*(PSHEX.vkln(4)+1)+PSHEX.vkln(4)+1,(j-1)*(PSHEX.vkln(4)+1)+1:(j-1)*(PSHEX.vkln(4)+1)+PSHEX.vkln(4)+1)=MBF;

end

end

%Расчет клеше

V = length(XY_FP);

ist.Z = Params.A*Params.L/(Params.A-Params.a); % Определение фокусной плоскости

delenie=Params.NumT; %дискретность расчета плоскостей АФ

h1=waitbar(0,'Klesh');

for i=1:V

ist.X = XY_FP(i,1);

ist.Y = XY_FP(i,2);

% s = Shadow_Hex(ist, Params, XY_Det,Params.ind);

s = Shadow_Hex_New_Cir(ist); %поиск площади пересечения тени от ячейки КК с ячейками ПЧД

K(i,:) = s;

waitbar(i/V, h1, ['Kleshe ', num2str(i),' of ',num2str(V)]);

end

close(h1)

name=['massi&kleshe\',num2str(Params.R),'kleshe'];

save(name,'K') %сохранение клеше

Z = ist.Z-Params.L+delenie:delenie:2*ist.Z-Params.L;

LZ = length(Z);

massiv_tochkek_AF = zeros(V, LZ);

h=waitbar(1);

for i=1:LZ %расчет для каждой плоскости

IST.Z=Z(i);

for perem_xy=1:V; %для каждой открытой ячейки КК

IST.X = XY_FP(perem_xy,1);

IST.Y = XY_FP(perem_xy,2);

% S = Shadow_Hex(IST, Params, XY_Det,Params.ind);

S = Shadow_Hex_New_Cir(IST);

Y = K*S'; %свертка с клеше

X2 = RESTORE_HEX(Y, PSHEX); %восстановление

massiv_tochkek_AF(perem_xy,i)=max(abs(X2)); %получение АФ

end;

waitbar(i/LZ, h, ['Counted ', num2str(i),' plane of ',num2str(LZ)]);

end;

close(h)

Расчет площади пересечения тени КК с ячейками ПЧД

function S = Shadow_Hex_New_Cir(IST)

b=sqrt(3)/2;

global Params;

va = Params.a;

vA = Params.A;

vL = Params.L;

S_sot2 = pi*vA^2*3/4;

f2 = (vA*vL/(vA-va))^2;

x = IST.X;

y = IST.Y;

z = IST.Z;

az = z/(z-vL);

AA= va*az*b;

aa=b*vA;

global XY_Det X_d Y_d ;

V = length(XY_Det);

% d1=0.5*numel(Params.MKK);

% X_d = repmat(XY_Det(:,1), 1, d1);

% Y_d = repmat(XY_Det(:,2), 1, d1);

ff = ones(V,1);

X_SH = (x+(Params.MKK(:,1)-x)*az)';

Y_SH = (y+(Params.MKK(:,2)-y)*az)';

X_sh = ff*X_SH;

Y_sh = ff*Y_SH;

DDD = (X_sh-X_d).^2+(Y_sh-Y_d).^2; %критерий отбора

[I, J] = find(DDD<(AA+aa)^2); %номера ячеек ПЧД, которые могут попасть в тень

op = length(I);

A = max(AA,aa);

a = min(AA,aa);

a2 = a*a;

A2 = A*A;

s = zeros(1, V);

for index = 1:op %по всем номерам пригодных ячеек ПЧД

j = J(index);

k = I(index);

d = sqrt(DDD(k,j));

%расчет площади пересечения двух кругов

if A>=(a+d)

s(k)=s(k)+a2*pi;

else

x1=(d*d-a2+A2)/2/d;

ug1=2*acos(x1/A);

ug2=2*acos((d-x1)/a);

y1 = x1*tan(ug1/2);

S1=A2*ug1/2-x1*y1;

S2=a2*ug2/2-(d-x1)*y1;

s(k)=s(k)+S1+S2;

end;

end

S = s*f2/z/z/S_sot2; %нормировка на единицу (1=ячейка полностью лежит в тени)

Список сокращений

АФ - аппаратная функция

ДКТ - двумерная кодовая таблица

ИКСИ - интегрально-кодовые системы измерений

КК - кодирующий коллиматор

МНР - метод направленного расхождения

МСС - метод скорейшего спуска

МФП - метод фокусных плоскостей

ПСП - псевдослучайная последовательность

ПЧД - позиционно-чувствительный детектор

РПСП - расширенная псевдослучайная последовательность

РТП - расширенная троичная последовательность

СКО - среднеквадратичное отклонение

СЛАУ - система линейных алгебраических уравнений

ТП - троичная последовательность

Заключение

Проведённые исследования большого числа двоичных и троичных кодирующих коллиматоров, показали необходимость в систематизации и всестороннем анализе всех возможных кодирующих коллиматоров. В этой работе описана база данных всех характеристик глубинных аппаратных функций, параметров сфокусированных и улучшенных сфокусированных изображений. Обнаружено, что лучшими аппаратными функциями обладают гексагональные кодирующие коллиматоры, построенные на основе расширенных троичных последовательностей с маленьким средним пропусканием и гексагональные кодирующие коллиматоры, построенные на основе классических и ассоциированных псевдослучайных последовательностей со средним пропусканием около 0.5. Полученная формула предполагаемого положения ложных пиков сходится с численным исследованием аппаратных функций. Обнаружено, что «семейства» гексагональных кодирующих коллиматоров, построенных из одной и той же последовательности обладают схожими свойствами и одинаковым положением ложных пиков.

Реализованные итерационные алгоритмы для гексагональных кодирующих коллиматоров показали медленную сходимость, что показывает необходимость в новых итерационных алгоритмах, позволяющих проводить более быстро и точно восстановление пространственного излучения объёмного источника.

Список литературы

1. Dicke R.H. Scatter-hole cameras for X-rays and gamma rays// The Astrophysical J. 1968. V. 153. P. L101-L106.

2. Федоров Г.А., С.А. Терещенко. Вычислительная эмиссионная томография. М.: Энергоатомиздат, 1990. - 184 с.

3. Barrett H.H. Fresnel zone plate in nuclear medinice // J. of Nuclear Medicine, 1972. - V. 13, No. 6. - P. 382-385.

4. Caroly J, Stephen J.B., Di Cocco G. et al. Coded aperture imaging in X-and Gamma-ray astronomy”. Space Science Reviews, 1987. v. 45. № 3/4. P. 349-403.

5. Могильнер А.И., Сальников О.А., Тимохин Л.А. Корреляционный метод измерения энергетических спектров ядерных частиц по времени пролета // ПТЭ, 1966. №2. с. 22.

6. Wilhelmi G., Gompf F. Binary sequences and error analysis for psedo-statistical neutron modulators with different duty cycles // Nucl. Instr. Meth., 1970. v. 81. № 1. p. 36.

7. Федоров Г.А. Радиационная интроскопия. Кодирование информации и оптимизация эксперимента. М.: Энергоатомиздат, 1982.

8. P.T. Durrant, M. Dallimore, I.D. Jupp, D. Ramsden. The application of pinhole and coded aperture imaging in the nuclear environment. // Nucl. Instr. and Meth. A, 1999, vol. 422. p. 667-671.

9. O. Gal, M. Gmar, O. Ivanov, F. Laine, F. Lamadie, C. Le Goaller, C. Mahe, E. Manach. Development of a portable gamma camera with coded aperture // Nucl. Instrum. and Meth. in Phys. Res. Section А, 2006. V. 563. Issue 1. P. 233-237.

10. Федоров Г.А., Терещенко С.А. Интегрально-кодовые системы регистрации ионизирующих излучений. 3. Аппаратные функции. //Измерительная техника, 1997, №2, с.44-50; Fedorov G.A., Tereshchenko S.A. // Measurement Techniques, 1997. V. 40. No. 2. P. 164.

11. Федоров Г.А., Терещенко С.А. Интегрально-кодовые системы регистрации ионизирующих излучений. 1. Коды и кодирующие устройства // Измерительная техника, 1995. № 11. с. 49-54; Fedorov G. A., Tereshchenko S.A. // Measurement Techniques, 1995. V. 38. No. 11. P. 1287.11. 126.

12. Федоров Г.А., Терещенко С.А. Интегрально-кодовые системы регистрации ионизирующих излучений. Использование итерационных алгоритмов реконструкции изображений в методе фокусных плоскостей // Измерительная техника, 2001. № 4. С. 57-60; Fedorov G.A., Tereshchenko S.A. // Measurement Techniques, 2001. V. 44. No. 4. P. 422.

13. Г.А. Федоров, С.А. Терещенко. Расширенные троичные последовательности и двумерные кодирующие коллиматоры на их основе. Измерительная техника, 2008, №8, с.62-68.

14. Федоров Г.А., Терещенко С.А. Использование обобщённых одномерных последовательностей для построения двумерных кодов и кодирующих устройств в интегральнокодовых системах измерений. Измерительная техника, 1999, №2, с.55-63.

15. Федоров Г.А., Терещенко С.А. Интегрально-кодовые системы регистрации ионизирующих излучений. Планирование и оптимизация эксперимента // Измерительная техника, 1996. № 9. с. 50-58; Fedorov G.A., Tereshchenko S.A. // Measurement Techniques, 1996. V. 39. No. 9. P. 965.

16. Федоров Г.А., Терещенко С.А. Расширенные псевдослучайные последовательности и двумерные кодирующие коллиматоры на их основе // Измерительная техника, 2007. №6. C. 66; Fedorov G.A., Tereshchenko S.A. // Measurement Techniques, 2007. V. 50. № 6. P. 681.

17. Brown C. Multiplex imaging with multiple-pinhole cameras //J. of Applied Physics, 1974, Vol.45, No.4, p.1806-1811.

18. Cook W.R., Finger M., Prince T.A., Stone E.C. Gamma-ray imaging with a rotating hexagonal uniformly redundant array // IEEE Tr. on Nuclear Science, 1984. V. NS-31. No. 1. P. 771-775

19. Федоров Г.А., Терещенко С.А., Антаков М.А., Бурнаевский И.С. Аппаратные функции интегрально-кодовых систем с многопинхольными гексагональными кодирующими коллиматорами // Измерительная техника (в печати)

20. Gottesman S.R., Fenimore E.E. New family of binary arrays for coded aperture imaging // Applied Optics, 1989. V. 28. No. 20. p. 4344-4352.

21. Р.Л. Блейк, А.Дж. Бурек, Е. Фенимор, Р. Поттер.Мультиплексорная камера-обскура для рентгенографического исследования Солнца/ Приборы для научных исследований, 1974, 45, №4, с.41-45.

22. Best linear decoding of random mask images / T.W. Woods, M.P. Ekstrom, T.M. Palmieri, R.E. Twogood.- IEEE Trans. Nuclear Science, 1975, NS-22, N2, p. 379-383.

23. Moffet A.T. Minimum-redundancy linear arrays. IEEE Trans. Antennas and Propagation, 1968, AP-16,N2, p.172-175.

24. Golay M.J.E. Point arrays having compact, nonredundant autocorrelations.- Journal of the Optical Society of America, 1970, 61, p.272-273

25. Федоров Г.А., Терещенко С.А., Антаков М.А. Оптимизация интегральнокодовых систем измерений, построенных на основе расширенных псевдослучайных последовательностей. //Измерительная техника, 2010, №3, с.47.

26. Сороко Л.М. Мультиплексные системы измерений в физике. М.: Атомиздат, 1980.

27. Венгринович В.Л. Золотарев С.А. Итерационные методы томографии, Минск: Белорусская наука, 2009. - 227 с.

28. Бахвалов Н.С. Численные методы.- М.: Наука, 1973.-632 с.

29. Бут Э.Д. Численные методы.- М.: ГИФМЛ, 1959.- 240 с.

30. Терещенко С.А., Антаков М.А., Бурнаевский И.С. Аппаратные функции гексагональных конфигураций на основе нерасширенных двоичных псевдослучайных последовательностей (тезисы). - Сборник трудов конференции БИОМЕДСИСТЕМЫ 2010 - С. 370-376.

31. Федоров Г.А., Терещенко С.А., Антаков М.А., Бурнаевский И.С. Аппаратные функции гексагональных многопинхольных кодирующих коллиматоров (тезисы). - Сборник аннотаций научной сессии НИЯУ МИФИ - 2011 - С. 87.

32. Бурнаевский И.С. Круговая аппроксимация в расчётах аппаратных функций гексагональных конфигураций многопинхольных кодирующих коллиматоров, построенных на основе двоичных псевдослучайных последовательностей (тезисы). - 18-я Всероссийская межвузовская научно-техническая конференция студентов и аспирантов «Микроэлектроника и информатика-2011»: Тезисы докладов. - М.: МИЭТ, 2011 - С. 255

33. Бурнаевский И.С. Семейства аппаратных функций гексагональных кодирующих коллиматоров по типу псевдослучайной последовательности (тезисы). - Тезисы докладов участников III окружной научно-технической конференции молодых ученых и специалистов - 2011 - С. 7

34. Бурнаевский И.С. Программа визуализации параметров кодирующих коллиматоров (тезисы). - Исследования и разработки молодых ученых, студентов и аспирантов в области электроники и приборостроения: Тезисы докладов. 2011 - С. 6

35. G.A. Fedorov, S.A. Tereshchenko, M.A. Antakov, I.S. Burnaevskiy Investigations of Point-Spread Functions for Hexagonal Multi-pinhole Coded Apertures, - Proceeding of the 7th Russian Bavarian Conference - 2011 - P. 54-57

36. Терещенко С.А., Бурнаевский И.С. Влияние на свойства гексагональных кодирующих коллиматоров циклических перестановок элемнтов в псевдослучайной последовательности (тезисы), - Восемнадцатая международная научно-техническая конференция студентов и аспирантов: Тезисы докладов. Том 1. - 2012 - С. 287-288

37. Бурнаевский И.С. Аппаратные функции интегрально-кодовых систем измерений с многопинхольными гексагональными кодирующими коллиматорами, построенными на основе троичных последовательностей (тезисы), - 19-я Всероссийская межвузовская научно-техническая конференция студентов и аспирантов «Микроэлектроника и информатика-2012»: Тезисы докладов. - М.: МИЭТ, 2011 - С. 239.


Подобные документы

  • Технические системы сбора телеметрической информации и охраны стационарных и подвижных объектов, методы обеспечения целостности информации. Разработка алгоритма и схемы работы кодирующего устройства. Расчет технико-экономической эффективности проекта.

    дипломная работа [3,8 M], добавлен 28.06.2011

  • Разработка структурной схемы системы передачи данных. Конструирование кодирующего устройства для формирования сверточного кода, представление его функциональной схемы. Оценка вероятности правильного приема сообщения, закодированного рекуррентным кодом.

    практическая работа [367,6 K], добавлен 01.12.2010

  • Особенности помехоустойчивого кодирования. Основные виды избыточных кодов, их декодирующие свойства. Разработка восьмиразрядного кодирующего устройства на основе образующего многочлена с исправлением однократной ошибки. Скорость выдачи информации.

    реферат [230,9 K], добавлен 17.11.2013

  • Принцип работы кодирующего и декодирующего устройства циклического кода. Определение объема передаваемой информации. Нахождение емкости и построение диаграммы. Расчет надежностных показателей основного и обходного каналов. Выбор магистрали по карте.

    курсовая работа [769,9 K], добавлен 06.05.2015

  • Сущность кода Хэмминга. Схемы кодирующего устройства на четыре информационных разряда и декодера. Определение числа проверочных разрядов. Построение корректирующего кода Хэмминга с исправлением одиночной ошибки при десяти информационных разрядах.

    курсовая работа [1,1 M], добавлен 10.01.2013

  • Модель частичного описания дискретного канала (модель Л. Пуртова). Определение параметров циклического кода и порождающего полинома. Построение кодирующего и декодирующего устройства. Расчет характеристик для основного и обходного канала передачи данных.

    курсовая работа [1,8 M], добавлен 11.03.2015

  • Исследование и специфика использования инверсного кода и Хемминга. Структурная схема устройства передачи данных, его компоненты и принцип работы. Моделирование датчика температуры, а также кодирующего и декодирующего устройства для инверсного кода.

    курсовая работа [530,1 K], добавлен 30.01.2016

  • Помехоустойчивость как одна из важнейших характеристик современных систем передачи информации. Основные особенности построения биортогонального двоичного кода на базе матрицы Адамара. Анализ и характеристика схемы функционального кодирующего устройства.

    контрольная работа [853,8 K], добавлен 06.01.2013

  • Обработка результатов измерений искажений. Программное (теоретическое) кодирование. Разработка схемы кодирующего устройств. Значения вероятностей ошибок в комбинации. Введение избыточности для повышения верности при передаче дискретной информации.

    курсовая работа [423,6 K], добавлен 28.05.2013

  • Изучение практического применения связи новых свойств взаимных многочленов циклического кода со структурой кодового полинома и его весом. Рассмотрение схемы построение генераторов М-последовательности на основе регистров сдвига с обратными связями.

    реферат [136,4 K], добавлен 09.02.2010

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.