IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Oct 3, 2008, 11:00:18 AM (18 years ago)
Author:
eugene
Message:

converted pmMoments elements to Mxx, Myy, Mxxx, etc; using pmPeak position in most cases instead of moments

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/psphot/src/models/pmModel_STRAIL.c

    r15843 r19881  
    478478{
    479479    pmMoments *Smoments = source->moments;
     480    pmPeak    *peak    = source->peak;
    480481    psF32     *params  = model->params->data.F32;
    481482
     
    484485    psEllipseMoments moments;
    485486
    486     moments.x2 = PS_SQR(Smoments->Sx);
    487     moments.y2 = PS_SQR(Smoments->Sy);
    488     moments.xy = Smoments->Sxy;
     487    moments.x2 = Smoments->Mxx;
     488    moments.y2 = Smoments->Myy;
     489    moments.xy = Smoments->Mxy;
    489490    //sometimes these moment inputs are zero...why?
    490491
     
    494495    shape = psEllipseAxesToShape(axes);
    495496
    496     params[0] = Smoments->Sky;
    497     params[1] = Smoments->Peak - Smoments->Sky;
    498     params[2] = Smoments->x;
    499     params[3] = Smoments->y;
     497    params[PM_PAR_SKY]  = Smoments->Sky;
     498    params[PM_PAR_I0]   = peak->flux;
     499    params[PM_PAR_XPOS] = peak->xf;
     500    params[PM_PAR_YPOS] = peak->yf;
     501
    500502    params[7] = 2 * axes.major;
    501503    params[8] = axes.theta;
Note: See TracChangeset for help on using the changeset viewer.