A simple start would be to get any two points:
dx = abs(point2gpsx - point1gpsx)/abs(point2imagex - point1imagex)
So by using a simple ratio, find how long a single pixel is in the x direction in GPS coordinates. Repeat for y direction. Then you can calculate the edges of the image in GPS coordinates, then use a meshgrid to get the coordinates for the rest.