Quick question.
get_coords_meshgrid computes the angular distances between all the points on a given grid and its centre. Since these are angular separations and assuming a small-angle approximation, shouldn't the expression for the distance along the R.A. direction,
x = (np.linspace(0.5, -0.5 + 1./float(ncol), ncol, dtype=real_dtype)) * ncol * dxy
include a corrective cos(Dec0), with Dec0 given by the declination of the grid centre?
Quick question.
get_coords_meshgridcomputes the angular distances between all the points on a given grid and its centre. Since these are angular separations and assuming a small-angle approximation, shouldn't the expression for the distance along the R.A. direction,include a corrective
cos(Dec0), withDec0given by the declination of the grid centre?