AstroToolkit 1.0.0__py3-none-any.whl

AstroToolkit/Surveys/Gaia.py

@@ -0,0 +1,168 @@
+import pandas as pd
+import astropy.coordinates as coord
+from astroquery.vizier import Vizier
+import astropy.units as u
+
+catalogue_names=['I/345/gaia2','I/355/gaiadr3']
+
+row_limit=Vizier.ROW_LIMIT = -1
+
+def renameHeadersDR3(data):
+	#Manually rename all columns to their original column names. This one line is 6849 characters long.
+	data.rename(columns={'RA_ICRS':'ra','DE_ICRS':'dec','Source':'source_id','e_RA_ICRS':'ra_error','e_DE_ICRS':'dec_error','Plx':'parallax','e_Plx':'parallax_error','PM':'pm','pmRA':'pmra','e_pmRA':'pmra_error','pmDE':'pmdec','e_pmDE':'pmdec_error','RUWE':'ruwe','FG':'phot_g_mean_flux','e_FG':'phot_g_mean_flux_error','Gmag':'phot_g_mean_mag','FBP':'phot_bp_mean_flux','e_FBP':'phot_bp_mean_flux_error','BPmag':'phot_bp_mean_mag','FRP':'phot_rp_mean_flux','e_FRP':'phot_rp_mean_flux_error','RPmag':'phot_rp_mean_mag','BP-RP':'bp_rp','RV':'radial_velocity','e_RV':'radial_velocity_error','Vbroad':'vbroad','GRVSmag':'grvs_mag','QSO':'in_qso_candidates','Gal':'in_galaxy_candidates','NSS':'in_galaxy_candidates','XPcont':'has_xp_continuous','XPsamp':'has_xp_sampled','RVS':'has_rvs','EpochPh':'has_epoch_photometry','EpochRV':'has_epoch_photometry','MCMCGSP':'has_mcmc_gspphot','MCMCMSC':'has_mcmc_msc','And':'in_andromeda_survey','Teff':'teff_gspphot','logg':'logg_gspphot','__Fe_H_':'mh_gspphot','Dist':'distance_gspphot','A0':'distance_gspphot','HIP':'hip_original_ext_source_id','PS1':'ps1_original_ext_source_id','SDSS13':'sdss13_ext_source_id','SKYM2':'skym2_original_ext_source_id','TYC2':'tyc2_original_ext_source_id','URAT1':'urat1_original_ext_source_id','ALLWISE':'allwise_original_ext_source_id','APASS9':'apass9_original_ext_source_id','GSC23':'gsc23_original_ext_source_id','RAVE5':'gsc23_original_ext_source_id','_2MASS':'twomass_original_ext_source_id','RAVE6':'rave6_original_ext_source_id','RAJ2000':'ra2000','DEJ2000':'dec2000','DR3Name':'designation','SolID':'solution_id','RandomI':'random_index','RPlx':'parallax_over_error','RADEcor':'ra_dec_corr','RAPlxcor':'ra_parallax_corr','RApmRAcor':'ra_pmra_corr','RApmDEcor':'ra_pmdec_corr','DEPlxcor':'dec_parallax_corr','DEpmRAcor':'dec_pmra_corr','DEpmDEcor':'dec_pmdec_corr','PlxpmRAcor':'parallax_pmra_corr','PlxpmDEcor':'parallax_pmdec_corr','pmRApmDEcor':'pmra_pmdec_corr','NAL':'astrometric_n_obs_al','NAC':'astrometric_n_obs_ac','NgAL':'astrometric_n_good_obs_al','NbAL':'astrometric_n_bad_obs_al','gofAL':'astrometric_gof_al','chi2AL':'astrometric_chi2_al','epsi':'astrometric_chi2_al','sepsi':'astrometric_chi2_al','Solved':'astrometric_params_solved','APF':'astrometric_primary_flag','nueff':'nu_eff_used_in_astrometry','pscol':'pseudocolour','e_pscol':'pseudocolour_error','RApscolCorr':'ra_pseudocolour_corr','DEpscolCorr':'dec_pseudocolour_corr','PlxpscolCorr':'parallax_pseudocolour_corr','pmRApscolCorr':'pmra_pseudocolour_corr','pmDEpscolCorr':'pmdec_pseudocolour_corr','MatchObsA':'astrometric_matched_transits','Nper':'visibility_periods_used','amax':'astrometric_sigma5d_max','MatchObs':'matched_transits','NewMatchObs':'new_matched_transits','MatchObsrm':'matched_transits_removed','IPDgofha':'ipd_gof_harmonic_amplitude','IPDgofhp':'ipd_gof_harmonic_phase','IPDfmp':'ipd_frac_multi_peak','IPDfow':'ipd_frac_odd_win','SDSk1':'scan_direction_strength_k1','SDSk2':'scan_direction_strength_k2','SDSk3':'scan_direction_strength_k3','SDSk4':'scan_direction_strength_k4','SDMk1':'scan_direction_mean_k1','SDMk2':'scan_direction_mean_k2','SDMk3':'scan_direction_mean_k3','SDMk4':'scan_direction_mean_k4','Dup':'duplicated_source','o_Gmag':'phot_g_n_obs','RFG':'phot_g_mean_flux_over_error','e_Gmag':'phot_g_mean_mag_error','o_BPmag':'phot_bp_n_obs','RFBP':'phot_bp_mean_flux_over_error','e_BPmag':'phot_bp_mean_mag_error','o_RPmag':'phot_rp_n_obs','RFRP':'phot_rp_mean_flux_over_error','e_RPmag':'phot_rp_mean_mag_error','E_BP_RP_':'phot_bp_rp_excess_factor','NBPcont':'phot_bp_n_contaminated_transits','NBPblend':'phot_bp_n_blended_transits','NRPcont':'phot_rp_n_contaminated_transits','NRPblend':'phot_rp_n_blended_transits','Mode':'phot_proc_mode','BP-G':'bp_g','G-RP':'g_rp','n_RV':'rv_method_used','o_RV':'rv_nb_transits','o_RVd':'rv_nb_deblended_transits','RVNper':'rv_visibility_periods_used','RVS_N':'rv_expected_sig_to_noise','RVgof':'rv_renormalised_gof','RVchi2':'rv_chisq_pvalue','RVTdur':'rv_time_duration','RVamp':'rv_amplitude_robust','RVtempTeff':'rv_template_teff','RVtemplogg':'rv_template_logg','RVtemp_Fe_H_':'rv_template_fe_h','Vatmparam':'rv_atm_param_origin','e_Vbroad':'vbroad_error','o_Vbroad':'vbroad_nb_transits','e_GRVSmag':'grvs_mag_error','o_GRVSmag':'grvs_mag_nb_transits','RVSS_N':'rvs_spec_sig_to_noise','GLON':'l','GLAT':'b','ELON':'ecl_lon','ELAT':'ecl_lat','PQSO':'classprob_dsc_combmod_quasar','PGal':'classprob_dsc_combmod_galaxy','PSS':'classprob_dsc_combmod_star','b_Teff':'teff_gspphot_lower','B_Teff':'teff_gspphot_upper','b_logg':'logg_gspphot_lower','B_logg':'logg_gspphot_upper','b__Fe_H_':'mh_gspphot_lower','B__Fe_H_':'mh_gspphot_upper','b_Dist':'distance_gspphot_lower','B_Dist':'distance_gspphot_upper','b_A0':'azero_gspphot_lower','B_A0':'azero_gspphot_upper','AG':'ag_bspphot','b_AG':'ag_gspphot_lower','B_AG':'ag_gspphot_upper','E_BP-RP_':'ebpminrp_gspphot','b_E_BP-RP_':'ebpminrp_gspphot_lower','B_E_BP-RP_':'ebpminrp_gspphot_upper','Lib':'libname_gspphot','dHIP':'hip_angular_distance','nHIP':'hip_number_of_neighbours','f_HIP':'hip_xm_flag','PS1coid':'clean_panstarrs1_oid','dPS1':'ps1_angular_distance','nPS1':'ps1_number_of_neighbours','mPS1':'ps1_number_of_mates','f_PS1':'ps1_xm_flag','SDSS13coid':'clean_sdssdr13_oid','dSDSS13':'sdss13_angular_distance','nSDSS13':'sdss13_number_of_neighbours','mSDSS13':'sdss13_number_of_mates','f_SDSS13':'sdss13_xm_flag','dSKYM2':'skym2_angular_distance','nSKYM2':'skym2_number_of_neighbours','mSKYM2':'skym2_number_of_mates','f_SKYM2':'skym2_xm_flag','dTYC2':'tyc2_angular_distance','f_TYC2':'tyc2_xm_flag','TYC2moid':'tycho2tdsc_merge_oid','nTYC2':'tyc2_number_of_neighbours','dURAT1':'urat1_angular_distance','f_URAT1':'urat1_xm_flag','URAT1oid':'urat1_oid','nURAT1':'urat1_number_of_neighbours','mURAT1':'urat1_number_of_mates','dAllWISE':'allwise_angular_distance','f_AllWISE':'allwise_xm_flag','AllWISEoid':'allwise_oid','nAllWISE':'allwise_number_of_neighbours','mAllWISE':'allwise_number_of_mates','APASS9coid':'clean_apassdr9_oid','dAPASS9':'apass9_angular_distance','nAPASS9':'apass9_number_of_neighbours','mAPASS9':'apass9_number_of_mates','f_APASS9':'apass9_xm_flag','dGSC23':'gsc23_angular_distance','f_GSC23':'gsc23_xm_flag','GSC23coid':'clean_gsc23_oid','nGSC23':'gsc23_number_of_neighbours','mGSC23':'gsc23_number_of_mates','dRAVE5':'rave5_angular_distance','f_RAVE5':'rave5_xm_flag','RAVE5coid':'clean_ravedr5_oid','nRAVE5':'rave5_number_of_neighbours','d2MASS':'twomass_angular_distance','f_2MASS':'twomass_xm_flag','_2MASScoid':'clean_tmass_psc_xsc_oid','n2MASS':'twomass_number_of_neighbours','m2MASS':'twomass_number_of_mates','dRAVE6':'rave6_angular_distance','f_RAVE6':'rave6_xm_flag','RAVE6oid':'ravedr6_oid','nRAVE6':'rave6_number_of_neighbours','e_RAJ2000':'ra2000_error','e_DEJ2000':'dec2000_error','RADEcorJ2000':'ra2000_de2000_corr'}, inplace=True)
+	return data
+	
+#Returns all Gaia data from a given data release for an object, given its Gaia designation
+def GaiaQueryDesignation(Designation,catalogue='dr3'):
+	#Raises an exception if multiple designations are given, added since GaiaGetDesignation can return a tuple which then breaks if immediately passed back into this function without parsing first.
+	if isinstance(Designation,list):
+		raise Exception('multiple designations given')
+	#'GAIA DR...' prefix in object designations is not present in returned tables, but is required for Vizier query. Therefore need to create two separate parameters, one with the prefix (Search Designation) and one without (Designation)
+	if isinstance(Designation,str):
+		Designation=Designation.upper()
+	catalogue=catalogue.lower()
+	if catalogue!='dr2' and catalogue!='dr3':
+		raise Exception('invalid catalogue')
+		exit()
+
+	if 'GAIA DR3 ' in str(Designation) or 'GAIA DR2 ' in str(Designation):
+		SearchDesignation=Designation
+		if 'GAIA DR3' in str(Designation):
+			catalogue='dr3'
+		elif 'GAIA DR2' in str(Designation):
+			catalogue='dr2'
+		Designation=Designation[9:]
+	else:
+		if catalogue=='dr2':
+			SearchDesignation='GAIA DR2 '+str(Designation)
+		elif catalogue=='dr3':
+			SearchDesignation='GAIA DR3 '+str(Designation)
+
+	data=[]
+
+	#Sends Vizier query for given designation. Returns multiple objects within a radius, but as you are searching with a GAIA designation it is assumed that you already know the exact object that want and so all others are removed.
+	v=Vizier(columns=['**'],column_filters={"Source":"=="+str(Designation)},row_limit=row_limit)
+	
+	if catalogue=='dr2':
+		data.append(v.query_object(SearchDesignation,catalog=catalogue_names[0]))
+		if len(data[0])==0:
+			print('[Gaia: GaiaQueryDesignation] Error: no Gaia data found for given designation')
+			return None
+		data=data[0][0].to_pandas()
+		data=renameHeadersDR3(data)
+		return data
+	elif catalogue=='dr3':
+		data.append(v.query_object(SearchDesignation,catalog=catalogue_names[1],))
+		if len(data[0])==0:
+			print('[Gaia: GaiaQueryDesignation] Error: no Gaia data found for given designation')
+			return None
+		data=data[0][0].to_pandas()
+		data=renameHeadersDR3(data)
+		return data
+
+#Returns all Gaia data from a given data release for an object, given its coordinates
+def GaiaQueryCoords(ObjRa,ObjDec,radius=3,catalogue='dr3'):
+	catalogue=catalogue.lower()
+	
+	if not isinstance(radius,int) and not isinstance(radius,float):
+		raise Exception('search radius must be a float or an integer')
+	
+	if not isinstance(ObjRa,float):
+		try:
+			ObjRa=float(ObjRa)
+		except:
+			raise Exception('object RA must be a float or an integer')
+	if not isinstance(ObjDec,float):
+		try:
+			ObjDec=float(ObjDec)
+		except:
+			raise Exception('object DEC must be a float or an integer')
+	
+	if catalogue!='dr2' and catalogue!='dr3':
+		raise Exception('invalid catalogue')
+	
+	data=[]
+	
+	#Get all columns
+	v=Vizier(columns=['**'],row_limit=row_limit)
+	
+	#Sends Vizier query, can return multiple objects in given search radius
+	if catalogue=='dr2':
+		data.append(v.query_region(coord.SkyCoord(ra=ObjRa,dec=ObjDec,unit=(u.deg,u.deg),frame='icrs'),width=radius*u.arcsec,catalog=catalogue_names[0]))
+		if len(data[0])==0:
+			print('[Gaia: GaiaQueryCoords] Error: no Gaia data found at given coordinates')
+			return None
+		data=data[0][0].to_pandas()
+		data=data.reset_index(drop=True)
+		data=renameHeadersDR3(data)
+		return data
+	if catalogue=='dr3':
+		data.append(v.query_region(coord.SkyCoord(ra=ObjRa,dec=ObjDec,unit=(u.deg,u.deg),frame='icrs'),width=radius*u.arcsec,catalog=catalogue_names[1]))
+		if len(data[0])==0:
+			print('[Gaia: GaiaQueryCoords] Error: no Gaia data found at given coordinates')
+			return None
+		data=data[0][0].to_pandas()
+		data=data.reset_index(drop=True)
+		data=renameHeadersDR3(data)
+		return data
+
+#Returns an objects Gaia designation given its coordinates. Designations are given the prefix 'GAIA DR2' or 'GAIA DR3' depending on the catalogue given, as this is the format used by Vizier.
+def GaiaGetDesignation(ObjRa,ObjDec,radius=3,catalogue='dr3'):
+	data=GaiaQueryCoords(ObjRa,ObjDec,radius,catalogue)
+	
+	catalogue=catalogue.lower()
+	
+	Designation=data['source_id']
+	#Returns designation if only one object is found
+	if len(data)==1:
+		if catalogue=='dr2':
+			Designation='GAIA DR2 '+str(Designation[0])
+		elif catalogue=='dr3':
+			Designation='GAIA DR3 '+str(Designation[0])
+	#Returns designation as a tuple if multiple objects are found
+	if len(data)>1:
+		print('[Gaia: GaiaGetDesignation] Note: '+len(data), 'objects detected within search radius.')
+		Designation=Designation.to_list()
+	if not isinstance(data,pd.DataFrame):
+		return None
+	return Designation
+
+#Returns an object's coordinates in a given data release given a Gaia designation, returns them as tuple [RA,DEC]
+def GaiaGetCoords(Designation,catalogue='dr3'):
+	data=GaiaQueryDesignation(Designation,catalogue)
+	if not isinstance(data,pd.DataFrame) or data.empty:
+		return None
+	
+	#If used correctly this should never trigger as each object should have a unique Gaia designation
+	if len(data)>1:
+		raise Exception('too many objects returned by search')
+	
+	objRa=data['ra']
+	objDec=data['dec']
+	Coords=[objRa[0],objDec[0]]
+	
+	return Coords
+	
+def GaiaGetPhotometryDesignation(Designation,catalogue='dr3'):
+	data=GaiaQueryDesignation(Designation,catalogue)
+	if not isinstance(data,pd.DataFrame) or data.empty:
+		return None
+	
+	photometry=data[['ra','dec','source_id','phot_g_mean_mag','phot_g_mean_mag_error','phot_bp_mean_mag','phot_bp_mean_mag_error','phot_rp_mean_mag','phot_rp_mean_mag_error']].copy()
+	
+	return photometry
+	
+def GaiaGetPhotometryCoords(ObjRa,ObjDec,radius=3):
+	data=GaiaQueryCoords(ObjRa,ObjDec,radius=radius)
+	
+	if not isinstance(data,pd.DataFrame) or data.empty:
+		return None
+	
+	photometry=data[['ra','dec','source_id','phot_g_mean_mag','phot_g_mean_mag_error','phot_bp_mean_mag','phot_bp_mean_mag_error','phot_rp_mean_mag','phot_rp_mean_mag_error']].copy()
+
+	return photometry
+
+if __name__=='__main__':
+	print(GaiaQueryDesignation(6050296829033196032))
+	
+	

