Цели. Цель работы состояла в создании эффективного итерационного алгоритма для томографической реконструкции объектов с большими объемами исходных данных. В отличие от сверточного алгоритма проецирования, широко используемого в коммерческих промышленных и медицинских томографах, алгебраические итерационные методы реконструкции используют значительные объемы памяти и характеризуются большими временными затратами на реконструкцию. В то же время итерационные методы позволяют решать более широкий круг диагностических задач, где требуется большая точность реконструкции, а также в случаях использования ограниченного объема данных при малоракурсной съемке или съемке с ограниченным угловым диапазоном.Методы. Особенностью созданного алгоритма является использование полярной системы координат, в которой проекционные системные матрицы инвариантны по отношению к вращению объекта. Это дает возможность значительно сократить объемы памяти для хранения проекционных матриц и использовать для реконструкции графические процессоры. В отличие от простой полярной системы координат, используемой ранее, нами была использована система координат с дихотомическим делением поля реконструкции, что позволяет обеспечить инвариантность к вращениям и в тоже время достаточно равномерное распределение пространственного разрешения по полю реконструкции.Результаты. Был разработан алгоритм реконструкции, основанный на использовании парциальных системных матриц, соответствующих дихотомическому делению поля изображения на парциальные кольцевые области реконструкции. С использованием цифровых фантомов Шеппа – Логана и Де Фриза были исследованы особенности работы предложенного алгоритма реконструкции и показана его применимость для решения томографических задач.Выводы. Предложенный алгоритм дает возможность реализовать алгебраическую реконструкцию изображения с использованием стандартных библиотек для работы с разреженными матрицами на базе настольных компьютеров с графическими процессорами.