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