64 lines
2.0 KiB
Python
64 lines
2.0 KiB
Python
from django.core.management.base import BaseCommand, CommandError
|
|
|
|
from datetime import date
|
|
import datetime
|
|
from django.utils import timezone
|
|
import astropy
|
|
from astropy.io import ascii
|
|
import pandas as pd
|
|
import pymysql
|
|
from sqlalchemy import create_engine
|
|
|
|
from heasarc.tdat import tDat
|
|
from heasarc.models import HeasarcTable, TableColumn, HeasarcObjectClass, NSIDE_SOURCES, ORDER, HeasarcXrayMaster, HeasarcXMMSSC, HeasarcBase
|
|
from srgcat.models import ArtCat, ArtSource
|
|
from astropy.table import Table
|
|
from astropy_healpix import HEALPix, neighbours
|
|
from astropy.coordinates import SkyCoord # High-level coordinates
|
|
from astropy.coordinates import ICRS, Galactic, FK4, FK5 # Low-level frames
|
|
from astropy.coordinates import Angle, Latitude, Longitude # Angles
|
|
import astropy.units as u
|
|
|
|
from astropy.config import reload_config
|
|
from astropy.config import get_config_dir
|
|
|
|
class Command(BaseCommand):
|
|
help = 'Initiates data dase'
|
|
|
|
# def add_arguments(self, parser):
|
|
# parser.add_argument('poll_id', nargs='+', type=int)
|
|
|
|
|
|
def handle(self, *args, **options):
|
|
|
|
print(get_config_dir())
|
|
return
|
|
|
|
hp = HEALPix(nside=NSIDE_SOURCES, order=ORDER, frame=FK5())
|
|
|
|
try:
|
|
src = ArtSource.objects.get(name='NRTA J175524.3-260347')
|
|
except:
|
|
print('Not Found')
|
|
return
|
|
|
|
print('SRC:',src,src.healpix)
|
|
|
|
try:
|
|
refs = HeasarcXMMSSC.objects.filter(name__contains="J175525.1-260403")
|
|
except:
|
|
print("no XMMSSC found")
|
|
print("***")
|
|
#print("NRTA J175524.3-260347",src.healpix,src.ra, src.dec)
|
|
around = neighbours(src.healpix, NSIDE_SOURCES, order=ORDER)
|
|
for index in around:
|
|
sc=hp.healpix_to_skycoord(index)
|
|
print(index,sc.ra.deg, sc.dec.deg)
|
|
|
|
for ref in refs:
|
|
print("ref ***")
|
|
print(ref.name,ref.healpix,src.ra, src.dec)
|
|
|
|
|
|
self.stdout.write(self.style.SUCCESS('Done'))
|