6.9E-7, 6.789E-7, 6.682E-7, 6.58E-7, 6.482E-7, 6.388E-7, 6.298E-7, 6.212E-7, 6.129E-7, 6.049E-7, 5.972E-7, 5.897E-7, 5.826E-7, 5.757E-7, 5.69E-7, 5.626E-7, 5.564E-7, 5.503E-7, 5.445E-7, 5.389E-7, 5.334E-7, 5.281E-7, 5.23E-7, 5.18E-7, 5.132E-7, 5.085E-7, 5.039E-7, 4.995E-7, 4.952E-7, 4.91E-7, 4.869E-7, 4.829E-7, 4.791E-7, 4.753E-7, 4.717E-7, 4.681E-7, 4.646E-7, 4.612E-7, 4.579E-7, 4.547E-7, 4.515E-7, 4.484E-7, 4.454E-7, 4.425E-7, 4.396E-7, 4.368E-7, 4.34E-7, 4.314E-7, 4.287E-7, 4.262E-7, 4.236E-7, 4.212E-7, 4.188E-7, 4.164E-7, 4.141E-7, 4.118E-7, 4.096E-7, 4.074E-7, 4.053E-7, 4.032E-7, 4.011E-7, 3.991E-7, 3.972E-7, 3.952E-7, 3.933E-7, 3.915E-7, 3.896E-7, 3.878E-7, 3.86E-7, 3.843E-7, 3.826E-7, 3.809E-7, 3.793E-7, 3.777E-7, 3.761E-7, 3.745E-7, 3.73E-7, 3.714E-7, 3.7E-7, 3.685E-7, 3.671E-7, 3.656E-7, 3.642E-7, 3.629E-7, 3.615E-7, 3.602E-7, 3.589E-7, 3.576E-7, 3.563E-7, 3.551E-7, 3.538E-7, 3.526E-7, 3.514E-7, 3.502E-7, 3.491E-7, 3.479E-7, 3.468E-7, 3.457E-7, 3.446E-7, 3.435E-7, 3.424E-7, 3.414E-7, 3.403E-7, 3.393E-7, 3.383E-7, 3.373E-7, 3.363E-7, 3.354E-7, 3.344E-7, 3.335E-7, 3.325E-7, 3.316E-7, 3.307E-7, 3.298E-7, 3.289E-7, 3.28E-7, 3.272E-7, 3.263E-7, 3.255E-7, 3.247E-7, 3.238E-7, 3.23E-7, 3.222E-7, 3.214E-7, 3.206E-7, 3.199E-7, 3.191E-7, 3.184E-7, 3.176E-7, 3.169E-7, 3.161E-7, 3.154E-7, 3.147E-7, 3.14E-7, 3.133E-7, 3.126E-7 Digitized First Byurakan Survey (DFBS) Extracted Spectra The First Byurakan Survey (FBS) is the largest and the first systematic objective prism survey of the extragalactic sky. It covers 17,000 sq.deg. in the Northern sky together with a high galactic latitudes region in the Southern sky. The FBS has been carried out by B.E. Markarian, V.A. Lipovetski and J.A. Stepanian in 1965-1980 with the Byurakan Observatory 102/132/213 cm (40"/52"/84") Schmidt telescope using 1.5 deg. prism. Each FBS plate contains low-dispersion spectra of some 15,000-20,000 objects; the whole survey consists of about 20,000,000 objects. 30 and sp_class='OK' Since the table needs to be sequentially scanned for this, it will take a minute or so. Combine with an object selection (see below) or other criteria if possible. Get Average Spectra ------------------- You cannot currently use the ADQL aggregate function AVG with arrays (which should be fixed at some time in the future). Meanwhile, you can work around this with a clumsy construction like this (this query will give you average spectra by magnitude bin; don't run it just for fun, it'll take a while):: select round(magb) as bin, avg(flux[1]) as col1, avg(flux[2]) as col2, avg(flux[3]) as col3, avg(flux[4]) as col4, avg(flux[5]) as col5, avg(flux[6]) as col6, avg(flux[7]) as col7, avg(flux[8]) as col8, avg(flux[9]) as col9, avg(flux[10]) as col10, avg(flux[11]) as col11, avg(flux[12]) as col12, avg(flux[13]) as col13, avg(flux[14]) as col14, avg(flux[15]) as col15, avg(flux[16]) as col16, avg(flux[17]) as col17, avg(flux[18]) as col18, avg(flux[19]) as col19, avg(flux[20]) as col20, avg(flux[21]) as col21, avg(flux[22]) as col22, avg(flux[23]) as col23, avg(flux[24]) as col24, avg(flux[25]) as col25, avg(flux[26]) as col26, avg(flux[27]) as col27, avg(flux[28]) as col28, avg(flux[29]) as col29, avg(flux[30]) as col30, avg(flux[31]) as col31, avg(flux[32]) as col32, avg(flux[33]) as col33, avg(flux[34]) as col34, avg(flux[35]) as col35, avg(flux[36]) as col36, avg(flux[37]) as col37, avg(flux[38]) as col38, avg(flux[39]) as col39, avg(flux[40]) as col40, avg(flux[41]) as col41, avg(flux[42]) as col42, avg(flux[43]) as col43, avg(flux[44]) as col44, avg(flux[45]) as col45, avg(flux[46]) as col46, avg(flux[47]) as col47, avg(flux[48]) as col48, avg(flux[49]) as col49, avg(flux[50]) as col50, avg(flux[51]) as col51, avg(flux[52]) as col52, avg(flux[53]) as col53, avg(flux[54]) as col54, avg(flux[55]) as col55, avg(flux[56]) as col56, avg(flux[57]) as col57, avg(flux[58]) as col58, avg(flux[59]) as col59 from dfbsspec.spectra where sp_class='OK' group by bin To map ``col`` to wavelenghts, see the contents of (any) ``spectral`` column. Build Templates --------------- To compute an average spectrum for a class of objects, we suggest to pull positions of such objects from SIMBAD and then fetch the associate spectra from this database. Since the response function of the photographic plates had a strong magnitude dependence, restrict the objects to a small magnitude range, for instance:: select otype, ra, dec, flux from basic join flux on (oid=oidref) where otype='HS*' and dec>-15 and filter='G' and flux between 12.5 and 13.5 (to be executed on SIMBAD's TAP service, see also `SIMBAD object types`_). .. _SIMBAD object types: http://simbad.u-strasbg.fr/simbad/sim-display?data=otypes With the resulting table, go do this service and execute a query like:: SELECT specid, spectral, flux FROM dfbsspec.spectra AS db JOIN TAP_UPLOAD.t1 AS tc ON DISTANCE(tc.ra, tc.dec, db.ra, db.dec)<5./3600. WHERE sp_class='OK' (adjust t1 according to your client's rules; in TOPCAT, that's t plus the table number from the control window). History of this Data Collection =============================== The original aim of the First Byurakan Survey was to search for galaxies with UV excess (:bibcode:`1986ApJS...62..751M`, Markarian et al. 1989,1997- catalogue No. VII/172 at CDS). Successively, the amount of spectral information contained in the plates allowed the development of several other projects concerning the spectral classification of Seyfert Galaxies (Weedman and Kachikian 1971), the first definition of starburst galaxies (Weedman 1977 ), the discovery and investigation of blue stellar objects (Abrahamian and Mickaelian, 1996, Mickaelian et al 2001, 2002, CDS catalogue No II/223) and a survey for late-type stars (Gigoyan et al. 2002). All these results were obtained by eye inspection of the plates performed with the aid of a microscope at the Byurakan Observatory. The number and classes of new objects discovered FBS made clear the need of open access to FBS for the entire astronomical community. ]]> 2017-11-28T10:50:00Z 100 spectroscopy surveys active-galactic-nuclei Markarian, B.E. Byurakan 1m Schmidt Byurakan Astrophysical Observatory BAO 2007A&A...464.1177M Research Archive Optical //products#table//scs#q3cindex//scs#pgs-pos-index Raw metadata for the spectra, to be combined with image metadata like date_obs and friends for a complete spectrum descriptions. This also contains spectral and flux points in array-valued columns. specid
make_ssa_view make_tap_view # The spectra are in one table per source plate, and # the dump has the values of one such plate in one line. # So, I look for such lines and then parse from there with gzip.open(self.sourceToken) as f: for line in f: if line.startswith(b"INSERT INTO `fbs"): for row in parseSpectraFromDump(line.decode("ascii")): yield row "\schema.spectra" "\rdId/%s-%s"%(@plate, @objectid[5:]) makeAbsoluteURL( "\rdId/sdl/dlget?ID="+urllib.parse.quote("\rdId/%s-%s"%( @plate, @objectid[5:]))) "application/x-votable+xml" "image/png" makeAbsoluteURL("\rdId/preview/qp/%s-%s"%( @plate, @objectid[5:])) 20000 [float(item or 'nan') for item in @fluxes.split("#")[int(@redBord)+1:]] float(@ra)/24.*360 float(@dec) @plate+"-"+@objectid[5:] [\spectralBins] try: @lam_min = bins[len(@flux)-1] except IndexError: @lam_min = bins[-1] pgsphere.SPoint.fromDegrees(@ra, @dec) 'ivo://\getConfig{ivoa}{authority}/~?\rdId/'+@specid len(@flux) A view providing standard SSA metadata for DBFS metadata in \schema.spectrapresent//ssap#mixc//ssap#simpleCoverage//obscore#publishSSAPMIXCsdlssa_pubDID
39165 44374 2.879e-19 6.355e-19 Metadata for the plates making up the Byurakan spectral surveys, obtained from the WFPDB.
make_ssa_view make_tap_view from urllib.parse import urlencode from gavo import votable from gavo.stc import jYearToDateTime, dateTimeToMJD f = utils.urlopenRemote(self.sourceToken, data=urlencode({ "LANG": "ADQL", "QUERY": "select object as plateid, epoch, exptime, emulsion," " raj2000, dej2000" " from wfpdb.main" " where object is not null and object!=''" " and instr_id='BYU102A'" " and method='objective prism'"})) data, metadata = votable.load(f) for row in metadata.iterDicts(data): # HACK: duplicates in WFPDB. See README if row["plateid"]=="NPS": continue if row["plateid"]=="FBS 0966" and int(row["epoch"])==1974: row["plateid"] = "FBS 0966a" if row["plateid"]=="FBS 0326" and row["epoch"]>1971.05: row["plateid"] = "FBS 0326a" if row["plateid"]=="FBS 0449" and row["epoch"]>1971.38: row["plateid"] = "FBS 0449a" row["epoch"] = dateTimeToMJD(jYearToDateTime(row["epoch"])) row["plateid"] = row["plateid"].replace("FBS ", "fbs") row["ra_center"] = row["raj2000"] row["dec_center"] = row["dej2000"] yield row @epoch-@exptime/43200 @epoch+@exptime/43200 This table contains basic metadata as well as the spectra from the Digital First Byurakan Survey (DFBS). 1 CREATE VIEW \curtable AS ( SELECT \colNames FROM ( SELECT a.*, b.*, ('{\spectralBins}'::float[])[1:array_length(a.flux, 1)] AS spectral, '\getConfig{web}{serverURL}/dfbsspec/q/cutout/qp/' || specid as cutout_link FROM \schema.raw_spectra as a JOIN \schema.platemeta as b ON (plate=plateid)) AS t)
//ssap#sdm-instance
with base.getTableConn() as conn: res = list(conn.query( "select spectral, flux from dfbsspec.spectra" " where accref=%(accref)s", {"accref": self.sourceToken["accref"]}))[0] for lam, flx in zip(*res): yield {"spectral": lam, "flux": flx} DFBS Datalink Service data "\rdId#ssa" "\rdId#build_sdm_data" if descriptor.pubDID is None: return specId = descriptor.pubDID.split("/")[-1] yield descriptor.makeLink( '\getConfig{web}{serverURL}/dfbsspec/q/cutout/qp/'+specId, contentType="application/fits", contentLength=200*200*16, description="Cutout of the objective prism spectrum this" " was generated from.") Digitized First Byurakan Survey Browser Interface DFBS SSAP Digitized First Byurakan Survey (DFBS) Spectra Query Service pointed MAXREC=1 archival query DFBS spectra preview maker specid from gavo.helpers.processing import SpectralPreviewMaker from gavo import svcs with base.getTableConn() as conn: res = list(conn.query("SELECT spectral, flux" " FROM \schema.spectra" " WHERE specid=%(specid)s", inputTable.args)) if not res: raise svcs.UnknownURI("No such spectrum known here") return ("image/png", SpectralPreviewMaker.get2DPlot( list(zip(res[0][0], res[0][1])), linear=True)) DFBS spectra plate cutout service specid from gavo import svcs with base.getTableConn() as conn: res = list(conn.query( "SELECT plate, pos, px_length, pos_ang" " FROM \schema.raw_spectra" " WHERE specid=%(specid)s", inputTable.args)) if not res: raise svcs.UnknownURI("No such spectrum known here") plate, pos, px_length, pos_ang = res[0] x, y = pos.x/utils.DEG, pos.y/utils.DEG, pos_ang = pos_ang*utils.DEG w = 0.01 # degrees h = 0.05 # degrees h0 = 0.01 # degrees raise svcs.Found( "dfbs/q/dl/dlget?ID=plate/%s" "&POLYGON=%s%%20%s%%20%s%%20%s%%20%s%%20%s%%20%s%%20%s"%( plate, x-h0*math.cos(pos_ang)-w*math.sin(pos_ang), y+h0*math.sin(pos_ang)-w*math.cos(pos_ang), x-h0*math.cos(pos_ang)+w*math.sin(pos_ang), y+h0*math.sin(pos_ang)+w*math.cos(pos_ang), x+h*math.cos(pos_ang)-w*math.sin(pos_ang), y-h*math.sin(pos_ang)-w*math.cos(pos_ang), x+h*math.cos(pos_ang)+w*math.sin(pos_ang), y-h*math.sin(pos_ang)+w*math.cos(pos_ang))) getssa/ssap.xml row = self.getFirstVOTableRow() self.assertEqual(row['ssa_location'].asSODA(), '211.0680249982 44.899394444') self.assertEqual(row["ssa_dstitle"], 'DFBSJ140416.32+445357.8 spectrum from fbs0017') self.assertEqual(row["ssa_bandpass"], "IIAF bkd") self.assertEqual(row["ssa_length"], 112) self.assertAlmostEqual(row["ssa_specstart"], 3.316E-7) self.assertAlmostEqual(row["ssa_dateObs"], 40407, places=5) /getproduct/dfbsspec/q/fbs0017-141024.11+445238.4?preview=True self.assertHasStrings("PNG", "IDATx") sdl/dlget self.assertHasStrings('utype="spec:Spectrum"', 'name="citation"', 'value="UNCALIBRATED"') rows = self.getVOTableRows() self.assertAlmostEqual(rows[-1]["spectral"], 3.316e-07) self.assertAlmostEqual(rows[-1]["flux"], -5.02) sdl/dlmeta for row in self.getVOTableRows(): if row["semantics"]=="#progenitor": break else: self.assertFalse("No #progenitor link found.") self.assertEqual(row["content_type"], "application/fits") self.assertTrue(row["access_url"].endswith( "dfbsspec/q/cutout/qp/fbs0017-140416.32+445357.8")) /tap/sync rows = self.getVOTableRows() self.assertEqual(len(rows), 1) row = rows[0] self.assertAlmostEqual(row["lam_min"], 3.2060000e-07) self.assertAlmostEqual(row["dec"], 45.606705555556) self.assertEqual(row["emulsion"], "IIAF bkd") self.assertAlmostEqual(row["flux"][2], 16.26, 5) self.assertAlmostEqual(row["spectral"][2], 6.682e-7) self.assertEqual(len(row["spectral"]), len(row["flux"])) cutout/qp/fbs0017-140416.32+445357.8 self.assertHasStrings("DATEORIG= '1969-07-05'", "NAXIS2 = 70", b"\\x06\\xd5\\x06\\x49\\x06\\x49") /tap/sync row = self.getFirstVOTableRow(rejectExtras=False) self.assertEqual(len(row["ss"]), 130) allNULL = row["ss"][-6:] for val in allNULL: self.assertTrue(val is None)