I have several images with overlaping areas. I try to mosaic them using rasterio.merge library. But rasterio.merge funtion replaces those overlapping pixels with "nan" values. I then tried both gdalwarp and gdalmerge. They both do the same as rasterio.merge. Do you have any suggestion for fixing this issue?
Hi Mohammad, I found the reason for this problem is that when doing mosaicing overlapping images, those libraries take the value of the top image. One solution could be to creat a vrt image and then add a pixel function to it, as shown here https://lists.osgeo.org/pipermail/gdal-dev/2016-September/045134.html. Pixel function is written using numpy, basically it sums all images and divides each pixel by the number of overlapping images for that particular pixel.
And finally, transform it to raster using gdal_translate and gdal python option set to 'yes':
gdal_translate --config GDAL_VRT_ENABLE_PYTHON YES raster.vrt raster.tif
- Badges
- Science topic
- Similar topics
- Biological Science
- Histology
- Fixatives