CREATE TEMP FUNCTION
ArcSecondDistance(p1 GEOGRAPHY, p2 GEOGRAPHY,
d FLOAT64) AS (ST_DISTANCE(p1,p2) < d * 30.8874796235);
SELECT
source_id_mf,
point
FROM
`bigquery-public-data.wise_all_sky_data_release.mep_wise`
WHERE
ArcSecondDistance(point,
ST_GEOGPOINT(201.5, -2.6),60)
AND ST_CONTAINS(
ST_GEOGFROMTEXT('Polygon((201.00 -3.10,201.00 -2.10,202.00 -2.10,202.00 -3.10, 201.00 -3.10))'),point)