I am using python and owslib to request a WMTS server for tiles in a given region using EPSG:4326 CRS.
I assemble the tiles with PIL by simply pasting them next to each other.
The image I got seems to have projection problem. Especially when I compare it with the same region I got from a WMS server with the same CRS.
The first image is the WMTS built by rendering tiles: Don't mind the black border. Also the scale is wrong here. Each layer shows the same distorsion. The band is centered on (48.387943, -4.398273) with a zoom level is 17.
What am I doing wrong ? Is there something else to do that simply pasting the tiles ?