Transform pixel coordinates (in FITS file) to equatorial Say I have a FITS file with an observed field of stars, each one with x,y coordinates given in pixels. Assume the image is aligned so that the x axis is the right ascension, and the y axis the declination.

I'm looking for a tutorial to transform the pixel coordinates of each star to equatorial, even if the transformation is approximate.

I'v been looking for a while but found no tutorial/book/site that enumerates the steps needed to produce this transformation.
Assuming that each (square) pixel has the same angular scale (not a given if the field of view is large) of $\theta$ degrees/pixel.

Then the declination (in degrees): $\delta \simeq \delta_0 + \theta y$ 

The right ascension (in degrees): $\alpha \simeq \alpha_0 - \theta x/\cos \delta$

where $\alpha_0$ and $\delta_0$ are the RA and Dec at $x=0$, $y=0$. [The minus sign is there because right ascension increases towards the left of a sky image.]

This approximatimation becomes poor as the field of view gets larger.

A slightly more complex, but accurate, approach is described in http://gtn.sonoma.edu/data_reduction/astrometry.php