Как назначаются трехзначные коды обсерватории Центра малых планет МАС? Почему некоторые буквы так популярны?

Ответы на широту, долготу и высоту обсерватории Ккойрама; где такие вещи можно посмотреть? укажите на https://www.minorplanetcenter.net/iau/lists/ObsCodesF.html , где в настоящее время перечислены 2198 обсерваторий.

Первый символ из трех символьных кодов имеет 36 возможных буквенно-цифровых значений. а вторые два - только цифры. Это позволяет получить 36*10*10 = 3600 возможных кодов.

Я нарисовал те, которые используются в настоящее время, и картина интересна. Почему некоторые буквенные коды используются полностью (все 100 двузначных суффиксов заполнены), а другие полностью пусты? Алфавитных ассоциаций не заметил.

Y = от 0 до 9 — это цифры, а от 10 до 36 — от A до Z: Затенение означает, что код используется.

коды

По какой-то причине я не могу перестать сравнивать данные со знаменитым сообщением Аресибо , крошечным 1-битным растровым изображением, отправленным в космос, содержащим ДНК, аминокислоты, нашу солнечную систему и другие вкусности:

Изображение Аресибо

import numpy as np
import matplotlib.pyplot as plt
# blob is the unformatted lines from https://www.minorplanetcenter.net/iau/lists/ObsCodes.html
threes = [line[:3] for line in blob.splitlines()]
key = '0123456789ABCDEFGHIJKLMNOPQRSTUVWXYZ'
a = [[x] for x in key]
for t in threes:
    a[key.find(t[0])].append(t)
b = np.zeros((36, 100))
for i, thing in enumerate(a):
    c = [int(x[1:]) for x in thing[1:]]
    for d in c:
        b[i, d] = 1
plt.imshow(-b, interpolation='nearest', cmap='gray')
plt.gca().set_aspect(2)
plt.show()
pairs = []
todegs = 180/np.pi
for line in blob.splitlines():
    try:
        cos, sin = [float(x) for x in (line[13:21], line[21:30])]
        lat = todegs * np.arctan2(sin, cos)
        lon = float(line[4:13])
        pairs.append([lon, lat])
    except:
        pass
print(len(pairs), len(blob.splitlines()), float(len(pairs))/len(blob.splitlines()) )
lon, lat = np.array(list(zip(*pairs)))
lon = np.mod(lon+180, 360) - 180
plt.figure()
plt.plot(lon, lat, 'ok', markersize=1)
plt.ylim(-90, 90)
plt.xlim(-180, 180)
plt.gca().set_aspect('equal')
plt.show()
Я подозреваю, что только Гарет Уильямс из MPC знает наверняка, но я предполагаю, что последние сгруппированы по долготе. Когда я подал заявку на получение трех кодов площадок для наших метровых телескопов в CTIO в Чили, мы получили W85-W87. DECam на той же горе позже получил W84, и когда мы подали заявку на коды для более новых 0,4-метровых телескопов (в разное время мы получили W89 и W79. Аналогично на Гавайях, мы получили T03 и T04, более поздние коды были T05-T14).
@astrosnapper, спасибо! Я проведу корреляционный анализ. Хотя широта недоступна напрямую в этом списке, возможно, я могу использовать коэффициенты параллакса. Могу ли я использовать коэффициенты параллакса для обсерваторий в качестве прокси для широты, используя arctan2 (sin, cos)?
По какой-то причине «Если бы все три были буквами, было бы 46656 возможных кодов» было удалено из моего вопроса, поэтому я добавлю его сюда в качестве комментария.

Ответы (1)

Я спросил в центре малых планет, как были определены коды, и ответ был

Исторически коды обсерваторий присваивались по возрастанию долготы к востоку (от нулевого меридиана): 360 градусов делились на числа. Когда трехзначных числовых кодов было недостаточно, буквы плюс две цифры снова использовались в полосах к востоку. Некоторые полосы уже заполнены (например, H**). Новый код обсерватории присваивается в доступном диапазоне букв и цифр (в зависимости от ее долготы).

Как назначаются трехзначные коды обсерватории Центра малых планет МАС? Почему некоторые буквы так популярны?
СпросиСетьПолитика конфиденциальности1y, 0v