ImageMath.py 7.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263
  1. #
  2. # The Python Imaging Library
  3. # $Id$
  4. #
  5. # a simple math add-on for the Python Imaging Library
  6. #
  7. # History:
  8. # 1999-02-15 fl Original PIL Plus release
  9. # 2005-05-05 fl Simplified and cleaned up for PIL 1.1.6
  10. # 2005-09-12 fl Fixed int() and float() for Python 2.4.1
  11. #
  12. # Copyright (c) 1999-2005 by Secret Labs AB
  13. # Copyright (c) 2005 by Fredrik Lundh
  14. #
  15. # See the README file for information on usage and redistribution.
  16. #
  17. import builtins
  18. from . import Image, _imagingmath
  19. def _isconstant(v):
  20. return isinstance(v, (int, float))
  21. class _Operand:
  22. """Wraps an image operand, providing standard operators"""
  23. def __init__(self, im):
  24. self.im = im
  25. def __fixup(self, im1):
  26. # convert image to suitable mode
  27. if isinstance(im1, _Operand):
  28. # argument was an image.
  29. if im1.im.mode in ("1", "L"):
  30. return im1.im.convert("I")
  31. elif im1.im.mode in ("I", "F"):
  32. return im1.im
  33. else:
  34. msg = f"unsupported mode: {im1.im.mode}"
  35. raise ValueError(msg)
  36. else:
  37. # argument was a constant
  38. if _isconstant(im1) and self.im.mode in ("1", "L", "I"):
  39. return Image.new("I", self.im.size, im1)
  40. else:
  41. return Image.new("F", self.im.size, im1)
  42. def apply(self, op, im1, im2=None, mode=None):
  43. im1 = self.__fixup(im1)
  44. if im2 is None:
  45. # unary operation
  46. out = Image.new(mode or im1.mode, im1.size, None)
  47. im1.load()
  48. try:
  49. op = getattr(_imagingmath, op + "_" + im1.mode)
  50. except AttributeError as e:
  51. msg = f"bad operand type for '{op}'"
  52. raise TypeError(msg) from e
  53. _imagingmath.unop(op, out.im.id, im1.im.id)
  54. else:
  55. # binary operation
  56. im2 = self.__fixup(im2)
  57. if im1.mode != im2.mode:
  58. # convert both arguments to floating point
  59. if im1.mode != "F":
  60. im1 = im1.convert("F")
  61. if im2.mode != "F":
  62. im2 = im2.convert("F")
  63. if im1.size != im2.size:
  64. # crop both arguments to a common size
  65. size = (min(im1.size[0], im2.size[0]), min(im1.size[1], im2.size[1]))
  66. if im1.size != size:
  67. im1 = im1.crop((0, 0) + size)
  68. if im2.size != size:
  69. im2 = im2.crop((0, 0) + size)
  70. out = Image.new(mode or im1.mode, im1.size, None)
  71. im1.load()
  72. im2.load()
  73. try:
  74. op = getattr(_imagingmath, op + "_" + im1.mode)
  75. except AttributeError as e:
  76. msg = f"bad operand type for '{op}'"
  77. raise TypeError(msg) from e
  78. _imagingmath.binop(op, out.im.id, im1.im.id, im2.im.id)
  79. return _Operand(out)
  80. # unary operators
  81. def __bool__(self):
  82. # an image is "true" if it contains at least one non-zero pixel
  83. return self.im.getbbox() is not None
  84. def __abs__(self):
  85. return self.apply("abs", self)
  86. def __pos__(self):
  87. return self
  88. def __neg__(self):
  89. return self.apply("neg", self)
  90. # binary operators
  91. def __add__(self, other):
  92. return self.apply("add", self, other)
  93. def __radd__(self, other):
  94. return self.apply("add", other, self)
  95. def __sub__(self, other):
  96. return self.apply("sub", self, other)
  97. def __rsub__(self, other):
  98. return self.apply("sub", other, self)
  99. def __mul__(self, other):
  100. return self.apply("mul", self, other)
  101. def __rmul__(self, other):
  102. return self.apply("mul", other, self)
  103. def __truediv__(self, other):
  104. return self.apply("div", self, other)
  105. def __rtruediv__(self, other):
  106. return self.apply("div", other, self)
  107. def __mod__(self, other):
  108. return self.apply("mod", self, other)
  109. def __rmod__(self, other):
  110. return self.apply("mod", other, self)
  111. def __pow__(self, other):
  112. return self.apply("pow", self, other)
  113. def __rpow__(self, other):
  114. return self.apply("pow", other, self)
  115. # bitwise
  116. def __invert__(self):
  117. return self.apply("invert", self)
  118. def __and__(self, other):
  119. return self.apply("and", self, other)
  120. def __rand__(self, other):
  121. return self.apply("and", other, self)
  122. def __or__(self, other):
  123. return self.apply("or", self, other)
  124. def __ror__(self, other):
  125. return self.apply("or", other, self)
  126. def __xor__(self, other):
  127. return self.apply("xor", self, other)
  128. def __rxor__(self, other):
  129. return self.apply("xor", other, self)
  130. def __lshift__(self, other):
  131. return self.apply("lshift", self, other)
  132. def __rshift__(self, other):
  133. return self.apply("rshift", self, other)
  134. # logical
  135. def __eq__(self, other):
  136. return self.apply("eq", self, other)
  137. def __ne__(self, other):
  138. return self.apply("ne", self, other)
  139. def __lt__(self, other):
  140. return self.apply("lt", self, other)
  141. def __le__(self, other):
  142. return self.apply("le", self, other)
  143. def __gt__(self, other):
  144. return self.apply("gt", self, other)
  145. def __ge__(self, other):
  146. return self.apply("ge", self, other)
  147. # conversions
  148. def imagemath_int(self):
  149. return _Operand(self.im.convert("I"))
  150. def imagemath_float(self):
  151. return _Operand(self.im.convert("F"))
  152. # logical
  153. def imagemath_equal(self, other):
  154. return self.apply("eq", self, other, mode="I")
  155. def imagemath_notequal(self, other):
  156. return self.apply("ne", self, other, mode="I")
  157. def imagemath_min(self, other):
  158. return self.apply("min", self, other)
  159. def imagemath_max(self, other):
  160. return self.apply("max", self, other)
  161. def imagemath_convert(self, mode):
  162. return _Operand(self.im.convert(mode))
  163. ops = {}
  164. for k, v in list(globals().items()):
  165. if k[:10] == "imagemath_":
  166. ops[k[10:]] = v
  167. def eval(expression, _dict={}, **kw):
  168. """
  169. Evaluates an image expression.
  170. :param expression: A string containing a Python-style expression.
  171. :param options: Values to add to the evaluation context. You
  172. can either use a dictionary, or one or more keyword
  173. arguments.
  174. :return: The evaluated expression. This is usually an image object, but can
  175. also be an integer, a floating point value, or a pixel tuple,
  176. depending on the expression.
  177. """
  178. # build execution namespace
  179. args = ops.copy()
  180. args.update(_dict)
  181. args.update(kw)
  182. for k, v in list(args.items()):
  183. if hasattr(v, "im"):
  184. args[k] = _Operand(v)
  185. compiled_code = compile(expression, "<string>", "eval")
  186. def scan(code):
  187. for const in code.co_consts:
  188. if type(const) is type(compiled_code):
  189. scan(const)
  190. for name in code.co_names:
  191. if name not in args and name != "abs":
  192. msg = f"'{name}' not allowed"
  193. raise ValueError(msg)
  194. scan(compiled_code)
  195. out = builtins.eval(expression, {"__builtins": {"abs": abs}}, args)
  196. try:
  197. return out.im
  198. except AttributeError:
  199. return out