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         default:
 134                 Debug("Adding %s -> %s", m_objects.bounds[object_id].Str().c_str(), TransformToQuadChild(m_objects.bounds[object_id], type).Str().c_str());
 135                 m_objects.bounds.push_back(TransformToQuadChild(m_objects.bounds[object_id], type));
 136                 m_objects.types.push_back(m_objects.types[object_id]);
 137                 m_objects.data_indices.push_back(m_objects.data_indices[object_id]);
 138                 return 1;
 139         }
 140         return 0;
 141 }
 142 QuadTreeIndex Document::GenQuadChild(QuadTreeIndex parent, QuadTreeNodeChildren type)
 143 {
 144         QuadTreeIndex new_index = m_quadtree.nodes.size();
 145         m_quadtree.nodes.push_back(QuadTreeNode{QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, parent, type, 0, 0});
 146
 147         m_quadtree.nodes[new_index].object_begin = m_objects.bounds.size();
 148         for (unsigned i = m_quadtree.nodes[parent].object_begin; i < m_quadtree.nodes[parent].object_end; ++i)
 149         {
 150                 if (IntersectsQuadChild(m_objects.bounds[i], type))
 151                 {
 152                         m_count += ClipObjectToQuadChild(i, type);
 153                 }
 154         }
 155         m_quadtree.nodes[new_index].object_end = m_objects.bounds.size();
 156         switch (type)
 157         {
 158                 case QTC_TOP_LEFT:
 159                         m_quadtree.nodes[parent].top_left = new_index;
 160                         break;
 161                 case QTC_TOP_RIGHT:
 162                         m_quadtree.nodes[parent].top_right = new_index;
 163                         break;
 164                 case QTC_BOTTOM_LEFT:
 165                         m_quadtree.nodes[parent].bottom_left = new_index;
 166                         break;
 167                 case QTC_BOTTOM_RIGHT:
 168                         m_quadtree.nodes[parent].bottom_right = new_index;
 169                         break;
 170                 default:
 171                         Fatal("Tried to add a QuadTree child of invalid type!");
 172         }
 173         return new_index;
 174 }
 175
 176 // Reparent a quadtree node, making it the "type" child of a new node.
 177 QuadTreeIndex Document::GenQuadParent(QuadTreeIndex child, QuadTreeNodeChildren type)
 178 {
 179         QuadTreeIndex new_index = m_quadtree.nodes.size();
 180         m_quadtree.nodes.push_back(QuadTreeNode{QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, QUADTREE_EMPTY, -1, QTC_UNKNOWN, 0, 0});
 181
 182         m_quadtree.nodes[new_index].object_begin = m_objects.bounds.size();
 183         for (unsigned i = m_quadtree.nodes[child].object_begin; i < m_quadtree.nodes[child].object_end; ++i)
 184         {
 185                 m_objects.bounds.push_back(TransformFromQuadChild(m_objects.bounds[i], type));
 186                 m_objects.types.push_back(m_objects.types[i]);
 187                 m_objects.data_indices.push_back(m_objects.data_indices[i]);
 188                 m_count++;
 189         }
 190         m_quadtree.nodes[new_index].object_end = m_objects.bounds.size();
 191         switch (type)
 192         {
 193                 case QTC_TOP_LEFT:
 194                         m_quadtree.nodes[new_index].top_left = child;
 195                         break;
 196                 case QTC_TOP_RIGHT:
 197                         m_quadtree.nodes[new_index].top_right = child;
 198                         break;
 199                 case QTC_BOTTOM_LEFT:
 200                         m_quadtree.nodes[new_index].bottom_left = child;
 201                         break;
 202                 case QTC_BOTTOM_RIGHT:
 203                         m_quadtree.nodes[new_index].bottom_right = child;
 204                         break;
 205                 default:
 206                         Fatal("Tried to add a QuadTree child of invalid type!");
 207         }
 208         return new_index;
 209 }
 210
 211 #endif
 212
 213 void Document::Load(const string & filename)
 214 {
 215         m_objects.bounds.clear();
 216         m_count = 0;
 217         if (filename == "")
 218         {
 219                 Debug("Loaded empty document.");
 220                 return;
 221         }
 222         Debug("Loading document from file \"%s\"", filename.c_str());
 223         FILE * file = fopen(filename.c_str(), "r");
 224         if (file == NULL)
 225                 Fatal("Couldn't open file \"%s\"", filename.c_str(), strerror(errno));
 226
 227         size_t read;
 228
 229         DocChunkTypes chunk_type;
 230         uint32_t chunk_size;
 231         while (ReadChunkHeader(file, chunk_type, chunk_size))
 232         {
 233                 switch(chunk_type)
 234                 {
 235                 case CT_NUMOBJS:
 236                         read = fread(&m_count, sizeof(m_count), 1, file);
 237                         if (read != 1)
 238                                 Fatal("Failed to read number of objects!");
 239                         Debug("Number of objects: %u", ObjectCount());
 240                         break;
 241                 case CT_OBJTYPES:
 242                         Debug("Object types...");
 243                         LoadStructVector<ObjectType>(file, chunk_size/sizeof(ObjectType), m_objects.types);
 244                         break;
 245                 case CT_OBJBOUNDS:
 246                         Debug("Object bounds...");
 247                         LoadStructVector<Rect>(file, chunk_size/sizeof(Rect), m_objects.bounds);
 248                         break;
 249                 case CT_OBJINDICES:
 250                         Debug("Object data indices...");
 251                         LoadStructVector<unsigned>(file, chunk_size/sizeof(unsigned), m_objects.data_indices);
 252                         break;
 253                 case CT_OBJBEZIERS:
 254                         Debug("Bezier data...");
 255                         LoadStructVector<Bezier>(file, chunk_size/sizeof(Bezier), m_objects.beziers);
 256                         break;
 257
 258                 case CT_OBJGROUPS:
 259                         Debug("Group data...");
 260                         Warn("Not handled because lazy");
 261                         break;
 262                 }
 263         }
 264         Debug("Successfully loaded %u objects from \"%s\"", ObjectCount(), filename.c_str());
 265 #ifndef QUADTREE_DISABLED
 266         if (m_quadtree.root_id == QUADTREE_EMPTY)
 267         {
 268                 GenBaseQuadtree();
 269         }
 270 #endif
 271 }
 272
 273 unsigned Document::AddGroup(unsigned start_index, unsigned end_index)
 274 {
 275         //TODO: Set bounds rect?
 276         unsigned result = Add(GROUP, Rect(0,0,1,1),0);
 277         m_objects.groups[m_count-1].first = start_index;
 278         m_objects.groups[m_count-1].second = end_index;
 279         return result;
 280 }
 281
 282 unsigned Document::AddBezier(const Bezier & bezier)
 283 {
 284         unsigned index = AddBezierData(bezier);
 285         return Add(BEZIER, bezier.SolveBounds(), index);
 286 }
 287
 288 unsigned Document::Add(ObjectType type, const Rect & bounds, unsigned data_index)
 289 {
 290         m_objects.types.push_back(type);
 291         m_objects.bounds.push_back(bounds);
 292         m_objects.data_indices.push_back(data_index);
 293         m_objects.groups.push_back(pair<unsigned, unsigned>(data_index, data_index));
 294         return (m_count++); // Why can't we just use the size of types or something?
 295 }
 296
 297 unsigned Document::AddBezierData(const Bezier & bezier)
 298 {
 299         m_objects.beziers.push_back(bezier);
 300         return m_objects.beziers.size()-1;
 301 }
 302
 303
 304 void Document::DebugDumpObjects()
 305 {
 306         Debug("Objects for Document %p are:", this);
 307         for (unsigned id = 0; id < ObjectCount(); ++id)
 308         {
 309                 Debug("%u. \tType: %u\tBounds: %s", id, m_objects.types[id], m_objects.bounds[id].Str().c_str());
 310         }
 311 }
 312
 313 bool Document::operator==(const Document & equ) const
 314 {
 315         return (ObjectCount() == equ.ObjectCount()
 316                 && memcmp(m_objects.bounds.data(), equ.m_objects.bounds.data(), ObjectCount() * sizeof(Rect)) == 0
 317                 && memcmp(m_objects.data_indices.data(), equ.m_objects.data_indices.data(), ObjectCount() * sizeof(unsigned)) == 0
 318                 && memcmp(m_objects.beziers.data(), equ.m_objects.beziers.data(), m_objects.beziers.size() * sizeof(Bezier)) == 0);
 319 }
 320
 321
 322
 323 // Behold my amazing tokenizing abilities
 324 static string & GetToken(const string & d, string & token, unsigned & i, const string & delims = "()[],{}<>;:=")
 325 {
 326         token.clear();
 327         while (i < d.size() && iswspace(d[i]))
 328         {
 329                 ++i;
 330         }
 331
 332         while (i < d.size())
 333         {
 334                 if (iswspace(d[i]) || strchr(delims.c_str(),d[i]) != NULL)
 335                 {
 336                         if (token.size() == 0 && !iswspace(d[i]))
 337                         {
 338                                 token += d[i++];
 339                         }
 340                         break;
 341                 }
 342                 token += d[i++];
 343         }
 344         //Debug("Got token \"%s\"", token.c_str());
 345         return token;
 346 }
 347
 348 static void GetXYPair(const string & d, Real & x, Real & y, unsigned & i,const string & delims = "()[],{}<>;:=")
 349 {
 350         string token("");
 351         while (GetToken(d, token, i, delims) == ",");
 352         x = strtod(token.c_str(),NULL);
 353         if (GetToken(d, token, i, delims) != ",")
 354         {
 355                 Fatal("Expected \",\" seperating x,y pair");
 356         }
 357         y = strtod(GetToken(d, token, i, delims).c_str(),NULL);
 358 }
 359
 360 static void TransformXYPair(Real & x, Real & y, const SVGMatrix & transform)
 361 {
 362         Real x0(x);
 363         x = transform.a * x + transform.c * y + transform.e;
 364         y = transform.b * x0 + transform.d * y + transform.f;
 365 }
 366
 367 void Document::ParseSVGTransform(const string & s, SVGMatrix & transform)
 368 {
 369         Debug("Parsing transform %s", s.c_str());
 370         string token;
 371         string command;
 372         unsigned i = 0;
 373
 374         while (i < s.size())
 375         {
 376                 GetToken(s, command, i);
 377                 if (command == "," || command == "" || command == ":")
 378                 {
 379                         if (i < s.size())
 380                                 GetToken(s, command, i);
 381                         else
 382                                 return;
 383                 }
 384                 Debug("Token is \"%s\"", command.c_str());
 385
 386                 SVGMatrix delta = {1,0,0,0,1,0};
 387
 388
 389                 assert(GetToken(s,token, i) == "(");
 390                 if (command == "translate")
 391                 {
 392                         GetXYPair(s, delta.e, delta.f, i);
 393                         assert(GetToken(s,token, i) == ")");
 394                 }
 395                 else if (command == "matrix")
 396                 {
 397                         GetXYPair(s, delta.a, delta.b,i);
 398                         GetXYPair(s, delta.c, delta.d,i);
 399                         GetXYPair(s, delta.e, delta.f,i);
 400                         assert(GetToken(s,token, i) == ")");
 401                 }
 402                 else if (command == "scale")
 403                 {
 404                         delta.a = (strtod(GetToken(s,token,i).c_str(), NULL));
 405                         GetToken(s, token, i);
 406                         if (token == ",")
 407                         {
 408                                 delta.d = (strtod(GetToken(s,token,i).c_str(), NULL));
 409                                 assert(GetToken(s, token, i) == ")");
 410                         }
 411                         else
 412                         {
 413                                 delta.d = delta.a;
 414                                 assert(token == ")");
 415                         }
 416
 417                 }
 418                 else
 419                 {
 420                         Warn("Unrecognised transform \"%s\", using identity", command.c_str());
 421                 }
 422
 423                 Debug("Old transform is {%f,%f,%f,%f,%f,%f}", transform.a, transform.b, transform.c, transform.d,transform.e,transform.f);
 424                 Debug("Delta transform is {%f,%f,%f,%f,%f,%f}", delta.a, delta.b, delta.c, delta.d,delta.e,delta.f);
 425
 426                 SVGMatrix old(transform);
 427                 transform.a = old.a * delta.a + old.c * delta.b;
 428                 transform.c = old.a * delta.c + old.c * delta.d;
 429                 transform.e = old.a * delta.e + old.c * delta.f + old.e;
 430
 431                 transform.b = old.b * delta.a + old.d * delta.b;
 432                 transform.d = old.b * delta.c + old.d * delta.d;
 433                 transform.f = old.b * delta.e + old.d * delta.f + old.f;
 434
 435                 Debug("New transform is {%f,%f,%f,%f,%f,%f}", transform.a, transform.b, transform.c, transform.d,transform.e,transform.f);
 436         }
 437 }
 438
 439 void Document::ParseSVGNode(pugi::xml_node & root, SVGMatrix & parent_transform)
 440 {
 441         Debug("Parse node <%s>", root.name());
 442
 443
 444         for (pugi::xml_node child = root.first_child(); child; child = child.next_sibling())
 445         {
 446                 SVGMatrix transform(parent_transform);
 447                 pugi::xml_attribute attrib_trans = child.attribute("transform");
 448                 if (!attrib_trans.empty())
 449                 {
 450                         ParseSVGTransform(attrib_trans.as_string(), transform);
 451                 }
 452
 453                 if (strcmp(child.name(), "svg") == 0 || strcmp(child.name(),"g") == 0
 454                         || strcmp(child.name(), "group") == 0)
 455                 {
 456
 457                         ParseSVGNode(child, transform);
 458                         continue;
 459                 }
 460                 else if (strcmp(child.name(), "path") == 0)
 461                 {
 462                         string d = child.attribute("d").as_string();
 463                         Debug("Path data attribute is \"%s\"", d.c_str());
 464                         pair<unsigned, unsigned> range = ParseSVGPathData(d, transform);
 465                         AddGroup(range.first, range.second);
 466
 467                 }
 468                 else if (strcmp(child.name(), "line") == 0)
 469                 {
 470                         Real x0(child.attribute("x1").as_float());
 471                         Real y0(child.attribute("y1").as_float());
 472                         Real x1(child.attribute("x2").as_float());
 473                         Real y1(child.attribute("y2").as_float());
 474                         TransformXYPair(x0,y0,transform);
 475                         TransformXYPair(x1,y1,transform);
 476                         AddBezier(Bezier(x0,y0,x1,y1,x1,y1,x1,y1));
 477                 }
 478                 else if (strcmp(child.name(), "rect") == 0)
 479                 {
 480                         Real coords[4];
 481                         const char * attrib_names[] = {"x", "y", "width", "height"};
 482                         for (size_t i = 0; i < 4; ++i)
 483                                 coords[i] = child.attribute(attrib_names[i]).as_float();
 484
 485                         Real x2(coords[0]+coords[2]);
 486                         Real y2(coords[1]+coords[3]);
 487                         TransformXYPair(coords[0],coords[1],transform); // x, y, transform
 488                         TransformXYPair(x2,y2,transform);
 489                         coords[2] = x2 - coords[0];
 490                         coords[3] = y2 - coords[1];
 491
 492                         bool outline = !(child.attribute("fill") && strcmp(child.attribute("fill").as_string(),"none") != 0);
 493                         Add(outline?RECT_OUTLINE:RECT_FILLED, Rect(coords[0], coords[1], coords[2], coords[3]),0);
 494                 }
 495                 else if (strcmp(child.name(), "circle") == 0)
 496                 {
 497                         Real cx = child.attribute("cx").as_float();
 498                         Real cy = child.attribute("cy").as_float();
 499                         Real r = child.attribute("r").as_float();
 500
 501                         Real x = (cx - r);
 502                         Real y = (cy - r);
 503                         TransformXYPair(x,y,transform);
 504                         Real w = Real(2)*r*transform.a; // width scales
 505                         Real h = Real(2)*r*transform.d; // height scales
 506
 507
 508                         Rect rect(x,y,w,h);
 509                         Add(CIRCLE_FILLED, rect,0);
 510                         Debug("Added Circle %s", rect.Str().c_str());
 511                 }
 512                 else if (strcmp(child.name(), "text") == 0)
 513                 {
 514                         Real x = child.attribute("x").as_float();
 515                         Real y = child.attribute("y").as_float();
 516                         TransformXYPair(x,y,transform);
 517                         Debug("Add text \"%s\"", child.child_value());
 518                         AddText(child.child_value(), 0.05, x, y);
 519                 }
 520         }
 521 }
 522
 523 /**
 524  * Load an SVG into a rectangle
 525  */
 526 void Document::LoadSVG(const string & filename, const Rect & bounds)
 527 {
 528         using namespace pugi;
 529
 530         xml_document doc_xml;
 531         ifstream input(filename.c_str(), ios_base::in);
 532         xml_parse_result result = doc_xml.load(input);
 533
 534         if (!result)
 535                 Fatal("Couldn't load \"%s\" - %s", filename.c_str(), result.description());
 536
 537         Debug("Loaded XML - %s", result.description());
 538
 539         input.close();
 540                                                 // a c e, b d f
 541         SVGMatrix transform = {bounds.w, 0,bounds.x, 0,bounds.h,bounds.y};
 542         ParseSVGNode(doc_xml, transform);
 543 }
 544
 545
 546
 547 // Fear the wrath of the tokenizing svg data
 548 // Seriously this isn't really very DOM-like at all is it?
 549 pair<unsigned, unsigned> Document::ParseSVGPathData(const string & d, const SVGMatrix & transform)
 550 {
 551         Real x[4] = {0,0,0,0};
 552         Real y[4] = {0,0,0,0};
 553
 554         string token("");
 555         string command("m");
 556
 557         Real x0(0);
 558         Real y0(0);
 559
 560         unsigned i = 0;
 561         unsigned prev_i = 0;
 562
 563         bool start = false;
 564
 565
 566         static string delims("()[],{}<>;:=LlHhVvmMqQzZcC");
 567
 568         pair<unsigned, unsigned> range(m_count, m_count);
 569
 570         while (i < d.size() && GetToken(d, token, i, delims).size() > 0)
 571         {
 572                 if (isalpha(token[0]))
 573                         command = token;
 574                 else
 575                 {
 576                         i = prev_i; // hax
 577                         if(command == "")
 578                                 command = "L";
 579                 }
 580
 581                 bool relative = islower(command[0]);
 582
 583                 if (command == "m" || command == "M")
 584                 {
 585                         //Debug("Construct moveto command");
 586                         Real dx = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
 587                         assert(GetToken(d,token,i,delims) == ",");
 588                         Real dy = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
 589
 590                         x[0] = (relative) ? x[0] + dx : dx;
 591                         y[0] = (relative) ? y[0] + dy : dy;
 592
 593                         //Debug("mmoveto %f,%f", Float(x[0]),Float(y[0]));
 594                         command = (command == "m") ? "l" : "L";
 595                 }
 596                 else if (command == "c" || command == "C" || command == "q" || command == "Q")
 597                 {
 598                         //Debug("Construct curveto command");
 599                         Real dx = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
 600                         assert(GetToken(d,token,i,delims) == ",");
 601                         Real dy = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
 602
 603                         x[1] = (relative) ? x[0] + dx : dx;
 604                         y[1] = (relative) ? y[0] + dy : dy;
 605
 606                         dx = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
 607                         assert(GetToken(d,token,i,delims) == ",");
 608                         dy = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
 609
 610                         x[2] = (relative) ? x[0] + dx : dx;
 611                         y[2] = (relative) ? y[0] + dy : dy;
 612
 613                         if (command != "q" && command != "Q")
 614                         {
 615                                 dx = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
 616                                 assert(GetToken(d,token,i,delims) == ",");
 617                                 dy = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
 618                                 x[3] = (relative) ? x[0] + dx : dx;
 619                                 y[3] = (relative) ? y[0] + dy : dy;
 620                         }
 621                         else
 622                         {
 623                                 x[3] = x[2];
 624                                 y[3] = y[2];
 625                                 Real old_x1(x[1]), old_y1(y[1]);
 626                                 x[1] = x[0] + Real(2) * (old_x1 - x[0])/ Real(3);
 627                                 y[1] = y[0] + Real(2) * (old_y1 - y[0])/ Real(3);
 628                                 x[2] = x[3] + Real(2) * (old_x1 - x[3])/ Real(3);
 629                                 y[2] = y[3] + Real(2) * (old_y1 - y[3])/ Real(3);
 630                         }
 631
 632                         Real x3(x[3]);
 633                         Real y3(y[3]);
 634                         for (int j = 0; j < 4; ++j)
 635                                 TransformXYPair(x[j],y[j], transform);
 636
 637                         range.second = AddBezier(Bezier(x[0],y[0],x[1],y[1],x[2],y[2],x[3],y[3]));
 638
 639                         //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]));
 640
 641                         x[0] = x3;
 642                         y[0] = y3;
 643
 644
 645                 }
 646                 else if (command == "l" || command == "L" || command == "h" || command == "H" || command == "v" || command == "V")
 647                 {
 648                         Debug("Construct lineto command, relative %d", relative);
 649
 650                         Real dx = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
 651                         Real dy;
 652                         if (command == "l" || command == "L")
 653                         {
 654                                 assert(GetToken(d,token,i,delims) == ",");
 655                                 dy = Real(strtod(GetToken(d,token,i,delims).c_str(),NULL));
 656                         }
 657                         else if (command == "v" || command == "V")
 658                         {
 659                                 swap(dx,dy);
 660                         }
 661
 662                         x[1] = (relative) ? x[0] + dx : dx;
 663                         y[1] = (relative) ? y[0] + dy : dy;
 664                         if (command == "v" || command == "V")
 665                         {
 666                                 x[1] = x[0];
 667                         }
 668                         else if (command == "h" || command == "H")
 669                         {
 670                                 y[1] = y[0];
 671                         }
 672
 673                         Real x1(x[1]);
 674                         Real y1(y[1]);
 675
 676                         TransformXYPair(x[0],y[0],transform);
 677                         TransformXYPair(x[1],y[1],transform);
 678
 679
 680                         range.second = AddBezier(Bezier(x[0],y[0],x[1],y[1],x[1],y[1],x[1],y[1]));
 681
 682                         //Debug("[%u] lineto %f,%f %f,%f", index, Float(x[0]),Float(y[0]),Float(x[1]),Float(y[1]));
 683
 684                         x[0] = x1;
 685                         y[0] = y1;
 686
 687                 }
 688                 else if (command == "z" || command == "Z")
 689                 {
 690                         //Debug("Construct returnto command");
 691                         x[1] = x0;
 692                         y[1] = y0;
 693                         x[2] = x0;
 694                         y[2] = y0;
 695                         x[3] = x0;
 696                         y[3] = y0;
 697
 698                         Real x3(x[3]);
 699                         Real y3(y[3]);
 700                         for (int j = 0; j < 4; ++j)
 701                                 TransformXYPair(x[j],y[j], transform);
 702
 703                         range.second = AddBezier(Bezier(x[0],y[0],x[1],y[1],x[2],y[2],x[3],y[3]));
 704                         //Debug("[%u] returnto %f,%f %f,%f", index, Float(x[0]),Float(y[0]),Float(x[1]),Float(y[1]));
 705
 706                         x[0] = x3;
 707                         y[0] = y3;
 708                         command = "m";
 709                 }
 710                 else
 711                 {
 712                         Warn("Unrecognised command \"%s\", set to \"m\"", command.c_str());
 713                         command = "m";
 714                 }
 715
 716                 if (!start)
 717                 {
 718                         x0 = x[0];
 719                         y0 = y[0];
 720                         start = true;
 721                 }
 722                 prev_i = i;
 723         }
 724         return range;
 725 }
 726
 727 void Document::SetFont(const string & font_filename)
 728 {
 729         if (m_font_data != NULL)
 730         {
 731                 free(m_font_data);
 732         }
 733
 734         FILE *font_file = fopen("DejaVuSansMono.ttf", "rb");
 735         fseek(font_file, 0, SEEK_END);
 736         size_t font_file_size = ftell(font_file);
 737         fseek(font_file, 0, SEEK_SET);
 738         m_font_data = (unsigned char*)malloc(font_file_size);
 739         size_t read = fread(m_font_data, 1, font_file_size, font_file);
 740         if (read != font_file_size)
 741         {
 742                 Fatal("Failed to read font data from \"%s\" - Read %u bytes expected %u - %s", font_filename.c_str(), read, font_file_size, strerror(errno));
 743         }
 744         fclose(font_file);
 745         stbtt_InitFont(&m_font, m_font_data, 0);
 746 }
 747
 748 void Document::AddText(const string & text, Real scale, Real x, Real y)
 749 {
 750         if (m_font_data == NULL)
 751         {
 752                 Warn("No font loaded");
 753                 return;
 754         }
 755
 756         float font_scale = stbtt_ScaleForPixelHeight(&m_font, scale);
 757         Real x0(x);
 758         //Real y0(y);
 759         int ascent = 0, descent = 0, line_gap = 0;
 760         stbtt_GetFontVMetrics(&m_font, &ascent, &descent, &line_gap);
 761         Real y_advance = Real(font_scale) * Real(ascent - descent + line_gap);
 762         for (unsigned i = 0; i < text.size(); ++i)
 763         {
 764                 if (text[i] == '\n')
 765                 {
 766                         y += y_advance;
 767                         x = x0;
 768                 }
 769                 if (!isprint(text[i]))
 770                         continue;
 771
 772                 int advance_width = 0, left_side_bearing = 0, kerning = 0;
 773                 stbtt_GetCodepointHMetrics(&m_font, text[i], &advance_width, &left_side_bearing);
 774                 if (i > 1)
 775                 {
 776                         kerning = stbtt_GetCodepointKernAdvance(&m_font, text[i-1], text[i]);
 777                 }
 778                 x += Real(font_scale) * Real(left_side_bearing + kerning);
 779                 AddFontGlyphAtPoint(&m_font, text[i], font_scale, x, y);
 780                 x += Real(font_scale) * Real(advance_width);
 781         }
 782 }
 783
 784 void Document::AddFontGlyphAtPoint(stbtt_fontinfo *font, int character, Real scale, Real x, Real y)
 785 {
 786         int glyph_index = stbtt_FindGlyphIndex(font, character);
 787
 788         // Check if there is actully a glyph to render.
 789         if (stbtt_IsGlyphEmpty(font, glyph_index))
 790         {
 791                 return;
 792         }
 793
 794         stbtt_vertex *instructions;
 795         int num_instructions = stbtt_GetGlyphShape(font, glyph_index, &instructions);
 796
 797         Real current_x(0), current_y(0);
 798
 799         for (int i = 0; i < num_instructions; ++i)
 800         {
 801                 // TTF uses 16-bit signed ints for coordinates:
 802                 // with the y-axis inverted compared to us.
 803                 // Convert and scale any data.
 804                 Real inst_x = Real(instructions[i].x)*scale;
 805                 Real inst_y = Real(instructions[i].y)*-scale;
 806                 Real inst_cx = Real(instructions[i].cx)*scale;
 807                 Real inst_cy = Real(instructions[i].cy)*-scale;
 808                 Real old_x(current_x), old_y(current_y);
 809                 current_x = inst_x;
 810                 current_y = inst_y;
 811                 unsigned bezier_index;
 812                 switch(instructions[i].type)
 813                 {
 814                 // Move To
 815                 case STBTT_vmove:
 816                         break;
 817                 // Line To
 818                 case STBTT_vline:
 819                         bezier_index = AddBezierData(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));
 820                         Add(BEZIER,Rect(0,0,1,1),bezier_index);
 821                         break;
 822                 // Quadratic Bezier To:
 823                 case STBTT_vcurve:
 824                         // Quadratic -> Cubic:
 825                         // - Endpoints are the same.
 826                         // - cubic1 = quad0+(2/3)*(quad1-quad0)
 827                         // - cubic2 = quad2+(2/3)*(quad1-quad2)
 828                         bezier_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,
 829                                                 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));
 830                         break;
 831                 }
 832         }
 833
 834         stbtt_FreeShape(font, instructions);
 835 }