CellProfiler-nightly 5.0.0.dev318__py3-none-any.whl → 5.0.0.dev324__py3-none-any.whl

cellprofiler/_version.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '5.0.0.dev318'
-__version_tuple__ = version_tuple = (5, 0, 0, 'dev318')
+__version__ = version = '5.0.0.dev324'
+__version_tuple__ = version_tuple = (5, 0, 0, 'dev324')
 
 __commit_id__ = commit_id = None

cellprofiler/modules/crop.py

@@ -38,14 +38,13 @@ help for more information on saving cropping shapes.
 """
 
 import logging
-import math
-
 import centrosome.filter
 import matplotlib.axes
 import matplotlib.cm
 import matplotlib.figure
 import matplotlib.patches
 import numpy
+from cellprofiler_core.constants.measurement import GROUP_INDEX
 from cellprofiler_core.image import Image
 from cellprofiler_core.module import Module
 from cellprofiler_core.preferences import get_primary_outline_color
@@ -57,39 +56,11 @@ from cellprofiler_core.setting.subscriber import ImageSubscriber
 from cellprofiler_core.setting.subscriber import LabelSubscriber
 from cellprofiler_core.setting.text import CropImageName
 from cellprofiler_core.setting.text import Integer
-from cellprofiler_core.utilities.image import crop_image
-
+from cellprofiler_library.functions.image_processing import get_ellipse_cropping, get_rectangle_cropping
+from cellprofiler_library.modules._crop import crop, get_measurements
+from cellprofiler_library.opts.crop import RemovalMethod, Measurement, Shape, CroppingMethod, CroppingPattern, Limits, Ellipse, Rectangle
 LOGGER = logging.getLogger(__name__)
 
-SH_RECTANGLE = "Rectangle"
-SH_ELLIPSE = "Ellipse"
-SH_IMAGE = "Image"
-SH_OBJECTS = "Objects"
-SH_CROPPING = "Previous cropping"
-CM_COORDINATES = "Coordinates"
-CM_MOUSE = "Mouse"
-IO_INDIVIDUALLY = "Every"
-IO_FIRST = "First"
-RM_NO = "No"
-RM_EDGES = "Edges"
-RM_ALL = "All"
-
-# Doesn't seem to like importing defs from cellprofiler.gui.moduleview so define here
-ABSOLUTE = "Absolute"
-FROM_EDGE = "From edge"
-
-EL_XCENTER = "xcenter"
-EL_YCENTER = "ycenter"
-EL_XRADIUS = "xradius"
-EL_YRADIUS = "yradius"
-
-RE_LEFT = "left"
-RE_TOP = "top"
-RE_RIGHT = "right"
-RE_BOTTOM = "bottom"
-
-FF_AREA_RETAINED = "Crop_AreaRetainedAfterCropping_%s"
-FF_ORIGINAL_AREA = "Crop_OriginalImageArea_%s"
 
 OFF_IMAGE_NAME = 0
 OFF_CROPPED_IMAGE_NAME = 1
@@ -130,8 +101,8 @@ class Crop(Module):
 
         self.shape = Choice(
             text="Select the cropping shape",
-            choices=[SH_RECTANGLE, SH_ELLIPSE, SH_IMAGE, SH_OBJECTS, SH_CROPPING],
-            value=SH_RECTANGLE,
+            choices=[Shape.RECTANGLE.value, Shape.ELLIPSE.value, Shape.IMAGE.value, Shape.OBJECTS.value, Shape.CROPPING.value],
+            value=Shape.RECTANGLE.value,
             doc="""\
 Choose the shape into which you would like to crop:
 
@@ -161,19 +132,19 @@ Choose the shape into which you would like to crop:
    cropping that was used to generate whichever image you choose.
 """.format(
                 **{
-                    "SH_RECTANGLE": SH_RECTANGLE,
-                    "SH_ELLIPSE": SH_ELLIPSE,
-                    "SH_IMAGE": SH_IMAGE,
-                    "SH_OBJECTS": SH_OBJECTS,
-                    "SH_CROPPING": SH_CROPPING,
+                    "SH_RECTANGLE": Shape.RECTANGLE.value,
+                    "SH_ELLIPSE": Shape.ELLIPSE.value,
+                    "SH_IMAGE": Shape.IMAGE.value,
+                    "SH_OBJECTS": Shape.OBJECTS.value,
+                    "SH_CROPPING": Shape.CROPPING.value,
                 }
             ),
         )
 
         self.crop_method = Choice(
             text="Select the cropping method",
-            choices=[CM_COORDINATES, CM_MOUSE],
-            value=CM_COORDINATES,
+            choices=[CroppingMethod.COORDINATES.value, CroppingMethod.MOUSE.value],
+            value=CroppingMethod.COORDINATES.value,
             doc="""\
 Choose whether you would like to crop by typing in pixel coordinates or
 clicking with the mouse.
