Add Reflection dataFilter

src/reflectionReconstruction/dataFilter.cpp

@@ -0,0 +1,26 @@
+#include "dataFilter.h"
+#include "Function1D.h"
+#include "Function2D.h"
+#include "Function3D.h"
+
+#define PI 3.141592653589793238462
+using namespace Aurora;
+namespace Recon {
+Matrix filterReflectionData(const Matrix &aVReceiverPositionBlock,
+                            const Matrix &aVSenderPositionBlock,
+                            const Matrix &aVSenderNormalBlock,
+                            int aConstrictReflectionAngles,
+                            double aAngleLowerLimit, double aAngleUpperLimit) {
+    auto reflectData = ones(1, aVReceiverPositionBlock.getDimSize(1));
+    if (aConstrictReflectionAngles == 1) {
+        auto positionDifs = aVReceiverPositionBlock - aVSenderPositionBlock;
+        auto inbetweenAngle = acos(dot(aVSenderNormalBlock, positionDifs) /
+                                vecnorm(positionDifs, Norm2, 1));
+        auto flags = Aurora::auroraNot((inbetweenAngle < (PI * aAngleLowerLimit / 180)) +
+                    (inbetweenAngle > (PI * aAngleUpperLimit / 180)));
+        reflectData = reflectData * flags;
+    }
+    return reflectData;
+}
+
+} // namespace Recon

src/reflectionReconstruction/dataFilter.h

@@ -0,0 +1,15 @@
+#ifndef __REF_DATAFILTER_H__
+#define __REF_DATAFILTER_H__
+
+#include "Matrix.h"
+namespace Recon {
+    Aurora::Matrix filterReflectionData(const Aurora::Matrix &aVReceiverPositionBlock,
+                                    const Aurora::Matrix &aVSenderPositionBlock,
+                                    const Aurora::Matrix &aVSenderNormalBlock,
+                                    int aConstrictReflectionAngles,
+                                    double aAngleLowerLimit,
+                                    double aAngleUpperLimit);
+}
+
+
+#endif // __DATAFILTER_H__