Mosaicking satellite datasets with GDAL VRT

February 11, 2019 | category 👉 GDAL, image processing, python, remote sensing

Share on:


Mosaicking of remotely sensed datasets are essential requirement when our area of interest covering two or more of image footprints.

As mentioned in my previous post, GDAL VRT can help us to save lots of space, and very helpful to speed-up your processing work-flow.

GDAL is like our swiss army knife, small and versatile. so here are tips on how to use DAL VRT to help you in your projects either it big or small.

1. gdavbuildvrt for mosaicing satellite images.

If your area of interest covering two satellite’s scene, you can run gdalbuildvrt from command line.

gdalbuildvrt mosaic.vrt scene1.tif scene2.tif

In case you have more images to mosaic, you can list all the images into *.txt file.

gdalbuildvrt -input_file_list image_scene_list.txt mosaic.vrt

2. Merge satellite images into RGB or composite.

It is possible to create RGB image, by using different combination of bands. Let say we have landsat 7 dataset which have Band1, Band2, Band3, Band4, etc… in our dataset and we need to merge them to create a new composite image. Remember the combination will be Red, Blue, and Green. So, when we run gdalbuildvrt it will be like this

gdalbuildvrt RGB_Image.vrt Band3.tif Band2.tif Band1.tif

To create composite images, we add -separate into the command line.

gdalbuildvrt -separate Composite_image.vrt Band1.tif Band2.tif band3.tif Band4.tif Band5.tif

or, just make a list of your images into *.txt file

gdalbuildvrt -separate -input_file_list image_scene_list.txt Composite_image.vrt

If you want to learn more about gdalbuildvrt you can visit here and here.

⤎ back to posts