@mikita87 another way to do this is if you do not care about exact precision cosi corr may provide, is to use ipfind/ipmatch from ames stereo pipeline to define ground control points to then use GDAL warp to align the images, roughly approximating the auto geo-registration tool in arcgis pro. I’ve implemented it here: https://github.com/AndrewAnnex/asap_stereo/blob/8e2f38b7144d0e972320a7ca9e3a0287c2b3afca/asap_stereo/asap.py#L1752, all you would need to do is select one image as a base reference image and then align your N many hirise images to it. The “im_feeling_lucky” function will do this all in one command. There will still be some distortion between the images from viewing angles so you may have to orthorectify all of your data against a high resolution DEM first, then use ASAP’s method.