I made a point shapefile (using PySHP in Python) that represent objects detected through AI.<br> This shapefile has no CRS and each point is X,Y relative to the image it was processed from (actually a TIFF converted to JPG)<br> ex : point1 = (200,300) ; point2 = (600, 500) [![point shapefile][1]][1] I have the original TIFF and I can extract its information using some Python library (Rasterio):<br> tiff shape : (7579, 13887)<br> EPSG:4326<br> BoundingBox(left=-73.3363245134031, bottom=45.18278294324435, right=-73.33487807770662, top=45.18357239348942)<br> [![enter image description here][2]][2] What is the correct process to build the point shapefile in order for its extent to be the same as the original TIFF?<br> -Should I modify the EPSG and boundaries of the shapefile so it recover the TIFF?<br> -Should I reproject each point while creating the SHP using pySHP? It seems this is the legitimate method, but how can I recalculate each point coordinates knowing only the EPSG and the boundaries of my TIFF?<br> [1]: https://i.sstatic.net/ufYln.png [2]: https://i.sstatic.net/ohn0p.png