webknossos.geometry.mag

View Source
import re
from functools import total_ordering
from math import log2
from typing import Any, List

import numpy as np


@total_ordering
class Mag(object):
    def __init__(self, mag: Any):
        self.mag: List[int] = []

        if isinstance(mag, int):
            self.mag = [mag] * 3
        elif isinstance(mag, list):
            self.mag = mag
        elif isinstance(mag, tuple):
            self.mag = [mag_d for mag_d in mag]
        elif isinstance(mag, str):
            if re.match(r"^\d+$", mag) is not None:
                self.mag = [int(mag)] * 3
            elif re.match(r"^\d+-\d+-\d+$", mag) is not None:
                self.mag = [int(m) for m in mag.split("-")]
        elif isinstance(mag, Mag):
            self.mag = mag.mag
        elif isinstance(mag, np.ndarray):
            assert mag.shape == (3,)
            self.mag = list(mag)

        if self.mag is None or len(self.mag) != 3:
            raise ValueError(
                "magnification must be int or a vector3 of ints or a string shaped like e.g. 2-2-1"
            )

        for m in self.mag:
            assert log2(m) % 1 == 0, "magnification needs to be power of 2."

    def __lt__(self, other: Any) -> bool:
        return max(self.mag) < (max(Mag(other).to_array()))

    def __le__(self, other: Any) -> bool:
        return max(self.mag) <= (max(Mag(other).to_array()))

    def __eq__(self, other: Any) -> bool:
        return all(m1 == m2 for m1, m2 in zip(self.mag, Mag(other).mag))

    def __str__(self) -> str:
        return self.to_layer_name()

    def __expr__(self) -> str:
        return f"Mag({self.to_layer_name()})"

    def to_layer_name(self) -> str:
        x, y, z = self.mag
        if x == y and y == z:
            return str(x)
        else:
            return self.to_long_layer_name()

    def to_long_layer_name(self) -> str:
        x, y, z = self.mag
        return "{}-{}-{}".format(x, y, z)

    def to_array(self) -> List[int]:
        return self.mag

    def scaled_by(self, factor: int) -> "Mag":
        return Mag([mag * factor for mag in self.mag])

    def scale_by(self, factor: int) -> None:
        self.mag = [mag * factor for mag in self.mag]

    def divided(self, coord: List[int]) -> List[int]:
        return [c // m for c, m in zip(coord, self.mag)]

    def divide_by(self, d: int) -> None:
        self.mag = [mag // d for mag in self.mag]

    def divided_by(self, d: int) -> "Mag":
        return Mag([mag // d for mag in self.mag])

    def as_np(self) -> np.ndarray:
        return np.array(self.mag)

    def __hash__(self) -> int:
        return hash(tuple(self.mag))
#   class Mag:
View Source
class Mag(object):
    def __init__(self, mag: Any):
        self.mag: List[int] = []

        if isinstance(mag, int):
            self.mag = [mag] * 3
        elif isinstance(mag, list):
            self.mag = mag
        elif isinstance(mag, tuple):
            self.mag = [mag_d for mag_d in mag]
        elif isinstance(mag, str):
            if re.match(r"^\d+$", mag) is not None:
                self.mag = [int(mag)] * 3
            elif re.match(r"^\d+-\d+-\d+$", mag) is not None:
                self.mag = [int(m) for m in mag.split("-")]
        elif isinstance(mag, Mag):
            self.mag = mag.mag
        elif isinstance(mag, np.ndarray):
            assert mag.shape == (3,)
            self.mag = list(mag)

        if self.mag is None or len(self.mag) != 3:
            raise ValueError(
                "magnification must be int or a vector3 of ints or a string shaped like e.g. 2-2-1"
            )

        for m in self.mag:
            assert log2(m) % 1 == 0, "magnification needs to be power of 2."

    def __lt__(self, other: Any) -> bool:
        return max(self.mag) < (max(Mag(other).to_array()))

    def __le__(self, other: Any) -> bool:
        return max(self.mag) <= (max(Mag(other).to_array()))

    def __eq__(self, other: Any) -> bool:
        return all(m1 == m2 for m1, m2 in zip(self.mag, Mag(other).mag))

    def __str__(self) -> str:
        return self.to_layer_name()

    def __expr__(self) -> str:
        return f"Mag({self.to_layer_name()})"

    def to_layer_name(self) -> str:
        x, y, z = self.mag
        if x == y and y == z:
            return str(x)
        else:
            return self.to_long_layer_name()

    def to_long_layer_name(self) -> str:
        x, y, z = self.mag
        return "{}-{}-{}".format(x, y, z)

    def to_array(self) -> List[int]:
        return self.mag

    def scaled_by(self, factor: int) -> "Mag":
        return Mag([mag * factor for mag in self.mag])

    def scale_by(self, factor: int) -> None:
        self.mag = [mag * factor for mag in self.mag]

    def divided(self, coord: List[int]) -> List[int]:
        return [c // m for c, m in zip(coord, self.mag)]

    def divide_by(self, d: int) -> None:
        self.mag = [mag // d for mag in self.mag]

    def divided_by(self, d: int) -> "Mag":
        return Mag([mag // d for mag in self.mag])

    def as_np(self) -> np.ndarray:
        return np.array(self.mag)

    def __hash__(self) -> int:
        return hash(tuple(self.mag))
#   Mag(mag: Any)
View Source
    def __init__(self, mag: Any):
        self.mag: List[int] = []

        if isinstance(mag, int):
            self.mag = [mag] * 3
        elif isinstance(mag, list):
            self.mag = mag
        elif isinstance(mag, tuple):
            self.mag = [mag_d for mag_d in mag]
        elif isinstance(mag, str):
            if re.match(r"^\d+$", mag) is not None:
                self.mag = [int(mag)] * 3
            elif re.match(r"^\d+-\d+-\d+$", mag) is not None:
                self.mag = [int(m) for m in mag.split("-")]
        elif isinstance(mag, Mag):
            self.mag = mag.mag
        elif isinstance(mag, np.ndarray):
            assert mag.shape == (3,)
            self.mag = list(mag)

        if self.mag is None or len(self.mag) != 3:
            raise ValueError(
                "magnification must be int or a vector3 of ints or a string shaped like e.g. 2-2-1"
            )

        for m in self.mag:
            assert log2(m) % 1 == 0, "magnification needs to be power of 2."
#   def to_layer_name(self) -> str:
View Source
    def to_layer_name(self) -> str:
        x, y, z = self.mag
        if x == y and y == z:
            return str(x)
        else:
            return self.to_long_layer_name()
#   def to_long_layer_name(self) -> str:
View Source
    def to_long_layer_name(self) -> str:
        x, y, z = self.mag
        return "{}-{}-{}".format(x, y, z)
#   def to_array(self) -> List[int]:
View Source
    def to_array(self) -> List[int]:
        return self.mag
#   def scaled_by(self, factor: int) -> webknossos.geometry.mag.Mag:
View Source
    def scaled_by(self, factor: int) -> "Mag":
        return Mag([mag * factor for mag in self.mag])
#   def scale_by(self, factor: int) -> None:
View Source
    def scale_by(self, factor: int) -> None:
        self.mag = [mag * factor for mag in self.mag]
#   def divided(self, coord: List[int]) -> List[int]:
View Source
    def divided(self, coord: List[int]) -> List[int]:
        return [c // m for c, m in zip(coord, self.mag)]
#   def divide_by(self, d: int) -> None:
View Source
    def divide_by(self, d: int) -> None:
        self.mag = [mag // d for mag in self.mag]
#   def divided_by(self, d: int) -> webknossos.geometry.mag.Mag:
View Source
    def divided_by(self, d: int) -> "Mag":
        return Mag([mag // d for mag in self.mag])
#   def as_np(self) -> numpy.ndarray:
View Source
    def as_np(self) -> np.ndarray:
        return np.array(self.mag)