aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
Diffstat (limited to 'src')
-rw-r--r--src/argaze/AreaOfInterest/AOIFeatures.py31
1 files changed, 18 insertions, 13 deletions
diff --git a/src/argaze/AreaOfInterest/AOIFeatures.py b/src/argaze/AreaOfInterest/AOIFeatures.py
index c2714ca..25046ff 100644
--- a/src/argaze/AreaOfInterest/AOIFeatures.py
+++ b/src/argaze/AreaOfInterest/AOIFeatures.py
@@ -40,6 +40,10 @@ class AreaOfInterest(numpy.ndarray):
return numpy.array(points).view(AreaOfInterest)
+ def __init__(self, points: numpy.array = numpy.empty(0)):
+
+ self.__polygon = Polygon(points)
+
def __repr__(self):
"""String representation"""
@@ -81,8 +85,7 @@ class AreaOfInterest(numpy.ndarray):
# TODO: Use pygeos
N = 32
- points = [(math.cos(2 * math.pi / N * x) * radius + cx, math.sin(2 * math.pi / N * x) * radius + cy) for x
- in range(0, N + 1)]
+ points = [(math.cos(2 * math.pi / N * x) * radius + cx, math.sin(2 * math.pi / N * x) * radius + cy) for x in range(0, N + 1)]
return AreaOfInterest(points)
@@ -95,8 +98,7 @@ class AreaOfInterest(numpy.ndarray):
# TODO: Use pygeos
N = 32
- points = [(math.cos(2 * math.pi / N * x) * rx + cx, math.sin(2 * math.pi / N * x) * ry + cy) for x in
- range(0, N + 1)]
+ points = [(math.cos(2 * math.pi / N * x) * rx + cx, math.sin(2 * math.pi / N * x) * ry + cy) for x in range(0, N + 1)]
@property
def dimension(self) -> int:
@@ -135,7 +137,7 @@ class AreaOfInterest(numpy.ndarray):
@property
def area(self) -> float:
"""Area of the polygon defined by aoi's points."""
- return Polygon(self).area
+ return self.__polygon.area
@property
def bounding_box(self) -> numpy.array:
@@ -151,6 +153,12 @@ class AreaOfInterest(numpy.ndarray):
return numpy.array([(min_x, min_y), (max_x, min_y), (max_x, max_y), (min_x, max_y)])
+ @property
+ def is_valid(self) -> bool:
+ """Is area polygon valid?"""
+
+ return self.__polygon.is_valid
+
def clockwise(self) -> Self:
"""Get area points in clockwise order.
!!! warning
@@ -232,9 +240,7 @@ class AreaOfInterest(numpy.ndarray):
assert (self.dimension == 2)
- self_polygon = Polygon(self)
-
- if not self_polygon.is_valid:
+ if not self.__polygon.is_valid:
logging.warning('AreaOfInterest.circle_intersection: AOI polygon is not valid.')
@@ -244,14 +250,13 @@ class AreaOfInterest(numpy.ndarray):
args_circle = Point(center).buffer(radius)
- if self_polygon.intersects(args_circle):
+ if self.__polygon.intersects(args_circle):
- intersection = self_polygon.intersection(args_circle)
+ intersection = self.__polygon.intersection(args_circle)
- intersection_array = numpy.array([list(xy) for xy in intersection.exterior.coords[:]]).astype(
- numpy.float32).view(AreaOfInterest)
+ intersection_array = numpy.array([list(xy) for xy in intersection.exterior.coords[:]]).astype(numpy.float32).view(AreaOfInterest)
- return intersection_array, intersection.area / self_polygon.area, intersection.area / args_circle.area
+ return intersection_array, intersection.area / self.__polygon.area, intersection.area / args_circle.area
else: