Homec4science
Diffusion RADIANCE Photon Map (pmapmerge) Browse

RadiancePmap/pmaproi.c
2623f4232550pmapmerge

Last Change
Raw File

pmaproi.c
View Options

#ifndef lint
static const char RCSid[] = "$Id$";
#endif
/*
======================================================================
Regions of interest used by mkpmap to constrain location of photons
Roland Schregle (roland.schregle@{hslu.ch, gmail.com}
(c) Fraunhofer Institute for Solar Energy Systems,
supported by the German Research Foundation
(DFG LU-204/10-2, "Fassadenintegrierte Regelsysteme" (FARESYS))
(c) Lucerne University of Applied Sciences and Arts,
supported by the Swiss National Science Foundation
(SNSF #147053, "Daylight Redirecting Components",
SNSF #179067, "Light Fields for Spatio-Temporal Glare Assessment")
(c) Tokyo University of Science,
supported by the JSPS Grants-in-Aid for Scientific Research
(KAKENHI JP19KK0115, "Three-Dimensional Light Flow")
======================================================================
$Id$
*/
#include "pmaproi.h"
#include "otspecial.h"
#include "otypes.h"
/* Rectangular / spherical regions of interest */
unsigned pmapNumROI = 0;
PhotonMapROI *pmapROI = NULL;
/* Polyhedral regions of interest by modifier */
char *pmapROImodList [MAXSET + 1] = {NULL};
OBJECT polyROIset [MAXSET + 1] = {0};
void getPolyROIs (char **roiModList)
/* Find geometry comprising polyhedral regions of interest from modifiers in
* roiModList */
{
OBJECT i;
OBJREC *obj, *mat;
char **lp;
unsigned lastNumPolys = 0;
/* Init poly ROI set */
polyROIset [0] = 0;
if (!roiModList [0])
return;
for (lp = roiModList; *lp; lp++) {
lastNumPolys = polyROIset [0];
for (i = 0; i < nobjects; i++) {
obj = objptr(i);
mat = findmaterial(obj);
/* Check if object is geometry and its modifier is in list */
if (issurface(obj -> otype) && mat && !strcmp(mat -> oname, *lp)) {
if (polyROIset [0] >= PMAP_MAXROI)
error(USER, "too many ROI objects");
insertelem(polyROIset, i);
}
}
if (polyROIset [0] == lastNumPolys) {
/* No geometry found for this modifier */
sprintf(errmsg, "no ROI geometry found for modifier %s", *lp);
error(USER, errmsg);
}
}
if (!polyROIset [0])
/* Made redundant by check above? */
error(USER, "no ROI modifiers found");
}
int photonInROI (const RAY *photonRay)
/* Returns nonzero if photon defined by photonRay lies in any defined region
* of interest */
{
unsigned inROI = 0, i;
if (pmapNumROI && pmapROI) {
/* Check rectangular / spherical ROIs */
FVECT photonDist;
for (i = 0; i < pmapNumROI; i++) {
VSUB(photonDist, photonRay -> rop, pmapROI [i].pos);
/* NOTE: Size of spherical ROI is _squared_ */
inROI = (PMAP_ROI_ISSPHERE(pmapROI + i)
? DOT(photonDist, photonDist) <= pmapROI [i].siz [0]
: fabs(photonDist [0]) <= pmapROI [i].siz [0] &&
fabs(photonDist [1]) <= pmapROI [i].siz [1] &&
fabs(photonDist [2]) <= pmapROI [i].siz [2]
);
if (inROI)
/* Bail out early, skip remaining ROIs */
return 1;
}
}
if (polyROIset [0]) {
/* Check polyhedral ROIs; send ray and count number of crossed
* polygons in ROI set */
RAY testRay;
/* Init inclusion test ray from photon origin and incident direction */
rayorigin(&testRay, PRIMARY, NULL, NULL);
/* Offset origin along reversed photon direction in case we have a
* surface-bound (i.e. global, caustic) photon, thereby avoiding
* potential boundary issues. */
VSUM(testRay.rorg, photonRay -> rop, photonRay -> rdir, -2 * FTINY);
VCOPY(testRay.rdir, photonRay -> rdir);
testRay.rmax = FHUGE;
/* Follow inclusion test ray until it leaves the scene */
while (localhit(&testRay, &thescene)) {
/* Toggle odd/even number of crossings if intersected object is in
ROI set (note this works since inset() returns 0 or 1!) */
inROI ^= inset(polyROIset, testRay.robj);
/* Update ray origin with intersection for next iteration */
VCOPY(testRay.rorg, testRay.rop);
}
return inROI;
}
/* No ROIs defined; accept photon unconditionally */
return 1;
}
c4science · Help