161def query_disc(nside, ra, dec, max_rad, min_rad=0):
162 """Get the list of healpix indices within max_rad, min_rad given in radians
163 around ra,dec given in degrees
168 Resolution of the healpixels to search/return.
172 Declination in degrees
174 Max distance in radians to search nearby healpixels.
175 min_rad : `float`, optional
176 Minimum distance in radians to search healpixels. Default = 0.
178 ra = np.atleast_1d(ra)
179 dec = np.atleast_1d(dec)
181 max_rad_deg = np.rad2deg(max_rad)
184 [hpg.query_circle(nside, a, b, max_rad_deg, nest=
False)
185 for (a, b)
in zip(ra, dec)])
187 if min_rad > 0
and len(pixels) > 0:
189 min_rad2 = min_rad**2
191 dsq = np.sum((vecs - vec0)**2, axis=1)
192 match = dsq > min_rad2
193 pixels = pixels[match]