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