utils.py 9.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320
  1. #!/usr/bin/env python
  2. """
  3. Miscellaneous Routines.
  4. """
  5. import struct
  6. from sys import maxint as INF
  7. ## PNG Predictor
  8. ##
  9. def apply_png_predictor(pred, colors, columns, bitspercomponent, data):
  10. if bitspercomponent != 8:
  11. # unsupported
  12. raise ValueError(bitspercomponent)
  13. nbytes = colors*columns*bitspercomponent//8
  14. i = 0
  15. buf = ''
  16. line0 = '\x00' * columns
  17. for i in xrange(0, len(data), nbytes+1):
  18. ft = data[i]
  19. i += 1
  20. line1 = data[i:i+nbytes]
  21. line2 = ''
  22. if ft == '\x00':
  23. # PNG none
  24. line2 += line1
  25. elif ft == '\x01':
  26. # PNG sub (UNTESTED)
  27. c = 0
  28. for b in line1:
  29. c = (c+ord(b)) & 255
  30. line2 += chr(c)
  31. elif ft == '\x02':
  32. # PNG up
  33. for (a, b) in zip(line0, line1):
  34. c = (ord(a)+ord(b)) & 255
  35. line2 += chr(c)
  36. elif ft == '\x03':
  37. # PNG average (UNTESTED)
  38. c = 0
  39. for (a, b) in zip(line0, line1):
  40. c = ((c+ord(a)+ord(b))//2) & 255
  41. line2 += chr(c)
  42. else:
  43. # unsupported
  44. raise ValueError(ft)
  45. buf += line2
  46. line0 = line2
  47. return buf
  48. ## Matrix operations
  49. ##
  50. MATRIX_IDENTITY = (1, 0, 0, 1, 0, 0)
  51. def mult_matrix((a1, b1, c1, d1, e1, f1), (a0, b0, c0, d0, e0, f0)):
  52. """Returns the multiplication of two matrices."""
  53. return (a0*a1+c0*b1, b0*a1+d0*b1,
  54. a0*c1+c0*d1, b0*c1+d0*d1,
  55. a0*e1+c0*f1+e0, b0*e1+d0*f1+f0)
  56. def translate_matrix((a, b, c, d, e, f), (x, y)):
  57. """Translates a matrix by (x, y)."""
  58. return (a, b, c, d, x*a+y*c+e, x*b+y*d+f)
  59. def apply_matrix_pt((a, b, c, d, e, f), (x, y)):
  60. """Applies a matrix to a point."""
  61. return (a*x+c*y+e, b*x+d*y+f)
  62. def apply_matrix_norm((a, b, c, d, e, f), (p, q)):
  63. """Equivalent to apply_matrix_pt(M, (p,q)) - apply_matrix_pt(M, (0,0))"""
  64. return (a*p+c*q, b*p+d*q)
  65. ## Utility functions
  66. ##
  67. # isnumber
  68. def isnumber(x):
  69. return isinstance(x, (int, long, float))
  70. # uniq
  71. def uniq(objs):
  72. """Eliminates duplicated elements."""
  73. done = set()
  74. for obj in objs:
  75. if obj in done:
  76. continue
  77. done.add(obj)
  78. yield obj
  79. return
  80. # csort
  81. def csort(objs, key=lambda x: x):
  82. """Order-preserving sorting function."""
  83. idxs = dict((obj, i) for (i, obj) in enumerate(objs))
  84. return sorted(objs, key=lambda obj: (key(obj), idxs[obj]))
  85. # fsplit
  86. def fsplit(pred, objs):
  87. """Split a list into two classes according to the predicate."""
  88. t = []
  89. f = []
  90. for obj in objs:
  91. if pred(obj):
  92. t.append(obj)
  93. else:
  94. f.append(obj)
  95. return (t, f)
  96. # drange
  97. def drange(v0, v1, d):
  98. """Returns a discrete range."""
  99. assert v0 < v1
  100. return xrange(int(v0)//d, int(v1+d)//d)
  101. # get_bound
  102. def get_bound(pts):
  103. """Compute a minimal rectangle that covers all the points."""
  104. (x0, y0, x1, y1) = (INF, INF, -INF, -INF)
  105. for (x, y) in pts:
  106. x0 = min(x0, x)
  107. y0 = min(y0, y)
  108. x1 = max(x1, x)
  109. y1 = max(y1, y)
  110. return (x0, y0, x1, y1)
  111. # pick
  112. def pick(seq, func, maxobj=None):
  113. """Picks the object obj where func(obj) has the highest value."""
  114. maxscore = None
  115. for obj in seq:
  116. score = func(obj)
  117. if maxscore is None or maxscore < score:
  118. (maxscore, maxobj) = (score, obj)
  119. return maxobj
  120. # choplist
  121. def choplist(n, seq):
  122. """Groups every n elements of the list."""
  123. r = []
  124. for x in seq:
  125. r.append(x)
  126. if len(r) == n:
  127. yield tuple(r)
  128. r = []
  129. return
  130. # nunpack
  131. def nunpack(s, default=0):
  132. """Unpacks 1 to 4 byte integers (big endian)."""
  133. l = len(s)
  134. if not l:
  135. return default
  136. elif l == 1:
  137. return ord(s)
  138. elif l == 2:
  139. return struct.unpack('>H', s)[0]
  140. elif l == 3:
  141. return struct.unpack('>L', '\x00'+s)[0]
  142. elif l == 4:
  143. return struct.unpack('>L', s)[0]
  144. else:
  145. raise TypeError('invalid length: %d' % l)
  146. # decode_text
  147. PDFDocEncoding = ''.join(unichr(x) for x in (
  148. 0x0000, 0x0001, 0x0002, 0x0003, 0x0004, 0x0005, 0x0006, 0x0007,
  149. 0x0008, 0x0009, 0x000a, 0x000b, 0x000c, 0x000d, 0x000e, 0x000f,
  150. 0x0010, 0x0011, 0x0012, 0x0013, 0x0014, 0x0015, 0x0017, 0x0017,
  151. 0x02d8, 0x02c7, 0x02c6, 0x02d9, 0x02dd, 0x02db, 0x02da, 0x02dc,
  152. 0x0020, 0x0021, 0x0022, 0x0023, 0x0024, 0x0025, 0x0026, 0x0027,
  153. 0x0028, 0x0029, 0x002a, 0x002b, 0x002c, 0x002d, 0x002e, 0x002f,
  154. 0x0030, 0x0031, 0x0032, 0x0033, 0x0034, 0x0035, 0x0036, 0x0037,
  155. 0x0038, 0x0039, 0x003a, 0x003b, 0x003c, 0x003d, 0x003e, 0x003f,
  156. 0x0040, 0x0041, 0x0042, 0x0043, 0x0044, 0x0045, 0x0046, 0x0047,
  157. 0x0048, 0x0049, 0x004a, 0x004b, 0x004c, 0x004d, 0x004e, 0x004f,
  158. 0x0050, 0x0051, 0x0052, 0x0053, 0x0054, 0x0055, 0x0056, 0x0057,
  159. 0x0058, 0x0059, 0x005a, 0x005b, 0x005c, 0x005d, 0x005e, 0x005f,
  160. 0x0060, 0x0061, 0x0062, 0x0063, 0x0064, 0x0065, 0x0066, 0x0067,
  161. 0x0068, 0x0069, 0x006a, 0x006b, 0x006c, 0x006d, 0x006e, 0x006f,
  162. 0x0070, 0x0071, 0x0072, 0x0073, 0x0074, 0x0075, 0x0076, 0x0077,
  163. 0x0078, 0x0079, 0x007a, 0x007b, 0x007c, 0x007d, 0x007e, 0x0000,
  164. 0x2022, 0x2020, 0x2021, 0x2026, 0x2014, 0x2013, 0x0192, 0x2044,
  165. 0x2039, 0x203a, 0x2212, 0x2030, 0x201e, 0x201c, 0x201d, 0x2018,
  166. 0x2019, 0x201a, 0x2122, 0xfb01, 0xfb02, 0x0141, 0x0152, 0x0160,
  167. 0x0178, 0x017d, 0x0131, 0x0142, 0x0153, 0x0161, 0x017e, 0x0000,
  168. 0x20ac, 0x00a1, 0x00a2, 0x00a3, 0x00a4, 0x00a5, 0x00a6, 0x00a7,
  169. 0x00a8, 0x00a9, 0x00aa, 0x00ab, 0x00ac, 0x0000, 0x00ae, 0x00af,
  170. 0x00b0, 0x00b1, 0x00b2, 0x00b3, 0x00b4, 0x00b5, 0x00b6, 0x00b7,
  171. 0x00b8, 0x00b9, 0x00ba, 0x00bb, 0x00bc, 0x00bd, 0x00be, 0x00bf,
  172. 0x00c0, 0x00c1, 0x00c2, 0x00c3, 0x00c4, 0x00c5, 0x00c6, 0x00c7,
  173. 0x00c8, 0x00c9, 0x00ca, 0x00cb, 0x00cc, 0x00cd, 0x00ce, 0x00cf,
  174. 0x00d0, 0x00d1, 0x00d2, 0x00d3, 0x00d4, 0x00d5, 0x00d6, 0x00d7,
  175. 0x00d8, 0x00d9, 0x00da, 0x00db, 0x00dc, 0x00dd, 0x00de, 0x00df,
  176. 0x00e0, 0x00e1, 0x00e2, 0x00e3, 0x00e4, 0x00e5, 0x00e6, 0x00e7,
  177. 0x00e8, 0x00e9, 0x00ea, 0x00eb, 0x00ec, 0x00ed, 0x00ee, 0x00ef,
  178. 0x00f0, 0x00f1, 0x00f2, 0x00f3, 0x00f4, 0x00f5, 0x00f6, 0x00f7,
  179. 0x00f8, 0x00f9, 0x00fa, 0x00fb, 0x00fc, 0x00fd, 0x00fe, 0x00ff,
  180. ))
  181. def decode_text(s):
  182. """Decodes a PDFDocEncoding string to Unicode."""
  183. if s.startswith('\xfe\xff'):
  184. return unicode(s[2:], 'utf-16be', 'ignore')
  185. else:
  186. return ''.join(PDFDocEncoding[ord(c)] for c in s)
  187. # enc
  188. def enc(x, codec='ascii'):
  189. """Encodes a string for SGML/XML/HTML"""
  190. x = x.replace('&', '&amp;').replace('>', '&gt;').replace('<', '&lt;').replace('"', '&quot;')
  191. return x.encode(codec, 'xmlcharrefreplace')
  192. def bbox2str((x0, y0, x1, y1)):
  193. return '%.3f,%.3f,%.3f,%.3f' % (x0, y0, x1, y1)
  194. def matrix2str((a, b, c, d, e, f)):
  195. return '[%.2f,%.2f,%.2f,%.2f, (%.2f,%.2f)]' % (a, b, c, d, e, f)
  196. ## Plane
  197. ##
  198. ## A set-like data structure for objects placed on a plane.
  199. ## Can efficiently find objects in a certain rectangular area.
  200. ## It maintains two parallel lists of objects, each of
  201. ## which is sorted by its x or y coordinate.
  202. ##
  203. class Plane(object):
  204. def __init__(self, bbox, gridsize=50):
  205. self._objs = set()
  206. self._grid = {}
  207. self.gridsize = gridsize
  208. (self.x0, self.y0, self.x1, self.y1) = bbox
  209. return
  210. def __repr__(self):
  211. return ('<Plane objs=%r>' % list(self))
  212. def __iter__(self):
  213. return iter(self._objs)
  214. def __len__(self):
  215. return len(self._objs)
  216. def __contains__(self, obj):
  217. return obj in self._objs
  218. def _getrange(self, (x0, y0, x1, y1)):
  219. if (x1 <= self.x0 or self.x1 <= x0 or
  220. y1 <= self.y0 or self.y1 <= y0): return
  221. x0 = max(self.x0, x0)
  222. y0 = max(self.y0, y0)
  223. x1 = min(self.x1, x1)
  224. y1 = min(self.y1, y1)
  225. for y in drange(y0, y1, self.gridsize):
  226. for x in drange(x0, x1, self.gridsize):
  227. yield (x, y)
  228. return
  229. # extend(objs)
  230. def extend(self, objs):
  231. for obj in objs:
  232. self.add(obj)
  233. return
  234. # add(obj): place an object.
  235. def add(self, obj):
  236. for k in self._getrange((obj.x0, obj.y0, obj.x1, obj.y1)):
  237. if k not in self._grid:
  238. r = []
  239. self._grid[k] = r
  240. else:
  241. r = self._grid[k]
  242. r.append(obj)
  243. self._objs.add(obj)
  244. return
  245. # remove(obj): displace an object.
  246. def remove(self, obj):
  247. for k in self._getrange((obj.x0, obj.y0, obj.x1, obj.y1)):
  248. try:
  249. self._grid[k].remove(obj)
  250. except (KeyError, ValueError):
  251. pass
  252. self._objs.remove(obj)
  253. return
  254. # find(): finds objects that are in a certain area.
  255. def find(self, (x0, y0, x1, y1)):
  256. done = set()
  257. for k in self._getrange((x0, y0, x1, y1)):
  258. if k not in self._grid:
  259. continue
  260. for obj in self._grid[k]:
  261. if obj in done:
  262. continue
  263. done.add(obj)
  264. if (obj.x1 <= x0 or x1 <= obj.x0 or
  265. obj.y1 <= y0 or y1 <= obj.y0):
  266. continue
  267. yield obj
  268. return