Разброс Pyplot на стереографической проекции базовой картыi c - PullRequest
Новая
       136

Разброс Pyplot на стереографической проекции базовой картыi c

0 голосов
/ 09 мая 2020

Мне нужны указатели, пожалуйста ...

У меня есть звездный каталог, содержащий координаты прямого восхождения и склонения 1,3 миллиона звезд, а также ряд других атрибутов. Единственные другие атрибуты, которые меня беспокоят на данный момент, - это величина и созвездие.

Этот (быстрый и грязный) код считывается в каталоге, группирует его по созвездию, а затем для каждой группы создает диаграмму рассеяния. Я буду использовать тот же механизм для группировки по другим атрибутам и построения графиков. Он сохранит график, но пока я использую plt.show () для проверки вывода. Он также назначает размер маркера величине каждой звезды - ярче = больше, общепринятое соглашение на звездных картах. 

import matplotlib.pyplot as plt
import numpy as np
from astropy.io import ascii
table = ascii.read('xxxx.ecsv')
tablegr = table.group_by('Short Constellation')
constellations = tablegr.groups.keys
sizeTable = [150, 110, 70, 50, 30, 20, 12, 7, 4, 2, 1, 0.75, 0.5, 0.5]
for i in range(0, 88):
    sizes = []
    const = tablegr.groups[i]
    mag = np.array(const['magnitude'])
    print(mag)
    for row in mag:
        sizes.append(sizeTable[int(row) + 1])
    print(sizes)
    fig = plt.figure()
    ax = plt.scatter(const['RA'], const['Dec'], marker='.', s=sizes, c='k')
    plt.title(constellations[i])
    plt.show()

Я хочу построить стереографическую проекцию c с линиями сетки, но я Я изо всех сил пытаюсь понять базовую карту.

Я знаю, что мне нужно настроить проекцию примерно так: 

map = Basemap(projection='stere', width, height, lat_ts, lat_0, lon_0, celestial=True, rsphere=180/math.pi)

Где Celestial = True заставляет астрономические соглашения о долготе.

Итак, я установил атрибуты базовой карты на основе набора данных (max / min / central RA / De c et c), и я получил график, который не кажется стерографическим c.

Если я добавлю рисованные параллели и рисованные меридианы, я получу линии сетки, но никакие данные не будут нанесены на график

Мне явно не хватает чего-то в отношении присоединения атрибутов проекции к границам данных - Любые идеи ?

Спасибо.

JF

...