---+ Example Atmosphere and related code provided by Friends
Code for Subtracting the atmosphere using PCA, from Tom Crawford (SPT) 2006.
This code was given to Tom from Sunil (Caltech) who uses it with Bolocam. It DOES NOT run on IDL, but instead uses the open source version called Octave. See tom's comments about porting it to IDL.
-------- Original Message --------
Subject: more atmo stuff
Date: Fri, 21 Apr 2006 19:23:44 -0500 (CDT)
From: Tom Crawford
Reply-To: Tom Crawford
To: Matt Dobbs
Hi Matt -
Here's the PCA routine and some stuff you'll need to run the array-mode
projector I sent you in the other e-mail. The PCA routine makes an
external call to an octave routine called "exteigen" (which I've attached)
to diagonalize the covariance matrix and find its eigenvectors and
eigenvalues. This is only because octave is much faster than IDL at
linear algebra -- it may be easier to start out doing this the slow way in
IDL, and you can do this by uncommenting some lines in pca_clean_tc.pro
and commenting out some others (it should be clear which ones; if not,
just let me know).
Have fun!
Tom
Code for simulating, measureing, and subtracting the atmosphere using the Triangular Fourier method, from Tom Crawford (SPT) 2004.
My array-mode-projection method is buried inside that
"spt_scan_sim_paredown_temp.pro" routine. Unfortunately, it doesn't exist
in any much cleaner state.
-------- Original Message --------
Subject: atmo code
Date: Fri, 08 Oct 2004 18:24:33 -0500 (CDT)
From: Tom Crawford
Reply-To: Tom Crawford
To: Matt Dobbs
Hi Matt -
I've attached the code I use to make the screen of atmospheric
fluctuations that get blown past / scanned by the detector array. I've
also attached the huge, awful routine that does the observation and the
atmospheric subtraction, because buried somewhere in there is how I
incorporate the atmosphere screen into sky observations. The atmospheric
screen creator is pretty simple and should make sense if you've read the
Lay & Halverson paper, but the other piece of kludged-up spaghetti code
may be impenetrable, in which case feel free to ask questions.
Enjoy.
--
Tom Crawford
University of Chicago
Department of Astronomy & Astrophysics
773-834-2103
This topic: APEX_SZ
> WebHome >
2006_AtmosphereSubtraction > 2006_AtmosphereCode
Topic revision: r1 - 2006-05-11 - MattDobbs