Sdss sample spectra
From MediaWiki
SDSS spektri
Rebinovanje spektara
Napisala sam proceduru koja rebinuje u ln skali spektre u zadatim granicama. U svom najprostijem obliku nalazi se u File:Ppxf example ngc3325 sdssdr8.tgz kao myrebin.pro. Daje identicne rezultate Kapelarijevoj log_rebin proceduri. Treba otpakovati tgz, snimiti u ppxf dir i pokrenuti ppxf_example_ngc3325_sdssdr8.tgz. Inace, galaksija u primeru je NGC 3325. Prosto, treba odrediti nove granice i izvrsiti linearnu interpolaciju izmedju novih i starih tacaka. Posto se skala menja iz linearne u logaritamsku, nove logaritamske granice moraju se vratiti na linearne:
lamRange = sxpar(head,'CRVAL1') + [0d,sxpar(head,'CDELT1')*(sxpar(head,'NAXIS1')-1d)] lambda=exp(range(alog(lamRange[0]),alog(lamRange[1]),n+1)) step=(lambda[n-1] - lambda[0])/(n-1) lim = (lamRange + [-0.5d,0.5d]*step ; Prvi i poslednji pixel moraju isto da se interp pa se produzavaju za pola newLambda = range(lim[0],lim[1],n+1) logLim = alog(lamRange + [-0.5d,0.5d]*step) ili se skracuje za velicinu velScale: if n_elements(velScale) gt 0 then begin ; Velocity scale is set by user logScale = velScale/c m = floor((logLim[1]-logLim[0])/logScale) ; Number of output pixels logLim[1] = logLim[0] + m*logScale endif else $ velScale = (logLim[1]-logLim[0])/m*c ; Only for output newLimits = exp(range(logLim[0],logLim[1],m+1)) ; Ovim logaritamske granice vracamo na normalne