emdbva 0.0.1.dev135__tar.gz → 0.0.1.dev136__tar.gz

{emdbva-0.0.1.dev135/emdbva.egg-info → emdbva-0.0.1.dev136}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: emdbva
-Version: 0.0.1.dev135
+Version: 0.0.1.dev136
 Summary: CryoEM validation toolkit
 Home-page: https://test.pypi.org/project/va/
 Author: Zhe Wang

{emdbva-0.0.1.dev135 → emdbva-0.0.1.dev136/emdbva.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: emdbva
-Version: 0.0.1.dev135
+Version: 0.0.1.dev136
 Summary: CryoEM validation toolkit
 Home-page: https://test.pypi.org/project/va/
 Author: Zhe Wang

{emdbva-0.0.1.dev135 → emdbva-0.0.1.dev136}/va/metrics/projections.py

@@ -183,7 +183,10 @@ class Projections:
     def _green_percentage(self, image_obj):
         """
         image_obj can be a file path, PIL Image, or numpy array.
-        Returns green percentage as a float.
+        Returns a dictionary with:
+          - percentage: total green percentage in the whole image
+          - diff_vertical: left-half green percentage minus right-half green percentage
+          - diff_horizontal: top-half green percentage minus bottom-half green percentage
         """
         if isinstance(image_obj, str):
             img = Image.open(image_obj).convert("RGB")
@@ -193,10 +196,43 @@ class Projections:
             img = Image.fromarray(image_obj).convert("RGB")
 
         img_array = np.array(img)
-        green_pixels = np.sum(np.all(img_array == [0, 138, 0], axis=-1))
-        total_pixels = img_array.shape[0] * img_array.shape[1]
-        proportion_green = (green_pixels / total_pixels) * 100 if total_pixels else 0.0
-        return round(proportion_green, 2)
+        green_mask = np.all(img_array == [0, 138, 0], axis=-1)
+        height, width = green_mask.shape
+        total_pixels = height * width
+        if total_pixels == 0:
+            return {
+                'percentage': 0.0,
+                'diff_vertical': 0.0,
+                'diff_horizontal': 0.0,
+            }
+
+        green_pixels = np.sum(green_mask)
+        proportion_green = (green_pixels / total_pixels) * 100
+
+        mid_vertical = width // 2
+        left_total = height * mid_vertical
+        right_total = height * (width - mid_vertical)
+        green_pixels_left_half = np.sum(green_mask[:, :mid_vertical])
+        green_pixels_right_half = np.sum(green_mask[:, mid_vertical:])
+        proportion_green_left_half = (green_pixels_left_half / left_total) * 100 if left_total else 0.0
+        proportion_green_right_half = (green_pixels_right_half / right_total) * 100 if right_total else 0.0
+        diff_vertical = proportion_green_left_half - proportion_green_right_half
+
+        mid_horizontal = height // 2
+        top_total = mid_horizontal * width
+        bottom_total = (height - mid_horizontal) * width
+        green_pixels_top_half = np.sum(green_mask[:mid_horizontal, :])
+        green_pixels_bottom_half = np.sum(green_mask[mid_horizontal:, :])
+        proportion_green_top_half = (green_pixels_top_half / top_total) * 100 if top_total else 0.0
+        proportion_green_bottom_half = (green_pixels_bottom_half / bottom_total) * 100 if bottom_total else 0.0
+        diff_horizontal = proportion_green_top_half - proportion_green_bottom_half
+
+        return {
+            'percentage': round(proportion_green, 2),
+            'diff_vertical': round(diff_vertical, 2),
+            'diff_horizontal': round(diff_horizontal, 2),
+        }
+
 
     def orthogonal_projections(self, mapin=None, workdir=None, type=None, label=''):
         map, workdir = self.mapincheck(mapin, workdir)

{emdbva-0.0.1.dev135 → emdbva-0.0.1.dev136}/va/version.py

@@ -18,5 +18,5 @@ under the License.
 
 """
 
-__version__ = '0.0.1.dev135'
+__version__ = '0.0.1.dev136'
 __em_statistics_version__ = '202505.v01'