Skip to content
Snippets Groups Projects
Commit e072616a authored by Roberto Susino's avatar Roberto Susino
Browse files

Fix headers in FITS extensions

parent 7a4420be
Branches
Tags Version-2.0.4
No related merge requests found
Show whitespace changes
Inline Side-by-side
metis_l2_prep_uv.pro
+ 26
− 8
View file @ e072616a
...@@ -189,10 +189,25 @@ pro metis_l2_prep_uv ...@@ -189,10 +189,25 @@ pro metis_l2_prep_uv
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null base_header = primary_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count' sxdelpar, base_header, 'EXTEND'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count' sxdelpar, base_header, 'DATASUM'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name' sxdelpar, base_header, 'CHECKSUM'
sxdelpar, base_header, 'COMMENT'
sxdelpar, base_header, 'HISTORY'
extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Pixel quality'
fxaddpar, extension_header, 'BUNIT', 'None'
fxaddpar, extension_header, 'DATAMIN', min(quality_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(quality_matrix, /nan)
fxaddpar, extension_header, 'COMMENT', 'Quality matrix values:'
fxaddpar, extension_header, 'COMMENT', ' NaN = saturated or null L0 pixel counts'
fxaddpar, extension_header, 'COMMENT', ' 0 = unreliable pixel value'
fxaddpar, extension_header, 'COMMENT', ' 1 = good pixel'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'UV') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'UV')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name, extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name, extension_header, /no_comment, /silent
...@@ -201,11 +216,14 @@ pro metis_l2_prep_uv ...@@ -201,11 +216,14 @@ pro metis_l2_prep_uv
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = data * sqrt(error) error_matrix = data * sqrt(error)
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Absolute error'
fxaddpar, extension_header, 'DATAMIN', min(error_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(error_matrix, /nan)
if not ref_detector then error_matrix = metis_rectify(error_matrix, 'UV') if not ref_detector then error_matrix = metis_rectify(error_matrix, 'UV')
fits_add_checksum, extension_header, float(error_matrix) fits_add_checksum, extension_header, float(error_matrix)
mwrfits, float(error_matrix), out_file_name, extension_header, /no_comment, /silent mwrfits, float(error_matrix), out_file_name, extension_header, /no_comment, /silent
......
metis_l2_prep_vl_generic.pro
+ 26
− 8
View file @ e072616a
...@@ -230,10 +230,25 @@ pro metis_l2_prep_vl_generic ...@@ -230,10 +230,25 @@ pro metis_l2_prep_vl_generic
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null base_header = primary_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count' sxdelpar, base_header, 'EXTEND'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count' sxdelpar, base_header, 'DATASUM'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name' sxdelpar, base_header, 'CHECKSUM'
sxdelpar, base_header, 'COMMENT'
sxdelpar, base_header, 'HISTORY'
extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Pixel quality'
fxaddpar, extension_header, 'BUNIT', 'None'
fxaddpar, extension_header, 'DATAMIN', min(quality_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(quality_matrix, /nan)
fxaddpar, extension_header, 'COMMENT', 'Quality matrix values:'
fxaddpar, extension_header, 'COMMENT', ' NaN = saturated or null L0 pixel counts'
fxaddpar, extension_header, 'COMMENT', ' 0 = unreliable pixel value'
fxaddpar, extension_header, 'COMMENT', ' 1 = good pixel'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name, extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name, extension_header, /no_comment, /silent
...@@ -242,11 +257,14 @@ pro metis_l2_prep_vl_generic ...@@ -242,11 +257,14 @@ pro metis_l2_prep_vl_generic
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = data * sqrt(error) error_matrix = data * sqrt(error)
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Absolute error'
fxaddpar, extension_header, 'DATAMIN', min(error_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(error_matrix, /nan)
if not ref_detector then data = metis_rectify(error_matrix, 'VL') if not ref_detector then data = metis_rectify(error_matrix, 'VL')
fits_add_checksum, extension_header, float(error_matrix) fits_add_checksum, extension_header, float(error_matrix)
mwrfits, float(error_matrix), out_file_name, extension_header, /no_comment, /silent mwrfits, float(error_matrix), out_file_name, extension_header, /no_comment, /silent
......
metis_l2_prep_vl_polariz.pro
+ 114
− 44
View file @ e072616a
...@@ -403,10 +403,25 @@ pro metis_l2_prep_vl_polariz ...@@ -403,10 +403,25 @@ pro metis_l2_prep_vl_polariz
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null base_header = primary_pb_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count' sxdelpar, base_header, 'EXTEND'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count' sxdelpar, base_header, 'DATASUM'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name' sxdelpar, base_header, 'CHECKSUM'
sxdelpar, base_header, 'COMMENT'
sxdelpar, base_header, 'HISTORY'
extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Pixel quality'
fxaddpar, extension_header, 'BUNIT', 'None'
fxaddpar, extension_header, 'DATAMIN', min(quality_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(quality_matrix, /nan)
fxaddpar, extension_header, 'COMMENT', 'Quality matrix values:'
fxaddpar, extension_header, 'COMMENT', ' NaN = saturated or null L0 pixel counts'
fxaddpar, extension_header, 'COMMENT', ' 0 = unreliable pixel value'
fxaddpar, extension_header, 'COMMENT', ' 1 = good pixel'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name[0], extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name[0], extension_header, /no_comment, /silent
...@@ -415,11 +430,14 @@ pro metis_l2_prep_vl_polariz ...@@ -415,11 +430,14 @@ pro metis_l2_prep_vl_polariz
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = intarr(header.naxis1, header.naxis2) error_matrix = intarr(header.naxis1, header.naxis2)
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Absolute error'
fxaddpar, extension_header, 'DATAMIN', min(error_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(error_matrix, /nan)
if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL') if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')
fits_add_checksum, extension_header, error_matrix fits_add_checksum, extension_header, error_matrix
mwrfits, float(error_matrix), out_file_name[0], extension_header, /no_comment, /silent mwrfits, float(error_matrix), out_file_name[0], extension_header, /no_comment, /silent
...@@ -457,10 +475,25 @@ pro metis_l2_prep_vl_polariz ...@@ -457,10 +475,25 @@ pro metis_l2_prep_vl_polariz
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null base_header = primary_tb_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count' sxdelpar, base_header, 'EXTEND'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count' sxdelpar, base_header, 'DATASUM'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name' sxdelpar, base_header, 'CHECKSUM'
sxdelpar, base_header, 'COMMENT'
sxdelpar, base_header, 'HISTORY'
extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Pixel quality'
fxaddpar, extension_header, 'BUNIT', 'None'
fxaddpar, extension_header, 'DATAMIN', min(quality_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(quality_matrix, /nan)
fxaddpar, extension_header, 'COMMENT', 'Quality matrix values:'
fxaddpar, extension_header, 'COMMENT', ' NaN = saturated or null L0 pixel counts'
fxaddpar, extension_header, 'COMMENT', ' 0 = unreliable pixel value'
fxaddpar, extension_header, 'COMMENT', ' 1 = good pixel'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name[1], extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name[1], extension_header, /no_comment, /silent
...@@ -469,11 +502,14 @@ pro metis_l2_prep_vl_polariz ...@@ -469,11 +502,14 @@ pro metis_l2_prep_vl_polariz
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = intarr(header.naxis1, header.naxis2) error_matrix = intarr(header.naxis1, header.naxis2)
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Absolute error'
fxaddpar, extension_header, 'DATAMIN', min(error_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(error_matrix, /nan)
if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL') if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')
fits_add_checksum, extension_header, error_matrix fits_add_checksum, extension_header, error_matrix
mwrfits, float(error_matrix), out_file_name[1], extension_header, /no_comment, /silent mwrfits, float(error_matrix), out_file_name[1], extension_header, /no_comment, /silent
...@@ -511,10 +547,25 @@ pro metis_l2_prep_vl_polariz ...@@ -511,10 +547,25 @@ pro metis_l2_prep_vl_polariz
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null base_header = primary_polangle_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count' sxdelpar, base_header, 'EXTEND'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count' sxdelpar, base_header, 'DATASUM'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name' sxdelpar, base_header, 'CHECKSUM'
sxdelpar, base_header, 'COMMENT'
sxdelpar, base_header, 'HISTORY'
extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Pixel quality'
fxaddpar, extension_header, 'BUNIT', 'None'
fxaddpar, extension_header, 'DATAMIN', min(quality_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(quality_matrix, /nan)
fxaddpar, extension_header, 'COMMENT', 'Quality matrix values:'
fxaddpar, extension_header, 'COMMENT', ' NaN = saturated or null L0 pixel counts'
fxaddpar, extension_header, 'COMMENT', ' 0 = unreliable pixel value'
fxaddpar, extension_header, 'COMMENT', ' 1 = good pixel'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name[2], extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name[2], extension_header, /no_comment, /silent
...@@ -523,11 +574,14 @@ pro metis_l2_prep_vl_polariz ...@@ -523,11 +574,14 @@ pro metis_l2_prep_vl_polariz
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = intarr(header.naxis1, header.naxis2) error_matrix = intarr(header.naxis1, header.naxis2)
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Absolute error'
fxaddpar, extension_header, 'DATAMIN', min(error_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(error_matrix, /nan)
if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL') if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')
fits_add_checksum, extension_header, error_matrix fits_add_checksum, extension_header, error_matrix
mwrfits, float(error_matrix), out_file_name[2], extension_header, /no_comment, /silent mwrfits, float(error_matrix), out_file_name[2], extension_header, /no_comment, /silent
...@@ -587,15 +641,15 @@ pro metis_l2_prep_vl_polariz ...@@ -587,15 +641,15 @@ pro metis_l2_prep_vl_polariz
; add the extension with the stokes q image ; add the extension with the stokes q image
extension_header = !null extension_header = primary_stokes_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Stokes Q', 'extension name' fxaddpar, extension_header, 'EXTNAME', 'Stokes Q', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Stokes Q' fxaddpar, extension_header, 'BTYPE', 'Stokes Q'
fxaddpar, extension_header, 'BUNIT', cal_pack.vl_channel.cal_units fxaddpar, extension_header, 'BUNIT', cal_pack.vl_channel.cal_units
fxaddpar, extension_header, 'DATAMIN', min(q, /nan) fxaddpar, extension_header, 'DATAMIN', min(q, /nan)
fxaddpar, extension_header, 'DATAMAX', max(q, /nan) fxaddpar, extension_header, 'DATAMAX', max(q, /nan)
sxdelpar, extension_header, 'EXTEND'
if not ref_detector then q = metis_rectify(q, 'VL') if not ref_detector then q = metis_rectify(q, 'VL')
fits_add_checksum, extension_header, q fits_add_checksum, extension_header, q
mwrfits, float(q), out_file_name[3], extension_header, /no_comment, /silent mwrfits, float(q), out_file_name[3], extension_header, /no_comment, /silent
...@@ -604,15 +658,15 @@ pro metis_l2_prep_vl_polariz ...@@ -604,15 +658,15 @@ pro metis_l2_prep_vl_polariz
; add the extension with the stokes u image ; add the extension with the stokes u image
extension_header = !null extension_header = primary_stokes_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count' fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count' fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Stokes U', 'extension name' fxaddpar, extension_header, 'EXTNAME', 'Stokes U', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Stokes U' fxaddpar, extension_header, 'BTYPE', 'Stokes U'
fxaddpar, extension_header, 'BUNIT', cal_pack.vl_channel.cal_units fxaddpar, extension_header, 'BUNIT', cal_pack.vl_channel.cal_units
fxaddpar, extension_header, 'DATAMIN', min(u, /nan) fxaddpar, extension_header, 'DATAMIN', min(u, /nan)
fxaddpar, extension_header, 'DATAMAX', max(u, /nan) fxaddpar, extension_header, 'DATAMAX', max(u, /nan)
sxdelpar, extension_header, 'EXTEND'
if not ref_detector then u = metis_rectify(u, 'VL') if not ref_detector then u = metis_rectify(u, 'VL')
fits_add_checksum, extension_header, u fits_add_checksum, extension_header, u
mwrfits, float(u), out_file_name[3], extension_header, /no_comment, /silent mwrfits, float(u), out_file_name[3], extension_header, /no_comment, /silent
...@@ -621,11 +675,25 @@ pro metis_l2_prep_vl_polariz ...@@ -621,11 +675,25 @@ pro metis_l2_prep_vl_polariz
; add the extension with the quality matrix ; add the extension with the quality matrix
extension_header = !null base_header = primary_stokes_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count' sxdelpar, base_header, 'EXTEND'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count' sxdelpar, base_header, 'DATASUM'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name' sxdelpar, base_header, 'CHECKSUM'
sxdelpar, base_header, 'COMMENT'
sxdelpar, base_header, 'HISTORY'
extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Quality matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Pixel quality'
fxaddpar, extension_header, 'BUNIT', 'None'
fxaddpar, extension_header, 'DATAMIN', min(quality_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(quality_matrix, /nan)
fxaddpar, extension_header, 'COMMENT', 'Quality matrix values:'
fxaddpar, extension_header, 'COMMENT', ' NaN = saturated or null L0 pixel counts'
fxaddpar, extension_header, 'COMMENT', ' 0 = unreliable pixel value'
fxaddpar, extension_header, 'COMMENT', ' 1 = good pixel'
if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL') if not ref_detector then quality_matrix = metis_rectify(quality_matrix, 'VL')
fits_add_checksum, extension_header, quality_matrix fits_add_checksum, extension_header, quality_matrix
mwrfits, float(quality_matrix), out_file_name[3], extension_header, /no_comment, /silent mwrfits, float(quality_matrix), out_file_name[3], extension_header, /no_comment, /silent
...@@ -634,12 +702,14 @@ pro metis_l2_prep_vl_polariz ...@@ -634,12 +702,14 @@ pro metis_l2_prep_vl_polariz
; add the extension with the error matrix ; add the extension with the error matrix
extension_header = !null extension_header = base_header
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name'
error_matrix = intarr(header.naxis1, header.naxis2) error_matrix = intarr(header.naxis1, header.naxis2)
fxaddpar, extension_header, 'PCOUNT', 0, 'parameter count', before = 'LONGSTRN'
fxaddpar, extension_header, 'GCOUNT', 1, 'group count', before = 'LONGSTRN'
fxaddpar, extension_header, 'EXTNAME', 'Error matrix', 'extension name', before = 'LONGSTRN'
fxaddpar, extension_header, 'BTYPE', 'Absolute error'
fxaddpar, extension_header, 'DATAMIN', min(error_matrix, /nan)
fxaddpar, extension_header, 'DATAMAX', max(error_matrix, /nan)
if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL') if not ref_detector then error_matrix = metis_rectify(error_matrix, 'VL')
fits_add_checksum, extension_header, error_matrix fits_add_checksum, extension_header, error_matrix
mwrfits, float(error_matrix), out_file_name[3], extension_header, /no_comment, /silent mwrfits, float(error_matrix), out_file_name[3], extension_header, /no_comment, /silent
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment