OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

#3876 Add reference to Shoelace formula for counter clockwise polygon test.

Browse Source
This commit is contained in:
Gaute Lindkvist 2018-12-21 09:44:11 +01:00
parent 9fdec77c94
commit b1b278fffc
1 changed files with 6 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
Fwk/AppFwk/CommonCode/cafContourLines.cpp
View File

@ -281,14 +281,15 @@ std::vector<caf::ContourLines::ClosedPolygons> caf::ContourLines::create(const s
closedPolygonDeque.push_back(closedPolygonDeque.front()); closedPolygonDeque.push_back(closedPolygonDeque.front());
} }
// Make sure it is counter clockwise // Make sure it is counter clockwise. Use Shoelace formula to calculate signed area.
double sum = 0.0; // https://en.wikipedia.org/wiki/Shoelace_formula
double signedArea = 0.0;
for (size_t j = 0; j < closedPolygonDeque.size() - 1; ++j) for (size_t j = 0; j < closedPolygonDeque.size() - 1; ++j)
{ {
sum += (closedPolygonDeque[j + 1].x() - closedPolygonDeque[j].x()) * signedArea += (closedPolygonDeque[j + 1].x() - closedPolygonDeque[j].x()) *
(closedPolygonDeque[j + 1].y() + closedPolygonDeque[j].y()); (closedPolygonDeque[j + 1].y() + closedPolygonDeque[j].y());
} }
if (sum < 0.0) if (signedArea < 0.0)
{ {
closedPolygonsPerLevel[i].emplace_back(closedPolygonDeque.rbegin(), closedPolygonDeque.rend()); closedPolygonsPerLevel[i].emplace_back(closedPolygonDeque.rbegin(), closedPolygonDeque.rend());
} }