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