Welcome to iraf.net Friday, March 29 2024 @ 10:28 AM GMT
Edward Olszewski |
05/24/2006 04:43PM (Read 4025 times)
|
|
|
Status: offline
Registered: 02/13/2006
Posts: 6
|
I was reducing some MMT Megacam data. I'd forgotten SAO's warning
that "there's a bug in merge amps in ccproc that adds a blank column",
and i didn't use their "megamerge" routine. So indeed I have an mef with
36 extensions of data, all of which are 1025x2304, where col 513 is
a colum of zeros. It should be 1024x2304. I don't have a backup copy of the "before merging" data,
so unless there's a clever way to get rif of col 513 in each extension I'll
have to completely re-reduce the data. Thoughts?I just looked at my 90prime data and the merge files are 4096x4096, ie,
no extra column. So it must be something special or weird about
megacam data that brings out this "bug". ed olszewski
|
|
|
|
fitz |
05/24/2006 04:43PM
|
|
|
Status: offline
Registered: 09/30/2005
Posts: 4040
|
Hi Ed,I'll let Frank comment on the "merge amps bug", but if what you want is remove a column from the middle of 36 extensions in a MEF then try something like:[code:1:301a6ceceb]
list = "myfiles.txt"
while (fscan (list, s1) != EOF) { # loop over file list
for (i=1; i <= 36; i=i+1) { # loop over extensions
imcopy (s1//"["//i//"][514:1025,*]",
s1//"["//i//"][513:1024,*]") # shift columns
imcopy (s1//"["//i//"][1:1024*]",
s1//"[i//",overwrite+]") # trim size
}
}
[/code:1:301a6ceceb]The file list should contain the ".fits" extension in the filenames or else be sure your 'imtype' is set for fits, but basically this just copied the section past the blank column to shift it over, and then overwrites the extension with a trimmed section so you end up with a 1024x2304 array.Note, I [b:301a6ceceb][u:301a6ceceb][color=red:301a6ceceb]strongly[/color:301a6ceceb][/u:301a6ceceb][/b:301a6ceceb] suggest you work on a backup copy of the data first and tweak the above as needed (i.e. I didn't really test it before posting). This could easily be made into a procedure script if you'll need it again. Hope this helps.Cheers,
-Mike
|
|
|
|
| |
|
Content generated in: 0.08 seconds |
|