Home > matutils > beamwidth.m

beamwidth

PURPOSE

width=beamwidth(ad,beam,lev)

SYNOPSIS

function width=beamwidth(ad,beam,lev)

DESCRIPTION

 width=beamwidth(ad,beam,lev)

 Measure the width of a 2D input beam
 by taking a simple 1D cut across it

 Optional arg lev defaults to 0.5 giving FWHM

CROSS-REFERENCE INFORMATION

This function calls:

This function is called by:

SOURCE CODE

0001 function width=beamwidth(ad,beam,lev)
0002 % width=beamwidth(ad,beam,lev)
0003 %
0004 % Measure the width of a 2D input beam
0005 % by taking a simple 1D cut across it
0006 %
0007 % Optional arg lev defaults to 0.5 giving FWHM
0008 
0009 if(~exist('lev'))
0010   lev=0.5;
0011 end
0012 
0013 c=ad.N_pix/2+1;
0014 beam1d=beam(c,c:end)/beam(c,c);
0015 width=interp1(beam1d,ad.t_val_deg(c:end),lev)*2;
0016 
0017 return