IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
May 8, 2025, 4:29:52 PM (15 months ago)
Author:
eugene
Message:

merge changes from eam_branches/ipp-20230313

Location:
trunk/Ohana
Files:
19 edited

Legend:

Unmodified
Added
Removed

  • trunk/Ohana

  • trunk/Ohana/src/opihi/cmd.astro/cdensify.c

    r41515 r42821  
    11# include "data.h"
     2# define dCOS(A)   ((double) cos ((double)RAD_DEG*A))
    23
    34# define CHECKVAL(ARG) if (!isfinite(ARG)) { gprint (GP_ERR, "illegal value for %s: %f\n", #ARG, ARG); return (FALSE); }
    … …  
    116117  PutCoords (&newcoords, &bf[0].header);
    117118 
     119  // use the mask to prevent double-counting
     120  ALLOCATE_PTR (mask, char, Nx*Ny);
     121
    118122  float scalescale = scale*scale;
    119123  float scale2 = (scale + 1.0) * (scale + 1.0);
    … …  
    142146    coords.crval1 = rn;
    143147    coords.crval2 = *d;
     148   
     149    // XXX do not oversample by more than a factor of 10
     150    float dXn = dX * MAX(fabs(dCOS(*d)), 0.1);
    144151
    145152    float F = 1.0;
    146153    if (vv) { F = isFloatScale ? Fs[i] : Is[i]; }
     154
     155    // reset the mask so we do not double-count
     156    memset (mask, 0, Nx*Ny);
    147157
    148158    switch (PSFTYPE) {
    … …  
    162172        break;
    163173      case IS_SQUARE:
    164         for (ix = -scale; ix <= scale; ix += dX) {
     174        for (ix = -scale; ix <= scale; ix += dXn) {
    165175          for (iy = -scale; iy <= scale; iy += dY) {
    166176            double rp, dp;
    … …  
    184194        break;
    185195      case IS_CIRCLE:
    186         for (ix = -scale; ix <= scale; ix += dX) {
    187           for (iy = -scale; iy <= scale; iy += dY) {
     196        for (ix = -scale - dXn; ix <= scale + dXn; ix += dXn) {
     197          for (iy = -scale - dY; iy <= scale + dY; iy += dY) {
    188198            float r2 = ix*ix + iy*iy;
    189199            double rp, dp;
    … …  
    199209            if (Xb < 0) continue;
    200210            if (Yb < 0) continue;
     211            if (mask[Xb + Yb*Nx]) continue;
     212            mask[Xb + Yb*Nx] = 1;
    201213            if (vv) {
    202214              val[Xb + Yb*Nx] += Normalize ? fCircle*F : F;
    … …  
    208220        break;
    209221      case IS_GAUSS:
    210         for (ix = -3.0*scale; ix <= 3.0*scale; ix += dX) {
     222        for (ix = -3.0*scale; ix <= 3.0*scale; ix += dXn) {
    211223          for (iy = -3.0*scale; iy <= 3.0*scale; iy += dY) {
    212224            float r2 = ix*ix + iy*iy;
    … …  
    232244    }
    233245  }
     246
     247  free (mask);
    234248  return (TRUE);
    235249}

  • trunk/Ohana/src/opihi/cmd.astro/imfit-fgauss-pol.c

    r42078 r42821  
    5151  set_variable ("Zpk",  par[5]);
    5252  set_variable ("Sg",   par[6]);
     53
     54  /*
     55  set_variable ("dXg",   dpar[0]);
     56  set_variable ("dYg",   dpar[1]);
     57  set_variable ("dSXg",  NAN);
     58  set_variable ("dSYg",  NAN);
     59  set_variable ("dSXYg", NAN);
     60  set_variable ("dZpk",  dpar[5]);
     61  set_variable ("dSg",   dpar[6]);
     62  */
    5363}
    5464

  • trunk/Ohana/src/opihi/cmd.astro/imfit-fgauss.c

    r42078 r42821  
    3939  set_variable ("Zpk",  par[5]);
    4040  set_variable ("Sg",   par[6]);
     41
     42/*
     43  set_variable ("dXg",   dpar[0]);
     44  set_variable ("dYg",   dpar[1]);
     45  set_variable ("dSXg",  2.35 / dpar[2]);
     46  set_variable ("dSYg",  2.35 / dpar[3]);
     47  set_variable ("dSXYg", dpar[4]);
     48  set_variable ("dZpk",  dpar[5]);
     49  set_variable ("dSg",   dpar[6]);
     50*/
    4151}
    4252

  • trunk/Ohana/src/opihi/cmd.astro/imfit-pgauss-psf.c

    r39457 r42821  
    3535  set_variable ("Zpk",  par[2]);
    3636  set_variable ("Sg",   par[3]);
     37
     38/*
     39  set_variable ("dXg",   dpar[0]);
     40  set_variable ("dYg",   dpar[1]);
     41  set_variable ("dZpk",  dpar[2]);
     42  set_variable ("dSg",   dpar[3]);
     43*/
    3744}
    3845

  • trunk/Ohana/src/opihi/cmd.astro/imfit-pgauss.c

    r39457 r42821  
    3737  set_variable ("Zpk",  par[5]);
    3838  set_variable ("Sg",   par[6]);
     39
     40/*
     41  set_variable ("dXg",   dpar[0]);
     42  set_variable ("dYg",   dpar[1]);
     43  set_variable ("dSXg",  2.35 / dpar[2]);
     44  set_variable ("dSYg",  2.35 / dpar[3]);
     45  set_variable ("dSXYg", dpar[4]);
     46  set_variable ("dZpk",  dpar[5]);
     47  set_variable ("dSg",   dpar[6]);
     48*/
    3949}
    4050

  • trunk/Ohana/src/opihi/cmd.astro/imfit-q2gauss.c

    r42078 r42821  
    4141  set_variable ("Sg",   par[6]);
    4242  set_variable ("Sr", par[7]);
     43
     44/*
     45  set_variable ("dXg",   dpar[0]);
     46  set_variable ("dYg",   dpar[1]);
     47  set_variable ("dSXg",  2.35 / dpar[2]);
     48  set_variable ("dSYg",  2.35 / dpar[3]);
     49  set_variable ("dSXYg", dpar[4]);
     50  set_variable ("dZpk",  dpar[5]);
     51  set_variable ("dSg",   dpar[6]);
     52  set_variable ("dSr",   dpar[7]);
     53*/
    4354}
    4455

  • trunk/Ohana/src/opihi/cmd.astro/imfit-qfgauss.c

    r39457 r42821  
    4141  set_variable ("Zpk",  par[5]);
    4242  set_variable ("Sg",   par[6]);
     43
     44/*
     45  set_variable ("dXg",   dpar[0]);
     46  set_variable ("dYg",   dpar[1]);
     47  set_variable ("dSXg",  2.35 / dpar[2]);
     48  set_variable ("dSYg",  2.35 / dpar[3]);
     49  set_variable ("dSXYg", dpar[4]);
     50  set_variable ("dZpk",  dpar[5]);
     51  set_variable ("dSg",   dpar[6]);
     52*/
    4353}
    4454

  • trunk/Ohana/src/opihi/cmd.astro/imfit-qgauss-psf.c

    r39457 r42821  
    3838  set_variable ("Zpk", par[2]);
    3939  set_variable ("Sg",  par[3]);
     40
     41  /*
     42  set_variable ("dXg",  dpar[0]);
     43  set_variable ("dYg",  dpar[1]);
     44  set_variable ("dZpk", dpar[2]);
     45  set_variable ("dSg",  dpar[3]);
     46  */
    4047}
    4148

  • trunk/Ohana/src/opihi/cmd.astro/imfit-qgauss.c

    r42078 r42821  
    4141  set_variable ("Sg",   par[6]);
    4242  set_variable ("Sr",   par[7]);
     43
     44/*
     45  set_variable ("dXg",   dpar[0]);
     46  set_variable ("dYg",   dpar[1]);
     47  set_variable ("dSXg",  2.35 / dpar[2]);
     48  set_variable ("dSYg",  2.35 / dpar[3]);
     49  set_variable ("dSXYg", dpar[4]);
     50  set_variable ("dZpk",  dpar[5]);
     51  set_variable ("dSg",   dpar[6]);
     52  set_variable ("dSr",   dpar[7]);
     53*/
    4354}
    4455

  • trunk/Ohana/src/opihi/cmd.astro/imfit-qrgauss.c

    r42078 r42821  
    4242  set_variable ("Sg",   par[6]);
    4343  set_variable ("Npow", par[7]);
     44
     45/*
     46  set_variable ("dXg",   dpar[0]);
     47  set_variable ("dYg",   dpar[1]);
     48  set_variable ("dSXg",  2.35 / dpar[2]);
     49  set_variable ("dSYg",  2.35 / dpar[3]);
     50  set_variable ("dSXYg", dpar[4]);
     51  set_variable ("dZpk",  dpar[5]);
     52  set_variable ("dSg",   dpar[6]);
     53  set_variable ("dNpow",   dpar[7]);
     54*/
    4455}
    4556

  • trunk/Ohana/src/opihi/cmd.astro/imfit-r2gauss.c

    r42078 r42821  
    8585  set_variable ("SYf", 2.35 * sqrt(2.0) / par[8]);
    8686  set_variable ("SXYf", par[9]);
     87
     88/*
     89  set_variable ("dXg",   dpar[0]);
     90  set_variable ("dYg",   dpar[1]);
     91  set_variable ("dSXg",  2.35 * sqrt(2.0) / dpar[2]);
     92  set_variable ("dSYg",  2.35 * sqrt(2.0) / dpar[3]);
     93  set_variable ("dSXYg", dpar[4]);
     94  set_variable ("dZpk",  dpar[5]);
     95  set_variable ("dSg",   dpar[6]);
     96  set_variable ("dSXf", 2.35 * sqrt(2.0) / dpar[7]);
     97  set_variable ("dSYf", 2.35 * sqrt(2.0) / dpar[8]);
     98  set_variable ("dSXYf", dpar[9]);
     99*/
    87100}

  • trunk/Ohana/src/opihi/cmd.astro/imfit-rgauss.c

    r42078 r42821  
    4242  set_variable ("Sg",   par[6]);
    4343  set_variable ("Sr",   par[7]);
     44
     45/*
     46  set_variable ("dXg",   dpar[0]);
     47  set_variable ("dYg",   dpar[1]);
     48  set_variable ("dSXg",  2.35 / dpar[2]);
     49  set_variable ("dSYg",  2.35 / dpar[3]);
     50  set_variable ("dSXYg", dpar[4]);
     51  set_variable ("dZpk",  dpar[5]);
     52  set_variable ("dSg",   dpar[6]);
     53  set_variable ("dSr",   dpar[7]);
     54*/
    4455}
    4556

  • trunk/Ohana/src/opihi/cmd.astro/imfit-sgauss.c

    r39457 r42821  
    4949  set_variable ("SYf", 2.35 / par[8]);
    5050  set_variable ("SXYf", par[9]);
     51
     52/*
     53  set_variable ("dXg",   dpar[0]);
     54  set_variable ("dYg",   dpar[1]);
     55  set_variable ("dSXg",  2.35 / dpar[2]);
     56  set_variable ("dSYg",  2.35 / dpar[3]);
     57  set_variable ("dSXYg", dpar[4]);
     58  set_variable ("dZpk",  dpar[5]);
     59  set_variable ("dSg",   dpar[6]);
     60  set_variable ("dSXf", 2.35 / dpar[7]);
     61  set_variable ("dSYf", 2.35 / dpar[8]);
     62  set_variable ("dSXYf", dpar[9]);
     63*/
    5164}
    5265

  • trunk/Ohana/src/opihi/cmd.astro/imfit-tgauss.c

    r39457 r42821  
    8080  set_variable ("SYf",  2.35 * sqrt(2.0) / par[8]);
    8181  set_variable ("SXYf", par[9]);
     82
     83/*
     84  set_variable ("dXg",   dpar[0]);
     85  set_variable ("dYg",   dpar[1]);
     86  set_variable ("dSXg",  2.35 * sqrt(2.0) / dpar[2]);
     87  set_variable ("dSYg",  2.35 * sqrt(2.0) / dpar[3]);
     88  set_variable ("dSXYg", dpar[4]);
     89  set_variable ("dZpk",  dpar[5]);
     90  set_variable ("dSg",   dpar[6]);
     91  set_variable ("dSXf",  2.35 * sqrt(2.0) / dpar[7]);
     92  set_variable ("dSYf",  2.35 * sqrt(2.0) / dpar[8]);
     93  set_variable ("dSXYf", dpar[9]);
     94*/
    8295}

  • trunk/Ohana/src/opihi/cmd.astro/imfit-vgauss.c

    r39457 r42821  
    6969  set_variable ("SXf", par[7]);
    7070  set_variable ("SYf", par[8]);
     71
     72  /*
     73  set_variable ("dXg",   dpar[0]);
     74  set_variable ("dYg",   dpar[1]);
     75  set_variable ("dSXg",  2.35 / dpar[2]);
     76  set_variable ("dSYg",  2.35 / dpar[3]);
     77  set_variable ("dSXYg", dpar[4]);
     78  set_variable ("dZpk",  dpar[5]);
     79  set_variable ("dSg",   dpar[6]);
     80  set_variable ("dSXf",  dpar[7]);
     81  set_variable ("dSYf",  dpar[8]);
     82  */
    7183}

  • trunk/Ohana/src/opihi/cmd.astro/imfit.c

    r42078 r42821  
    1313  }
    1414
     15  Buffer *var = NULL;
     16  if ((N = get_argument (argc, argv, "-var-image"))) {
     17    remove_argument (N, &argc, argv);
     18    var = SelectBuffer (argv[N], OLDBUFFER, TRUE);
     19    if (!var) {
     20      gprint (GP_ERR, "unknown buffer for variance %s\n", argv[N]);
     21      FREE (Save);
     22      return (FALSE);
     23    }
     24    remove_argument (N, &argc, argv);
     25  }
     26
    1527  int Insert = FALSE;
    1628  if ((N = get_argument (argc, argv, "-insert"))) {
    … …  
    1830    Insert = TRUE;
    1931    if (Save) { gprint (GP_ERR, "-save and -insert are mutually exclusive\n"); free (Save); return (FALSE); }
     32    if (var)  { gprint (GP_ERR, "-save and -var-image are mutually exclusive\n"); return (FALSE); }
    2033  }
    2134
    … …  
    4053    Gain = atof(argv[N]);
    4154    remove_argument (N, &argc, argv);
     55    if (var) { gprint (GP_ERR, "warning: -var-image selected, -gain will have no effect\n"); }
    4256  }
    4357
    … …  
    4862    RDnoise = atof(argv[N]);
    4963    remove_argument (N, &argc, argv);
     64    if (var) { gprint (GP_ERR, "warning: -var-image selected, -rdnoise will have no effect\n"); }
    5065  }
    5166
    … …  
    140155    if (j + sy >= Ny) continue;
    141156    float *V = (float *)(buf[0].matrix.buffer) + (j+sy)*buf[0].matrix.Naxis[0] + sx;
     157    float *dV = (var) ? ((float *)(var[0].matrix.buffer) + (j+sy)*var[0].matrix.Naxis[0] + sx) : NULL;
    142158    for (int i = 0; i < dX; i++) {
    143159      if (i + sx < 0) continue;
    … …  
    145161      if (*V > SatThreshold) goto next; // skip pixels above threshold
    146162      if (!isfinite(*V)) goto next; // skip nan pixels
    147       dz[N] = (SQ(RDnoise) + MAX(0.0, *V/Gain)); // treat negative pixels as pure read noise
     163      if (dV) {
     164        dz[N] = *dV;
     165      } else {
     166        dz[N] = (SQ(RDnoise) + MAX(0.0, *V/Gain)); // treat negative pixels as pure read noise
     167      }
    148168      if (dz[N] <= 0) goto next;
    149169      dz[N] = 1.0 / dz[N];
    … …  
    154174    next:
    155175      V++;
     176      if (dV) { dV++; }
    156177    }
    157178  }

  • trunk/Ohana/src/opihi/cmd.astro/mkgauss.c

    r41341 r42821  
    1111  Buffer *buf;
    1212
     13  // if TRUE, integrated flux is Flux, else Io is Flux
    1314  int Normalize = FALSE;
    1415  if ((N = get_argument (argc, argv, "-norm"))) {
    … …  
    1617    remove_argument (N, &argc, argv);
    1718  }   
     19
     20  // Io or integrated flux (depending on value of Normalize)
     21  float Flux = 1.0;
     22  if ((N = get_argument (argc, argv, "-flux"))) {
     23    remove_argument (N, &argc, argv);
     24    Flux = atof(argv[N]);
     25    remove_argument (N, &argc, argv);
     26  }   
     27
    1828
    1929  // this should be Nx/2, Ny/2 if not set
    … …  
    2939  if ((argc < 3) || (argc > 5)) {
    3040    gprint (GP_ERR, "USAGE: mkgauss (buffer) (sigma) [[sy/sx] angle]\n");
     41    gprint (GP_ERR, " -flux flux : integral or peak is flux (default is 1.0)\n");
     42    gprint (GP_ERR, " -norm : integral is flux (else peak)\n");
     43    gprint (GP_ERR, " -c X Y : place center at X,Y\n");
    3144    return (FALSE);
    3245  }
    … …  
    6679  /* f = exp (-r), r = (x^2 / 2Sx) + (y^2 / 2Sy) + Sxy*x*y */
    6780
    68   double Io = Normalize ? 1.0 / (2.0 * M_PI * Sig_x * Sig_y) : 1.0;
     81  double Io = Normalize ? Flux / (2.0 * M_PI * Sig_x * Sig_y) : Flux;
    6982
    7083  in = (float *) buf[0].matrix.buffer;

  • trunk/Ohana/src/opihi/cmd.astro/region.c

    r42078 r42821  
    160160    if (!strcasecmp (argv[CtypeArg], "PAR")) { strcpy (graphmode.coords.ctype, "DEC--PAR"); goto got_ctype; }
    161161    if (!strcasecmp (argv[CtypeArg], "MOL")) { strcpy (graphmode.coords.ctype, "DEC--MOL"); goto got_ctype; }
     162    if (!strcasecmp (argv[CtypeArg], "LIN")) { strcpy (graphmode.coords.ctype, "DEC--LIN"); goto got_ctype; }
     163    if (!strcasecmp (argv[CtypeArg], "CAR")) { strcpy (graphmode.coords.ctype, "DEC--CAT"); goto got_ctype; }
    162164    gprint (GP_ERR, "ERROR: invalid projection type %s\n", argv[CtypeArg]);
    163     gprint (GP_ERR, "allowed values: TAN, SIN, ARC, STG, ZEA, AIT, GLS, PAR, MOL\n");
     165    gprint (GP_ERR, "allowed values: TAN, SIN, ARC, STG, ZEA, AIT, GLS, PAR, MOL, LIN, CAR\n");
    164166    return FALSE;
    165167  }
Note: See TracChangeset for help on using the changeset viewer.