Sam & Hsin-Fang Chiang: this looks good (having recently added the same feature to my plotting scripts, I certainly appreciate how useful it can be!). I have one comment, though. The figure you show is for 3 widely separated visits, so it looks quite nice. However, for a large number of visits/tighter dithering pattern, labelling every CCD will result in significant clutter. May I suggest adding the ability to add CCD labels to a restricted number of visits? This could be done in many ways:
- only add labels to one visit (if the data all come from the same instrument, this will likely provide all the info you need to id the ccds on the non-labeled visits)
- only add labels if the visit is separated from any other labeled visit by some minimum amount
If you do opt to label a restricted number of visits, I also like emphasizing (darker shading) of the border for just those visits, e.g. by adding the line
pyplot.fill(ra, dec, fill=False, alpha=0.4, color=color, edgecolor=color)
when adding the CCD label.
The above is all aesthetics, so feel free to disregard if you prefer the all-or-nothing approach.
Finally, is it worth specifying in the command line help message that by "ccd numbers", you are getting what is returned by getSerial() (i.e. you could just add the word "serial"). There is often confusion about what is meant by ccd blank (where blank == number vs. name vs. id vs. serial...).