FITS files of Fields contain the science data in SDSS images, one per band (or filter).
Using PyFITS, you can select the image data for any region within a 1489 x 2048 pixel FITS image. Five times (one for each band/filter).
But how to select the (five) regions so they are each centered on the same (RA, Dec) point (within half a pixel)? ???
I'm guessing you can work out the position, in pixel coordinates, of a particular (RA, Dec) point using the Header fields "RA" ("1st row - Right ascension of telescope boresigh"), "Dec" ("1st row - Declination of telescope boresight (d"), "SPA" ("1st row - Camera col position angle wrt north ("), etc ... but how? ??? ???
I just posted this in the Galaxy Zoo forum. I'll be interested to see how the responses differ ...