Locally optimal Delaunay triangulations are constructed to improve previous image approximation schemes. Our construction relies on a local optimization procedure, termed exchange. The efficient implementation of the exchange algorithm is addressed, and its complexity is discussed. The good performance of our improved image approximation is illustrated by numerical comparisons.