91 lines
4.0 KiB
Python
91 lines
4.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
|
|
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
|
|
from srglib.utils import load_skymap_sources_for_artsurvey
|
|
|
|
import glob
|
|
import os
|
|
import logging
|
|
|
|
|
|
class Command(BaseCommand):
|
|
help = 'Initiates data dase'
|
|
|
|
def handle(self, *args, **options):
|
|
logging.getLogger().setLevel(logging.DEBUG)
|
|
|
|
#path='/export/django/srg/data/Rodion_GC/18Apr2022'
|
|
#load_skymap_sources_for_artsurvey("{}/src_list_0.fits".format(path),'GC.E0',force=True)
|
|
#return
|
|
|
|
#path='/export/django/srg/data/Andy_L20'
|
|
#load_skymap_sources_for_artsurvey("{}/cat_lp20.fits".format(path),'L20.E0',force=True)
|
|
#load_skymap_sources_for_artsurvey("{}/lp20_20220307.fits".format(path),'L20.E0',force=True)
|
|
|
|
#path='/export/django/srg/data/Galplane_l20'
|
|
#load_skymap_sources_for_artsurvey("{}/src_list_0.fits".format(path),'L20rb.E0',force=True)
|
|
|
|
#path='/export/django/srg/data/Galplane_l20'
|
|
#load_skymap_sources_for_artsurvey("{}/src_list_0.fits".format(path),'L20rb_ml.E0',force=True)
|
|
|
|
#path='/srg/a1/work/oper/staff/kate/skymap/catalog_artdays/'
|
|
#load_skymap_sources_for_artsurvey("{}/22y_cat_0.fits".format(path),'22y.E0',force=True)
|
|
#load_skymap_sources_for_artsurvey("{}/22y_cat_5.fits".format(path),'22y.E5',force=True)
|
|
|
|
path='/srg/a1/work/kate/skymap_artdays/catalog/'
|
|
load_skymap_sources_for_artsurvey("{}/22y_cat_5.fits".format(path),'22y.E5',force=True)
|
|
|
|
|
|
#load_skymap_sources_for_artsurvey("{}/scan_tot_cat.fits".format(path),'L20.E0',force=True)
|
|
#load_skymap_sources_for_artsurvey("{}/scan_0_cat.fits".format(path),'L20p1.E0',force=True)
|
|
#load_skymap_sources_for_artsurvey("{}/scan_1_cat.fits".format(path),'L20p2.E0',force=True)
|
|
#load_skymap_sources_for_artsurvey("{}/scan_2_cat.fits".format(path),'L20p3.E0',force=True)
|
|
#load_skymap_sources_for_artsurvey("{}/scan_3_cat.fits".format(path),'L20p4.E0',force=True)
|
|
|
|
|
|
force=True
|
|
ebands=['0','1','2','3','4']
|
|
#fields=['040','025','020','015','010','005','000','355','350','345']
|
|
#fields=['070','075','080','085','090',]
|
|
#fields=['270','275','280','285','290','295','300','305','310','215','320','325']
|
|
#fields=['315']
|
|
fields=['330','335','340']
|
|
path='/srg/a1/work/rodion/SRG/art-xc/data/galmap/field_data/'
|
|
"""
|
|
for field in fields:
|
|
for eband in ebands:
|
|
load_skymap_sources_for_artsurvey("{}/{}/1g/cat_1g_{}.fits".format(path,field,eband),
|
|
"GP{}.E{}".format(field,eband),
|
|
force=force)
|
|
"""
|
|
self.stdout.write(self.style.SUCCESS('Done'))
|