[daip] [!7632]: AIPS - IMFIT
Lynn Matthews
do-not-reply at nrao.edu
Fri Nov 20 13:38:39 EST 2015
Lynn Matthews updated #7632
---------------------------
IMFIT
-----
Ticket ID: 7632
URL: https://help.nrao.edu/staff/index.php?/Tickets/Ticket/View/7632
Full Name: Lynn Matthews
Email: lmatthew at haystack.mit.edu
Creator: User
Department: AIPS Data Reduction
Staff (Owner): -- Unassigned --
Type: Issue
Status: Open
Priority: Default
SLA: NRAO E2E
Template Group: Default
Created: 20 November 2015 06:38 PM
Updated: 20 November 2015 06:38 PM
Due: 24 November 2015 06:38 PM (4d 0h 0m)
Resolution Due: 30 November 2015 06:38 PM (10d 0h 0m)
I am experiencing some unexpected behavior with IMFIT. Regardless of the input image, if the box around the source is fairly
tight, for a single-component Gaussian fit reports messages like the following:
.
.
.
IMFIT1: MORE POINTS IN CALL TO FUNCTION THAN WE HAVE
IMFIT1: MORE POINTS IN CALL TO FUNCTION THAN WE HAVE
IMFIT1: STRANGE VALUE FOR COMP= 1 PARAM= 2 VALUE= 5.115E+02
IMFIT1: RIDICULOUS VALUE FOR SOME PARAMETER AT NITER= 15
IMFIT1: CHECK INPUT MODEL. MODEL VALUE IS CRAZY
IMFIT1: returns adverbs to AIPS
IMFIT1: Purports to die of UNNATURAL causes
JMFIT has no problem with the same inputs (appended below). If I make the box larger for IMFIT, it will do only a single iteration, and report something like the following:
IMFIT1: IMPROPER INPUT PARAMETERS. NITER= 1
IMFIT1: ****** RMS of fit= 1.1690E-05 per point ******
IMFIT1: Begin tweaking each parameter
IMFIT1: MORE POINTS IN CALL TO FUNCTION THAN WE HAVE
IMFIT1: MORE POINTS IN CALL TO FUNCTION THAN WE HAVE
.
.
.
IMFIT1: ****** New RMS of fit= 2.4454E-06 per point ******
IMFIT1: Writing ST file version 3 attached to input image
IMFIT1:
IMFIT1: ********* Solution from IMFIT *********************************
IMFIT1: Error estimate based on rms 4.711E-06 JY/BEAM
IMFIT1:
IMFIT1: Component 1 - Gaussian
IMFIT1: Peak intensity = 1.0160E-04 +/- 4.71E-06 JY/BEAM ( 21.57)
.
.
.
IMFIT1: returns adverbs to AIPS
IMFIT1: Appears to have ended successfully
Lynn Matthews
AIPS 1: IMFIT: Task to fit gaussian models to an image
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME 'O CETI AP-FX' Input image name (name)
AIPS 1: INCLASS 'ICL001' Input image name (class)
AIPS 1: INSEQ 4 Input image name (seq. #)
AIPS 1: INDISK 2 Input image disk drive #
AIPS 1: BLC 119 120 Bottom Left corner of fit
AIPS 1: 1 1 1 1
AIPS 1: 1
AIPS 1: TRC 140 139 Top Right corner of fit
AIPS 1: 1 1 1 1
AIPS 1: 1
AIPS 1: OUTNAME ' ' Output image name (name)
AIPS 1: OUTCLASS ' ' Output image name (class)
AIPS 1: OUTSEQ 0 Output image name (seq. #)
AIPS 1: OUTDISK 0 Output image disk drive #
AIPS 1: NGAUSS 1 Number of components
AIPS 1: CTYPE *all 0 Model types; one for each
AIPS 1: component
AIPS 1: 0->1=Gaussian
AIPS 1: 2=Zero level; 3=Zero+slope
AIPS 1: 4=Zero+slope+curve;
AIPS 1: 5=see HELP IMFIT
AIPS 1: Guess of model parameters
AIPS 1: GMAX *all 0 Peak of component (JY)
AIPS 1: 0-> Use maximum value
AIPS 1: GPOS *all 0 (X,Y) position (pixels)
AIPS 1: 0-> Use position of max
AIPS 1: GWIDTH *all 0 (BMAJ, BMIN, PA) of comp.
AIPS 1: (pixels,pixels,deg)
AIPS 1: 0->Use clean beam
AIPS 1: DOMAX 1.804E-05 1 Solve for GMAX? >0 -> yes
AIPS 1: 1 1 returns the error
AIPS 1: DOPOS 0.01285 0.00887 Solve for GPOS? >0 -> yes
AIPS 1: 1 1 returns the error
AIPS 1: 1 1 1 1
AIPS 1: DOWIDTH 0.0304 0.02066 Solve for GWIDTH? >0-> yes
AIPS 1: 0.3956 1 returns the error
AIPS 1: 1 1 1 1
AIPS 1: 1 1 1 1
AIPS 1: BWSMEAR 0 Bandwidth smearing corr.
AIPS 1: RADIUS 0 Radius for finding RMS
AIPS 1: NITER 20 Maximum # of iteration
AIPS 1: 0->200*NGAUSS
AIPS 1: Solve for model parameters?
AIPS 1: DOPRINT 0 <0 -> Print maps and
AIPS 1: solutions on line printer
AIPS 1: >0 -> print on terminal
AIPS 1: =0 -> maps not printed
AIPS 1: FITOUT ' '
AIPS 1: Disk file to save fit info
AIPS 1: DOOUTPUT -1 >0 -> 1 Catalog residual map
AIPS 1: OFFSET 0 Cutoff level. 0-> None
AIPS 1: DOMODEL -1 > 0 => put solutions in a CC
AIPS 1: file with input image
AIPS 1: OUTVERS 0 MF table version number
AIPS 1: -1 => none, 0 => new
AIPS 1: STVERS 0 STars output file version
AIPS 1: -1 => none, 0 => new
AIPS 1: PBPARM *all 0 Primary beam parameters:
AIPS 1: (1) level to believe - <= 0
AIPS 1: means do not apply a primary
AIPS 1: beam (2) > 0 use (3)-(7)
AIPS 1: EFACTOR 0 Scale width sigmas in the
AIPS 1: deconvolution tests
------------------------------------------------------
Staff CP: https://help.nrao.edu/staff
More information about the Daip
mailing list