Optical methods in experimental mechanics are important because their results are accurate and they may be used for both full field interpretation and analysis of the local rapid variation of the stresses produced by the stress concentrators. Researchers conceived several graphical, analytical and numerical methods for the experimental data reduction. The paper presents an original computer method employed to compute the analytic functions of the isostatics, using the pattern of isoclinics of a photoelastic model or coating. The resulting software instrument may be included in hybrid models consisting of analytical, numerical and experimental studies. The computer-based integration of the results of these studies offers a higher level of understanding of the phenomena. A thorough examination of the sources of inaccuracy of this computer based numerical method was done and the conclusions were tested using the original computer code which implements the algorithm.