@@ -190,18 +161,18 @@ clicking with the mouse.
    are in the interior of the region you wish to retain.
 """.format(
                 **{
-                    "CM_COORDINATES": CM_COORDINATES,
-                    "SH_ELLIPSE": SH_ELLIPSE,
-                    "SH_RECTANGLE": SH_RECTANGLE,
-                    "CM_MOUSE": CM_MOUSE,
+                    "CM_COORDINATES": CroppingMethod.COORDINATES.value,
+                    "SH_ELLIPSE": Shape.ELLIPSE.value,
+                    "SH_RECTANGLE": Shape.RECTANGLE.value,
+                    "CM_MOUSE": CroppingMethod.MOUSE.value,
                 }
             ),
         )
 
         self.individual_or_once = Choice(
             text="Apply which cycle's cropping pattern?",
-            choices=[IO_INDIVIDUALLY, IO_FIRST],
-            value=IO_INDIVIDUALLY,
+            choices=[CroppingPattern.INDIVIDUALLY.value, CroppingPattern.FIRST.value],
+            value=CroppingPattern.INDIVIDUALLY.value,
             doc="""\
 Specify how a given cropping pattern should be applied to other image cycles:
 
@@ -210,7 +181,7 @@ Specify how a given cropping pattern should be applied to other image cycles:
    intended to function as a template in some fashion.
 -  *{IO_INDIVIDUALLY}:* Every image cycle is cropped individually.
 """.format(
-                **{"IO_FIRST": IO_FIRST, "IO_INDIVIDUALLY": IO_INDIVIDUALLY}
+                **{"IO_FIRST": CroppingPattern.FIRST.value, "IO_INDIVIDUALLY": CroppingPattern.INDIVIDUALLY.value}
             ),
         )
 
@@ -232,9 +203,9 @@ Specify the left and right positions for the bounding rectangle by selecting one
    of the image’s original size.
 """.format(
                 **{
-                    "SH_RECTANGLE": SH_RECTANGLE,
-                    "ABSOLUTE": ABSOLUTE,
-                    "FROM_EDGE": FROM_EDGE,
+                    "SH_RECTANGLE": Shape.RECTANGLE.value,
+                    "ABSOLUTE": Limits.ABSOLUTE.value,
+                    "FROM_EDGE": Limits.FROM_EDGE.value,
                 }
             ),
         )
@@ -255,9 +226,9 @@ Specify the top and bottom positions for the bounding rectangle by selecting one
    from the edges of your images irrespective of their size.
 """.format(
                 **{
-                    "SH_RECTANGLE": SH_RECTANGLE,
-                    "ABSOLUTE": ABSOLUTE,
-                    "FROM_EDGE": FROM_EDGE,
+                    "SH_RECTANGLE": Shape.RECTANGLE.value,
+                    "ABSOLUTE": Limits.ABSOLUTE.value,
+                    "FROM_EDGE": Limits.FROM_EDGE.value,
                 }
             ),
         )
@@ -270,7 +241,7 @@ Specify the top and bottom positions for the bounding rectangle by selecting one
 
 Specify the center pixel position of the ellipse.
 """.format(
-                **{"SH_ELLIPSE": SH_ELLIPSE}
+                **{"SH_ELLIPSE": Shape.ELLIPSE.value}
             ),
         )
 
@@ -282,7 +253,7 @@ Specify the center pixel position of the ellipse.
 
 Specify the radius of the ellipse in the X direction.
 """.format(
-                **{"SH_ELLIPSE": SH_ELLIPSE}
+                **{"SH_ELLIPSE": Shape.ELLIPSE.value}
             ),
         )
 
@@ -294,7 +265,7 @@ Specify the radius of the ellipse in the X direction.
 
 Specify the radius of the ellipse in the Y direction.
 """.format(
-                **{"SH_ELLIPSE": SH_ELLIPSE}
+                **{"SH_ELLIPSE": Shape.ELLIPSE.value}
             ),
         )
 
@@ -306,7 +277,7 @@ Specify the radius of the ellipse in the Y direction.
 
 Select the image to be use as a cropping mask.
 """.format(
-                **{"SH_IMAGE": SH_IMAGE}
+                **{"SH_IMAGE": Shape.IMAGE.value}
             ),
         )
 
@@ -318,7 +289,7 @@ Select the image to be use as a cropping mask.
 
 Select the image associated with the cropping mask that you want to use.
 """.format(
-                **{"SH_CROPPING": SH_CROPPING}
+                **{"SH_CROPPING": Shape.CROPPING.value}
             ),
         )
 
