;-------------------------------------------------------------------
; function : M A K E _ C O H
;
; written by Stéphane Guillaso
;-------------------------------------------------------------------
; description :
;  	calculate the complex coherency between 2 images
;
; input :
;		im1 --> image 1 (complex data)
;		im2 --> image 2 (complex data)
;		fe --> line array containing the flatearth phase information
;				 (do fe=0 if you don't need this)
;		coef --> coeficient for smooth function
;
; output :
;		complex coherency
;
; call :
;		coh = make_coh(im1,im2,fl,coef)
;-------------------------------------------------------------------
; attention :
;		this function needs the make_interf function
;-------------------------------------------------------------------
; This software has been released under the terms of the GNU Public
; license. See http://www.gnu.org/copyleft/gpl.html for details.
;------------------------------------------------------------------

function make_coh,im1,im2,fl,coef
	return,smooth(make_interf(im1,im2,fe=fl),coef) / sqrt((abs(smooth(im1*conj(im1),coef))*abs(smooth(im2*conj(im2),coef))))
end