This work presents a new technique for improving stress analysis in photoelastic plane models. This technique is based on the determination of coefficients of approximated stress functions through the collection of isochromatic fringe data, and by the applications of boundary conditions to photoelastic plane models. The technique has as advantages the accurate determination of stress concentration factors and the stress without the need for collecting isoclinic fringe data. Another great advantage, is the immediate use of this method coupled with an image analyzer device for automatic extraction of Kt. Excellent results were found in problems where stress concentration or concentrated loading were present. The coefficients of the stress functions were determined by a BASIC program, using a Newton-Raphson technique coupled to a least square algorithm. Many example were tested to evaluate its performance. Among them are: concentrated force at a point of a straight boundary; infinite and finite plate with circular holes subjected to uniaxial tensile loads; circular disk under compressive diametral forces; bar subjected to pure tensile and pure bending loads; and a double U-notched bar under pure bending.