Resampling techniques are commonly required in digital image processing systems. Many times the classical interpolation functions are used, i.e., nearest-neighbour interpolation and bilinear interpolation, which are prone to the introduction of undesirable artifacts due to aliasing such as moire patterns. This paper presents a novel approach which minimizes the loss of information, in a least-squares sense, while resampling between orthogonal and hexagonal lattices. Making use of an extension of 2D splines to hexagonal lattices, the proper reconstruction function is derived. Experimental results for a printing application demonstrate the feasibility of the proposed method and are compared against the classical techniques.