Coordinate mesh grid
lucadimascolo opened this issue · comments
Luca Di Mascolo commented
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?