014a7b779c4807ab8426dfebfc686aa03db98a2a
1 #!/usr/bin/env python
2 """
3 Copyright (C) 2005 Aaron Spike, aaron@ekips.org
5 This program is free software; you can redistribute it and/or modify
6 it under the terms of the GNU General Public License as published by
7 the Free Software Foundation; either version 2 of the License, or
8 (at your option) any later version.
10 This program is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 GNU General Public License for more details.
15 You should have received a copy of the GNU General Public License
16 along with this program; if not, write to the Free Software
17 Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
19 Perspective approach & math by Dmitry Platonov, shadowjack@mail.ru, 2006
20 """
21 import sys, inkex, os, re, simplepath, cubicsuperpath
22 from ffgeom import *
23 try:
24 from numpy import *
25 from numpy.linalg import *
26 except:
27 inkex.debug("Failed to import the numpy or numpy.linalg modules. These modules are required by this extension. Please install them and try again.")
28 sys.exit()
30 uuconv = {'in':90.0, 'pt':1.25, 'px':1, 'mm':3.5433070866, 'cm':35.433070866, 'pc':15.0}
31 def unittouu(string):
32 unit = re.compile('(%s)$' % '|'.join(uuconv.keys()))
33 param = re.compile(r'(([-+]?[0-9]+(\.[0-9]*)?|[-+]?\.[0-9]+)([eE][-+]?[0-9]+)?)')
35 p = param.match(string)
36 u = unit.search(string)
37 if p:
38 retval = float(p.string[p.start():p.end()])
39 else:
40 retval = 0.0
41 if u:
42 try:
43 return retval * uuconv[u.string[u.start():u.end()]]
44 except KeyError:
45 pass
46 return retval
48 class Project(inkex.Effect):
49 def __init__(self):
50 inkex.Effect.__init__(self)
51 def effect(self):
52 if len(self.options.ids) < 2:
53 inkex.debug("Requires two selected paths. The second must be exactly four nodes long.")
54 sys.exit()
56 #obj is selected second
57 obj = self.selected[self.options.ids[0]]
58 envelope = self.selected[self.options.ids[1]]
59 if (obj.tagName == 'path' or obj.tagName == 'g') and envelope.tagName == 'path':
60 path = cubicsuperpath.parsePath(envelope.attributes.getNamedItem('d').value)
61 dp = zeros((4,2), dtype=float64)
62 for i in range(4):
63 dp[i][0] = path[0][i][1][0]
64 dp[i][1] = path[0][i][1][1]
66 #query inkscape about the bounding box of obj
67 q = {'x':0,'y':0,'width':0,'height':0}
68 file = self.args[-1]
69 id = self.options.ids[0]
70 for query in q.keys():
71 _,f,err = os.popen3("inkscape --query-%s --query-id=%s %s" % (query,id,file))
72 q[query] = float(f.read())
73 f.close()
74 err.close()
75 sp = array([[q['x'], q['y']+q['height']],[q['x'], q['y']],[q['x']+q['width'], q['y']],[q['x']+q['width'], q['y']+q['height']]], dtype=float64)
77 solmatrix = zeros((8,8), dtype=float64)
78 free_term = zeros((8), dtype=float64)
79 for i in (0,1,2,3):
80 solmatrix[i][0] = sp[i][0]
81 solmatrix[i][1] = sp[i][1]
82 solmatrix[i][2] = 1
83 solmatrix[i][6] = -dp[i][0]*sp[i][0]
84 solmatrix[i][7] = -dp[i][0]*sp[i][1]
85 solmatrix[i+4][3] = sp[i][0]
86 solmatrix[i+4][4] = sp[i][1]
87 solmatrix[i+4][5] = 1
88 solmatrix[i+4][6] = -dp[i][1]*sp[i][0]
89 solmatrix[i+4][7] = -dp[i][1]*sp[i][1]
90 free_term[i] = dp[i][0]
91 free_term[i+4] = dp[i][1]
93 res = solve(solmatrix, free_term)
94 projmatrix = array([[res[0],res[1],res[2]],[res[3],res[4],res[5]],[res[6],res[7],1.0]],dtype=float64)
95 if obj.tagName == "path":
96 self.process_path(obj,projmatrix)
97 if obj.tagName == "g":
98 self.process_group(obj,projmatrix)
101 def process_group(self,group,m):
102 for node in group.childNodes:
103 if node.nodeType==node.ELEMENT_NODE:
104 if node.tagName == 'path':
105 self.process_path(node,m)
106 if node.tagName == 'g':
107 self.process_group(node,m)
110 def process_path(self,path,m):
111 d = path.attributes.getNamedItem('d')
112 p = cubicsuperpath.parsePath(d.value)
113 for subs in p:
114 for csp in subs:
115 csp[0] = self.project_point(csp[0],m)
116 csp[1] = self.project_point(csp[1],m)
117 csp[2] = self.project_point(csp[2],m)
118 d.value = cubicsuperpath.formatPath(p)
122 def project_point(self,p,m):
123 x = p[0]
124 y = p[1]
125 return [(x*m[0][0] + y*m[0][1] + m[0][2])/(x*m[2][0]+y*m[2][1]+m[2][2]),(x*m[1][0] + y*m[1][1] + m[1][2])/(x*m[2][0]+y*m[2][1]+m[2][2])]
127 e = Project()
128 e.affect()