**
** ROUTINE: taSegmentResolve
**
**
** Resolve all objects within a segment, where the segment geometry
** is fully specified by the input TA_GEOMETRY structure.
** databases. All objects whose flag OK_RUN is set and which lie within
** the range of mu are considered to belong to the segment; if there are
** abutting PRIMARY segments along the mu border then we must
** guarantee a unique and complete set of resolved obects in the border
** region (of radius matchRadius) where the segments meet.
** This is done by removing from the segment any object with a matching
** object which belongs to an abutting PRIMARY segment, and adding to the
** segment any object outside the mu limits but within
** matchRadius which has a matching object an
** abutting PRIMARY segments but which does not
** belong to that segment. All objects belonging to a segment have
** their RESOLVED flag set, as well as the PSEGMENT flag if the segment
** is a PRIMARY segment. Additionally, if
** the object belongs to the first field in the segment, then its
** FIRST_FIELD flag is set to distinguish objects lying in
** fields shared by two segments.
**
** If the object lies
** within the primary range of great circle latitude for the scanline,
** then its OK_SCANLINE flag is set. A resolution process similar to
** that performed in mu is also done here to guarantee a set of unique
** and complete detections in the border region between PRIMARY segments
** in adjacent scanlines. For objects whose OK_SCANLINE flag is set,
** if its lies within the primary range of survey latitude for the
** stripe, then its OK_STRIPE flag is set. For objects whose OK_STRIPE
** flag is set, if it lies within the primary range of survey longitude
** for the stripe then its PRIMARY flag is set. Again, a resolution
** process
** is performed to guarantee a set of unique and complete detections
** around the border region between stripes (no latitude cut is applied
** to stripes in the southern survey). If the PRIMARY flag is not
** set, then the SECONDARY flag is (the PRIMARY or SECONDARY flag is set
** for objects of class STAR, GALAXY, UNK, and SKY only).
**
** For PRIMARY segments, it is assumed that there are no overlapping
** PRIMARY segments within the same scanline already defined in the
** database.
**
** It is assumed that segments within a given chunk are processed in
** order of increasing scanline number; thus, we need to check
** matching objects in adjacent segments only along the lower nu
** border, and not along the upper nu border. However, we want to
** know if this object could match an object along the upper nu border
** or either mu border to facilitate processing an entire chunk.
** To accomodate that, the routine returns 1 if the object is matchable
** (OK_RUN, STAR || GALAXY || UNKNOWN), else 0.
**
**
**
**============================================================================
*/
int
taSegmentResolve(TA_CALIB_OBJ *obj /* MOD: calibrated object */,
int firstField, /* IN: 1=object in first field
* in segments, 0=not */
const TA_GEOMETRY *geometry, /* IN: survey geometry */
const TA_BORDER_OBJ *muBorder[2],
/* IN: objects from previous/next segment in border region
* with this field */
int nMuBorder[2], /* IN: number of border objects in
* previous/next segment */
const TA_BORDER_OBJ *nuBorder, /* IN: objects from previous
* scanline in border region
* with this field */
int nNuBorder, /* IN: number of border objects in
* previous scanline */
const TA_BORDER_OBJ *etaBorder[2],
/* IN: objects from previous/next chunk in border region
* with this chunk */
int nEtaBorder[2] /* IN: number of border objects in
* previous/next chunk */)
{
int i, sgu, inSurvey, ok, match;
double mu, nu, matchRadius;
/* Skip objects which have already been resolved */
shAssert(! (obj->status & AR_OBJECT_STATUS_RESOLVED));
/* Skip objects which are not OK_RUN. */
if (! (obj->status & AR_OBJECT_STATUS_OK_RUN)) return 0;
/* Only objects of class STAR, GALAXY, UNK, and SKY will be classed
* as either SECONDARY or PRIMARY. Only objects of class STAR,
* GALAXY, or UNK are matched.
*/
sgu = (obj->objc_type == AR_OBJECT_TYPE_STAR ||
obj->objc_type == AR_OBJECT_TYPE_GALAXY ||
obj->objc_type == AR_OBJECT_TYPE_UNK);
inSurvey = (sgu || obj->objc_type == AR_OBJECT_TYPE_SKY);
/* Get great circle coordinates. */
atEqToGC(obj->ra, obj->dec, &mu, &nu, geometry->stripeNode,
geometry->stripeIncl);
if (mu < geometry->muLimits[0] - 30.) {
mu += 360.0;
}
/* Does this fall within the segment's limits in mu?
* If it abuts another primary segment, then we need to guarantee
* complete and unique detections within a region of radius
* "matchRadius" around the border between the two. Only objects of class
* STAR, GALAXY, or UNK are matched, so we need check borders for those.
*/
matchRadius = geometry->matchRadius * at_asec2Deg;
ok = (mu >= geometry->muLimits[0] && mu < geometry->muLimits[1]);
if (sgu)
for (i = 0; i < 2; i++)
if (nMuBorder[i] > 0)
if (mu > geometry->muLimits[i] - matchRadius &&
mu < geometry->muLimits[i] + matchRadius)
{
shAssert(muBorder[i] != NULL);
match = taBorderMatch(obj, muBorder[i], nMuBorder[i], 0,
matchRadius, NULL);
if (match >= 0)
{
if (muBorder[i][match].status & AR_OBJECT_STATUS_PSEGMENT)
ok = 0;
else
ok = 1;
}
break;
}
/* If outside mu range, then not in the segment. */
if (! ok) return 0;
/* This object belongs to this segment. Flag accordingly. */
obj->status |= (AR_OBJECT_STATUS_RESOLVED | AR_OBJECT_STATUS_PSEGMENT);
if (firstField) obj->status |= AR_OBJECT_STATUS_FIRST_FIELD;
/* Is this object in primary nu range for this scanline?
* Within a given chunk, we require that segments be processed in order
* of increasing scanline number, so we need check only for matches along
* the lower nu border.
*/
ok = (nu >= geometry->nuLimits[0] && nu < geometry->nuLimits[1]);
if (sgu && nNuBorder > 0)
if (nu > geometry->nuLimits[0] - matchRadius &&
nu < geometry->nuLimits[0] + matchRadius)
{
match = taBorderMatch(obj, nuBorder, nNuBorder, 0, matchRadius, NULL);
if (match >= 0)
{
if (nuBorder[match].status & AR_OBJECT_STATUS_OK_SCANLINE)
ok = 0;
else
ok = 1;
}
}
if (! ok)
{
if (inSurvey) obj->status |= AR_OBJECT_STATUS_SECONDARY;
return sgu;
}
obj->status |= AR_OBJECT_STATUS_OK_SCANLINE;
/* Is this object in primary eta range for this stave? */
if (geometry->stripe > 45)
ok = 1;
else
{
ok = (obj->eta >= geometry->etaLimits[0] &&
obj->eta < geometry->etaLimits[1]);
if (sgu)
for (i = 0; i < 2; i++)
if (nEtaBorder[i] > 0)
if (obj->eta > geometry->etaLimits[i] - matchRadius &&
obj->eta < geometry->etaLimits[i] + matchRadius)
{
match = taBorderMatch(obj, etaBorder[i], nEtaBorder[i], 0,
matchRadius, NULL);
if (match >= 0)
{
if (etaBorder[i][match].status&AR_OBJECT_STATUS_OK_STRIPE)
ok = 0;
else
ok = 1;
}
break;
}
}
if (! ok)
{
if (inSurvey) obj->status |= AR_OBJECT_STATUS_SECONDARY;
return sgu;
}
obj->status |= AR_OBJECT_STATUS_OK_STRIPE;
/* Is this object in the lambda range for the stripe? */
if (inSurvey)
{
if (obj->lambda >= geometry->lambdaLimits[0] &&
obj->lambda < geometry->lambdaLimits[1])
obj->status |= AR_OBJECT_STATUS_PRIMARY;
else
obj->status |= AR_OBJECT_STATUS_SECONDARY;
}
return sgu;
}