import pyfits import os,sys import time #source pour Python 2.7 # recuperer Python 2.7 ici: http://www.python.org/download/releases/2.7.6/ # attention, il faut garder la coherence entre python et les package. C'est 32 bit ou 63 bit au choix. # aller ici pour recuperer pyfits: http://www.lfd.uci.edu/~gohlke/pythonlibs/#pyfits # return [min,max] entre les longeur d onde xmin,xmax def foundMinMax(data,start,step,xmin,xmax): ixmin=int((xmin-start)/step) ixmax=int((xmax-start)/step) n=len(data) if (ixmin<0): ixmin=0 if (ixmax<0): ixmax=0 if (ixmin>=n): ixmin=0 if (ixmax>=n): ixmax=n return [min(data[ixmin:ixmax]),max(data[ixmin:ixmax])] def WriteDatFile(data,start,step): f=open(PathFileTmpDat,'w') wave=start for intensity in data: l=str(wave)+" "+str(intensity)+"\n" f.write(l) wave+=step f.close() def WritePngFile(DstFile,title,xmin,xmax,ymin,ymax): print("DstFile: "+DstFile) print("Title: "+title) f=open(PathGnuplot+'gnuplotcmd.txt','w') l1="call 'C:\\\\gnuplot\\\\std.gnu' '"+PathFileTmpDat+"' '"+title+"' '" l1+=str(xmin)+":"+str(xmax)+"' '"+str(ymin)+":"+str(ymax)+"' '800,400' '"+DstFile+"'\n" l2="pause -1" f.write(l1) # f.write(l2) f.close() cmd=PathGnuplot+"binary\\pgnuplot.exe C:\gnuplot\gnuplotcmd.txt" print ("CMD = "+cmd) sortie=os.system(cmd) print(sortie) # le programme commence ici #le repertoire contenant les fichiers .fits au format bess est dans le chemin courant et s apelle web-xxxxx # verifie le nombre d'argument if len(sys.argv)==2: filename=sys.argv[1] print "nom du repertoire: "+filename else: print "nombre d'argument incorrect" print "utiliser un arguments 'nom de repertoire'" exit(1) StatusDic={\ 'BD+64 106':['Already publish
not in Bess','emission line signaled in Catalogue of Stellar Spectral Classifications (Skiff, 2009-2014)
a high resolution spectrum R=50k was acquired by C.Buil'],\ 'HR 8107':['already in Bess as HD201836','discovered in 1989'],\ 'BD+28 4177':['','a high resolution spectrum R=50k was acquired by C.Buil'],\ 'BD+60 133':['','a high resolution spectrum R=50kwas acquired by C.Buil'],\ 'BD+60 2523':['','a high resolution spectrum R=50k was acquired by C.Buil'],\ 'HD 237126':['','a high resolution spectrum R=50k was acquired by C.Buil'],\ } AltNameDic={'HR 8107':'HD201836'} BasePath=sys.path[0].replace("\\","\\\\")+"\\\\" dbSourcePath=BasePath+"web-"+filename PathWeb=BasePath+"web-"+filename+"\\\\" PathFileTmpDat=BasePath+"@.dat" PathGnuplot="C:\\gnuplot\\" WebPageFileName="auto-"+filename+".html" WebTitle="Spectroscopic observation of "+filename dirList= os.listdir(dbSourcePath) dirList=sorted(dirList) #efface les fichier .png present for file in dirList: if (file.endswith(".png")): print "remove " + file os.remove(dbSourcePath+"\\"+file) #la generation du code html commence ici f=open(PathWeb+WebPageFileName,'w') f.write("\n\nBe search, Be candidate\n") f.write("""""") f.write("\n\n") f.write("

Be search survey

\n") f.write("Last update at "+time.strftime("%d/%m/%Y %H:%M")+" from fits files. ") f.write("Contact: thierry dot lemoult arobase gmail dot com
\n") f.write("

Introduction

\n") f.write('This page in a first attenpt to make a synthesis of a new Field: Be Search
') f.write(' ARAS Spectroscopy sub forum dedicated to Be Search
\n') f.write("

Be candidate list

