clipping aerial photographic with POLYGONS

[zoometh]

👋 @AhmedMAMahmoud, @nsheldrick, @yourvick

In order to prepare the learning/training dataset we need to use spatial extensions of CVNS recorded in these SHP to clip a high definition aerial photography.

• complete dataset, download: https://github.com/eamena-project/eamena-gee/raw/main/data/cvns-polygons.zip
• only 3 CVNS, download: https://github.com/eamena-project/eamena-gee/raw/main/data/cvns-polygons-n3.zip

Apart from GEE, we see two ways to do it:

1. Create a QGIS atlas with 1 page = 1 CVNS, see this QGIS project, download: https://github.com/eamena-project/eamena-gee/raw/main/ML-CVNS/cvns-polygons.qgz
2. use gdal_translate, like:

gdal_translate -of VRT -projwin 1774709 6166823 1868433 6112638 "<GDAL_WMS><Service name='TMS'><ServerUrl>https://tile.openstreetmap.org/${z}/${x}/${y}.png</ServerUrl></Service><DataWindow><UpperLeftX>-20037508.34</UpperLeftX><UpperLeftY>20037508.34</UpperLeftY><LowerRightX>20037508.34</LowerRightX><LowerRightY>-20037508.34</LowerRightY><TileLevel>18</TileLevel><TileCountX>1</TileCountX><TileCountY>1</TileCountY><YOrigin>top</YOrigin></DataWindow><Projection>EPSG:3857</Projection><BlockSizeX>256</BlockSizeX><BlockSizeY>256</BlockSizeY><BandsCount>3</BandsCount><UserAgent>QGIS</UserAgent><Cache /></GDAL_WMS>" /home/archesadmin/vienna.vrt

directly copied from the web

My QGIS (solution 1) is associated to my Python and the latter is very buggy (conflict of versions). I only start to go through gdal_translate (solution 2)