Astropy: ValueError: один из входных данных для search_around_sky — это скаляр в search_around_sky.

В этом смысле это «скаляр», потому что это всего лишь одна координата (которую математически мы могли бы не рассматривать как скаляр, но поскольку SkyCoord может содержать массив координат, это скаляр SkyCoord). Я действительно думаю, что немного странно, что функция search_around_sky не может принять скалярную координату, потому что в принципе она может обнаружить этот случай (на самом деле она уже делает это и вызывает конкретное исключение, которое вы получаете), но вместо того, чтобы вызвать исключение, она может попытаться чтобы «продвинуть» ваш аргумент в однокоординатный массив. Но я также частично понимаю аргументацию.

Согласно этим документам, идея состоит в том, что вы должны иметь возможность вызвать catalog.search_around_sky(some_coords, ...) и получить пару массивов индексов как в some_coords, так и в catalog, где совпадающие координаты находятся в каждом из соответствующих массивов координат. Более того, каждый из этих массивов индексов должен быть одинакового размера и формы, поскольку совпадения должны быть найдены в каждом из массивов координат.

Поэтому, если вы вводите «координаты» — это просто одна координата, а не массив, логика этой функции не имеет никакого смысла. Это также указано в документах:

Однако обратите внимание, что эта двойная индексация означает, что search_around_* не работает, если одна из координат является скаляром, потому что возвращаемый индекс не имеет смысла для скаляра:

>>> scalarc = SkyCoord(1*u.deg, 2*u.deg)  
>>> idxscalarc, idxcatalog, d2d, d3d = catalog.search_around_sky(scalarc, 1*u.deg)  
# THIS DOESN'T ACTUALLY WORK  
>>> scalarc[idxscalarc]  
IndexError: 0-d arrays can't be indexed

В результате (и в любом случае, поскольку алгоритм search_around_* неэффективен в скалярном случае) лучший подход для этого сценария — вместо этого использовать методы separation*:

 >>> d2d = scalarc.separation(catalog)  
 >>> catalogmsk = d2d < 1*u.deg  
 >>> d3d = scalarc.separation_3d(catalog)  
 >>> catalog3dmsk = d3d < 1*u.kpc  

В этом примере d2d представляет собой массив двумерных расстояний между scalarc и координатами в catalog. Тогда catalogmsk — это логическая маска тех элементов в catalog, где расстояние было меньше 1 градуса. Затем вы можете использовать эту маску для извлечения этих элементов именно из catalog, например catalog[catalogmsk]. Или вы можете полностью пропустить назначение маски промежуточной переменной и просто сделать, например, catalog[d2d < 1*u.deg].

Если вы посмотрите на реализация search_around_sky вы можете видеть, что он использует SkyCoord.separation внутри, но все еще немного сложнее, потому что он должен использовать методы поиска, чтобы эффективно находить совпадающие пары координат.

Если у вас есть только одна координата, которую вы хотите найти в пределах некоторого расстояния от координат в каталоге, гораздо проще и понятнее использовать метод SkyCoord.separation напрямую, как описано выше.

Проблема была с c=SkyCoord(ra=x*u.deg,dec=y*u.deg, frame=FK5). Который считался (не знаю почему) скалярным.

search_around_sky принимает координаты массива.

Так что он отлично работает, когда вместо c=SkyCoord(ra=x*u.deg,dec=y*u.deg, frame=FK5)
я использовал: c=SkyCoord(ra=[x]*u.deg,dec=[y]*u.deg, frame=FK5)