AstroToolkit/Surveys/PanSTARRS.py

@@ -0,0 +1,221 @@
+from astropy.table import Table
+from astropy.io import fits
+import astropy
+from astropy.visualization import PercentileInterval, AsinhStretch
+from astropy.wcs import WCS
+import PIL
+from PIL import Image
+from io import BytesIO
+import numpy as np
+import requests
+import re
+import warnings
+from astropy.utils.exceptions import AstropyWarning
+import pandas as pd
+from bokeh.plotting import figure
+from bokeh.models import Range1d
+
+#Ignore warning about deprecated header
+warnings.simplefilter('ignore', category=AstropyWarning)
+
+def get_image(ra,dec,size=30,band='g'):
+	# Construct URL for fetching a list of image urls
+	service='https://ps1images.stsci.edu/cgi-bin/ps1filenames.py'
+	url=f'{service}?ra={ra}&dec={dec}&band={band}'
+	
+	# Try to read the resulting table
+	try:
+		table=Table.read(url, format='ascii')
+	except:
+		print('[PanSTARRS: get_image_data] Error: experiencing issues with Pan-STARRS.')
+		return None, None, None
+	
+	if not len(table)>0:
+		print('[PanSTARRS: get_image_data] Error: no PanSTARRS images found using given fields.')
+		return None, None, None
+	
+	# Finds any missing filters in returned data and removes them from the 'band' string to signify this
+	shortnames=[]
+	for i in range(0,len(table)):
+		shortnames.append(table['shortname'][i])
+	
+	missing_filters=[]
+	if len(table)!=len(band):
+		for i in range(0,len(band)):
+			filterCheck='.'+band[i]+'.'
+			if any(filterCheck in shortname for shortname in shortnames):
+				pass
+			else:
+				print('[PanSTARRS: get_image_data] Note:  no PanSTARRS data in '+band[i]+' band')
+				missing_filters.append(band[i])
+	
+	for i in range(0,len(missing_filters)):
+		band=band.replace(missing_filters[i],'')
+
+	# Constructs the URL used to fetch the images
+	url=(f'https://ps1images.stsci.edu/cgi-bin/fitscut.cgi?ra={ra}&dec={dec}&size={size}&format=fits')
+	
+	filterList=['yzirg'.find(x) for x in table['filter']]
+	table=table[np.argsort(filterList)]
+	
+	urlbase=url+'&red='
+	url=[]
+	for filename in table['filename']:
+		url.append(urlbase+filename)
+
+	# Grab images using constructed URL
+	try:
+		r=requests.get(url[0],timeout=15)
+		if r.status_code!=200:
+			print('[PanSTARRS: getFitsImage] Error:  experiencing issues with Pan-STARRS.')
+			return None, None, None
+	except:
+		print('[PanSTARRS: get_image_data] Error: experiencing issues with Pan-STARRS.')
+		return None, None, None
+	
+	# Read and format image data
+	try:
+		fh=fits.open(BytesIO(r.content))
+		fitsImage=fh[0].data
+		image_header=fh[0].header
+		fitsImage[np.isnan(fitsImage)]=0.0
+		transform=AsinhStretch()+PercentileInterval(95)
+		image=transform(fitsImage)
+
+		return image, image_header, band
+	except:
+		print('[PanSTARRS: get_image_data] Error: no PanSTARRS images found using given fields. 2')
+		return None, None, None
+
+'''
+Function for returning the wcs (transformation) and mjd of an image
+'''
+def get_info(ra,dec,size=30,band='g'):
+	# Verify inputs
+	ra,dec,size,sizeAS,band=input_management(ra=ra,dec=dec,size=size,band=band)
+	# Get image header
+	image_data=get_image(ra=ra,dec=dec,size=size,band=band)[1]
+	
+	if not isinstance(image_data,astropy.io.fits.header.Header):
+		return None, None
+	
+	wcs=WCS(image_data)
+	mjd=image_data['MJD-OBS']
+	
+	return wcs, mjd
+
+'''
+Verifies all inputs and transforms them into the units required by queries
+'''
+def input_management(ra,dec,size=30,band='g'):
+	if len(band)!=len(set(band)):
+		raise Exception('duplicate detected in filter selection')
+	
+	if not isinstance(ra,float):
+		try:
+			ra=float(ra)
+		except:
+			raise Exception('object RA must be a float or an integer')
+	if not isinstance(dec,float):
+		try:
+			dec=float(dec)
+		except:
+			raise Exception('object DEC must be a float or an integer')
+		
+	band=band.lower()
+	
+	if not isinstance(size,int):
+		try:
+			size=int(size)
+		except:
+			raise Exception('size must take an integer value')
+	
+	# Adjusts size to unit required by PanSTARRS (0.25"/pixel)
+	sizeAS=size
+	size=size*4
+	
+	if size>6000:
+		raise Exception('maximum size is 25 arcmin')
+	
+	if re.match('^[grizy]+$', band):
+			pass
+	else:
+		raise Exception('invalid bands, supported bands: [g,r,i,z,y]')
+	
+	return ra,dec,size,sizeAS,band
+
+def get_plot(ra,dec,size=30,band='g'):
+	# Verify and transform inputs
+	ra,dec,size,sizeAS,band=input_management(ra=ra,dec=dec,size=size,band=band)
+
+	# Create axes and plot image
+	plot=figure(width=400,height=400,title=f'PanSTARRS Image ({sizeAS}")')
+	# Fixes scaling issue
+	plot.min_border=75
+
+	plot.xgrid.grid_line_color=None
+	plot.ygrid.grid_line_color=None	
+
+	image,image_data,_=get_image(ra=ra,dec=dec,size=size,band=band)
+	
+	if not isinstance(image,np.ndarray) or not isinstance(image_data,astropy.io.fits.header.Header):
+		return None
+
+	# Get x and y limit (in pixels) from image header
+	xlim=image_data['NAXIS1']
+	ylim=image_data['NAXIS2']
+	
+	# Get points on axes in world coordinates (deg)
+	wcs=WCS(image_data)
+	
+	# Get array of pixels in x and y
+	x_points=np.arange(start=0,stop=xlim+1,step=1)
+	y_points=np.arange(start=0,stop=ylim+1,step=1)
+
+	# Apply transformation to coordinates
+	coords=wcs.all_pix2world(x_points,y_points,1)
+	x_points,y_points=coords[0],coords[1]
+	
+	# Get image size in degrees
+	x_range=max(x_points)-min(x_points)
+	y_range=max(y_points)-min(y_points)
+
+	plot.x_range=Range1d(max(x_points),min(x_points))
+	plot.y_range=Range1d(min(y_points),max(y_points))
+
+	plot.image(image=[image],x=x_points[0],y=y_points[0],dw=x_range,dh=y_range,palette='Greys256',level="image",origin='bottom_right',anchor='bottom_right')
+	return plot
+
+def PanSTARRSQueryCoords(ObjRa,ObjDec,radius=3):
+	radius=radius/3600
+
+	#Had to use JSON format to get this to work, the CSV format used is not easily formattable by pandas
+	url=f'https://gsss.stsci.edu/webservices/vo/CatalogSearch.aspx?CAT=PS1dr2OBJECTS&RA={ObjRa}&DEC={ObjDec}&SR={radius}&FORMAT=json'
+	
+	try:
+		r=requests.get(url)
+	except:
+		print('[PanSTARRS: PanSTARRSQueryCoords] Error: experiencing issues with Pan-STARRS')
+		return None
+	
+	if r.status_code!=200:
+		print('[PanSTARRS: PanSTARRSQueryCoords] Error: experiencing issues with Pan-STARRS')
+		return None
+
+	data=BytesIO(r.content)
+	try:
+		data=pd.read_json(data)
+	except:
+		print('[PanSTARRS: PanSTARRSQueryCoords] Note: No data found using given coordinates')
+	
+	return data
+	
+def PanSTARRSGetPhotometryCoords(ObjRa,ObjDec,radius=3):
+	data=PanSTARRSQueryCoords(ObjRa,ObjDec,radius=radius)
+	
+	if not isinstance(data,pd.DataFrame) or data.empty:
+		print('[PanSTARRS: PanSTARRSGetPhotometryCoords] Note: No PanSTARRS data at given coordinates.')
+		return None
+	
+	photometry=data[['ra','dec','objID','gMeanPSFMag','gMeanPSFMagErr','rMeanPSFMag','rMeanPSFMagErr','iMeanPSFMag','iMeanPSFMagErr','zMeanPSFMag','zMeanPSFMagErr','yMeanPSFMag','yMeanPSFMagErr']].copy()
+	return photometry

