gdal_translate "http://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer?f=json" wms.xml -of WMS
gdalwarp wms.xml -te -1e6 -1e6 1e6 1e6 -ts 512 512 -t_srs "+proj=laea +lon_0=147 +lat_0=-42" out.tif -of COG
library(terra)
plotRGB(terra::rast("out.tif"))