Infinite Quadtree precision now works with some B├ęziers
[ipdf/code.git] / src / document.cpp
1 #include "document.h"
2 #include "bezier.h"
3 #include <cstdio>
4 #include <fstream>
5
6 #include "../contrib/pugixml-1.4/src/pugixml.cpp"
7
8 #include "stb_truetype.h"
9
10 using namespace IPDF;
11 using namespace std;
12
13 //TODO: Make this work for variable sized Reals
14
15 // Loads an std::vector<T> of size num_elements from a file.
16 template<typename T>
17 static void LoadStructVector(FILE *src_file, size_t num_elems, std::vector<T>& dest)
18 {
19         size_t structsread = 0;
20         dest.resize(num_elems);
21         structsread = fread(dest.data(), sizeof(T), num_elems, src_file);
22         if (structsread != num_elems)
23                 Fatal("Only read %u structs (expected %u)!", structsread, num_elems);
24 }
25
26 // Saves an std::vector<T> to a file. Size must be saves separately.
27 template<typename T>
28 static void SaveStructVector(FILE *dst_file, std::vector<T>& src)
29 {
30         size_t written = 0;
31         written = fwrite(src.data(), sizeof(T), src.size(), dst_file);
32         if (written != src.size())
33                 Fatal("Only wrote %u structs (expected %u)!", written, src.size());
34 }
35
36 static void WriteChunkHeader(FILE *dst_file, DocChunkTypes type, uint32_t size)
37 {
38         size_t written = 0;
39         written = fwrite(&type, sizeof(type), 1, dst_file);
40         if (written != 1)
41                 Fatal("Could not write Chunk header! (ID)");
42         written = fwrite(&size, sizeof(size), 1, dst_file);
43         if (written != 1)
44                 Fatal("Could not write Chunk header (size)!");
45 }
46
47 static bool ReadChunkHeader(FILE *src_file, DocChunkTypes& type, uint32_t& size)
48 {
49         if (fread(&type, sizeof(DocChunkTypes), 1, src_file) != 1)
50                 return false;
51         if (fread(&size, sizeof(uint32_t), 1, src_file) != 1)
52                 return false;
53         return true;
54 }
55
56 void Document::Save(const string & filename)
57 {
58         Debug("Saving document to file \"%s\"...", filename.c_str());
59         FILE * file = fopen(filename.c_str(), "w");
60         if (file == NULL)
61                 Fatal("Couldn't open file \"%s\" - %s", filename.c_str(), strerror(errno));
62
63         size_t written;
64         Debug("Number of objects (%u)...", ObjectCount());
65         WriteChunkHeader(file, CT_NUMOBJS, sizeof(m_count));
66         written = fwrite(&m_count, sizeof(m_count), 1, file);
67         if (written != 1)
68                 Fatal("Failed to write number of objects!");
69
70         Debug("Object types...");
71         WriteChunkHeader(file, CT_OBJTYPES, m_objects.types.size() * sizeof(ObjectType));
72         SaveStructVector<ObjectType>(file, m_objects.types);
73
74         Debug("Object bounds...");
75         WriteChunkHeader(file, CT_OBJBOUNDS, m_objects.bounds.size() * sizeof(Rect));
76         SaveStructVector<Rect>(file, m_objects.bounds);
77
78         Debug("Object data indices...");
79         WriteChunkHeader(file, CT_OBJINDICES, m_objects.data_indices.size() * sizeof(unsigned));
80         SaveStructVector<unsigned>(file, m_objects.data_indices);
81         
82         Debug("Bezier data...");
83         WriteChunkHeader(file, CT_OBJBEZIERS, m_objects.beziers.size() * sizeof(uint8_t));
84         SaveStructVector<Bezier>(file, m_objects.beziers);
85
86         int err = fclose(file);
87         if (err != 0)
88                 Fatal("Failed to close file \"%s\" - %s", filename.c_str(), strerror(err));
89
90         Debug("Successfully saved %u objects to \"%s\"", ObjectCount(), filename.c_str());
91 }
92
93 #ifndef QUADTREE_DISABLED
94
95 void Document::GenBaseQuadtree()
96 {
97         m_quadtree.nodes.push_back(QuadTreeNode{QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, QTC_UNKNOWN, 0, ObjectCount()});
98         m_quadtree.root_id = 0;
99         GenQuadChild(0, QTC_TOP_LEFT);
100 }
101
102 int Document::ClipObjectToQuadChild(int object_id, QuadTreeNodeChildren type)
103 {
104         switch (m_objects.types[object_id])
105         {
106         case RECT_FILLED:
107         case RECT_OUTLINE:
108                 {
109                 Rect obj_bounds = TransformToQuadChild(m_objects.bounds[object_id], type);
110                 if (obj_bounds.x < 0)
111                 {
112                         obj_bounds.w += obj_bounds.x;
113                         obj_bounds.x = 0;
114                 }
115                 if (obj_bounds.y < 0)
116                 {
117                         obj_bounds.h += obj_bounds.y;
118                         obj_bounds.y = 0;
119                 }
120                 if (obj_bounds.x + obj_bounds.w > 1)
121                 {
122                         obj_bounds.w += (1 - (obj_bounds.x + obj_bounds.w));
123                 }
124                 if (obj_bounds.y + obj_bounds.h > 1)
125                 {
126                         obj_bounds.h += (1 - (obj_bounds.y + obj_bounds.h));
127                 }
128                 m_objects.bounds.push_back(obj_bounds);
129                 m_objects.types.push_back(m_objects.types[object_id]);
130                 m_objects.data_indices.push_back(m_objects.data_indices[object_id]);
131                 return 1;
132                 }
133         case BEZIER:
134                 {
135 //              Rect child_node_bounds = m_objects.bounds[object_id];
136                 Rect clip_bezier_bounds = TransformRectCoordinates(m_objects.bounds[object_id], TransformFromQuadChild(Rect{0,0,1,1}, type)); 
137                 std::vector<Bezier> new_curves = m_objects.beziers[m_objects.data_indices[object_id]].ClipToRectangle(clip_bezier_bounds);
138                 for (size_t i = 0; i < new_curves.size(); ++i)
139                 {
140                         Rect new_bounds = TransformToQuadChild(m_objects.bounds[object_id], type);
141                         Bezier new_curve_data = new_curves[i].ToAbsolute(TransformToQuadChild(m_objects.bounds[object_id],type));
142                         new_bounds = new_curve_data.SolveBounds();
143                         new_curve_data = new_curve_data.ToRelative(new_bounds);
144                         unsigned index = AddBezierData(new_curve_data);
145                         m_objects.bounds.push_back(new_bounds);
146                         m_objects.types.push_back(BEZIER);
147                         m_objects.data_indices.push_back(index);
148                 }
149                 return new_curves.size();
150                 }
151         default:
152                 Debug("Adding %s -> %s", m_objects.bounds[object_id].Str().c_str(), TransformToQuadChild(m_objects.bounds[object_id], type).Str().c_str());
153                 m_objects.bounds.push_back(TransformToQuadChild(m_objects.bounds[object_id], type));
154                 m_objects.types.push_back(m_objects.types[object_id]);
155                 m_objects.data_indices.push_back(m_objects.data_indices[object_id]);
156                 return 1;
157         }
158         return 0;
159 }
160 QuadTreeIndex Document::GenQuadChild(QuadTreeIndex parent, QuadTreeNodeChildren type)
161 {
162         QuadTreeIndex new_index = m_quadtree.nodes.size();
163         m_quadtree.nodes.push_back(QuadTreeNode{QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, parent, type, 0, 0});
164
165         m_quadtree.nodes[new_index].object_begin = m_objects.bounds.size();
166         for (unsigned i = m_quadtree.nodes[parent].object_begin; i < m_quadtree.nodes[parent].object_end; ++i)
167         {
168                 if (IntersectsQuadChild(m_objects.bounds[i], type))
169                 {
170                         m_count += ClipObjectToQuadChild(i, type);
171                 }
172         }
173         m_quadtree.nodes[new_index].object_end = m_objects.bounds.size();
174         switch (type)
175         {
176                 case QTC_TOP_LEFT:
177                         m_quadtree.nodes[parent].top_left = new_index;
178                         break;
179                 case QTC_TOP_RIGHT:
180                         m_quadtree.nodes[parent].top_right = new_index;
181                         break;
182                 case QTC_BOTTOM_LEFT:
183                         m_quadtree.nodes[parent].bottom_left = new_index;
184                         break;
185                 case QTC_BOTTOM_RIGHT:
186                         m_quadtree.nodes[parent].bottom_right = new_index;
187                         break;
188                 default:
189                         Fatal("Tried to add a QuadTree child of invalid type!");
190         }
191         return new_index;
192 }
193
194 // Reparent a quadtree node, making it the "type" child of a new node.
195 QuadTreeIndex Document::GenQuadParent(QuadTreeIndex child, QuadTreeNodeChildren type)
196 {
197         QuadTreeIndex new_index = m_quadtree.nodes.size();
198         m_quadtree.nodes.push_back(QuadTreeNode{QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, -1, QTC_UNKNOWN, 0, 0});
199
200         m_quadtree.nodes[new_index].object_begin = m_objects.bounds.size();
201         for (unsigned i = m_quadtree.nodes[child].object_begin; i < m_quadtree.nodes[child].object_end; ++i)
202         {
203                 m_objects.bounds.push_back(TransformFromQuadChild(m_objects.bounds[i], type));
204                 m_objects.types.push_back(m_objects.types[i]);
205                 m_objects.data_indices.push_back(m_objects.data_indices[i]);
206                 m_count++;
207         }
208         m_quadtree.nodes[new_index].object_end = m_objects.bounds.size();
209         switch (type)
210         {
211                 case QTC_TOP_LEFT:
212                         m_quadtree.nodes[new_index].top_left = child;
213                         break;
214                 case QTC_TOP_RIGHT:
215                         m_quadtree.nodes[new_index].top_right = child;
216                         break;
217                 case QTC_BOTTOM_LEFT:
218                         m_quadtree.nodes[new_index].bottom_left = child;
219                         break;
220                 case QTC_BOTTOM_RIGHT:
221                         m_quadtree.nodes[new_index].bottom_right = child;
222                         break;
223                 default:
224                         Fatal("Tried to add a QuadTree child of invalid type!");
225         }
226         return new_index;
227 }
228
229 #endif
230
231 void Document::Load(const string & filename)
232 {
233         m_objects.bounds.clear();
234         m_count = 0;
235         if (filename == "")
236         {
237                 Debug("Loaded empty document.");
238                 return;
239         }
240         Debug("Loading document from file \"%s\"", filename.c_str());
241         FILE * file = fopen(filename.c_str(), "r");
242         if (file == NULL)
243                 Fatal("Couldn't open file \"%s\"", filename.c_str(), strerror(errno));
244
245         size_t read;
246
247         DocChunkTypes chunk_type;
248         uint32_t chunk_size;
249         while (ReadChunkHeader(file, chunk_type, chunk_size))
250         {
251                 switch(chunk_type)
252                 {
253                 case CT_NUMOBJS:
254                         read = fread(&m_count, sizeof(m_count), 1, file);
255                         if (read != 1)
256                                 Fatal("Failed to read number of objects!");
257                         Debug("Number of objects: %u", ObjectCount());
258                         break;
259                 case CT_OBJTYPES:
260                         Debug("Object types...");
261                         LoadStructVector<ObjectType>(file, chunk_size/sizeof(ObjectType), m_objects.types);
262                         break;
263                 case CT_OBJBOUNDS:
264                         Debug("Object bounds...");
265                         LoadStructVector<Rect>(file, chunk_size/sizeof(Rect), m_objects.bounds);
266                         break;
267                 case CT_OBJINDICES:
268                         Debug("Object data indices...");
269                         LoadStructVector<unsigned>(file, chunk_size/sizeof(unsigned), m_objects.data_indices);
270                         break;
271                 case CT_OBJBEZIERS:
272                         Debug("Bezier data...");
273                         LoadStructVector<Bezier>(file, chunk_size/sizeof(Bezier), m_objects.beziers);
274                         break;
275                         
276                 case CT_OBJPATHS:
277                         Debug("Path data...");
278                         Warn("Not handled because lazy");
279                         break;
280                 }
281         }
282         Debug("Successfully loaded %u objects from \"%s\"", ObjectCount(), filename.c_str());
283 #ifndef QUADTREE_DISABLED
284         if (m_quadtree.root_id == QUADTREE_EMPTY)
285         {
286                 GenBaseQuadtree();
287         }
288 #endif
289 }
290
291 unsigned Document::AddPath(unsigned start_index, unsigned end_index, const Colour & fill)
292 {
293         Path path(m_objects, start_index, end_index, fill);
294         unsigned data_index = AddPathData(path);
295         Rect bounds = path.SolveBounds(m_objects);
296         unsigned result = Add(PATH, bounds,data_index);
297         return result;
298 }
299
300 /**
301  * Add a Bezier using Absolute coords
302  */
303 unsigned Document::AddBezier(const Bezier & bezier)
304 {
305         Rect bounds = bezier.SolveBounds();
306         Bezier data = bezier.ToRelative(bounds); // Relative
307         if (data.ToAbsolute(bounds) != bezier)
308         {
309                 Error("%s != %s", data.ToAbsolute(Rect(0,0,1,1)).Str().c_str(),
310                         bezier.Str().c_str());
311                 Fatal("ToAbsolute on ToRelative does not give original Bezier");
312         }
313         unsigned index = AddBezierData(data);
314         return Add(BEZIER, bounds, index);
315 }
316
317 unsigned Document::Add(ObjectType type, const Rect & bounds, unsigned data_index)
318 {
319         m_objects.types.push_back(type);
320         m_objects.bounds.push_back(bounds);
321         m_objects.data_indices.push_back(data_index);
322         return (m_count++); // Why can't we just use the size of types or something?
323 }
324
325 unsigned Document::AddBezierData(const Bezier & bezier)
326 {
327         m_objects.beziers.push_back(bezier);
328         return m_objects.beziers.size()-1;
329 }
330
331 unsigned Document::AddPathData(const Path & path)
332 {
333         m_objects.paths.push_back(path);
334         return m_objects.paths.size()-1;
335 }
336
337 void Document::DebugDumpObjects()
338 {
339         Debug("Objects for Document %p are:", this);
340         for (unsigned id = 0; id < ObjectCount(); ++id)
341         {
342                 Debug("%u. \tType: %u\tBounds: %s", id, m_objects.types[id], m_objects.bounds[id].Str().c_str());
343         }
344 }
345
346 bool Document::operator==(const Document & equ) const
347 {
348         return (ObjectCount() == equ.ObjectCount() 
349                 && memcmp(m_objects.bounds.data(), equ.m_objects.bounds.data(), ObjectCount() * sizeof(Rect)) == 0
350                 && memcmp(m_objects.data_indices.data(), equ.m_objects.data_indices.data(), ObjectCount() * sizeof(unsigned)) == 0
351                 && memcmp(m_objects.beziers.data(), equ.m_objects.beziers.data(), m_objects.beziers.size() * sizeof(Bezier)) == 0);
352 }
353
354
355
356 // Behold my amazing tokenizing abilities
357 static string & GetToken(const string & d, string & token, unsigned & i, const string & delims = "()[],{}<>;:=")
358 {
359         token.clear();
360         while (i < d.size() && iswspace(d[i]))
361         {
362                 ++i;
363         }
364         
365         while (i < d.size())
366         {
367                 if (iswspace(d[i]) || strchr(delims.c_str(),d[i]) != NULL)
368                 {
369                         if (token.size() == 0 && !iswspace(d[i]))
370                         {
371                                 token += d[i++];
372                         }
373                         break;  
374                 }
375                 token += d[i++];
376         }
377         //Debug("Got token \"%s\"", token.c_str());
378         return token;
379 }
380
381 static void GetXYPair(const string & d, Real & x, Real & y, unsigned & i,const string & delims = "()[],{}<>;:=")
382 {
383         string token("");
384         while (GetToken(d, token, i, delims) == ",");
385         x = strtod(token.c_str(),NULL);
386         if (GetToken(d, token, i, delims) != ",")
387         {
388                 Fatal("Expected \",\" seperating x,y pair");
389         }
390         y = strtod(GetToken(d, token, i, delims).c_str(),NULL);
391 }
392
393 static bool GetKeyValuePair(const string & d, string & key, string & value, unsigned & i, const string & delims = "()[],{}<>;:=")
394 {
395         key = "";
396         string token;
397         while (GetToken(d, token, i, delims) == ":" || token == ";");
398         key = token;
399         if (GetToken(d, token, i, delims) != ":")
400         {
401                 Error("Expected \":\" seperating key:value pair");
402                 return false;
403         }
404         value = "";
405         GetToken(d, value, i, delims);
406         return true;
407 }
408
409 static void TransformXYPair(Real & x, Real & y, const SVGMatrix & transform)
410 {
411         Real x0(x);
412         x = transform.a * x + transform.c * y + transform.e;
413         y = transform.b * x0 + transform.d * y + transform.f;
414 }
415
416 void Document::ParseSVGTransform(const string & s, SVGMatrix & transform)
417 {
418         //Debug("Parsing transform %s", s.c_str());
419         string token;
420         string command;
421         unsigned i = 0;
422         
423         while (i < s.size())
424         {
425                 GetToken(s, command, i);
426                 if (command == "," || command == "" || command == ":")
427                 {
428                         if (i < s.size())
429                                 GetToken(s, command, i);
430                         else
431                                 return;
432                 }
433                 //Debug("Token is \"%s\"", command.c_str());
434         
435                 SVGMatrix delta = {1,0,0,0,1,0};
436         
437         
438                 assert(GetToken(s,token, i) == "(");
439                 if (command == "translate")
440                 {
441                         GetXYPair(s, delta.e, delta.f, i);
442                         assert(GetToken(s,token, i) == ")");    
443                 }
444                 else if (command == "matrix")
445                 {
446                         GetXYPair(s, delta.a, delta.b,i);
447                         GetXYPair(s, delta.c, delta.d,i);
448                         GetXYPair(s, delta.e, delta.f,i);
449                         assert(GetToken(s,token, i) == ")");    
450                 }
451                 else if (command == "scale")
452                 {
453                         delta.a = (strtod(GetToken(s,token,i).c_str(), NULL));
454                         GetToken(s, token, i);
455                         if (token == ",")
456                         {
457                                 delta.d = (strtod(GetToken(s,token,i).c_str(), NULL));
458                                 assert(GetToken(s, token, i) == ")");
459                         }
460                         else
461                         {
462                                 delta.d = delta.a;
463                                 assert(token == ")");
464                         }
465                         
466                 }
467                 else
468                 {
469                         Warn("Unrecognised transform \"%s\", using identity", command.c_str());
470                 }
471         
472                 //Debug("Old transform is {%f,%f,%f,%f,%f,%f}", transform.a, transform.b, transform.c, transform.d,transform.e,transform.f);
473                 //Debug("Delta transform is {%f,%f,%f,%f,%f,%f}", delta.a, delta.b, delta.c, delta.d,delta.e,delta.f);
474         
475                 SVGMatrix old(transform);
476                 transform.a = old.a * delta.a + old.c * delta.b;
477                 transform.c = old.a * delta.c + old.c * delta.d;
478                 transform.e = old.a * delta.e + old.c * delta.f + old.e;
479         
480                 transform.b = old.b * delta.a + old.d * delta.b;
481                 transform.d = old.b * delta.c + old.d * delta.d;
482                 transform.f = old.b * delta.e + old.d * delta.f + old.f;
483         
484                 //Debug("New transform is {%f,%f,%f,%f,%f,%f}", transform.a, transform.b, transform.c, transform.d,transform.e,transform.f);
485         }
486 }
487
488 void Document::ParseSVGNode(pugi::xml_node & root, SVGMatrix & parent_transform)
489 {
490         //Debug("Parse node <%s>", root.name());
491
492                 
493         for (pugi::xml_node child = root.first_child(); child; child = child.next_sibling())
494         {
495                 SVGMatrix transform(parent_transform);  
496                 pugi::xml_attribute attrib_trans = child.attribute("transform");
497                 if (!attrib_trans.empty())
498                 {
499                         ParseSVGTransform(attrib_trans.as_string(), transform);
500                 }
501                 
502                 
503                 
504                 if (strcmp(child.name(), "svg") == 0 || strcmp(child.name(),"g") == 0
505                         || strcmp(child.name(), "group") == 0)
506                 {
507                         
508                         ParseSVGNode(child, transform);
509                         continue;
510                 }
511                 else if (strcmp(child.name(), "path") == 0)
512                 {
513                         string d = child.attribute("d").as_string();
514                         //Debug("Path data attribute is \"%s\"", d.c_str());
515                         bool closed = false;
516                         pair<unsigned, unsigned> range = ParseSVGPathData(d, transform, closed);
517                         if (closed)
518                         {
519                                 Colour c(0,0,0,0);
520                                 string colour_str("");
521                                 map<string, string> style;
522                                 if (child.attribute("style"))
523                                 {
524                                         ParseSVGStyleData(child.attribute("style").as_string(), style);
525                                 }
526                                 
527                                 // Determine shading colour
528                                 if (child.attribute("fill"))
529                                 {
530                                         colour_str = child.attribute("fill").as_string();
531                                 }
532                                 else if (style.find("fill") != style.end())
533                                 {
534                                         colour_str = style["fill"];
535                                 }
536                                 if (colour_str == "red")
537                                         c = {1,0,0,1};
538                                 else if (colour_str == "blue")
539                                         c = {0,0,1,1};
540                                 else if (colour_str == "green")
541                                         c = {0,1,0,1};
542                                 else if (colour_str == "black")
543                                         c = {0,0,0,1};
544                                 else if (colour_str == "white")
545                                         c = {1,1,1,1};
546                                 else if (colour_str.size() == 7 && colour_str[0] == '#')
547                                 {
548                                         Debug("Parse colour string: \"%s\"", colour_str.c_str());
549                                         char comp[2] = {colour_str[1], colour_str[2]};
550                                         c.r = Real(strtoul(comp, NULL, 16))/Real(255);
551                                         comp[0] = colour_str[3]; comp[1] = colour_str[4];
552                                         c.g = Real(strtoul(comp, NULL, 16))/Real(255);
553                                         comp[0] = colour_str[5]; comp[1] = colour_str[6];
554                                         c.b = Real(strtoul(comp, NULL, 16))/Real(255);
555                                         c.a = 1;
556                                         Debug("Colour is: %f, %f, %f, %f", Float(c.r), Float(c.g), Float(c.b), Float(c.a));
557                                 }
558                                 
559                                 // Determin shading alpha
560                                 if (child.attribute("fill-opacity"))
561                                 {
562                                         c.a = child.attribute("fill-opacity").as_float();
563                                 }
564                                 else if (style.find("fill-opacity") != style.end())
565                                 {
566                                         c.a = strtod(style["fill-opacity"].c_str(), NULL);
567                                 }
568                                 
569                                         Debug("fill-opacity is %f", Float(c.a));
570                                 AddPath(range.first, range.second, c);
571                         }
572                         
573                 }
574                 else if (strcmp(child.name(), "line") == 0)
575                 {
576                         Real x0(child.attribute("x1").as_float());
577                         Real y0(child.attribute("y1").as_float());
578                         Real x1(child.attribute("x2").as_float());
579                         Real y1(child.attribute("y2").as_float());
580                         TransformXYPair(x0,y0,transform);
581                         TransformXYPair(x1,y1,transform);
582                         AddBezier(Bezier(x0,y0,x1,y1,x1,y1,x1,y1));
583                 }
584                 else if (strcmp(child.name(), "rect") == 0)
585                 {
586                         Real coords[4];
587                         const char * attrib_names[] = {"x", "y", "width", "height"};
588                         for (size_t i = 0; i < 4; ++i)
589                                 coords[i] = child.attribute(attrib_names[i]).as_float();
590                         
591                         Real x2(coords[0]+coords[2]);
592                         Real y2(coords[1]+coords[3]);
593                         TransformXYPair(coords[0],coords[1],transform); // x, y, transform
594                         TransformXYPair(x2,y2,transform);
595                         coords[2] = x2 - coords[0];
596                         coords[3] = y2 - coords[1];
597                         
598                         bool outline = !(child.attribute("fill") && strcmp(child.attribute("fill").as_string(),"none") != 0);
599                         Add(outline?RECT_OUTLINE:RECT_FILLED, Rect(coords[0], coords[1], coords[2], coords[3]),0);
600                 }
601                 else if (strcmp(child.name(), "circle") == 0)
602                 {
603                         Real cx = child.attribute("cx").as_float();
604                         Real cy = child.attribute("cy").as_float();
605                         Real r = child.attribute("r").as_float();
606                         
607                         Real x = (cx - r);
608                         Real y = (cy - r);
609                         TransformXYPair(x,y,transform);
610                         Real w = Real(2)*r*transform.a; // width scales
611                         Real h = Real(2)*r*transform.d; // height scales
612                         
613                         
614                         Rect rect(x,y,w,h);
615                         Add(CIRCLE_FILLED, rect,0);
616                         Debug("Added Circle %s", rect.Str().c_str());                   
617                 }
618                 else if (strcmp(child.name(), "text") == 0)
619                 {
620                         Real x = child.attribute("x").as_float();
621                         Real y = child.attribute("y").as_float();
622                         TransformXYPair(x,y,transform);
623                         Debug("Add text \"%s\"", child.child_value());
624                         AddText(child.child_value(), 0.05, x, y);
625                 }
626         }
627 }
628
629 void Document::ParseSVGStyleData(const string & style, map<string, string> & results)
630 {
631         unsigned i = 0;
632         string key;
633         string value;
634         while (i < style.size() && GetKeyValuePair(style, key, value, i))
635         {
636                 results[key] = value;
637         }
638 }
639
640 /**
641  * Parse an SVG string into a rectangle
642  */
643 void Document::ParseSVG(const string & input, const Rect & bounds)
644 {
645         using namespace pugi;
646         
647         xml_document doc_xml;
648         xml_parse_result result = doc_xml.load(input.c_str());
649         
650         if (!result)
651                 Error("Couldn't parse SVG input - %s", result.description());
652                 
653         Debug("Loaded XML - %s", result.description());
654         SVGMatrix transform = {bounds.w, 0,bounds.x, 0,bounds.h,bounds.y};
655         ParseSVGNode(doc_xml, transform);
656 }
657
658 /**
659  * Load an SVG into a rectangle
660  */
661 void Document::LoadSVG(const string & filename, const Rect & bounds)
662 {
663         using namespace pugi;
664         
665         xml_document doc_xml;
666         ifstream input(filename.c_str(), ios_base::in);
667         xml_parse_result result = doc_xml.load(input);
668         
669         if (!result)
670                 Error("Couldn't load \"%s\" - %s", filename.c_str(), result.description());
671                 
672         Debug("Loaded XML - %s", result.description());
673         
674         input.close();
675                                                 // a c e, b d f
676         SVGMatrix transform = {bounds.w, 0,bounds.x, 0,bounds.h,bounds.y};
677         ParseSVGNode(doc_xml, transform);
678 }
679
680
681
682 // Fear the wrath of the tokenizing svg data
683 // Seriously this isn't really very DOM-like at all is it?
684 pair<unsigned, unsigned> Document::ParseSVGPathData(const string & d, const SVGMatrix & transform, bool & closed)
685 {
686         closed = false;
687         Real x[4] = {0,0,0,0};
688         Real y[4] = {0,0,0,0};
689         
690         string token("");
691         string command("m");
692         
693         Real x0(0);
694         Real y0(0);
695         
696         unsigned i = 0;
697         unsigned prev_i = 0;
698         
699         bool start = false;
700         
701
702         static string delims("()[],{}<>;:=LlHhVvmMqQzZcC");
703
704         pair<unsigned, unsigned> range(m_count, m_count);
705         
706         while (i < d.size() && GetToken(d, token, i, delims).size() > 0)
707         {
708                 if (isalpha(token[0]))
709                         command = token;
710                 else
711                 {
712                         i = prev_i; // hax
713                         if(command == "")
714                                 command = "L";
715                 }
716                 
717                 bool relative = islower(command[0]);
718                         
719                 if (command == "m" || command == "M")
720                 {
721                         //Debug("Construct moveto command");
722                         Real dx = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
723                         assert(GetToken(d,token,i,delims) == ",");
724                         Real dy = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
725                         
726                         x[0] = (relative) ? x[0] + dx : dx;
727                         y[0] = (relative) ? y[0] + dy : dy;
728                         
729                         x0 = x[0];
730                         y0 = y[0];
731                         //Debug("mmoveto %f,%f", Float(x[0]),Float(y[0]));
732                         command = (command == "m") ? "l" : "L";
733                 }
734                 else if (command == "c" || command == "C" || command == "q" || command == "Q")
735                 {
736                         //Debug("Construct curveto command");
737                         Real dx = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
738                         assert(GetToken(d,token,i,delims) == ",");
739                         Real dy = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
740                         
741                         x[1] = (relative) ? x[0] + dx : dx;
742                         y[1] = (relative) ? y[0] + dy : dy;
743                         
744                         dx = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
745                         assert(GetToken(d,token,i,delims) == ",");
746                         dy = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
747                         
748                         x[2] = (relative) ? x[0] + dx : dx;
749                         y[2] = (relative) ? y[0] + dy : dy;
750                         
751                         if (command != "q" && command != "Q")
752                         {
753                                 dx = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
754                                 assert(GetToken(d,token,i,delims) == ",");
755                                 dy = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
756                                 x[3] = (relative) ? x[0] + dx : dx;
757                                 y[3] = (relative) ? y[0] + dy : dy;
758                         }
759                         else
760                         {
761                                 x[3] = x[2];
762                                 y[3] = y[2];
763                                 Real old_x1(x[1]), old_y1(y[1]);
764                                 x[1] = x[0] + Real(2) * (old_x1 - x[0])/ Real(3);
765                                 y[1] = y[0] + Real(2) * (old_y1 - y[0])/ Real(3);
766                                 x[2] = x[3] + Real(2) * (old_x1 - x[3])/ Real(3);
767                                 y[2] = y[3] + Real(2) * (old_y1 - y[3])/ Real(3);
768                         }
769                         
770                         Real x3(x[3]);
771                         Real y3(y[3]);
772                         for (int j = 0; j < 4; ++j)
773                                 TransformXYPair(x[j],y[j], transform);
774
775                         range.second = AddBezier(Bezier(x[0],y[0],x[1],y[1],x[2],y[2],x[3],y[3]));
776                         
777                         //Debug("[%u] curveto %f,%f %f,%f %f,%f", index, Float(x[1]),Float(y[1]),Float(x[2]),Float(y[2]),Float(x[3]),Float(y[3]));
778                         
779                         x[0] = x3;
780                         y[0] = y3;
781
782                         
783                 }
784                 else if (command == "l" || command == "L" || command == "h" || command == "H" || command == "v" || command == "V")
785                 {
786                         //Debug("Construct lineto command, relative %d", relative);
787                 
788                         Real dx = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
789                         Real dy;
790                         if (command == "l" || command == "L")
791                         {
792                                 assert(GetToken(d,token,i,delims) == ",");
793                                 dy = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
794                         }
795                         else if (command == "v" || command == "V")
796                         {
797                                 swap(dx,dy);
798                         }
799                         
800                         x[1] = (relative) ? x[0] + dx : dx;
801                         y[1] = (relative) ? y[0] + dy : dy;
802                         if (command == "v" || command == "V")
803                         {
804                                 x[1] = x[0];
805                         }
806                         else if (command == "h" || command == "H")
807                         {
808                                 y[1] = y[0];
809                         }
810                         
811                         Real x1(x[1]);
812                         Real y1(y[1]);
813                         
814                         TransformXYPair(x[0],y[0],transform);
815                         TransformXYPair(x[1],y[1],transform);
816
817
818                         range.second = AddBezier(Bezier(x[0],y[0],x[1],y[1],x[1],y[1],x[1],y[1]));
819                         
820                         //Debug("[%u] lineto %f,%f %f,%f", index, Float(x[0]),Float(y[0]),Float(x[1]),Float(y[1]));
821                         
822                         x[0] = x1;
823                         y[0] = y1;
824
825                 }
826                 else if (command == "z" || command == "Z")
827                 {
828                         //Debug("Construct returnto command");
829                         x[1] = x0;
830                         y[1] = y0;
831                         x[2] = x0;
832                         y[2] = y0;
833                         x[3] = x0;
834                         y[3] = y0;
835                         
836                         Real x3(x[3]);
837                         Real y3(y[3]);
838                         for (int j = 0; j < 4; ++j)
839                                 TransformXYPair(x[j],y[j], transform);
840
841                         range.second = AddBezier(Bezier(x[0],y[0],x[1],y[1],x[2],y[2],x[3],y[3]));
842                         //Debug("[%u] returnto %f,%f %f,%f", index, Float(x[0]),Float(y[0]),Float(x[1]),Float(y[1]));
843                         
844                         x[0] = x3;
845                         y[0] = y3;
846                         command = "m";
847                         closed = true;
848                 }
849                 else
850                 {
851                         Warn("Unrecognised command \"%s\", set to \"m\"", command.c_str());
852                         command = "m";
853                 }
854                 
855                 if (!start)
856                 {
857                         x0 = x[0];
858                         y0 = y[0];
859                         start = true;
860                 }
861                 prev_i = i;
862         }
863         return range;
864 }
865
866 void Document::SetFont(const string & font_filename)
867 {
868         if (m_font_data != NULL)
869         {
870                 free(m_font_data);
871         }
872         
873         FILE *font_file = fopen(font_filename.c_str(), "rb");
874         fseek(font_file, 0, SEEK_END);
875         size_t font_file_size = ftell(font_file);
876         fseek(font_file, 0, SEEK_SET);
877         m_font_data = (unsigned char*)malloc(font_file_size);
878         size_t read = fread(m_font_data, 1, font_file_size, font_file);
879         if (read != font_file_size)
880         {
881                 Fatal("Failed to read font data from \"%s\" - Read %u bytes expected %u - %s", font_filename.c_str(), read, font_file_size, strerror(errno));
882         }
883         fclose(font_file);
884         stbtt_InitFont(&m_font, m_font_data, 0);
885 }
886
887 void Document::AddText(const string & text, Real scale, Real x, Real y)
888 {
889         if (m_font_data == NULL)
890         {
891                 Warn("No font loaded");
892                 return;
893         }
894                 
895         Real x0(x);
896         //Real y0(y);
897         int ascent = 0, descent = 0, line_gap = 0;
898         stbtt_GetFontVMetrics(&m_font, &ascent, &descent, &line_gap);
899         float font_scale = scale / (float)(ascent - descent);
900         Real y_advance = Real(font_scale) * Real(ascent - descent + line_gap);
901         for (unsigned i = 0; i < text.size(); ++i)
902         {
903                 if (text[i] == '\n')
904                 {
905                         y += y_advance;
906                         x = x0;
907                 }
908                 if (!isprint(text[i]))
909                         continue;
910                         
911                 int advance_width = 0, left_side_bearing = 0, kerning = 0;
912                 stbtt_GetCodepointHMetrics(&m_font, text[i], &advance_width, &left_side_bearing);
913                 if (i >= 1)
914                 {
915                         kerning = stbtt_GetCodepointKernAdvance(&m_font, text[i-1], text[i]);
916                 }
917                 x += Real(font_scale) * Real(kerning);
918                 AddFontGlyphAtPoint(&m_font, text[i], font_scale, x, y);
919                 x += Real(font_scale) * Real(advance_width);
920         }
921 }
922
923 void Document::AddFontGlyphAtPoint(stbtt_fontinfo *font, int character, Real scale, Real x, Real y)
924 {
925         int glyph_index = stbtt_FindGlyphIndex(font, character);
926
927         // Check if there is actully a glyph to render.
928         if (stbtt_IsGlyphEmpty(font, glyph_index))
929         {
930                 return;
931         }
932
933         stbtt_vertex *instructions;
934         int num_instructions = stbtt_GetGlyphShape(font, glyph_index, &instructions);
935
936         Real current_x(0), current_y(0);
937         unsigned start_index = m_count;
938         unsigned end_index = m_count;
939         for (int i = 0; i < num_instructions; ++i)
940         {
941                 // TTF uses 16-bit signed ints for coordinates:
942                 // with the y-axis inverted compared to us.
943                 // Convert and scale any data.
944                 Real inst_x = Real(instructions[i].x)*scale;
945                 Real inst_y = Real(instructions[i].y)*-scale;
946                 Real inst_cx = Real(instructions[i].cx)*scale;
947                 Real inst_cy = Real(instructions[i].cy)*-scale;
948                 Real old_x(current_x), old_y(current_y);
949                 current_x = inst_x;
950                 current_y = inst_y;
951                 
952                 switch(instructions[i].type)
953                 {
954                 // Move To
955                 case STBTT_vmove:
956                         break;
957                 // Line To
958                 case STBTT_vline:
959                         end_index = AddBezier(Bezier(old_x + x, old_y + y, old_x + x, old_y + y, current_x + x, current_y + y, current_x + x, current_y + y));
960                         break;
961                 // Quadratic Bezier To:
962                 case STBTT_vcurve:
963                         // Quadratic -> Cubic:
964                         // - Endpoints are the same.
965                         // - cubic1 = quad0+(2/3)*(quad1-quad0)
966                         // - cubic2 = quad2+(2/3)*(quad1-quad2)
967                         end_index = AddBezier(Bezier(old_x + x, old_y + y, old_x + Real(2)*(inst_cx-old_x)/Real(3) + x, old_y + Real(2)*(inst_cy-old_y)/Real(3) + y,
968                                                 current_x + Real(2)*(inst_cx-current_x)/Real(3) + x, current_y + Real(2)*(inst_cy-current_y)/Real(3) + y, current_x + x, current_y + y));
969                         break;
970                 }
971         }
972         
973         if (start_index < m_count && end_index < m_count)
974         {
975                 AddPath(start_index, end_index);
976         }
977         Debug("Added Glyph \"%c\" at %f %f, scale %f", (char)character, Float(x), Float(y), Float(scale));
978
979         stbtt_FreeShape(font, instructions);
980 }

UCC git Repository :: git.ucc.asn.au