iraf.net

---------- Forwarded message ----------
Date: Mon, 26 Sep 2005 14:16:37 -0700
From: Joshua Jacob Turner <jturner5@uoregon.edu>
To: Mike Fitzpatrick <fitz@tucana.tuc.noao.edu>
Subject: RE: @-filesHello Mike,sorry for the unclear message. I've made a script to do
cross-correlations, with surface fits and background subtractions and
the like, but want to run the same script on a lot of images. I
actually would like to compare all 50 in one cube to every single
image in another cube, eventually, but I mentioned the comparison of
one image to a cube for starters, if i could actually just do that, i
would be happy.my file[x,y.*] was really suppose to be file[*,*,*] as I am comparing
the whole image, and wanted to increment the 3rd spot in the image
section. i would think it would be an easy solution, putting a "for"
loop in the script, but for some reason I cant get it to work. I am
probably doing something silly...right now, my input is just file1 and file2, i guess i can turn the
input into an @-list maybe? i wasn't able to get that going either.
thanks,Josh

---------- Forwarded message ----------
Date: Mon, 26 Sep 2005 16:43:50 -0700 (MST)
From: Mike Fitzpatrick <fitz@tucana.tuc.noao.edu>
To: Joshua Jacob Turner <jturner5@uoregon.edu>
Cc: Mike Fitzpatrick <fitz@tucana.tuc.noao.edu>
Subject: RE: @-files
Hi Josh,
It sounds like your task is expecting a 2-D image plane and
so looping over each band of the image is your only choice. Try
something like for (i=1; i <= 50; i = i + 1) {
s1 = "cube1[*,*," // i // "]" # build "cube[*,*,<N>]" section
s2 = "cube2[*,*," // i // "]" compare (s1, s2)
}You can build the band sections on the fly using the string concatenation
operator as above and loop over each. If you have problems with your loop
send me a code fragment.
-Mike

---------- Forwarded message ----------
Date: Tue, 27 Sep 2005 10:36:39 -0700
From: Joshua Jacob Turner <jturner5@uoregon.edu>
To: Mike Fitzpatrick <fitz@tucana.tuc.noao.edu>
Subject: RE: @-filesabsolutely brilliant!! thanks mike, you are amazing. i actually was
about to send you a few variations of the code, and the different
errors it was giving me, but finally got it. this was EXTREMELY
helpful. for some reason it was giving me syntax errors with the "*"
but i just typed in the entire pixel range, so no biggy. thanks so
much.josh-------------------
>
> Hi Josh,
> It sounds like your task is expecting a 2-D image plane and
> so looping over each band of the image is your only choice. Try
> something like
>
> for (i=1; i <= 50; i = i + 1) {
> s1 = "cube1[*,*," // i // "]" # build "cube[*,*,<N>]"
section
> s2 = "cube2[*,*," // i // "]"
>
> compare (s1, s2)
> }
>
> You can build the band sections on the fly using the string
concatenation
> operator as above and loop over each. If you have problems with
your loop
> send me a code fragment.
>
>
> -Mike
>
>