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