use pixel_scale_matrix in celestial_scale util