This image was prepared from all or parts of 19 different Landsat-8 scenes captured in October of 2020 (from 5 different orbital paths). I used QGIS for the initial sizing and alignment of the data layers, and did the rest of the work (merging, composing to RGB, color + light correction, seamline feathering, etc.) in Affinity Photo. The relief shaded version was made in Blender using the BlenderGIS plugin. I used BlenderGIS to create a 3D mesh from 30m DEM data and then rendered it with a 30m version of the RGB satellite image image draped on top.