<div dir="ltr"><div dir="ltr">Hello.<div><br></div><div>I am trying to merge a few tiff images (each has 4 bands) in R and It was partly successful.</div><div><br></div><div>However a look (in QGIS) at the merged output revealed some sort of radiometric distortions in the output file.</div><div><br></div><div>I would appreciate help in this regard.</div><div><br></div><div>Please find below my source code in R.</div><div><br></div><div>I attached to this email also (for comparison) one of the image before merging and the same image section in the merger output.</div><div><br></div><div>Thanks - Adebola Hassan</div><div><br></div><div>#Rcode</div><div><div>#list all tif files in working directory</div><div>all_tiff_list = list.files(".", pattern = "tif", full.names = T)</div><div><br></div><div>#plot first img only</div><div>all_tiff_imgs1 = stack(all_tiff_list[1])</div><div>plot(all_tiff_imgs1)</div><div><br></div><div>#merge and plot first and second tiff file</div><div>all_tiff_imgs2 = stack(all_tiff_list[2])</div><div>tiff_merge_1_2 = merge(all_tiff_imgs1, all_tiff_imgs2)</div><div>plot(tiff_merge_1_2)</div><div><br></div><div>#merge and plot all tiff files</div><div>img_list = list()</div><div>for (i in 1:length(all_tiff_list)) {</div><div>  img_list[i] = stack(all_tiff_list[i])</div><div>}</div><div><br></div><div>img_list$filename <- 'test.tif'</div><div>img_list$overwrite <- T</div><div>all_img_merge = do.call(merge, img_list)</div><div>plot(all_img_merge)</div><div>ggRGB(all_img_merge, r = 1, g = 2, b = 3)</div></div><div><br></div><div>#image after merge</div><div><img src="cid:ii_jqrzjxhu0" alt="Merged_Image.PNG" width="452" height="279"><br></div><div>#before merge</div><div><img src="cid:ii_jqrzjxln1" alt="Single_Image.PNG" width="452" height="281"><br></div></div></div>