My scenario consist of 4 million geolocated images each measuring 1536x1536 pixels in size.
What is the most efficient way to pyramid these images?
My end goal is to serve these images using geoserver via WMS requests.
I have tried creating a virtual dataset of all images and putting them into gdal_retile.py. I have also tried merging each tile with a transparent "null" image that covers the globe then pyramid the result. Both methods seem like they will take forever and make me think there is a better way.
I also have the requirement of adding additional images in the future.
Is there an easy way to do this without having to regenerate everything again?
I am new to geoserver and gdal.