In this paper, we are proposing a novel motion correction algorithm for high-resolution OR-PAM imaging. Our algorithm combines a modified demons-based tracking approach with a newly developed multi-scale vascular feature matching method to track motion between adjacent B-scan images without needing any reference object. We first applied this algorithm to correct motion artifacts within one three-dimensional (3D) data segment of rat iris obtained with OR-PAM imaging. We then extended the application of this algorithm to correct motions to obtain vasculature imaging in the whole mouse back. In here, we stitched five adjacent 3D data segments (large field-of-view) obtained while changing the focus of OR-PAM differently for each subarea. The results showed that the motion artifacts of both large blood vessels and microvessels could be accurately corrected in both cases. Compared to the manually stitching method and the traditional SIFT algorithm, the algorithm proposed in this paper has better performance in stitching adjacent data segments. The high accuracy of the motion correction algorithm makes it valuable in OR-PAM for high-resolution imaging of large animals and for quantitative functional imaging.