AstroToolkit/Surveys/ROSAT.py

@@ -0,0 +1,51 @@
+import pandas as pd
+import astropy.coordinates as coord
+from astroquery.vizier import Vizier
+import astropy.units as u
+import math
+import warnings
+from astropy.utils.exceptions import AstropyWarning
+
+catalogue_name='J/A+A/588/A103'
+
+Vizier.ROW_LIMIT = -1
+
+warnings.simplefilter('ignore', category=AstropyWarning)
+
+def ROSATQueryCoords(ObjRa,ObjDec,SearchRadius=3):
+	data=[]
+	
+	#Filter inputs
+	if not isinstance(SearchRadius,int) and not isinstance(SearchRadius,float):
+		raise Exception('search radius must be a float or an integer')
+	
+	if not isinstance(ObjRa,float):
+		try:
+			ObjRa=float(ObjRa)
+		except:
+			raise Exception('object RA must be a float or an integer')
+	if not isinstance(ObjDec,float):
+		try:
+			ObjDec=float(ObjDec)
+		except:
+			raise Exception('object DEC must be a float or an integer')
+	
+	v=Vizier(columns=['**'])
+	
+	#Sends Vizier query, can return multiple objects in given search radius
+	data.append(v.query_region(coord.SkyCoord(ra=ObjRa,dec=ObjDec,unit=(u.deg,u.deg),frame='icrs'),width=SearchRadius*u.arcsec,catalog=catalogue_name))
+	if len(data[0])==0:
+		print('[ROSAT: ROSATQueryCoords] Error: no ROSAT data found at given coordinates')
+		return None
+	data=data[0][0].to_pandas()
+	data=data.reset_index(drop=True)
+	return data
+	
+def ROSATGetPhotometryCoords(ObjRa,ObjDec,SearchRadius=3):
+	data=ROSATQueryCoords(ObjRa,ObjDec,SearchRadius=SearchRadius)
+	if not isinstance(data,pd.DataFrame) or data.empty:
+		return None
+	
+	photometry=data[['RAJ2000','DEJ2000','VmagVV10','VTmag','BTmag','VmagBSC','VmagHMXB','VmagLMXB','VmagWD','VsphotWD']].copy()
+
+	return photometry

