Corrections for non-linearity

Corrections for non-linearity

The currently distributed version of the ORACDR does not correct for non-linearity. Hence a small script is provided here which will do the corrections. For applying the corrections, please copy the script to any directory and set the environment variable "ORAC_PRIMITIVE_DIR" to the path of the directory where the script is placed.

Click here to download the script.

Note that this correction applies only to the data collected using the ARC (formerly SDSU) controller commissioned with UIST in December 2006. For numbers specific to the old Edict system, please contact the instrument scientist.

For non-thermal data, the correction "y" is:
y = ( 1.7929E-14 * x³ ) - ( 1.9544E-10 * x² ) + ( 6.5558E-6 * x ) + 0.98657
for data where the counts are between 750 and 20000 ADUs.

For thermal data, the correction "y" is :
y = ( 3.6058E-27 * x⁶ ) - (4.1373E-22* x⁵) + ( 1.8632E-17 * x⁴ ) - ( 4.1209E-13 * x³ ) + ( 4.762E-9 * x² ) - ( 2.353E-5 * x ) + 1.0294
for data where the counts are between 890 and 36000 ADUs.

"x" is the observed counts on the raw data (in ADUs). "y" is the correction factor to be multiplied to "x" to derive the linearized images.

Note that for thermal images, the lower limit of 890 ADUs is adequate. We will always be above this level when we work in the thermal regime. However, for non-thermal images, a lower limit of 750 ADUs is not always sufficient. Improved corrections extending well below 750 ADUs will be available soon.

How do we tell ORACDR to use this script?