2D feature-based alignment of orthophotos using pc_align

99 views
Skip to first unread message

Benjamin Hemingway

unread,
Sep 14, 2020, 5:31:37 PM9/14/20
to Ames Stereo Pipeline Support

The documentation indicates that pc_align supports feature-based alignment (through hillshading DEMs). Similarly, I would like to align overlapping orthophotos using 2D feature-based alignment. Can this be accomplished in ASP?

Alexandrov, Oleg (ARC-TI)[KBR Wyle Services, LLC]

unread,
Sep 14, 2020, 5:58:09 PM9/14/20
to Benjamin Hemingway, Ames Stereo Pipeline Support
I guess it could. It would interpret each orthophoto as a DEM (if it is a geotiff) with each grayscale value a height, create a hillshaded image from it, find features in those, and align it that way.

Yes, inane this is. Ideally you would want features to be found in othophotos, and not in their hillshaded versions. You can look at how pc_align does its job. First it find matches among hillshaded images, then computes a transform based on them. You could in principle use the same calls to ipfind and ipmatch to find matches between the orthoimages themselves, and then pc_align can be told to use a given match-file you compute via the --match-file option. It would then apply whatever transform it finds onto one orthoimage to bring it on top of another one.

Not sure how it will work, but if you have no better tools you can give it a try. In principle, one could cook up something with python and opencv avoiding ASP altogether. Or maybe somebody already wrote something.



From: ames-stereo-pi...@googlegroups.com <ames-stereo-pi...@googlegroups.com> on behalf of Benjamin Hemingway <blhe...@asu.edu>
Sent: Monday, September 14, 2020 2:31 PM
To: Ames Stereo Pipeline Support <ames-stereo-pi...@googlegroups.com>
Subject: [EXTERNAL] 2D feature-based alignment of orthophotos using pc_align
 

The documentation indicates that pc_align supports feature-based alignment (through hillshading DEMs). Similarly, I would like to align overlapping orthophotos using 2D feature-based alignment. Can this be accomplished in ASP?

--
You received this message because you are subscribed to the Google Groups "Ames Stereo Pipeline Support" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ames-stereo-pipeline...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/ames-stereo-pipeline-support/b9838b75-8cd9-4069-9db1-08b4b2f79e1cn%40googlegroups.com.

Benjamin Hemingway

unread,
Sep 14, 2020, 6:56:11 PM9/14/20
to Alexandrov, Oleg (ARC-TI)[KBR Wyle Services, LLC], Ames Stereo Pipeline Support
Thank you, Oleg, for the helpful suggestions. I’ll will give that a shot. 

Best,
Ben

You received this message because you are subscribed to a topic in the Google Groups "Ames Stereo Pipeline Support" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/ames-stereo-pipeline-support/jM-wAfAyBOs/unsubscribe.
To unsubscribe from this group and all its topics, send an email to ames-stereo-pipeline...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/ames-stereo-pipeline-support/SA0PR09MB6393E7C8E9641EE604F8F3CF88230%40SA0PR09MB6393.namprd09.prod.outlook.com.

Matthieu Volat

unread,
Sep 15, 2020, 11:09:58 AM9/15/20
to Benjamin Hemingway, Alexandrov, Oleg (ARC-TI)[KBR Wyle Services, LLC], Ames Stereo Pipeline Support
Some time ago, I did this by running ipfind/ipmatch, then using a gdal script to read and insert the matches as GCP and generate a new affine transform I used to resample/rectify one image against another using gdalwarp. I did that with georeferenced images, I guess this could work with non georeferenced dataset by declaring a dummy georef matrix.

The concept kind of worked, my issue was as always to have a robust matching set.

--
Matthieu Volat
CNRS/Observatoire de Lyon


Shashank Bhushan

unread,
Sep 15, 2020, 12:12:14 PM9/15/20
to Matthieu Volat, Benjamin Hemingway, Alexandrov, Oleg (ARC-TI)[KBR Wyle Services, LLC], Ames Stereo Pipeline Support
Also, take a look at the package Arosics, which has wrapped up everything mentioned in this discussion above end to end for orthophotos :)
I find the package to be very nice.

Benjamin Hemingway

unread,
Sep 15, 2020, 6:09:06 PM9/15/20
to Shashank Bhushan, Matthieu Volat, Alexandrov, Oleg (ARC-TI)[KBR Wyle Services, LLC], Ames Stereo Pipeline Support
Thank you, Matthieu and Shashank. This is all really helpful. I should be able to accomplish what I am wanting to do. Thanks again.

-Ben

mikita belikau

unread,
Aug 28, 2021, 4:15:57 PM8/28/21
to Ames Stereo Pipeline Support
Guys, is there somewhere available source code for this task I can examine?
I want to do exactly the same thing - coregister several HiRISE JP2 (LBL) files. 
So far I have run ipfind/ipmatch which gave me *.match file, but the next steps are a bit fuzzy for me right now.
I understand that I need to perform affine transformation for this file and then somehow use this to create a gcp file and use it in gdal's translate/wrap commands. Is there any source code available in some public repo?

Oleg Alexandrov

unread,
Apr 6, 2022, 2:19:46 PM4/6/22
to Ames Stereo Pipeline Support
ASP added support for alignment of images. The documentation is at https://stereopipeline.readthedocs.io/en/latest/tools/image_align.html. The latest build having this functionality is at https://github.com/NeoGeographyToolkit/StereoPipeline/releases
Reply all
Reply to author
Forward
0 new messages