Mathematical research breakthroughs, in the domain of image processing filters, were recently accomplished by combining non-local means smoothing operations with polynomial regression techniques. We show how inherently noisy speckle interferograms can benefit from these methods, which have been shown to achieve near-optimal noise removal, while maximally preserving edge information corresponding to phase jump locations in optically captured phase fields.