diff options
Diffstat (limited to 'noao/imred/vtel/readss2.x')
| -rw-r--r-- | noao/imred/vtel/readss2.x | 174 |
1 files changed, 174 insertions, 0 deletions
diff --git a/noao/imred/vtel/readss2.x b/noao/imred/vtel/readss2.x new file mode 100644 index 00000000..71ae8758 --- /dev/null +++ b/noao/imred/vtel/readss2.x @@ -0,0 +1,174 @@ +include <mach.h> +include <imhdr.h> +include <fset.h> +include "vt.h" + +define WDSBRSTR 50 + +# READSS2 -- Read a type 2 sector scan from tape and format into 3 iraf images. +# Type two sector scans consist of three images with 16 bits per +# pixel. The three images are 1. velocity (16 bits) 2. select (16 bits) and +# 3. brightness (16 bits). The images are 512 pixels high. + +procedure readss2 (inputfd, filenumber, brief, select, bright, velocity, hs) + +int inputfd # file descriptor for input (usually tape) +int filenumber # file number on tape +bool brief # short output file names +bool select # flag to make select image +bool bright # flag to make bright image +bool velocity # flag to make velocity image +int hs # header data structure pointer + +char velimage[SZ_FNAME] # velocity image +char selimage[SZ_FNAME] # select image +char britimage[SZ_FNAME] # brightness image +short u[SWTH_HIGH] +int date, hour, minute, seconds, i, j, num, lrs +pointer velim, selim, britim, velsrp, selsrp, britsrp + +int read() +pointer immap(), impl2s() +errchk immap, impl2s + +begin + # Calculate the time. Assemble the output image names. + hour = int(VT_HTIME(hs)/3600) + minute = int((VT_HTIME(hs) - hour * 3600)/60) + seconds = int(VT_HTIME(hs) - hour * 3600 - minute * 60) + if (brief) { + call sprintf (velimage[1], SZ_FNAME, "v%03d") + call pargi (filenumber) + call sprintf (selimage[1], SZ_FNAME, "s%03d") + call pargi (filenumber) + call sprintf (britimage[1], SZ_FNAME, "b%03d") + call pargi (filenumber) + } else { + call sprintf (velimage[1], SZ_FNAME, "v%02d_%02d%02d_%03d") + call pargi (VT_HDAY(hs)) # day of month + call pargi (hour) + call pargi (minute) + call pargi (filenumber) + call sprintf (selimage[1], SZ_FNAME, "s%02d_%02d%02d_%03d") + call pargi (VT_HDAY(hs)) # day of month + call pargi (hour) + call pargi (minute) + call pargi (filenumber) + call sprintf (britimage[1], SZ_FNAME, "b%02d_%02d%02d_%03d") + call pargi (VT_HDAY(hs)) # day of month + call pargi (hour) + call pargi (minute) + call pargi (filenumber) + } + + if (select) { + selim = immap (selimage, NEW_IMAGE, 0) + IM_NDIM(selim) = 2 + IM_LEN(selim,1) = SWTH_HIGH + IM_LEN(selim,2) = VT_HNUMCOLS(hs) + IM_PIXTYPE(selim) = TY_SHORT + call imaddi (selim, "obs_time", VT_HTIME(hs)) + date = VT_HMONTH(hs) * 10000 + VT_HDAY(hs) * 100 + VT_HYEAR(hs) + call imaddi (selim, "obs_date", date ) + call imaddi (selim, "wv_lngth", VT_HWVLNGTH(hs)) + call imaddi (selim, "obs_type", VT_HOBSTYPE(hs)) + call imaddi (selim, "av_intns", VT_HAVINTENS(hs)) + call imaddi (selim, "num_cols", VT_HNUMCOLS(hs)) + call imaddi (selim, "intg/pix", VT_HINTGPIX(hs)) + call imaddi (selim, "rep_time", VT_HREPTIME(hs)) + } + if (bright) { + britim = immap (britimage, NEW_IMAGE, 0) + IM_NDIM(britim) = 2 + IM_LEN(britim,1) = SWTH_HIGH + IM_LEN(britim,2) = VT_HNUMCOLS(hs) + IM_PIXTYPE(britim) = TY_SHORT + call imaddi (britim, "obs_time", VT_HTIME(hs)) + date = VT_HMONTH(hs) * 10000 + VT_HDAY(hs) * 100 + VT_HYEAR(hs) + call imaddi (britim, "obs_date", date ) + call imaddi (britim, "wv_lngth", VT_HWVLNGTH(hs)) + call imaddi (britim, "obs_type", VT_HOBSTYPE(hs)) + call imaddi (britim, "av_intns", VT_HAVINTENS(hs)) + call imaddi (britim, "num_cols", VT_HNUMCOLS(hs)) + call imaddi (britim, "intg/pix", VT_HINTGPIX(hs)) + call imaddi (britim, "rep_time", VT_HREPTIME(hs)) + } + if (velocity) { + velim = immap (velimage, NEW_IMAGE, 0) + IM_NDIM(velim) = 2 + IM_LEN(velim,1) = SWTH_HIGH + IM_LEN(velim,2) = VT_HNUMCOLS(hs) + IM_PIXTYPE(velim) = TY_SHORT + call imaddi (velim, "obs_time", VT_HTIME(hs)) + date = VT_HMONTH(hs) * 10000 + VT_HDAY(hs) * 100 + VT_HYEAR(hs) + call imaddi (velim, "obs_date", date ) + call imaddi (velim, "wv_lngth", VT_HWVLNGTH(hs)) + call imaddi (velim, "obs_type", VT_HOBSTYPE(hs)) + call imaddi (velim, "av_intns", VT_HAVINTENS(hs)) + call imaddi (velim, "num_cols", VT_HNUMCOLS(hs)) + call imaddi (velim, "intg/pix", VT_HINTGPIX(hs)) + call imaddi (velim, "rep_time", VT_HREPTIME(hs)) + } + + do j = 1, VT_HNUMCOLS(hs) { + if (select) + selsrp = impl2s (selim, j) + if (bright) + britsrp = impl2s (britim, j) + if (velocity) + velsrp = impl2s (velim, j) + + iferr (num = read (inputfd, u, SWTH_HIGH*SZB_SHORT/SZB_CHAR)) { + call fseti (inputfd, F_VALIDATE, lrs*SZB_SHORT/SZB_CHAR) + call eprintf ("Error on tape read.\n") + num = read (inputfd, u, SWTH_HIGH*SZB_SHORT/SZB_CHAR) + } + lrs = num + if (num < SWTH_HIGH*SZB_SHORT/SZB_CHAR) + call error (0, "eof encountered when reading file") + if (BYTE_SWAP2 == YES) + call bswap2 (u, 1, u, 1, SWTH_HIGH * SZB_SHORT) + + if (velocity) + do i = 1, 512 + Mems[velsrp+i-1] = u[i] + + iferr (num = read (inputfd, u, SWTH_HIGH*SZB_SHORT/SZB_CHAR)) { + call fseti (inputfd, F_VALIDATE, lrs*SZB_SHORT/SZB_CHAR) + call eprintf ("Error on tape read.\n") + num = read (inputfd, u, SWTH_HIGH*SZB_SHORT/SZB_CHAR) + } + lrs = num + if (num < SWTH_HIGH*SZB_SHORT/SZB_CHAR) + call error (0, "eof encountered when reading file") + if (BYTE_SWAP2 == YES) + call bswap2 (u, 1, u, 1, SWTH_HIGH * SZB_SHORT) + + if (select) + do i = 1, 512 + Mems[selsrp+i-1] = u[i] + + iferr (num = read (inputfd, u, SWTH_HIGH*SZB_SHORT/SZB_CHAR)) { + call fseti (inputfd, F_VALIDATE, lrs*SZB_SHORT/SZB_CHAR) + call eprintf ("Error on tape read.\n") + num = read (inputfd, u, SWTH_HIGH*SZB_SHORT/SZB_CHAR) + } + lrs = num + if (num < SWTH_HIGH*SZB_SHORT/SZB_CHAR) + call error (0, "eof encountered when reading file") + if (BYTE_SWAP2 == YES) + call bswap2 (u, 1, u, 1, SWTH_HIGH * SZB_SHORT) + + if (bright) + do i = 1, 512 + Mems[britsrp+i-1] = u[i] + } + + # Unmap images. + if (select) + call imunmap (selim) + if (velocity) + call imunmap (velim) + if (bright) + call imunmap (britim) +end |