1. 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 precisely and 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.