I have a Sentinel 2A scene (downloaded from Sentinel AWS storage in UTM projection) and coordinates (a polygon of lat/lng) of a specific plot of land. How can I extract a piece corresponding to that field (preferably in EPSG 3857, so that I could overlay it on Google Maps) a piece corresponding to my plot of interest (or to a given Google Maps tile), using either Python or CLI tools?
I have used
gdalwarp to transform the scene image to EPSG:3857, but how can I then transform the lat/lng coordinates of my plot's borders (or the given tile's borders) to the coordinates of the pixels I need to crop? Are there any ready-made utilities or libraries for that?