The paper describes an algorithm of photoelastic tomography and its application for residual stress measurement in glass articles of complicated shape. The algorithm is based on a linearized solution of the equations of integrated photoelasticity. The problem of tensor field tomography is decomposed into several problems of scalar field tomography for normal stress components of the stress tensor. The method is implemented with an automated polariscope supplied with a rotary stage. Software for the control of tomographic measurements and for the calculation of the stress fields has been elaborated. Several examples illustrate application of the method.