AstroToolkit/Surveys/SDSS.py

@@ -0,0 +1,78 @@
+from astroquery.sdss import SDSS
+from astropy import coordinates as coords
+from astropy import units as u
+from astroquery.vizier import Vizier
+import astropy.coordinates as coord
+import pandas as pd
+from bokeh.plotting import figure
+import astropy
+
+catalogue_name='V/154/sdss16'
+
+Vizier.ROW_LIMIT = -1
+
+def get_spectrum_data(ra,dec,radius):
+	pos=coords.SkyCoord(ra,dec,unit='deg')
+	radius=radius/3600*u.deg
+
+	data=SDSS.get_spectra(pos,radius)
+	if data!=None:
+		data=data[0][1]
+	else:
+		return None, None
+	
+	spectrum_data=data.data
+	
+	log_wavelength=spectrum_data['loglam']
+	wavelength=10**log_wavelength*u.AA
+	flux=spectrum_data['flux']*10**-17*u.Unit('erg cm-2 s-1 AA-1')
+
+	return wavelength,flux
+
+def get_plot(ra,dec,radius=3):
+	x,y=get_spectrum_data(ra,dec,radius)
+	
+	if not isinstance(x,astropy.units.quantity.Quantity) or not isinstance(y,astropy.units.quantity.Quantity):
+		return None
+
+	plot=figure(width=400,height=400,title="SDSS Spectrum",x_axis_label=r'\[\lambda\text{ }[\text{AA}]\]',y_axis_label=r"\[\text{flux [erg}\text{ cm }^{-2}\text{ s }^{-1}\text{AA}^{-1}]\]")
+	plot.line(x,y,color='black',line_width=1)
+	
+	return plot
+
+def get_data(ra,dec,radius=3):
+	data=[]
+	
+	if not isinstance(radius,int) and not isinstance(radius,float):
+		raise Exception('search radius must be a float or an integer')
+	
+	if not isinstance(ra,float):
+		try:
+			ra=float(ra)
+		except:
+			raise Exception('object RA must be a float or an integer')
+	if not isinstance(dec,float):
+		try:
+			dec=float(dec)
+		except:
+			raise Exception('object DEC must be a float or an integer')
+	
+	v=Vizier(columns=['**'])
+
+	#Sends Vizier query, can return multiple objects in given search radius
+	data.append(v.query_region(coord.SkyCoord(ra=ra,dec=dec,unit=(u.deg,u.deg),frame='icrs'),width=radius*u.arcsec,catalog=catalogue_name))
+	if len(data[0])==0:
+		print('[SDSS: get_data] Error: no SDSS data found at given coordinates')
+		return None
+	data=data[0][0].to_pandas()
+	data=data.reset_index(drop=True)
+	return data
+	
+def get_photometry(ra,dec,radius=3):
+	SDSSdata=get_data(ra,dec,radius=radius)
+	if not isinstance(SDSSdata,pd.DataFrame) or SDSSdata.empty:
+		return None
+	
+	photometry=SDSSdata[['RA_ICRS','DE_ICRS','objID','uPmag','e_uPmag','gPmag','e_gPmag','rPmag','e_rPmag','iPmag','e_iPmag','zPmag','e_zPmag']].copy()
+
+	return photometry