81 lines
2.6 KiB
Python
81 lines
2.6 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
|
|
from srgcat.models import ArtCat, ArtSource
|
|
from astropy.table import Table
|
|
from astropy_healpix import HEALPix
|
|
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
|
|
import logging
|
|
from astropy.io import fits
|
|
from django.db.models import Q
|
|
from datetime import datetime
|
|
from astropy.time import Time
|
|
from astropy.io import fits
|
|
|
|
from srgcat.models import SkyMaps, SrgDataDump, SkyMapSource
|
|
from monthplan.models import SurveyHealpixPlate
|
|
from srglib.utils import TZ_MSK, load_srg_data_dumps, make_source_name, load_simbad_sources
|
|
from srglib.utils import find_heasarc_counterparts, update_all_skymap_sources_counterparts, load_skymap_sources_dir
|
|
from srglib.utils import mark_new_skymap_sources_in_latest, clean_skymap_sources, update_allsky_missed
|
|
|
|
import glob
|
|
import os
|
|
import logging
|
|
|
|
|
|
class Command(BaseCommand):
|
|
help = 'Initiates data dase'
|
|
|
|
def handle(self, *args, **options):
|
|
logging.getLogger().setLevel(logging.INFO)
|
|
|
|
load_skymap_sources_dir('/srg/work/oper/staff/kate/skymap/catalog')
|
|
|
|
update_allsky_missed()
|
|
|
|
#print('mark_new_skymap_sources_in_latest')
|
|
#mark_new_skymap_sources_in_latest()
|
|
#print('update_allsky_missed')
|
|
#update_allsky_missed()
|
|
|
|
"""
|
|
skymap = SkyMaps.objects.get(SMAPNR=230147)
|
|
dump=SrgDataDump.objects.get(name='srg_20200310_163523_000')
|
|
print(skymap)
|
|
srcs=skymap.skymapsource_set.all().filter(dump=dump).order_by('-rate')
|
|
clean_skymap_sources(srcs)
|
|
"""
|
|
|
|
"""
|
|
skymap = SkyMaps.objects.get(SMAPNR=273024)
|
|
dump=SrgDataDump.objects.get(name='srg_20200311_164642_000')
|
|
print(skymap)
|
|
srcs=skymap.skymapsource_set.all().filter(dump=dump).order_by('-rate')
|
|
for s in srcs:
|
|
print(s)
|
|
"""
|
|
|
|
|
|
|
|
#load_simbad_sources(srcs,minrad=60,maxdist=120)
|
|
#clean_skymap_sources(srcs)
|
|
|
|
return
|
|
|
|
load_skymap_sources_file('/srg/work/oper/staff/kate/skymap/catalog/srg_20200302_165649_000.fits')
|
|
""" update_all_skymap_sources_counterparts() """
|
|
self.stdout.write(self.style.SUCCESS('Done'))
|