Hi,
Okay, so you have a set of stars with known RA,Dec, and a set of known times t_i and pixel positions x_i, y_i. From RA,Dec, t_i and your GPS location, you can compute alt_i, az_i.
Now you have x_i, y_i, alt_i, az_i. Choose some arbitrary time t_r. Take t_r, your GPS location, alt_i, and az_i, and compute ra_i, dec_i.
Now you have x_i, y_i, ra_i, dec_i -- a bunch of predicted RA,Dec positions that would project to the x_i,y_i locations if you observed at time t_r.
Feed that list of x_i, y_i, ra_i, dec_i into the Astrometry.net fit-wcs routine, or your favourite other fitting routine -- it's just a non-linear least squares fit, no decisions to make.
Or do I have something wrong?
cheers,
dustin