How to equalize an image via reference points?

Hi there,

My problem: I've got a generated image of the a landscape (see attachment). I've compared my image to GPS-reference data for several points (>100). There is a distortion. I try to apply something to compensate this distortion.

Approaches: Imwarp seems to be the way to do it, but I'm not sure how. B = imwarp(A, tform), but how do I get my points into tform (or the displacement field)? Basically, I've got a 2d-matrix with 4 columns: x_GPS, y_GPS, X_image, Y_image. To get from these singular points to a geometrical transformation, I need to interpolate over the whole field. I guess it should be something with optical flow, but I'm not sure about how to proceed.

Any suggestions are highly appreciated.

[Image Analyst]

Give the coordinates, and upload the image without the stars on it. Otherwise it's hard for us to try things.