##################################################################
# phil massey, Lowell Observatory
###################################################################
procedure e2vlincor (quad,inimage,outimage)
#
string  quad            {prompt= "quadrant"}
file	inimage	        {prompt= "Input image "}
file	outimage	{prompt= "Output image \n"}
#

begin
string 	inimg, outimg, q
string a
real alpha, beta1, beta2,b,c,d
string expr

q=quad
inimg=inimage
outimg=outimage

      alpha=1.000
      beta1=0.0
      beta2=0.0 
    if(q=='c1') {
      alpha=1.1150
      beta1=-0.0338
      beta2=0.0057
      }
      else if (q=='c2') { 
      alpha=1.0123
      beta1=-0.0123
      beta2=0.0017
      }
      else if (q=='c3') {
      alpha=1.0000
      beta1=-0.0099
      beta2=0.0014
      }
      else if (q=='c4') {
      alpha=1.0707
      beta1=-0.0163
      beta2=0.0027
      }
    expr="a*b*(1.+c*(a/32000.)+d*(a/32000.)**2)"
    imexpr (expr=expr,output=outimg,a=inimg,b=alpha,c=beta1,d=beta2,\
    dims="auto",intype="auto",outtype="auto",refim="auto")
    print(inimg)
    hedit(outimg,'Lincor','Linearity corrected by e2vlincor.cl',add+,del-,\
    ver-,show-,update+)
end