Commit 5c0777a9 authored by Roberto Susino's avatar Roberto Susino
Browse files

Optimize radiometric calibration procedure

parent aec75b24

metis_rad_cal.pro

+10 −5
Original line number Diff line number Diff line
@@ -12,7 +12,7 @@ function metis_rad_cal, data, header, cal_pack, polarimetric = polarimetric, nor 
		if keyword_set(polarimetric) then pmp_factor = 2. else pmp_factor = 1.



		ndit = header.ndit

	end

	endif



	if header.filter.contains('UV', /fold) then begin

		channel = cal_pack.uv_channel
@@ -82,12 +82,17 @@ function metis_rad_cal, data, header, cal_pack, polarimetric = polarimetric, nor 
		

		cal_factor = cal_factor/rad_factor

		

		cal_error = sqrt((radiometry.rad_response.error/radiometry.rad_response.value)^2 + (cal_pack.instrument.pupil_area.error/cal_pack.instrument.pupil_area.value)^2 + (channel.angular_pixel.error/channel.angular_pixel.value)^2 + (unit_error/unit_factor)^2)

		cal_error = sqrt((radiometry.rad_response.error/radiometry.rad_response.value)^2 $

			+ (cal_pack.instrument.pupil_area.error/cal_pack.instrument.pupil_area.value)^2 $

			+ (channel.angular_pixel.error/channel.angular_pixel.value)^2 $

			+ (unit_error/unit_factor)^2)

		

		if header.filter.contains('VL', /fold) then begin

			history = [history, '  cal. factor = ' + string(cal_factor, format = '(E8.2)') + ' ' + units + '/DN']



			journal, '  total cal. factor = ' + string(cal_factor, format = '(E8.2)') + ' ' + units + '/DN'

			journal, '  cal. factor error = ' + string(cal_error, format = '(E8.2)') + ' ' + units + '/DN'

		endif

	endelse

		

	; radiometric calibration