\n") # tri des fichiers suivant la date d'observation ListOrder=[] for file in dirList: if (file.endswith(".fits") or file.endswith(".fit")): hdulist = pyfits.open(dbSourcePath+'/'+file) prihdr = hdulist[0].header ListOrder=ListOrder+[(file,prihdr['OBJNAME']+prihdr['DATE-OBS'])] ListOrder=sorted(ListOrder,key=lambda ListOrder:ListOrder[1]) dirList=[a[0] for a in ListOrder] # fin du tri des fichiers # creer un tableau resumant les observations f.write("""\n""") f.write("\n") for file in dirList: hdulist = pyfits.open(dbSourcePath+'/'+file) prihdr = hdulist[0].header try: exptime2=prihdr['EXPTIME2'] except KeyError: try: exptime2=str(prihdr['EXPTIME'])+"s" except KeyError: exptime2='undefine' try: instrument=prihdr['BSS_INST'] except KeyError: instrument='' f.write("") # nom alternatif si definit try: altName=AltNameDic[prihdr['OBJNAME']] except KeyError: altName='' f.write("") f.write("") # commentaire et statut si definit try: status=StatusDic[prihdr['OBJNAME']][0] comment=StatusDic[prihdr['OBJNAME']][1] except KeyError: status='' comment='' f.write("") f.write("") f.write("\n") f.write("\n") f.write("\n") f.write("\n") f.write("\n") f.write("\n") f.write("\n") f.write("") f.write("
NameName2Observation Date in UTStatusCommentExposureInstrumentfits file
"+prihdr['OBJNAME']+""+altName+""+prihdr['DATE-OBS']+""+status+""+comment+""+exptime2+""+instrument+""+file+"
HD 189847HIP 98525Already publish
not in Bess		Robin , C.Buil realised HR
published here http://arxiv.org/abs/1409.4668		Waiting for Andrew spectrum
HD 181751HIP 951092014-09-21Robin realised HRWaiting for Andrew spectrum
HD192445HIP99667Waiting for Andrew spectrum

\n") f.write('

Biblio

\n') f.write('
  • The frequency of Be Stars, Jaschek, 1983 http://adsabs.harvard.edu/abs/1983A%26A...117..357J
    \n') f.write('
  • Classification of Be Stars, Jaschek, 1980 http://adsabs.harvard.edu/abs/1980A%26AS...42..103J
    \n') f.write('
  • To be completed...\n') f.write('\n') f.write('

    Complete scan List

    \n') f.write('
  • List currently used by Thierry Lemoult automated telescope, extracted from SIMBAD 6<B mag < 10.9\n') f.write('
  • List of non giant B stars ( V < 10) from Simbad http://www.astronomie-amateur.fr/Documents%20Spectro/B_Cibles.xlsm\n') f.write('

    List of already observed stars

    \n') f.write('
  • Thierry Lemoult automated telescope\n') f.write('

    Be Candidate spectrum plot

    \n') #genere les graphiques f.write("""\n""") i=0 #index unique pour les nom de fichier for file in dirList: # boucle sur les fichiers fits #recuperation des mots clef du fichier fits print "File name: "+file hdulist = pyfits.open(dbSourcePath+'/'+file) prihdr = hdulist[0].header name=prihdr['OBJNAME'] dateObs=prihdr['DATE-OBS'] try: Instrument=prihdr['BSS_INST'] except KeyError: Instrument='undefine' try: exptime2=prihdr['EXPTIME2'] except KeyError: try: exptime2=str(prihdr['EXPTIME'])+"s" except KeyError: exptime2='undefine' try: resolStr=str(prihdr['BSS_ITRP']) except KeyError: resolStr="non definie" start=prihdr['CRVAL1'] step=prihdr['CDELT1'] nstep=prihdr['NAXIS1'] print "longueur d onde [" + str(start) + " ... " +str(start+step*nstep)+"] "+prihdr['CUNIT1'] print "************" #genere resume dans la premiere colonne. f.write("\n\n\n") WriteDatFile(hdulist[0].data,start,step) titre=name+" "+dateObs[0:10]+" "+dateObs[11:22]+" "+Instrument+" "+exptime2 if (nstep*step>1000): # si on couvre une large bande spectral, on fait des zoom... #definition des bandes a tracer graph=[[start+step*nstep/2.0,step*nstep/2.0,""],[6560,90,"Region Ha"]] else: # sinon on affiche juste tout le spectre graph=[[start+step*nstep/2.0,step*nstep/2.0,""]] for zoom in graph: i+=1 lmin=zoom[0]-zoom[1] lmax=zoom[0]+zoom[1] regionTitle=zoom[2] if (lmin>=start*0.98 and lmax<=(start+step*nstep)*1.02): [Imin,Imax]=foundMinMax(hdulist[0].data,start,step,lmin,lmax) if (Imin<-1): Imin=0 pngFileName=name+"_"+dateObs[0:10]+"_"+dateObs[11:13]+"h"+dateObs[14:16]+"m_"+regionTitle+"_"+str(i)+".png" #genere image WritePngFile(PathWeb+pngFileName,titre+" "+regionTitle,lmin,lmax,Imin,Imax) #insere le lien dans le html f.write(" \n") f.write("\n\n\n") # ferme le fichier html f.write("
    \n") f.write("

    "+name+"

    \n") f.write("
    • Date: "+dateObs[0:10]+" a "+dateObs[11:13]+"h"+dateObs[14:16]+" TU
      • \n") f.write("
    • Exposure time: "+exptime2+"
      • \n") f.write("
    • Instrument:"+Instrument+"
      • \n") f.write("
    • Resolution: "+resolStr+"
      • \n") try: Site=prihdr['BSS_SITE'] except KeyError: Site='undefined' f.write("
    • Site: "+Site+"
    \n") f.write("

    Spectrum fits file link:
    "+file+"\n

    \n\n") f.close()