# HG changeset patch # User Alain Mazy # Date 1623942036 -7200 # Node ID d7ac540e0149063742ebcf045df888f3aa39393e # Parent 80e029f8ec91e9a78f40270de65b574a7ca232a9# Parent f6eaf617d8e8ab6ceb8b3db10a3860c621225717 merge diff -r f6eaf617d8e8 -r d7ac540e0149 OrthancStone/Sources/Toolbox/DicomStructureSet.cpp --- a/OrthancStone/Sources/Toolbox/DicomStructureSet.cpp Thu Jun 17 15:59:26 2021 +0200 +++ b/OrthancStone/Sources/Toolbox/DicomStructureSet.cpp Thu Jun 17 17:00:36 2021 +0200 @@ -342,8 +342,8 @@ // if prev* and cur* are on opposite sides of y, this means that the // segment intersects the plane. - if ((prevY < y && curY > y) || - (prevY > y && curY < y)) + if ((prevY <= y && curY >= y) || + (prevY >= y && curY <= y)) { double p = (curX * prevY - curY * prevX + y * (prevX - curX)) / (prevY - curY); xmin = std::min(xmin, p); @@ -413,8 +413,8 @@ double curX, curY; geometry_.ProjectPoint2(curX, curY, points_[i]); - if ((prevX < x && curX > x) || - (prevX > x && curX < x)) + if ((prevX <= x && curX >= x) || + (prevX >= x && curX <= x)) { double p = (curX * prevY - curY * prevX + x * (curY - prevY)) / (curX - prevX); ymin = std::min(ymin, p);