Вычислительный томограф

СОЮЗ СОВЕТСКИХСОЦИАЛИСТИЧЕСКИХРЕСПУБЛИК 801287847 594 А 61 В 60 ГОСУДАРСТВЕННЫЙ КОМИТЕТ СССР ПО ДЕЛАМ ИЗОБРЕТЕНИЙ И ОТКРЫТ К А ВТОРСКОМУ СВИДЕТЕЛЬСТВ матор 10, сдвигающий регистр 12, коммутатор 11, блок 13 свертки, блок 18 обратной проекции, видеоконтрольное устройство 27.Источник 2 и детекторное устройство 3 приводят в поступательное движение относительно неподвижного объекта 4, не меняя их взаимного расположения, затем поворачивают их на некоторый угол, не меняя взаимного расположения, и вновь задают им поступательное движение, и так далее до поворота на 180. Детекторное устройство 3 фиксирует количество принятого излучения. Введение новых элементов и образование новых связей между элементами устройства позволяют без увеличения аппара турных затрат обеспечить работу в режиме реального времени, при котором рекон- с струкция изображения осуществляется параллельно со сканированием и изображение ,р срезов поступает на дисплей одновременно %У ф с окончанием сканирования. 6 ил.(71) Специальное проекти технологическое бюро(54) ВЫЧИСЛИТЕЛЬНЫЙ ТОМОГРАФ (57) Изобретение предназначено для нейродиагностики и позволяет сократить время обследования за счет уменьшения количества операций умножения при свертке. Томограф содержит сканирующий блок 1 с источником 2 проникающего излучения и детекторным устройством 3, блок 5 предварительной обработки данных, накапливающий сумСАНИЕ ИЗОБРЕТЕНИЯИзобретение относится к медицинской технике и может быть использовано для нейродиагностики, а также для контроля промышленных изделий.Цель изобретения - сокращение времени обследования за счет уменьшения количества операций умножения при свертке.На фиг. 1 приведена структурная схема вычислительного томографа; на фиг. 2 - 4 - структурные схемы блока предварительной обработки данных, блока свертки и блока обратной проекции соответственно; на фиг, 5 и 6 - схематически представлены процессы формирования проекции и реконструкции изображения (обратной проекции).Вычислительный томограф (фиг. 1) содержит сканирующий блок 1, в состав которого входит источник 2 проникающего излучения и детекторное устройство 3, между которыми размешен исследуемый объект 4.Выход сканирующего блока 1, которым является выход детекторного устройства 3, соединен с входом блока 5 предварительной обработки сигналов - с входом интегратора 6 (фиг. 2).Блок 5 предварительной обработки сигналов содержит соединенные последовательно интегратор 6, аналого-цифровой преобразователь 7, буферный блок 8 памяти и логарифматор 9. Выход логарифматора 9, являющийся выходом блока 5 предварительной обработки сигналов, соединен с входом накапливающего сумматора 10 (фиг, 1) и с первым входом коммутатора 11. Выход накапливающего сумматора 10 соединен с входом сдвигающего регистра 12, выход которого соединен с вторым входом коммутатора 11, выход которого соединен с входом блока 13 свертки - с входом входного блока 14 памяти (фиг. 3).Выход блока 14 памяти соединен с первым входом устройства 15 умножения, второй вход которого соединен с выходом блока 16 памяти ядра свертки, Выход устройства 15 умножения соединен с входом накапливающего сумматора 17, выход которого, являющийся выходом блока 13 свертки, соединен с входом блока 18 обратной проекции, - с первым входом коммутатора 19 блока 20 памяти модифицированных проекций (фиг. 4).Блок 21 памяти констант, также входяший в состав блока 18 обратной проекции, соединен через накапливающий сумматор 22 с адресным входом 23 блока 20 памяти модифицированнь 1 х проекций. Выход блока 20 памяти модифицированных проекций соединен с первым входом накапливающего сумматора 24, выход которого соединен с входом выходного блока 25 памяти и вторым входом коммутатора 19, выход которого соединен с информационным входом 26 блока 20 памяти модифицированных проекций. 5 О 15 20 25 30 35 40 45 50 55 Выход выходного блока 25 памяти соединен с вторым входом сумматора 24, икроме того, является выходом блока 18обратной проекции в целом и соединяется с входом видеоконтрольного устройства 27.Устройство работает следующим образом,При работе томографа источник 2 проникающего излучения и детекторное устройство 3 сканирующего блока 1 приводят в поступательное движение относительно неподвижного исследуемого объекта4, При этом детекторное устройство 3 неменяет своего положения относительно источника 2,Затем источник 2 и детекторное устройство 3 поворачивают на некоторый угол(например, 1), также не меняя взаимногорасположения, и вновь совершают поступательное движение и т.д. до поворотана 180.При всех этих перемещениях детекторное устройство 3 регистрирует количествопринятого излучения и преобразует его вэлектрический ток пропорциональной величины.Интегратор 6 (фиг. 2) формирует на выходе электрическое напряжение, пропорциональное интегралу величины тока за некоторое время (например, 1 мс), соответствующее перемещению источника 2 и детекторного устройства 3 на отрезок, соответствующий размеру элемента изображения (например. 1 мм).Аналого-цифровой преобразователь 7преобразует напряжение с выхода интегратора 6 в параллельный двоичный код,поступающий в буферный блок 8 памяти.Таким образом, после окончания очередного прохода в буферном блоке 8 памятиформируется ряд чисел, соответствуюющихвеличине поглощения полос 28 элементарной толщины в исследуемом слое объекта4, сориентированных в направлении, перпендикулярном движению источника 2(фиг. 5).Проходя через логарифматор 9, каждое изэтих чисел преобразуется в величину (число), соответствующую интегральной суммекоэффициентов поглощения на соответствующей полосе. Ряд таких чисел Ео, ЕЕ,соответствующих т параллельным полосам,полностью покрывающим исследуемый слойобъекта 4, принято называть проекциейисследуемого слоя.Числа Ео, Е 1Е очередной проекции,сформированной в блоке 5 (фиг. 1) предварительной обработки сигналов, поступаютпоследовательно на первый вход коммутатора 11, С выхода коммутатора 11 этичисла поступают во входной блок 14 памятиблока 13 свертки и запоминаются в последовательных т ячейках блока 14.После этого организуется повторное чтение содержимого буферного блока 8 памятиВ блоке 18 обратной проекции эта проекция так же, как и нулевая, поступаетпоследовательно на вход 26 блока 20 памяти модифицированных проекций. Но в отличие от нулевой модифицированной проекции каждый член первой левой модифицированной проекции поступает не в одну,а в п последовательных ячеек памяти.Так же в блоке 13 свертки формируютсявсе остальные модифицированные проекции,а их запись в блок 20 осуществляется потому же принципу с тем, чтобы в ячейкахблока 20, отстоящих друг от друга на т адресов, хранились числа, подлежащие сложению при формировании результирующей модифицированной проекции. В рассмотренномпримере заполнение ячеек памяти блока 20следует произвести так:ячейка 0: 8 (0)ячейка 1: Я (1)ячейка 2: 8 (2)ячейка 239: 5(239)ячейка 240: Я(0)ячейка 241: 5(0)ячейка 242: 8,(2)ячейка 243: Я(2)ячейка 244: К (4)ячейка 245: 5(4)ячейка 458: 8(218)ячейка 459: К (219)ячейка 460: 0ячейка 461: 0ячейка 480: 0ячейка 481: 0ячейка 482: 0ячейка 499: 0ячейка 500: Ь(20)ячейка 501: 8 (20)ячейка 502: 5 (22)ячейка 503: 5 (22)ячейка 718: 8 (238)ячейка 719: Ь (238)ячейка 720: 5 (0)ячейка 721: 8(0)ячейка 722: 8(0)ячейка 723: 8 г (0)ячейка 724: 5(4)ячейка 725: 8(4)ячейка 726: 8 (4)ячейка 727: 8(4)ячейка 916: Я (196)ячейка 917: Ь(196)ячейка 918: 8 (196)ячейка 919: 8(196)ячейка 920: 0ячейка 921: 0ячейка 922: 0ячейка 960: 0ячейка 961: 0ячейка 999: О,ячейка 1000: Ь (40)ячейка 1001: 5 (40)ячейка 1002: 8(40)ячейка .1003: Я(40)ячейка 1196: Я (236)ячейка 1197: Я ( 236)ячейка 1198: У(236)ячейка 1199: % (236)Нули в соответствующие ячейки могут быть вписаны как предварительный очисткой всего блока 20, так и непосредственно в процессе записи проекций.Такой способ расположения информации в блоке 20 позволяет просто организовать вычисление результирующей проекции. Для этого из блока 21 памяти констант в накапливающий сумматор 22 записывается адрес нулевого члена нулевой модифицированной проекции. Считанное по этому адресу число Ь(0) поступает в накапливающий сумматор 24. Затем из блока 21 в сумматор 22 поступает приращение адреса т (например, 240) и складывается с адресом, хранившимся в сумматоре 22. Считанное по вычисленному в сумматоре 22 адресу число К (О) поступает в сумматор 24 и складывается с хранившимся там числом 8(0). Затем из блока 21 в сумматор 22 снова поступает приращение т и складывается с вычисленным ранее адресом. Считанное по новому адресу число (например, 0), складывается в сумматоре 24 с хранившейся там суммой.Работа блоков 20 и 21 и сумматоров 22 и 24 в этом режиме продолжается до тех пор, пока в .сумматоре 24 не сформируется нулевой член результирующей модифицированной проекции 8 о. Это число через входной коммутатор 19, который в этом режиме передает информацию с второго входа, поступает на информационный вход 26 блока 20. Одновременно из блока 2 считывается константа (1 х 1 - 1) гп (где 1 х 1 - число модифицированных проекций в блоке 20), например К=5, (М - 1) тп = - 960. 5 10 15 20 25 30 35 40 45 50 55 Сложение этого числа с адресом, хранящимся к этому моменту в сумматоре 22, дает в сумме адрес нулевого члена нуле вой модифицированной проекции. В ячейку с этим адресом производится запись числа о, после чего к адресу в сумматоре 22 прибавляется 1 и из блока 20 в сумматор 24 поступает число 8 (1). Аналогично описанному организуется последовательное чтение из блока 20 всех чисел, подлежащих сложению с 8 (1), и полученный в результате первый член результирующей модифицированной проекции 5 поступает на место числа 8(1) в блок 20.Точно так же вычисляются и заносятся в блок 20 остальные члены результирующей модифицированной проекции.После этого в блоке 18 осуществляется обратная проекция полученной результирующей модифицированной проекции на массив изображения, формируемый в выходном блоке 25 памяти. При этом в ячейки блока 25, каждая из которых соответствует элементарному участку изображения, поступают члены результируюгцей модифицированной(фиг. 2) и с выхода логарифматора 9 таже последовательность чисел Ео, Еь ",Епоступает в накапливаюший сумматор 1 О (фиг. 1). В накапливающем сумматоре 10 формируются суммы нескольких соседних членов проекции.Числа Е 1(0), Е 1(п 1), Е(2 п) каждое из которых определяется как сумма Е 1(1 п 1) = Екп + Е(л,:+;+ + Ец 0, 1,поступают в сдвигаюший регистр, где они уменьшаются в К раз за счет сдвига вправо на й разрядов, К= 2",С выхода сдвигающего регистра 12 по- следовательность Е 1(0) /К Е(п 1) /Кь Е(2 п 1) /Кь поступает на второй вход коммутатора 11, который в этот период находится в режиме передачи на выход информации с второго входа.С выхода коммутатора 11 эта последовательность также поступает во входной блок 14 памяти блока 13 свертки и запоминается в последуюших ячейках памяти блока 14. Затем снова повторяется чтение содержимого буферного блока 8 памяти (фиг. 2) и в той же последовательности происходят формирование второй дополнительной проекции Е(О)/К, Е(п)/К, Ег(2 п)/К,и ее запоминание в блоке 14 (фиг. 3) и т. д. до тех пор, пока все дополнительные проекции видаЕо(0)/КЕ(п;)/К;, Е,(2 п;)/К,не запишутся в блоке 14.В блоке 16 ядра свертки постоянно хранится ядро, состоящее из участков: 0 - 1. - р+1, 1. - р+2 1. - 1, 1.о, 1.ь, (.р- 1;1 (левая половина) - КМ( - р+1), К 1 МХ Х( - р+ 2)К,М( - р,);1 (правая половина) - К 1 М (р 1), КМХ Х (р 1+ 1)К 1 М(р. - 1);(левая половина) К,М ( - р,1),",К;М ( - р,);1 (правая половина) - К М (р;)К М (р,.1 --- 1).В блоке 13 свертки производится свертка каждой из проекций, хранящихся в блоке 14, с соответствуюшими участками ядра из блока 16.При этом сначала формируется нулевая модифицированная проекция Формирование нулевого члена этой последовательности осуществляется в следующем порядке 25 30 35 40 45 50 55 Из блока 14 последовательно считываются члены исходной проекции Ео, Е, Еи поступают на первый вход устройства15 умножения. Одновременно из блока16 последовательно считываются члены нулевого участка ядра1.о. 1 ь 1.и поступают на второй вход устройства 15умножения. С выхода устройства 15 умножения произведенияЕо 1-о, Е 1- - 1, Ь 1,поступают в накапливающий сумматор 17, в котором формируется сумма 8(0)= Ео 1.о+ +Е 11. + Е 1. +, С выхода накапливающего сумматора 17 нулевой член 8(0) нулевой модифицированной проекции поступает в блок 18 обратной проекции.Точно так же производится формирование первого члена 8 (1) нулевой модифицированной проекции с той только разницей, что из блока 16 последовательно считываются члены нулевого участка ядра, начиная с 1.11.ь 1.о, 1Так же последовательно формируются остальные члены нулевой модифицированной проекции. Поступая с выхода блока 13 свертки в блок 18 обратной проекции, члены этой проекции попадают на первый вход входного коммутатора 19. В этот период коммутатор 19 находится в режиме передачи на выход информации с первого входа, поэтому члены нулевой модифицированной проекции последовательно поступают на информационный вход 26 блока 20 памяти модифицированных проекций.Перед началом операции свертки из блока 21 памяти констант считывается в накапливающий сумматор 22 адрес, начиная с которого должна быть записана нулевая модифицированная проекция в блок 20. С выхода накапливаюшего сумматора 22 этот адрес поступает на адресный вход 23 блока 20. После записи числа 8 (0) по этому адресу код адреса в накапливающем сумматоре 22 увеличивается на единицу и после окончания вычисления числа 8(1) в блоке 13 свертки оно записывается по этому адресу в блок 20, Точно так же последовательно записываются в блок 20 все числа нулевой модифицированной проекции.Затем в блоке 13 свертки (фиг. 3) производится формирование первой левой модифицированной проекции. Для этого из блока 14 последовательно считываются члены первой дополнительной проекции, а из блока 16 - соответствующие члены первого левого участка ядра. С выхода устройства 15 умножения произведения соответствующих членов поступают в накапливаюший сумматор 17, на выходе которого последовательно формируются члены первой левой модифицированной проекции81(0), 81(п, ), Я(2 п), 1287847проекции, причем в каждую из ячеек поступает член проекции, соответствующий полосе, на которой находится центр соответствующего элементарного участка.В числовом видезона изображения, представляющая собой круг, разбита на элементарные квадраты (фиг. 6). В блоке 25 первая ячейка памяти соответствует крайнему слева квадрату в нижнем ряду изображения. Следующие ячейки 25 соответствуют следующим квадратам того же нижнего ряда. Затем идут ячейки следующего ряда слева направо и т.д. до крайней правой ячейки в верхнем ряду.В блоке 21 памяти констант хранится предварительно вычисленный адрес ячейки блока 20 памяти модифицированных проекций, в которой хранится число, подлежащее записи в первую ячейку изображения. Далее в блоке 21 хранится приращение адреса модифицированной проекции при переходе на один квадрат вправо. Последовательным прибавлением этой константы кадресу осуществляется выбор из блока 26всех чисел, подлежащих записи в ячейки блока 25, соответствующие элементарнымквадратам нижнего ряда. Затем из блока 21 извлекается константа приращения адреса при переходе от крайней правой ячейки первого ряда к крайней левой ячейке второго ряда и из блока 20 извлекается член модифицированной проекции, соответствующий полосе, на которой находится центр крайнего левого квадрата второго ряда,Таким образом производится заполнениевсех ячеек памяти блока 25 соответствую- шими числами.Каждая из остальных проекций исследуемого слоя, формируемая в блоке 5 предварительной обработки данных, проходит такую же процедуру обработки до формирования соответствующей результирующей модифицированной проекции.Затем производится обратная проекция, несколько отличная от описанной, Поскольку проекция расположена под другим углом. используется другой набор констант из блока 21 для определения адресов блока26 памяти модифицированных проекций 1020253540 Кроме того, члены результирующей мо 45дифицированной проекции не просто записываются в соответствующие ячейки блока 25,а складываются с содержимым этих ячеек. Для этого перед записью в очередную ячейку блока 25 производится чтение содержимого этой ячейки. Считанное число с выхода блока 25 поступает на второй вход сумматора 24 и складывается с поступающим на первый вход соответствующим членом модифицированной проекции.Сумма с выхода сумматора 24 поступает в блок 25 и записывается в ту же ячейку.После окончания обратной проекции всех модифицированных проекций блок 25 содержит реконструированное изображение исследуемого слоя, Затем эта информация из блока 25 передается в видеконтрольное устройство 27, которое воспроизводит изображение в виде полутоновой картины на экране.Предлагаемый вычислительный томограф позволяет без увеличения аппаратурных затрат обеспечить работу в режиме реального времени, при котором реконструкция изображения осуществляется параллельно со сканированием и изображение срезов поступает на дисплей одновременно с окончанием сканирования.Формула изобретенияВычислительный томограф, содержащий сканирующий блок, включающий источник проникающего излучения и детекторное устройство, блок предварительной обработки сигналов, подключенный к выходу детекторного устройства, блок свертки, блок обратной проекции, подключенный к выходу блока свертки, и видеоконтрольное устройство, подключенное к выходу блока обратной проекции, отличающийся тем, что, с целью сокрашения времени обследования за счет уменьшения количества операций умножения при свертке, введены накапливающий сумматор, сдвигающий регистр и коммутатор, причем выход блока предварительной обработки сигналов соеДинен с входом накапливающего сумматора и с первым входом коммутатора, второй вход которого соединен с выходом сдвигаюшего регистра, выход коммутатора соединен с входом блока свертки, а выход накапливающего сумматора соединен с входом сдвигающего регистра.1287847 гЯ Яг б ль ер СоставитТехред И.Тираж 617 го комитета С сква, Ж - 35, афическое пред едактор Л. Гратиллоаказ 7729/5ВНИИПИ Государственно113035, МоПроизводственно-полигр Н. Галамагаес Корректор Г. РешетникПодписноеСР по делам изобретений и открытийРаушская наб., д. 4/5риятие, г. Ужгород, ул. Проектная, 4