@@ -330,14 +301,14 @@ Select the image associated with the cropping mask that you want to use.
 
 Select the objects that are to be used as a cropping mask.
 """.format(
-                **{"SH_OBJECTS": SH_OBJECTS}
+                **{"SH_OBJECTS": Shape.OBJECTS.value}
             ),
         )
 
         self.remove_rows_and_columns = Choice(
             text="Remove empty rows and columns?",
-            choices=[RM_NO, RM_EDGES, RM_ALL],
-            value=RM_ALL,
+            choices=[RemovalMethod.NO.value, RemovalMethod.EDGES.value, RemovalMethod.ALL.value],
+            value=RemovalMethod.ALL.value,
             doc="""\
 Use this option to choose whether to remove rows and columns that lack
 objects:
@@ -349,7 +320,7 @@ objects:
 -  *{RM_ALL}:* Remove any row or column of all-blank pixels, even
    from the internal portion of the image.
 """.format(
-                **{"RM_NO": RM_NO, "RM_EDGES": RM_EDGES, "RM_ALL": RM_ALL}
+                **{"RM_NO": RemovalMethod.NO.value, "RM_EDGES": RemovalMethod.EDGES.value, "RM_ALL": RemovalMethod.ALL.value}
             ),
         )
 
@@ -373,22 +344,22 @@ objects:
 
     def visible_settings(self):
         result = [self.image_name, self.cropped_image_name, self.shape]
-        if self.shape.value in (SH_RECTANGLE, SH_ELLIPSE):
+        if self.shape.value in (Shape.RECTANGLE, Shape.ELLIPSE):
             result += [self.crop_method, self.individual_or_once]
-            if self.crop_method == CM_COORDINATES:
-                if self.shape == SH_RECTANGLE:
+            if self.crop_method.value == CroppingMethod.COORDINATES:
+                if self.shape.value == Shape.RECTANGLE:
                     result += [self.horizontal_limits, self.vertical_limits]
-                elif self.shape == SH_ELLIPSE:
+                elif self.shape.value == Shape.ELLIPSE:
                     result += [
                         self.ellipse_center,
                         self.ellipse_x_radius,
                         self.ellipse_y_radius,
                     ]
-        elif self.shape == SH_IMAGE:
+        elif self.shape.value == Shape.IMAGE:
             result += [self.image_mask_source]
-        elif self.shape == SH_CROPPING:
+        elif self.shape.value == Shape.CROPPING:
             result.append(self.cropping_mask_source)
-        elif self.shape == SH_OBJECTS:
+        elif self.shape.value == Shape.OBJECTS:
             result.append(self.objects_source)
         else:
             raise NotImplementedError("Unimplemented shape type: %s" % self.shape.value)
@@ -397,25 +368,25 @@ objects:
 
     def run(self, workspace):
         first_image_set = (
-            workspace.measurements.get_current_image_measurement("Group_Index") == 1
+            workspace.measurements.get_current_image_measurement(GROUP_INDEX) == 1
         )
         image_set_list = workspace.image_set_list
-        d = self.get_dictionary(image_set_list)
+        cache_dict = self.get_dictionary(image_set_list)
         orig_image = workspace.image_set.get_image(self.image_name.value)
         recalculate_flag = (
-            self.shape not in (SH_ELLIPSE, SH_RECTANGLE)
-            or self.individual_or_once == IO_INDIVIDUALLY
+            self.shape.value not in (Shape.ELLIPSE, Shape.RECTANGLE)
+            or self.individual_or_once.value == CroppingPattern.INDIVIDUALLY
             or first_image_set
             or workspace.pipeline.test_mode
         )
-        save_flag = self.individual_or_once == IO_FIRST and first_image_set
+        save_flag = self.individual_or_once.value == CroppingPattern.FIRST and first_image_set
         if not recalculate_flag:
-            if d[D_FIRST_CROPPING].shape != orig_image.pixel_data.shape[:2]:
+            if cache_dict[D_FIRST_CROPPING].shape != orig_image.pixel_data.shape[:2]:
                 recalculate_flag = True
                 LOGGER.warning(
                     """Image, "%s", size changed from %s to %s during cycle %d, recalculating""",
                     self.image_name.value,
-                    str(d[D_FIRST_CROPPING].shape),
+                    str(cache_dict[D_FIRST_CROPPING].shape),
                     str(orig_image.pixel_data.shape[:2]),
                     workspace.image_set.image_number,
                 )
@@ -423,59 +394,54 @@ objects:
         cropping = None
         masking_objects = None
         if not recalculate_flag:
-            cropping = d[D_FIRST_CROPPING]
-            mask = d[D_FIRST_CROPPING_MASK]
-        elif self.shape == SH_CROPPING:
+            cropping = cache_dict[D_FIRST_CROPPING]
+            mask = cache_dict[D_FIRST_CROPPING_MASK]
+        elif self.shape.value == Shape.CROPPING:
             cropping_image = workspace.image_set.get_image(
                 self.cropping_mask_source.value
             )
             cropping = cropping_image.crop_mask
-        elif self.shape == SH_IMAGE:
+        elif self.shape.value == Shape.IMAGE:
             source_image = workspace.image_set.get_image(
                 self.image_mask_source.value
             ).pixel_data
 
             cropping = source_image > 0
-        elif self.shape == SH_OBJECTS:
+        elif self.shape.value == Shape.OBJECTS:
             masking_objects = workspace.get_objects(self.objects_source.value)
             cropping = masking_objects.segmented > 0
-        elif self.crop_method == CM_MOUSE:
+        elif self.crop_method.value == CroppingMethod.MOUSE:
             cropping = self.ui_crop(workspace, orig_image)
-        elif self.shape == SH_ELLIPSE:
-            cropping = self.get_ellipse_cropping(workspace, orig_image)
-        elif self.shape == SH_RECTANGLE:
-            cropping = self.get_rectangle_cropping(workspace, orig_image)
-        if self.remove_rows_and_columns == RM_NO:
-            cropped_pixel_data = orig_image.pixel_data.copy()
-            if cropped_pixel_data.ndim == 3:
-                cropped_pixel_data[~cropping, :] = 0
-            else:
-                cropped_pixel_data[numpy.logical_not(cropping)] = 0
-            if mask is None:
-                mask = cropping
-            if orig_image.has_mask:
-                image_mask = mask & orig_image.mask
-            else:
-                image_mask = mask
-        else:
-            internal_cropping = self.remove_rows_and_columns == RM_ALL
-            cropped_pixel_data = crop_image(
-                orig_image.pixel_data, cropping, internal_cropping
+        elif self.shape.value == Shape.ELLIPSE:
+            cache_dict[Shape.ELLIPSE] = {
+                Ellipse.XCENTER: self.ellipse_center.x,   
+                Ellipse.YCENTER: self.ellipse_center.y,
+                Ellipse.XRADIUS: self.ellipse_x_radius.value,
+                Ellipse.YRADIUS: self.ellipse_y_radius.value,
+            }
+
+            cropping = get_ellipse_cropping(
+                orig_image.pixel_data,
+                (self.ellipse_center.x, self.ellipse_center.y),
+                (self.ellipse_x_radius.value, self.ellipse_y_radius.value)
             )
-            if mask is None:
-                mask = crop_image(cropping, cropping, internal_cropping)
-            if orig_image.has_mask:
-                image_mask = (
-                    crop_image(orig_image.mask, cropping, internal_cropping) & mask
-                )
-            else:
-                image_mask = mask
-
-            if cropped_pixel_data.ndim == 3:
-                cropped_pixel_data[~mask, :] = 0
-            else:
-                cropped_pixel_data[~mask] = 0
-        if self.shape == SH_OBJECTS:
+
+        elif self.shape.value == Shape.RECTANGLE:
+            h_min = self.horizontal_limits.min if not self.horizontal_limits.unbounded_min else None
+            h_max = self.horizontal_limits.max if not self.horizontal_limits.unbounded_max else None
+            v_min = self.vertical_limits.min if not self.vertical_limits.unbounded_min else None
+            v_max = self.vertical_limits.max if not self.vertical_limits.unbounded_max else None
+
+            cropping = get_rectangle_cropping(orig_image.pixel_data, (h_min, h_max, v_min, v_max), validate_boundaries=True)
+        else:
+            raise NotImplementedError(f"Cropping shape {self.shape.value} or crop method {self.crop_method} not supported.")
+
+        assert(cropping is not None)
+        assert(cropping.dtype == bool)
+
+        cropped_pixel_data, mask, image_mask = crop(orig_image.pixel_data, cropping, mask, orig_image.mask, self.remove_rows_and_columns.value)
+
+        if self.shape.value == Shape.OBJECTS:
             # Special handling for objects - masked objects instead of
             # mask and crop mask
             output_image = Image(
@@ -501,8 +467,8 @@ objects:
             )
 
         if save_flag:
-            d[D_FIRST_CROPPING_MASK] = mask
-            d[D_FIRST_CROPPING] = cropping
+            cache_dict[D_FIRST_CROPPING_MASK] = mask
+            cache_dict[D_FIRST_CROPPING] = cropping
         #
         # Save the image / cropping / mask
         #
@@ -510,13 +476,10 @@ objects:
         #
         # Save the old and new image sizes
         #
-        original_image_area = numpy.product(orig_image.pixel_data.shape[:2])
-        area_retained_after_cropping = numpy.sum(cropping)
-        feature = FF_AREA_RETAINED % self.cropped_image_name.value
         m = workspace.measurements
-        m.add_measurement("Image", feature, numpy.array([area_retained_after_cropping]))
-        feature = FF_ORIGINAL_AREA % self.cropped_image_name.value
-        m.add_measurement("Image", feature, numpy.array([original_image_area]))
+        for measurement in get_measurements(cropping, orig_image.pixel_data, self.cropped_image_name.value):
+            m.add_measurement("Image", measurement[1], numpy.array([measurement[2]]))
+
 
     def display(self, workspace, figure):
         orig_image_pixel_data = workspace.display_data.orig_image_pixel_data
@@ -536,20 +499,28 @@ objects:
         """Return information on the measurements made during cropping"""
         return [
             ("Image", x % self.cropped_image_name.value, "integer",)
-            for x in (FF_AREA_RETAINED, FF_ORIGINAL_AREA)
+            for x in (Measurement.AREA_RETAINED, Measurement.ORIGINAL_AREA)
         ]
 
     def ui_crop(self, workspace, orig_image):
         """Crop into a rectangle or ellipse, guided by UI"""
-        d = self.get_dictionary(workspace.image_set_list)
-        if (self.shape.value not in d) or self.individual_or_once == IO_INDIVIDUALLY:
-            d[self.shape.value] = workspace.interaction_request(
-                self, d.get(self.shape.value, None), orig_image.pixel_data
+        cache_dict = self.get_dictionary(workspace.image_set_list)
+        if (self.shape.value not in cache_dict) or self.individual_or_once.value == CroppingPattern.INDIVIDUALLY:
+            cache_dict[self.shape.value] = workspace.interaction_request(
+                self, cache_dict.get(self.shape.value, None), orig_image.pixel_data
             )
-        if self.shape == SH_ELLIPSE:
-            return self.apply_ellipse_cropping(workspace, orig_image)
+        if self.shape.value == Shape.ELLIPSE:
+            center = cache_dict[Shape.ELLIPSE][Ellipse.XCENTER], cache_dict[Shape.ELLIPSE][Ellipse.YCENTER]
+            radius = cache_dict[Shape.ELLIPSE][Ellipse.XRADIUS], cache_dict[Shape.ELLIPSE][Ellipse.YRADIUS]   
+            return get_ellipse_cropping(orig_image.pixel_data, center, radius)
         else:
-            return self.apply_rectangle_cropping(workspace, orig_image)
+            bounding_box = (
+                int(numpy.round(cache_dict[Shape.RECTANGLE][Rectangle.LEFT])),
+                int(numpy.round(cache_dict[Shape.RECTANGLE][Rectangle.RIGHT])),
+                int(numpy.round(cache_dict[Shape.RECTANGLE][Rectangle.TOP])),
+                int(numpy.round(cache_dict[Shape.RECTANGLE][Rectangle.BOTTOM])),
+            )
+            return get_rectangle_cropping(orig_image.pixel_data, bounding_box, validate_boundaries=True)
 
     def handle_interaction(self, current_shape, orig_image):
         from matplotlib.backends.backend_wxagg import FigureCanvasWxAgg
@@ -788,31 +759,31 @@ objects:
             def handles(self):
                 return [self.center_handle, self.radius_handle]
 
-        if self.shape == SH_ELLIPSE:
+        if self.shape.value == Shape.ELLIPSE:
             if current_shape is None:
                 current_shape = {
-                    EL_XCENTER: pixel_data.shape[1] / 2,
-                    EL_YCENTER: pixel_data.shape[0] / 2,
-                    EL_XRADIUS: pixel_data.shape[1] / 2,
-                    EL_YRADIUS: pixel_data.shape[0] / 2,
+                    Ellipse.XCENTER: pixel_data.shape[1] / 2,
+                    Ellipse.YCENTER: pixel_data.shape[0] / 2,
+                    Ellipse.XRADIUS: pixel_data.shape[1] / 2,
+                    Ellipse.YRADIUS: pixel_data.shape[0] / 2,
                 }
             ellipse = current_shape
             shape = CropEllipse(
-                (ellipse[EL_XCENTER], ellipse[EL_YCENTER]),
-                (ellipse[EL_XRADIUS], ellipse[EL_YRADIUS]),
+                (ellipse[Ellipse.XCENTER], ellipse[Ellipse.YCENTER]),
+                (ellipse[Ellipse.XRADIUS], ellipse[Ellipse.YRADIUS]),
             )
         else:
             if current_shape is None:
                 current_shape = {
-                    RE_LEFT: pixel_data.shape[1] / 4,
-                    RE_TOP: pixel_data.shape[0] / 4,
-                    RE_RIGHT: pixel_data.shape[1] * 3 / 4,
-                    RE_BOTTOM: pixel_data.shape[0] * 3 / 4,
+                    Rectangle.LEFT: pixel_data.shape[1] / 4,
+                    Rectangle.TOP: pixel_data.shape[0] / 4,
+                    Rectangle.RIGHT: pixel_data.shape[1] * 3 / 4,
+                    Rectangle.BOTTOM: pixel_data.shape[0] * 3 / 4,
                 }
             rectangle = current_shape
             shape = CropRectangle(
-                (rectangle[RE_LEFT], rectangle[RE_TOP]),
-                (rectangle[RE_RIGHT], rectangle[RE_BOTTOM]),
+                (rectangle[Rectangle.LEFT], rectangle[Rectangle.TOP]),
+                (rectangle[Rectangle.RIGHT], rectangle[Rectangle.BOTTOM]),
             )
         for patch in shape.patches:
             axes.add_artist(patch)
@@ -843,91 +814,21 @@ objects:
                 raise ValueError("Cancelled by user")
         finally:
             dialog_box.Destroy()
-        if self.shape == SH_RECTANGLE:
+        if self.shape.value == Shape.RECTANGLE:
             return {
-                RE_LEFT: shape.left,
-                RE_TOP: shape.top,
-                RE_RIGHT: shape.right,
-                RE_BOTTOM: shape.bottom,
+                Rectangle.LEFT: shape.left,
+                Rectangle.TOP: shape.top,
+                Rectangle.RIGHT: shape.right,
+                Rectangle.BOTTOM: shape.bottom,
             }
         else:
             return {
-                EL_XCENTER: shape.center_x,
-                EL_YCENTER: shape.center_y,
-                EL_XRADIUS: shape.width / 2,
-                EL_YRADIUS: shape.height / 2,
+                Ellipse.XCENTER: shape.center_x,
+                Ellipse.YCENTER: shape.center_y,
+                Ellipse.XRADIUS: shape.width / 2,
+                Ellipse.YRADIUS: shape.height / 2,
             }
 
-    def get_ellipse_cropping(self, workspace, orig_image):
-        """Crop into an ellipse using user-specified coordinates"""
-        x_center = self.ellipse_center.x
-        y_center = self.ellipse_center.y
-        x_radius = self.ellipse_x_radius.value
-        y_radius = self.ellipse_y_radius.value
-        d = self.get_dictionary(workspace.image_set_list)
-        d[SH_ELLIPSE] = {
-            EL_XCENTER: x_center,
-            EL_YCENTER: y_center,
-            EL_XRADIUS: x_radius,
-            EL_YRADIUS: y_radius,
-        }
-        return self.apply_ellipse_cropping(workspace, orig_image)
-
-    def apply_ellipse_cropping(self, workspace, orig_image):
-        d = self.get_dictionary(workspace.image_set_list)
-        ellipse = d[SH_ELLIPSE]
-        x_center, y_center, x_radius, y_radius = [
-            ellipse[x] for x in (EL_XCENTER, EL_YCENTER, EL_XRADIUS, EL_YRADIUS)
-        ]
-        pixel_data = orig_image.pixel_data
-        x_max = pixel_data.shape[1]
-        y_max = pixel_data.shape[0]
-        if x_radius > y_radius:
-            dist_x = math.sqrt(x_radius ** 2 - y_radius ** 2)
-            dist_y = 0
-            major_radius = x_radius
-        else:
-            dist_x = 0
-            dist_y = math.sqrt(y_radius ** 2 - x_radius ** 2)
-            major_radius = y_radius
-
-        focus_1_x, focus_1_y = (x_center - dist_x, y_center - dist_y)
-        focus_2_x, focus_2_y = (x_center + dist_x, y_center + dist_y)
-        y, x = numpy.mgrid[0:y_max, 0:x_max]
-        d1 = numpy.sqrt((x - focus_1_x) ** 2 + (y - focus_1_y) ** 2)
-        d2 = numpy.sqrt((x - focus_2_x) ** 2 + (y - focus_2_y) ** 2)
-        cropping = d1 + d2 <= major_radius * 2
-        return cropping
-
-    def get_rectangle_cropping(self, workspace, orig_image):
-        """Crop into a rectangle using user-specified coordinates"""
-        cropping = numpy.ones(orig_image.pixel_data.shape[:2], bool)
-        if not self.horizontal_limits.unbounded_min:
-            cropping[:, : self.horizontal_limits.min] = False
-        if not self.horizontal_limits.unbounded_max:
-            cropping[:, self.horizontal_limits.max :] = False
-        if not self.vertical_limits.unbounded_min:
-            cropping[: self.vertical_limits.min, :] = False
-        if not self.vertical_limits.unbounded_max:
-            cropping[self.vertical_limits.max :, :] = False
-        return cropping
-
-    def apply_rectangle_cropping(self, workspace, orig_image):
-        cropping = numpy.ones(orig_image.pixel_data.shape[:2], bool)
-        d = self.get_dictionary(workspace.image_set_list)
-        r = d[SH_RECTANGLE]
-        left, top, right, bottom = [
-            int(numpy.round(r[x])) for x in (RE_LEFT, RE_TOP, RE_RIGHT, RE_BOTTOM)
-        ]
-        if left > 0:
-            cropping[:, :left] = False
-        if right < cropping.shape[1]:
-            cropping[:, right:] = False
-        if top > 0:
-            cropping[:top, :] = False
-        if bottom < cropping.shape[0]:
-            cropping[bottom:, :] = False
-        return cropping
 
     def upgrade_settings(self, setting_values, variable_revision_number, module_name):
         if variable_revision_number == 1:
@@ -940,12 +841,12 @@ objects:
             # minor - "Cropping" changed to "Previous cropping"
             setting_values = list(setting_values)
             if setting_values[OFF_SHAPE] == "Cropping":
-                setting_values[OFF_SHAPE] = SH_CROPPING
+                setting_values[OFF_SHAPE] = Shape.CROPPING
             #
             # Individually changed to "every"
             #
             if setting_values[OFF_INDIVIDUAL_OR_ONCE] == "Individually":
-                setting_values[OFF_INDIVIDUAL_OR_ONCE] = IO_INDIVIDUALLY
+                setting_values[OFF_INDIVIDUAL_OR_ONCE] = CroppingPattern.INDIVIDUALLY
 
             setting_values = setting_values[:10] + setting_values[11:]
 

cellprofiler/modules/erodeimage.py

@@ -19,9 +19,8 @@ YES          YES          NO
 from cellprofiler_core.module import ImageProcessing
 from cellprofiler_core.setting import StructuringElement
 
-import cellprofiler.utilities.morphology
 from cellprofiler.modules._help import HELP_FOR_STREL
-
+from cellprofiler_library.modules._erodeimage import erode_image
 
 class ErodeImage(ImageProcessing):
     category = "Advanced"
@@ -48,6 +47,5 @@ class ErodeImage(ImageProcessing):
         return __settings__ + [self.structuring_element]
 
     def run(self, workspace):
-        self.function = cellprofiler.utilities.morphology.erosion
-
+        self.function = erode_image
         super(ErodeImage, self).run(workspace)

{cellprofiler_nightly-5.0.0.dev318.dist-info → cellprofiler_nightly-5.0.0.dev324.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: CellProfiler-nightly
-Version: 5.0.0.dev318
+Version: 5.0.0.dev324
 Summary: CellProfiler is a free open-source software designed to enable biologists without training in computer vision or programming to quantitatively measure phenotypes from thousands of images automatically.
 Author: Anne Carpenter, Thouis (Ray) Jones, Lee Kamentsky, Vebjorn Ljosa, David Logan, Mark Bray, Madison Swain-Bowden, Allen Goodman, Claire McQuinn, Alice Lucas, Callum Tromans-Coia
 Author-email: Beth Cimini <bcimini@broadinstitute.org>, David Stirling <dstirling@glencoesoftware.com>, Nodar Gogoberidze <ngogober@broadinstitute.org>

{cellprofiler_nightly-5.0.0.dev318.dist-info → cellprofiler_nightly-5.0.0.dev324.dist-info}/RECORD

@@ -1,6 +1,6 @@
 cellprofiler/__init__.py,sha256=AL2XeOBhIeYkBRyDd0QRgJan7j0DKjT1GD-RdzKvMUY,46
 cellprofiler/__main__.py,sha256=uy78oz5c6NBGRwDZkZ2Gl4HfhbJZQH6K1n54Qfl2OFc,38075
-cellprofiler/_version.py,sha256=VKwI68zHVxr7qKLpGyfYE5ws8C3HO1sHftyVZcFWLQk,721
+cellprofiler/_version.py,sha256=X_dsIS438JBC3vSFB9YEnaGPIPrvCqTB7mibbwXlkv0,721
 cellprofiler/knime_bridge.py,sha256=T6Op-KO79oULx92nGXRQ6lHsEcTutx1Uep1L4ZOKJgc,27767
 cellprofiler/misc.py,sha256=yqv873lP_mquxxkKcLgE_ZU4Hrc1trtuQ-NXLK2qQVc,553
 cellprofiler/data/examples/ExampleFly/ExampleFly.cppipe,sha256=JGZK9IZuYlOHOI4hi6a3DK36IahY69cfeLEd7eJ_rO0,15409
@@ -290,7 +290,7 @@ cellprofiler/modules/convertobjectstoimage.py,sha256=zvn2nlnF3-i6g16F0um1yNT-Y03
 cellprofiler/modules/correctilluminationapply.py,sha256=eEcXliF2xNkT36Bo7xkNNh7l2YGbSqkqpss4O7RbwaM,15862
 cellprofiler/modules/correctilluminationcalculate.py,sha256=wJMCy3W9zK-xPI8PUvh25MLtE_JSI3QuVVzdLjuTwD4,54200
 cellprofiler/modules/createbatchfiles.py,sha256=a8R5PpbaQYqF8OIvWBuC0uh-yhjzBaU4Uni4m3MRLok,19512
-cellprofiler/modules/crop.py,sha256=Lg1cQCIzO7o8wcwEgRlsccdt8mgiiugowcNEe44vT0k,36039
+cellprofiler/modules/crop.py,sha256=gnAHZNjPCNr-RLc7xQjZYViRpmKCTfJjuZ8MwHRiT1Y,34177
 cellprofiler/modules/definegrid.py,sha256=tjYcrzHu6TV7rkh1l2T0pNqcjBLwCfz9BCqTPWhQr_A,47160
 cellprofiler/modules/dilateimage.py,sha256=xcdsDdgvuplD69pcyJTOmVTxXFlwDWRNxapoQZHTiVc,1382
 cellprofiler/modules/dilateobjects.py,sha256=lnipjKkd6A3hhY3G1ttU9ZlhsTZp3oVU-QnU7_ayWbE,1828
@@ -302,7 +302,7 @@ cellprofiler/modules/displayscatterplot.py,sha256=sgB5Uf03mtUFMy3X6GqNNImhlKBWfZ
 cellprofiler/modules/editobjectsmanually.py,sha256=fkTFMTlrTfIiHNhhPzYsfFvZBSBDEIJgM7pZ7gDrMNA,21527
 cellprofiler/modules/enhanceedges.py,sha256=bivBs2OMNYvNnJ7tYrs9yheJP9JGIkgJDQfB2u837HE,11398
 cellprofiler/modules/enhanceorsuppressfeatures.py,sha256=vc_LLANt2gEsntqog9JwbzCVaaBQEM3ayvvY7CSJCF0,24062
-cellprofiler/modules/erodeimage.py,sha256=vdE2BjYq_wN5pKqNKdgs5xAINAK6vbQjB6rP_uLTwaY,1410
+cellprofiler/modules/erodeimage.py,sha256=GGO6EjgA_6A0pjrsATxcfs1a_QmZ1rS5Lec7Kkb21hw,1402
 cellprofiler/modules/erodeobjects.py,sha256=cDQR8y7XqSt-ssp1nAQp3wXdQk5OE1uT7GrhTgjfdlI,4381
 cellprofiler/modules/expandorshrinkobjects.py,sha256=fiFf4q2BplR4u6Cym-jdlOOdOKc6oZVtEkZE-O-Ansg,13911
 cellprofiler/modules/exporttodatabase.py,sha256=Tb5Hf329m96zUNaxPStQoTkr_MV62HGNBetZL_Y1iyk,217927
@@ -374,9 +374,9 @@ cellprofiler/modules/plugins/segmentationtemplatewithdependencies.py,sha256=Sh76
 cellprofiler/utilities/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 cellprofiler/utilities/morphology.py,sha256=8-81TrP8AmE3ETXIvQUKFD1vmKNBy2lfbc1QnM1eGIM,2685
 cellprofiler/utilities/rules.py,sha256=NoIHwFTA37zGvIP7vcB-aYeys0MDYVYxspfhLJe00OU,8790
-cellprofiler_nightly-5.0.0.dev318.dist-info/licenses/LICENSE,sha256=QLWaBS7kAioYx7PmJNXAMJaY8NODcFAag60YlUWuyz0,2276
-cellprofiler_nightly-5.0.0.dev318.dist-info/METADATA,sha256=x_cH5kn90YFPhDW1njHHwH4IHwMxcUKR5bfXgxSIKvQ,6063
-cellprofiler_nightly-5.0.0.dev318.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-cellprofiler_nightly-5.0.0.dev318.dist-info/entry_points.txt,sha256=MNDCjguFW3dKiS5Pcdu1NfWo4I0HHI3DekJLUJ4AKkY,60
-cellprofiler_nightly-5.0.0.dev318.dist-info/top_level.txt,sha256=bK7AacDeSj9qAmW8MGlO5wA79hDj6-ACt_mENUNKSIk,13
-cellprofiler_nightly-5.0.0.dev318.dist-info/RECORD,,
+cellprofiler_nightly-5.0.0.dev324.dist-info/licenses/LICENSE,sha256=QLWaBS7kAioYx7PmJNXAMJaY8NODcFAag60YlUWuyz0,2276
+cellprofiler_nightly-5.0.0.dev324.dist-info/METADATA,sha256=rcJmOQxLnbrpJc0xxW9H88YMgZrJgDMvTA-W-YDrcPE,6063
+cellprofiler_nightly-5.0.0.dev324.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+cellprofiler_nightly-5.0.0.dev324.dist-info/entry_points.txt,sha256=MNDCjguFW3dKiS5Pcdu1NfWo4I0HHI3DekJLUJ4AKkY,60
+cellprofiler_nightly-5.0.0.dev324.dist-info/top_level.txt,sha256=bK7AacDeSj9qAmW8MGlO5wA79hDj6-ACt_mENUNKSIk,13
+cellprofiler_nightly-5.0.0.dev324.dist-info/RECORD,,