A multiparameter high precision concurrent estimation method in image subpixel matching in which in an N-dimensional similarity space where N (N is an integer of four or greater) correspondence parameters between images are axes, said N correspondence parameters between images are estimated concurrently by using an N-dimensional similarity value between images obtained at discrete positions, comprising the steps of:
- determining a sub-sampling position where said N-dimensional similarity value between images is maximum or minimum on a line in parallel with a certain parameter axis, and determining an N-dimensional hyperplane that most approximates determined said sub-sampling position;
- determining N of said N-dimensional hyperplanes with respect to each parameter axis;
- determining an intersection point of N of said N-dimensional hyperplanes;
- and
- setting said intersection point as a sub-sampling grid estimation position for said correspondence parameter between images that gives a maximum value or a minimum value of N-dimensional similarity in said N-dimensional similarity space;
- and
- displaying an image deformed by said N correspondence parameters between images that were estimated.