<div dir="auto"><div style="font-family:sans-serif;font-size:12.8px" dir="auto"><div style="width:380px;margin:16px 0px"><div><div dir="ltr"><div>Hi,</div><div><br></div><div>I'm using gdal to reproject a set of web-mercator tiles to Quadralateralized-Spherified-Cube (QSC) for displaying data like LandSat on a 3D globe.</div><div><br></div><div>With some difficulty, I got this working. My work flow is:</div><div>1 - georeference the tiles using translate<br></div><div>2 - build a vrt from the georeferenced tiles.<br></div><div>3 - use gdalwarp to reproject (this yields 6 large tif files, one for each face of the cube).<br></div><div>4 - use gdal2tiles to dice up the resulting files into tiles.</div><div><br></div><div>My question is: how do I approach this with a very large tile set? (level 4 tiles result in almost 3GB tifs)</div><div><br></div><div>One idea is to use the -te parameter to limit the extent of the output of warp so that the output is in smaller chunks, then build another vrt to pass to gdal2tiles.</div><div>This approach would only work if warp uses the extent to limit the files that it reads, otherwise I'll run out of memory.</div><div><br></div><div>Is there are better way?<br></div><div><br></div><div>Thanks in advance for any help you can offer.</div><div style="color:rgb(136,136,136)"><div><br></div><div>Kurt<br></div></div></div></div></div><div style="height:0px"></div></div><br></div>