source: git/src/export.cc @ 27795245

/* export.cc
 * Export to CAD-like formats (DXF, Sketch, SVG, EPS) and also Compass PLT.
 */

/* Copyright (C) 1994-2004,2005,2006,2008,2010,2011,2012 Olly Betts
 * Copyright (C) 2004 John Pybus (SVG Output code)
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
 */

/* #define DEBUG_CAD3D */

#ifdef HAVE_CONFIG_H
#include <config.h>
#endif

#include "export.h"

#include "wx.h"
#include <wx/utils.h>
#include "exportfilter.h"
#include "hpgl.h"
#include "mainfrm.h"

#include <float.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>

#if defined(HAVE_GETPWUID) && !defined(__DJGPP__)
# include <pwd.h>
# include <sys/types.h>
# include <unistd.h>
#endif

#include "cmdline.h"
#include "debug.h"
#include "filename.h"
#include "hash.h"
#include "img.h"
#include "message.h"
#include "useful.h"

/* default values - can be overridden with --htext and --msize respectively */
#define TEXT_HEIGHT     0.6
#define MARKER_SIZE     0.8

#define GRID_SPACING    100

#define POINTS_PER_INCH 72.0
#define POINTS_PER_MM (POINTS_PER_INCH / MM_PER_INCH)

#define SQRT_2          1.41421356237309504880168872420969

#define LEGS 1
#define SURF 2
#define STNS 4
#define LABELS 8

static const char *layer_names[] = {
   NULL,
   "Legs",      // LEGS
   "Surface",   // SURF
   "Legs",      // LEGS|SURF
   "Stations",  // STNS
   NULL,
   NULL,
   NULL,
   "Labels"     // LABELS
};

/* default to DXF */
#define FMT_DEFAULT FMT_DXF

/* bounds */
static double min_x, min_y, min_z, max_x, max_y, max_z;

static double text_height; /* for station labels */
static double marker_size; /* for station markers */
static double grid; /* grid spacing (or 0 for no grid) */
static double scale = 500.0;
static double factor;
static const char *unit = "mm";
const double SVG_MARGIN = 5.0; // In units of "unit".

static const char *survey = NULL;

const int *
ExportFilter::passes() const
{
    static const int default_passes[] = { LEGS|SURF|STNS|LABELS, 0 };
    return default_passes;
}

class DXF : public ExportFilter {
  public:
    DXF() { }
    bool fopen(const char *fnm_out);
    void header(const char *);
    void line(const img_point *, const img_point *, bool, bool);
    void label(const img_point *, const char *, bool);
    void cross(const img_point *, bool);
    void footer();
};

bool
DXF::fopen(const char *fnm_out)
{
    // DXF gets written as text rather than binary.
    fh = ::fopen(fnm_out, "w");
    return (fh != NULL);
}

void
DXF::header(const char *)
{
   fprintf(fh, "0\nSECTION\n"
               "2\nHEADER\n");
   fprintf(fh, "9\n$EXTMIN\n"); /* lower left corner of drawing */
   fprintf(fh, "10\n%#-.6f\n", min_x); /* x */
   fprintf(fh, "20\n%#-.6f\n", min_y); /* y */
   fprintf(fh, "30\n%#-.6f\n", min_z); /* min z */
   fprintf(fh, "9\n$EXTMAX\n"); /* upper right corner of drawing */
   fprintf(fh, "10\n%#-.6f\n", max_x); /* x */
   fprintf(fh, "20\n%#-.6f\n", max_y); /* y */
   fprintf(fh, "30\n%#-.6f\n", max_z); /* max z */
   fprintf(fh, "9\n$PDMODE\n70\n3\n"); /* marker style as CROSS */
   fprintf(fh, "9\n$PDSIZE\n40\n%6.2f\n", marker_size); /* marker size */
   fprintf(fh, "0\nENDSEC\n");

   fprintf(fh, "0\nSECTION\n"
               "2\nTABLES\n");
   fprintf(fh, "0\nTABLE\n" /* Define CONTINUOUS and DASHED line types. */
               "2\nLTYPE\n"
               "70\n10\n"
               "0\nLTYPE\n"
               "2\nCONTINUOUS\n"
               "70\n64\n"
               "3\nContinuous\n"
               "72\n65\n"
               "73\n0\n"
               "40\n0.0\n"
               "0\nLTYPE\n"
               "2\nDASHED\n"
               "70\n64\n"
               "3\nDashed\n"
               "72\n65\n"
               "73\n2\n"
               "40\n2.5\n"
               "49\n1.25\n"
               "49\n-1.25\n"
               "0\nENDTAB\n");
   fprintf(fh, "0\nTABLE\n"
               "2\nLAYER\n");
   fprintf(fh, "70\n10\n"); /* max # off layers in this DXF file : 10 */
   /* First Layer: CentreLine */
   fprintf(fh, "0\nLAYER\n2\nCentreLine\n");
   fprintf(fh, "70\n64\n"); /* shows layer is referenced by entities */
   fprintf(fh, "62\n5\n"); /* color: kept the same used by SpeleoGen */
   fprintf(fh, "6\nCONTINUOUS\n"); /* linetype */
   /* Next Layer: Stations */
   fprintf(fh, "0\nLAYER\n2\nStations\n");
   fprintf(fh, "70\n64\n"); /* shows layer is referenced by entities */
   fprintf(fh, "62\n7\n"); /* color: kept the same used by SpeleoGen */
   fprintf(fh, "6\nCONTINUOUS\n"); /* linetype */
   /* Next Layer: Labels */
   fprintf(fh, "0\nLAYER\n2\nLabels\n");
   fprintf(fh, "70\n64\n"); /* shows layer is referenced by entities */
   fprintf(fh, "62\n7\n"); /* color: kept the same used by SpeleoGen */
   fprintf(fh, "6\nCONTINUOUS\n"); /* linetype */
   /* Next Layer: Surface */
   fprintf(fh, "0\nLAYER\n2\nSurface\n");
   fprintf(fh, "70\n64\n"); /* shows layer is referenced by entities */
   fprintf(fh, "62\n5\n"); /* color */
   fprintf(fh, "6\nDASHED\n"); /* linetype */
   /* Next Layer: SurfaceStations */
   fprintf(fh, "0\nLAYER\n2\nSurfaceStations\n");
   fprintf(fh, "70\n64\n"); /* shows layer is referenced by entities */
   fprintf(fh, "62\n7\n"); /* color */
   fprintf(fh, "6\nCONTINUOUS\n"); /* linetype */
   /* Next Layer: SurfaceLabels */
   fprintf(fh, "0\nLAYER\n2\nSurfaceLabels\n");
   fprintf(fh, "70\n64\n"); /* shows layer is referenced by entities */
   fprintf(fh, "62\n7\n"); /* color */
   fprintf(fh, "6\nCONTINUOUS\n"); /* linetype */
   if (grid > 0) {
      /* Next Layer: Grid */
      fprintf(fh, "0\nLAYER\n2\nGrid\n");
      fprintf(fh, "70\n64\n"); /* shows layer is referenced by entities */
      fprintf(fh, "62\n7\n"); /* color: kept the same used by SpeleoGen */
      fprintf(fh, "6\nCONTINUOUS\n"); /* linetype */
   }
   fprintf(fh, "0\nENDTAB\n"
               "0\nENDSEC\n");

   fprintf(fh, "0\nSECTION\n"
               "2\nENTITIES\n");

   if (grid > 0) {
      double x, y;
      x = floor(min_x / grid) * grid + grid;
      y = floor(min_y / grid) * grid + grid;
#ifdef DEBUG_CAD3D
      printf("x_min: %f  y_min: %f\n", x, y);
#endif
      while (x < max_x) {
         /* horizontal line */
         fprintf(fh, "0\nLINE\n");
         fprintf(fh, "8\nGrid\n"); /* Layer */
         fprintf(fh, "10\n%6.2f\n", x);
         fprintf(fh, "20\n%6.2f\n", min_y);
         fprintf(fh, "30\n0\n");
         fprintf(fh, "11\n%6.2f\n", x);
         fprintf(fh, "21\n%6.2f\n", max_y);
         fprintf(fh, "31\n0\n");
         x += grid;
      }
      while (y < max_y) {
         /* vertical line */
         fprintf(fh, "0\nLINE\n");
         fprintf(fh, "8\nGrid\n"); /* Layer */
         fprintf(fh, "10\n%6.2f\n", min_x);
         fprintf(fh, "20\n%6.2f\n", y);
         fprintf(fh, "30\n0\n");
         fprintf(fh, "11\n%6.2f\n", max_x);
         fprintf(fh, "21\n%6.2f\n", y);
         fprintf(fh, "31\n0\n");
         y += grid;
      }
   }
}

void
DXF::line(const img_point *p1, const img_point *p, bool fSurface, bool fPendingMove)
{
   fPendingMove = fPendingMove; /* unused */
   fprintf(fh, "0\nLINE\n");
   fprintf(fh, fSurface ? "8\nSurface\n" : "8\nCentreLine\n"); /* Layer */
   fprintf(fh, "10\n%6.2f\n", p1->x);
   fprintf(fh, "20\n%6.2f\n", p1->y);
   fprintf(fh, "30\n%6.2f\n", p1->z);
   fprintf(fh, "11\n%6.2f\n", p->x);
   fprintf(fh, "21\n%6.2f\n", p->y);
   fprintf(fh, "31\n%6.2f\n", p->z);
}

void
DXF::label(const img_point *p, const char *s, bool fSurface)
{
   /* write station label to dxf file */
   fprintf(fh, "0\nTEXT\n");
   fprintf(fh, fSurface ? "8\nSurfaceLabels\n" : "8\nLabels\n"); /* Layer */
   fprintf(fh, "10\n%6.2f\n", p->x);
   fprintf(fh, "20\n%6.2f\n", p->y);
   fprintf(fh, "30\n%6.2f\n", p->z);
   fprintf(fh, "40\n%6.2f\n", text_height);
   fprintf(fh, "1\n%s\n", s);
}

void
DXF::cross(const img_point *p, bool fSurface)
{
   /* write station marker to dxf file */
   fprintf(fh, "0\nPOINT\n");
   fprintf(fh, fSurface ? "8\nSurfaceStations\n" : "8\nStations\n"); /* Layer */
   fprintf(fh, "10\n%6.2f\n", p->x);
   fprintf(fh, "20\n%6.2f\n", p->y);
   fprintf(fh, "30\n%6.2f\n", p->z);
}

void
DXF::footer()
{
   fprintf(fh, "000\nENDSEC\n");
   fprintf(fh, "000\nEOF\n");
}

class Sketch : public ExportFilter {
  public:
    Sketch() { }
    const int * passes() const;
    void header(const char *);
    void start_pass(int layer);
    void line(const img_point *, const img_point *, bool, bool);
    void label(const img_point *, const char *, bool);
    void cross(const img_point *, bool);
    void footer();
};

const int *
Sketch::passes() const
{
    static const int sketch_passes[] = { LEGS|SURF, STNS, LABELS, 0 };
    return sketch_passes;
}

void
Sketch::header(const char *)
{
   fprintf(fh, "##Sketch 1 2\n"); /* Sketch file version */
   fprintf(fh, "document()\n");
   fprintf(fh, "layout((%.3f,%.3f),0)\n",
           (max_x - min_x) * factor, (max_y - min_y) * factor);
}

void
Sketch::start_pass(int layer)
{
   fprintf(fh, "layer('%s',1,1,0,0,(0,0,0))\n", layer_names[layer]);
}

void
Sketch::line(const img_point *p1, const img_point *p, bool fSurface, bool fPendingMove)
{
   fSurface = fSurface; /* unused */
   if (fPendingMove) {
       fprintf(fh, "b()\n");
       fprintf(fh, "bs(%.3f,%.3f,%.3f)\n", p1->x * factor, p1->y * factor, 0.0);
   }
   fprintf(fh, "bs(%.3f,%.3f,%.3f)\n", p->x * factor, p->y * factor, 0.0);
}

void
Sketch::label(const img_point *p, const char *s, bool fSurface)
{
   fSurface = fSurface; /* unused */
   fprintf(fh, "fp((0,0,0))\n");
   fprintf(fh, "le()\n");
   fprintf(fh, "Fn('Times-Roman')\n");
   fprintf(fh, "Fs(5)\n");
   fprintf(fh, "txt('");
   while (*s) {
      int ch = *s++;
      if (ch == '\'' || ch == '\\') PUTC('\\', fh);
      PUTC(ch, fh);
   }
   fprintf(fh, "',(%.3f,%.3f))\n", p->x * factor, p->y * factor);
}

void
Sketch::cross(const img_point *p, bool fSurface)
{
   fSurface = fSurface; /* unused */
   fprintf(fh, "b()\n");
   fprintf(fh, "bs(%.3f,%.3f,%.3f)\n",
           p->x * factor - MARKER_SIZE, p->y * factor - MARKER_SIZE, 0.0);
   fprintf(fh, "bs(%.3f,%.3f,%.3f)\n",
           p->x * factor + MARKER_SIZE, p->y * factor + MARKER_SIZE, 0.0);
   fprintf(fh, "bn()\n");
   fprintf(fh, "bs(%.3f,%.3f,%.3f)\n",
           p->x * factor + MARKER_SIZE, p->y * factor - MARKER_SIZE, 0.0);
   fprintf(fh, "bs(%.3f,%.3f,%.3f)\n",
           p->x * factor - MARKER_SIZE, p->y * factor + MARKER_SIZE, 0.0);
}

void
Sketch::footer(void)
{
   fprintf(fh, "guidelayer('Guide Lines',1,0,0,1,(0,0,1))\n");
   if (grid) {
      fprintf(fh, "grid((0,0,%.3f,%.3f),1,(0,0,1),'Grid')\n",
              grid * factor, grid * factor);
   }
}

typedef struct point {
   img_point p;
   const char *label;
   struct point *next;
} point;

#define HTAB_SIZE 0x2000

static point **htab;

static void
set_name(const img_point *p, const char *s)
{
   int hash;
   point *pt;
   union {
      char data[sizeof(int) * 3];
      int x[3];
   } u;

   u.x[0] = (int)(p->x * 100);
   u.x[1] = (int)(p->y * 100);
   u.x[2] = (int)(p->z * 100);
   hash = (hash_data(u.data, sizeof(int) * 3) & (HTAB_SIZE - 1));
   for (pt = htab[hash]; pt; pt = pt->next) {
      if (pt->p.x == p->x && pt->p.y == p->y && pt->p.z == p->z) {
         /* already got name for these coordinates */
         /* FIXME: what about multiple names for the same station? */
         return;
      }
   }

   pt = osnew(point);
   pt->label = osstrdup(s);
   pt->p = *p;
   pt->next = htab[hash];
   htab[hash] = pt;

   return;
}

static const char *
find_name(const img_point *p)
{
   int hash;
   point *pt;
   union {
      char data[sizeof(int) * 3];
      int x[3];
   } u;
   SVX_ASSERT(p);

   u.x[0] = (int)(p->x * 100);
   u.x[1] = (int)(p->y * 100);
   u.x[2] = (int)(p->z * 100);
   hash = (hash_data(u.data, sizeof(int) * 3) & (HTAB_SIZE - 1));
   for (pt = htab[hash]; pt; pt = pt->next) {
      if (pt->p.x == p->x && pt->p.y == p->y && pt->p.z == p->z)
         return pt->label;
   }
   return "?";
}

class SVG : public ExportFilter {
    bool to_close;
    bool close_g;

  public:
    SVG() : to_close(false), close_g(false) { }
    const int * passes() const;
    void header(const char *);
    void start_pass(int layer);
    void line(const img_point *, const img_point *, bool, bool);
    void label(const img_point *, const char *, bool);
    void cross(const img_point *, bool);
    void footer();
};

const int *
SVG::passes() const
{
    static const int svg_passes[] = { LEGS|SURF, LABELS, STNS, 0 };
    return svg_passes;
}

void
SVG::header(const char *)
{
   size_t i;
   htab = (point **)osmalloc(HTAB_SIZE * ossizeof(point *));
   for (i = 0; i < HTAB_SIZE; ++i) htab[i] = NULL;
   fprintf(fh, "<?xml version=\"1.0\" encoding=\"ISO-8859-1\"?>\n");
   double width = (max_x - min_x) * factor + SVG_MARGIN * 2;
   double height = (max_y - min_y) * factor + SVG_MARGIN * 2;
   fprintf(fh, "<svg version=\"1.1\" baseProfile=\"full\"\n"
               "xmlns=\"http://www.w3.org/2000/svg\"\n"
               "xmlns:xlink=\"http://www.w3.org/1999/xlink\"\n"
               "xmlns:ev=\"http://www.w3.org/2001/xml-events\"\n"
               "width=\"%.3f%s\" height=\"%.3f%s\"\n"
               "viewBox=\"0 0 %0.3f %0.3f\">\n",
           width, unit, height, unit, width, height);
   fprintf(fh, "<g transform=\"translate(%.3f %.3f)\">\n",
           SVG_MARGIN - min_x * factor, SVG_MARGIN + max_y * factor);
   to_close = false;
   close_g = false;
}

void
SVG::start_pass(int layer)
{
   if (to_close) {
      fprintf(fh, "\"/>\n");
      to_close = false;
   }
   if (close_g) {
      fprintf(fh, "</g>\n");
   }
   fprintf(fh, "<g id=\"%s\"", layer_names[layer]);
   if (layer & LEGS)
      fprintf(fh, " stroke=\"black\" fill=\"none\" stroke-width=\"0.4px\"");
   else if (layer & STNS)
      fprintf(fh, " stroke=\"black\" fill=\"none\" stroke-width=\"0.05px\"");
   else if (layer & LABELS)
      fprintf(fh, " font-size=\"%.3fem\"", text_height);
   fprintf(fh, ">\n");

   close_g = true;
}

void
SVG::line(const img_point *p1, const img_point *p, bool fSurface, bool fPendingMove)
{
   fSurface = fSurface; /* unused */
   if (fPendingMove) {
       if (to_close) {
           fprintf(fh, "\"/>\n");
       }
       fprintf(fh, "<path d=\"M%.3f %.3f", p1->x * factor, p1->y * -factor);
   }
   fprintf(fh, "L%.3f %.3f", p->x * factor, p->y * -factor);
   to_close = true;
}

void
SVG::label(const img_point *p, const char *s, bool fSurface)
{
   fSurface = fSurface; /* unused */
   fprintf(fh, "<text transform=\"translate(%.3f %.3f)\">",
           p->x * factor, p->y * -factor);
   fputs(s, fh);
   fputs("</text>\n", fh);
   set_name(p, s);
}

void
SVG::cross(const img_point *p, bool fSurface)
{
   fSurface = fSurface; /* unused */
   fprintf(fh, "<circle id=\"%s\" cx=\"%.3f\" cy=\"%.3f\" r=\"%.3f\"/>\n",
           find_name(p), p->x * factor, p->y * -factor, MARKER_SIZE * SQRT_2);
   fprintf(fh, "<path d=\"M%.3f %.3fL%.3f %.3fM%.3f %.3fL%.3f %.3f\"/>\n",
           p->x * factor - MARKER_SIZE, p->y * -factor - MARKER_SIZE,
           p->x * factor + MARKER_SIZE, p->y * -factor + MARKER_SIZE,
           p->x * factor + MARKER_SIZE, p->y * -factor - MARKER_SIZE,
           p->x * factor - MARKER_SIZE, p->y * -factor + MARKER_SIZE );
}

void
SVG::footer()
{
   if (to_close) {
      fprintf(fh, "\"/>\n");
      to_close = false;
   }
   if (close_g) {
      fprintf(fh, "</g>\n");
      close_g = false;
   }
   fprintf(fh, "</g>\n</svg>\n");
}

class PLT : public ExportFilter {
    string escaped;

    const char * find_name_plt(const img_point *p);

  public:
    PLT() { }
    const int * passes() const;
    void header(const char *);
    void line(const img_point *, const img_point *, bool, bool);
    void label(const img_point *, const char *, bool);
    void footer();
};

const int *
PLT::passes() const
{
    static const int plt_passes[] = { LABELS, LEGS|SURF, 0 };
    return plt_passes;
}

void
PLT::header(const char *title)
{
   size_t i;
   htab = (point **)osmalloc(HTAB_SIZE * ossizeof(point *));
   for (i = 0; i < HTAB_SIZE; ++i) htab[i] = NULL;
   /* Survex is E, N, Alt - PLT file is N, E, Alt */
   fprintf(fh, "Z %.3f %.3f %.3f %.3f %.3f %.3f\r\n",
           min_y / METRES_PER_FOOT, max_y / METRES_PER_FOOT,
           min_x / METRES_PER_FOOT, max_x / METRES_PER_FOOT,
           min_z / METRES_PER_FOOT, max_z / METRES_PER_FOOT);
   fprintf(fh, "N%s D 1 1 1 C%s\r\n", survey ? survey : "X",
           (title && title[0]) ? title : "X");
}

void
PLT::line(const img_point *p1, const img_point *p, bool fSurface, bool fPendingMove)
{
   fSurface = fSurface; /* unused */
   if (fPendingMove) {
       /* Survex is E, N, Alt - PLT file is N, E, Alt */
       fprintf(fh, "M %.3f %.3f %.3f ",
               p1->y / METRES_PER_FOOT, p1->x / METRES_PER_FOOT, p1->z / METRES_PER_FOOT);
       /* dummy passage dimensions are required to avoid compass bug */
       fprintf(fh, "S%s P -9 -9 -9 -9\r\n", find_name_plt(p1));
   }
   /* Survex is E, N, Alt - PLT file is N, E, Alt */
   fprintf(fh, "D %.3f %.3f %.3f ",
           p->y / METRES_PER_FOOT, p->x / METRES_PER_FOOT, p->z / METRES_PER_FOOT);
   /* dummy passage dimensions are required to avoid compass bug */
   fprintf(fh, "S%s P -9 -9 -9 -9\r\n", find_name_plt(p));
}

const char *
PLT::find_name_plt(const img_point *p)
{
    const char * s = find_name(p);
    escaped.resize(0);

    // PLT format can't handle spaces or control characters, so escape them
    // like in URLs (an arbitrary choice of escaping, but at least a familiar
    // one and % isn't likely to occur in station names).
    const char * q;
    for (q = s; *q; ++q) {
        unsigned char ch = *q;
        if (ch <= ' ' || ch == '%') {
            escaped.append(s, q - s);
            escaped += '%';
            escaped += "0123456789abcdef"[ch >> 4];
            escaped += "0123456789abcdef"[ch & 0x0f];
            s = q + 1;
        }
    }
    if (!escaped.empty()) {
        escaped.append(s, q - s);
        return escaped.c_str();
    }
    return s;
}

void
PLT::label(const img_point *p, const char *s, bool fSurface)
{
   fSurface = fSurface; /* unused */
   set_name(p, s);
}

void
PLT::footer(void)
{
   /* Survex is E, N, Alt - PLT file is N, E, Alt */
   fprintf(fh, "X %.3f %.3f %.3f %.3f %.3f %.3f\r\n",
           min_y / METRES_PER_FOOT, max_y / METRES_PER_FOOT,
           min_x / METRES_PER_FOOT, max_x / METRES_PER_FOOT,
           min_z / METRES_PER_FOOT, max_z / METRES_PER_FOOT);
   /* Yucky DOS "end of textfile" marker */
   PUTC('\x1a', fh);
}

class EPS : public ExportFilter {
  public:
    EPS() { }
    void header(const char *);
    void line(const img_point *, const img_point *, bool, bool);
    void label(const img_point *, const char *, bool);
    void cross(const img_point *, bool);
    void footer();
};

void
EPS::header(const char *title)
{
   const char * fontname_labels = "helvetica"; // FIXME
   int fontsize_labels = 10; // FIXME
   fputs("%!PS-Adobe-2.0 EPSF-1.2\n", fh);
   fputs("%%Creator: Survex "VERSION" EPS Output Filter\n", fh);

   if (title && title[0]) fprintf(fh, "%%%%Title: %s\n", title);

   char buf[64];
   time_t now = time(NULL);
   if (strftime(buf, sizeof(buf), "%Y-%m-%d %H:%M:%S %Z\n", localtime(&now))) {
      fputs("%%CreationDate: ", fh);
      fputs(buf, fh);
   }

   string name;
#if defined(HAVE_GETPWUID) && !defined(__DJGPP__)
   struct passwd * ent = getpwuid(getuid());
   if (ent && ent->pw_gecos[0]) name = ent->pw_gecos;
#endif
   if (name.empty()) {
       name = ::wxGetUserName().mb_str();
       if (name.empty()) {
           name = ::wxGetUserId().mb_str();
       }
   }
   if (!name.empty()) {
       fprintf(fh, "%%%%For: %s\n", name.c_str());
   }

   fprintf(fh, "%%%%BoundingBox: %d %d %d %d\n",
           int(floor(min_x * factor)), int(floor(min_y * factor)),
           int(ceil(max_x * factor)), int(ceil(max_y * factor)));
   fprintf(fh, "%%%%HiResBoundingBox: %.4f %.4f %.4f %.4f\n",
           min_x * factor, min_y * factor, max_x * factor, max_y * factor);
   fputs("%%LanguageLevel: 1\n"
         "%%PageOrder: Ascend\n"
         "%%Pages: 1\n"
         "%%Orientation: Portrait\n", fh);

   fprintf(fh, "%%%%DocumentFonts: %s\n", fontname_labels);

   fputs("%%EndComments\n"
         "%%Page 1 1\n"
         "save countdictstack mark\n", fh);

   /* this code adapted from a2ps */
   fputs("%%BeginResource: encoding ISO88591Encoding\n"
         "/ISO88591Encoding [\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs(
"/space /exclam /quotedbl /numbersign\n"
"/dollar /percent /ampersand /quoteright\n"
"/parenleft /parenright /asterisk /plus\n"
"/comma /minus /period /slash\n"
"/zero /one /two /three\n"
"/four /five /six /seven\n"
"/eight /nine /colon /semicolon\n"
"/less /equal /greater /question\n"
"/at /A /B /C /D /E /F /G\n"
"/H /I /J /K /L /M /N /O\n"
"/P /Q /R /S /T /U /V /W\n"
"/X /Y /Z /bracketleft\n"
"/backslash /bracketright /asciicircum /underscore\n"
"/quoteleft /a /b /c /d /e /f /g\n"
"/h /i /j /k /l /m /n /o\n"
"/p /q /r /s /t /u /v /w\n"
"/x /y /z /braceleft\n"
"/bar /braceright /asciitilde /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs("/.notdef /.notdef /.notdef /.notdef\n", fh);
   fputs(
"/space /exclamdown /cent /sterling\n"
"/currency /yen /brokenbar /section\n"
"/dieresis /copyright /ordfeminine /guillemotleft\n"
"/logicalnot /hyphen /registered /macron\n"
"/degree /plusminus /twosuperior /threesuperior\n"
"/acute /mu /paragraph /bullet\n"
"/cedilla /onesuperior /ordmasculine /guillemotright\n"
"/onequarter /onehalf /threequarters /questiondown\n"
"/Agrave /Aacute /Acircumflex /Atilde\n"
"/Adieresis /Aring /AE /Ccedilla\n"
"/Egrave /Eacute /Ecircumflex /Edieresis\n"
"/Igrave /Iacute /Icircumflex /Idieresis\n"
"/Eth /Ntilde /Ograve /Oacute\n"
"/Ocircumflex /Otilde /Odieresis /multiply\n"
"/Oslash /Ugrave /Uacute /Ucircumflex\n"
"/Udieresis /Yacute /Thorn /germandbls\n"
"/agrave /aacute /acircumflex /atilde\n"
"/adieresis /aring /ae /ccedilla\n"
"/egrave /eacute /ecircumflex /edieresis\n"
"/igrave /iacute /icircumflex /idieresis\n"
"/eth /ntilde /ograve /oacute\n"
"/ocircumflex /otilde /odieresis /divide\n"
"/oslash /ugrave /uacute /ucircumflex\n"
"/udieresis /yacute /thorn /ydieresis\n"
"] def\n"
"%%EndResource\n", fh);

   /* this code adapted from a2ps */
   fputs(
"/reencode {\n" /* def */
"dup length 5 add dict begin\n"
"{\n" /* forall */
"1 index /FID ne\n"
"{ def }{ pop pop } ifelse\n"
"} forall\n"
"/Encoding exch def\n"

/* Use the font's bounding box to determine the ascent, descent,
 * and overall height; don't forget that these values have to be
 * transformed using the font's matrix.
 * We use `load' because sometimes BBox is executable, sometimes not.
 * Since we need 4 numbers and not an array avoid BBox from being executed
 */
"/FontBBox load aload pop\n"
"FontMatrix transform /Ascent exch def pop\n"
"FontMatrix transform /Descent exch def pop\n"
"/FontHeight Ascent Descent sub def\n"

/* Define these in case they're not in the FontInfo (also, here
 * they're easier to get to.
 */
"/UnderlinePosition 1 def\n"
"/UnderlineThickness 1 def\n"

/* Get the underline position and thickness if they're defined. */
"currentdict /FontInfo known {\n"
"FontInfo\n"

"dup /UnderlinePosition known {\n"
"dup /UnderlinePosition get\n"
"0 exch FontMatrix transform exch pop\n"
"/UnderlinePosition exch def\n"
"} if\n"

"dup /UnderlineThickness known {\n"
"/UnderlineThickness get\n"
"0 exch FontMatrix transform exch pop\n"
"/UnderlineThickness exch def\n"
"} if\n"

"} if\n"
"currentdict\n"
"end\n"
"} bind def\n", fh);

   fprintf(fh, "/lab ISO88591Encoding /%s findfont reencode definefont pop\n",
           fontname_labels);

   fprintf(fh, "/lab findfont %d scalefont setfont\n", int(fontsize_labels));

   fprintf(fh, "0.1 setlinewidth\n");

#if 0
   /* C<digit> changes colour */
   /* FIXME: read from ini */
   {
      size_t i;
      for (i = 0; i < sizeof(colour) / sizeof(colour[0]); ++i) {
         fprintf(fh, "/C%u {stroke %.3f %.3f %.3f setrgbcolor} def\n", i,
                 (double)(colour[i] & 0xff0000) / 0xff0000,
                 (double)(colour[i] & 0xff00) / 0xff00,
                 (double)(colour[i] & 0xff) / 0xff);
      }
   }
   fputs("C0\n", fh);
#endif

   /* Postscript definition for drawing a cross */
   fprintf(fh, "/X {stroke moveto %.2f %.2f rmoveto %.2f %.2f rlineto "
           "%.2f 0 rmoveto %.2f %.2f rlineto %.2f %.2f rmoveto} def\n",
           -marker_size, -marker_size,  marker_size * 2, marker_size * 2,
           -marker_size * 2,  marker_size * 2, -marker_size * 2,
           -marker_size, marker_size );

   /* define some functions to keep file short */
   fputs("/M {stroke moveto} def\n"
         "/L {lineto} def\n"
/*       "/R {rlineto} def\n" */
         "/S {show} def\n", fh);

   fprintf(fh, "gsave %.8f dup scale\n", factor);
#if 0
   if (grid > 0) {
      double x, y;
      x = floor(min_x / grid) * grid + grid;
      y = floor(min_y / grid) * grid + grid;
#ifdef DEBUG_CAD3D
      printf("x_min: %f  y_min: %f\n", x, y);
#endif
      while (x < max_x) {
         /* horizontal line */
         fprintf(fh, "0\nLINE\n");
         fprintf(fh, "8\nGrid\n"); /* Layer */
         fprintf(fh, "10\n%6.2f\n", x);
         fprintf(fh, "20\n%6.2f\n", min_y);
         fprintf(fh, "30\n0\n");
         fprintf(fh, "11\n%6.2f\n", x);
         fprintf(fh, "21\n%6.2f\n", max_y);
         fprintf(fh, "31\n0\n");
         x += grid;
      }
      while (y < max_y) {
         /* vertical line */
         fprintf(fh, "0\nLINE\n");
         fprintf(fh, "8\nGrid\n"); /* Layer */
         fprintf(fh, "10\n%6.2f\n", min_x);
         fprintf(fh, "20\n%6.2f\n", y);
         fprintf(fh, "30\n0\n");
         fprintf(fh, "11\n%6.2f\n", max_x);
         fprintf(fh, "21\n%6.2f\n", y);
         fprintf(fh, "31\n0\n");
         y += grid;
      }
   }
#endif
}

void
EPS::line(const img_point *p1, const img_point *p, bool fSurface, bool fPendingMove)
{
   fSurface = fSurface; /* unused */
   if (fPendingMove) {
       fprintf(fh, "%.2f %.2f M\n", p1->x, p1->y);
   }
   fprintf(fh, "%.2f %.2f L\n", p->x, p->y);
}

void
EPS::label(const img_point *p, const char *s, bool fSurface)
{
   fprintf(fh, "%.2f %.2f M\n", p->x, p->y);
   PUTC('(', fh);
   while (*s) {
       unsigned char ch = *s++;
       switch (ch) {
           case '(': case ')': case '\\': /* need to escape these characters */
               PUTC('\\', fh);
               PUTC(ch, fh);
               break;
           default:
               PUTC(ch, fh);
               break;
       }
   }
   fputs(") S\n", fh);
}

void
EPS::cross(const img_point *p, bool fSurface)
{
   fSurface = fSurface; /* unused */
   fprintf(fh, "%.2f %.2f X\n", p->x, p->y);
}

void
EPS::footer(void)
{
   fputs("stroke showpage grestore\n"
         "%%Trailer\n"
         "cleartomark countdictstack exch sub { end } repeat restore\n"
         "%%EOF\n", fh);
}

typedef enum {
    FMT_DXF = 0, FMT_SVG, FMT_SKETCH, FMT_PLT, FMT_EPS, FMT_HPGL, FMT_ENDMARKER
} export_format;
static const char *extensions[] = { "dxf", "svg", "sk", "plt", "eps", "hpgl" };

bool
Export(const wxString &fnm_out, const wxString &title, const MainFrm * mainfrm,
       double pan, double tilt, bool labels, bool crosses, bool legs,
       bool surface)
{
   int fPendingMove = 0;
   img_point p, p1;
   double s = 0, c = 0;
   export_format format = FMT_DXF;
   const int *pass;
   bool elevation = (tilt == 0.0);

   /* Defaults */
   grid = 0;
   text_height = TEXT_HEIGHT;
   marker_size = MARKER_SIZE;
   // FIXME: allow these to be set from aven somehow!
#if 0
       case 'g': /* Grid */
         if (optarg) {
            grid = cmdline_double_arg();
         } else {
            grid = (double)GRID_SPACING;
         }
         break;
       case 'r': /* Reduction factor */
         scale = cmdline_double_arg();
         break;
       case 't': /* Text height */
         text_height = cmdline_double_arg();
         break;
       case 'm': /* Marker size */
         marker_size = cmdline_double_arg();
         break;
       case 's':
         survey = optarg;
         break;
#endif

   {
      size_t i;
      size_t len = fnm_out.length();
      for (i = 0; i < FMT_ENDMARKER; ++i) {
         size_t l = strlen(extensions[i]);
         if (len > l + 1 && fnm_out[len - l - 1] == FNM_SEP_EXT &&
             strcasecmp((const char *)fnm_out.mb_str() + len - l, extensions[i]) == 0) {
            format = export_format(i);
            break;
         }
      }
   }

   ExportFilter * filt;
   switch (format) {
       case FMT_DXF:
           filt = new DXF;
           break;
       case FMT_SKETCH:
           filt = new Sketch;
           factor = POINTS_PER_MM * 1000.0 / scale;
           break;
       case FMT_PLT:
           filt = new PLT;
           break;
       case FMT_SVG:
           filt = new SVG;
           factor = 1000.0 / scale;
           break;
       case FMT_EPS:
           filt = new EPS;
           factor = POINTS_PER_MM * 1000.0 / scale;
           break;
       case FMT_HPGL:
           filt = new HPGL;
           // factor = POINTS_PER_MM * 1000.0 / scale;
           break;
       default:
           return false;
   }

   // FIXME: This should really use fn_str() - currently we probably can't
   // save to a Unicode path on wxmsw.
   if (!filt->fopen(fnm_out.mb_str())) {
       delete filt;
       return false;
   }

   if (elevation) {
      s = sin(rad(pan));
      c = cos(rad(pan));
   }

   /* Get drawing corners */
   min_x = min_y = min_z = HUGE_VAL;
   max_x = max_y = max_z = -HUGE_VAL;
   list<traverse>::const_iterator trav = mainfrm->traverses_begin();
   list<traverse>::const_iterator tend = mainfrm->traverses_end();
   for ( ; trav != tend; ++trav) {
        vector<PointInfo>::const_iterator pos = trav->begin();
        vector<PointInfo>::const_iterator end = trav->end();
        for ( ; pos != end; ++pos) {
            p.x = pos->GetX();
            p.y = pos->GetY();
            p.z = pos->GetZ();

            if (elevation) {
                double xnew = p.x * c - p.y * s;
                double znew = - p.x * s - p.y * c;
                p.y = p.z;
                p.z = znew;
                p.x = xnew;
            }

            if (p.x < min_x) min_x = p.x;
            if (p.x > max_x) max_x = p.x;
            if (p.y < min_y) min_y = p.y;
            if (p.y > max_y) max_y = p.y;
            if (p.z < min_z) min_z = p.z;
            if (p.z > max_z) max_z = p.z;
        }
   }
   {
        list<LabelInfo*>::const_iterator pos = mainfrm->GetLabels();
        list<LabelInfo*>::const_iterator end = mainfrm->GetLabelsEnd();
        for ( ; pos != end; ++pos) {
            p.x = (*pos)->GetX();
            p.y = (*pos)->GetY();
            p.z = (*pos)->GetZ();

            if (elevation) {
                double xnew = p.x * c - p.y * s;
                double znew = - p.x * s - p.y * c;
                p.y = p.z;
                p.z = znew;
                p.x = xnew;
            }

            if (p.x < min_x) min_x = p.x;
            if (p.x > max_x) max_x = p.x;
            if (p.y < min_y) min_y = p.y;
            if (p.y > max_y) max_y = p.y;
            if (p.z < min_z) min_z = p.z;
            if (p.z > max_z) max_z = p.z;
        }
   }

   if (grid > 0) {
      min_x -= grid / 2;
      max_x += grid / 2;
      min_y -= grid / 2;
      max_y += grid / 2;
   }

   /* handle empty file gracefully */
   if (min_x > max_x) {
      min_x = min_y = min_z = 0;
      max_x = max_y = max_z = 0;
   }

   /* Header */
   filt->header(title.mb_str());

   p1.x = p1.y = p1.z = 0; /* avoid compiler warning */

   for (pass = filt->passes(); *pass; ++pass) {
      bool legs_this_pass = ((*pass & LEGS) && legs);
      bool surf_this_pass = ((*pass & SURF) && legs && surface);
      bool crosses_this_pass = ((*pass & STNS) && crosses);
      bool labels_this_pass = ((*pass & LABELS) && labels);
      if (!(legs_this_pass || surf_this_pass || crosses_this_pass ||
            labels_this_pass))
          continue;
      filt->start_pass(*pass);
      if (legs_this_pass) {
          trav = mainfrm->traverses_begin();
          tend = mainfrm->traverses_end();
          for ( ; trav != tend; ++trav) {
             assert(trav->size() > 1);
             vector<PointInfo>::const_iterator pos = trav->begin();
             vector<PointInfo>::const_iterator end = trav->end();
             for ( ; pos != end; ++pos) {
                 p.x = pos->GetX();
                 p.y = pos->GetY();
                 p.z = pos->GetZ();

                 if (format == FMT_SKETCH) {
                     p.x -= min_x;
                     p.y -= min_y;
                     p.z -= min_z;
                 } else if (format == FMT_HPGL) {
                     p.x -= (min_x + max_x) / 2;
                     p.y -= (min_y + max_y) / 2;
                     p.z -= (min_z + max_z) / 2;
                 }

                 if (elevation) {
                     double xnew = p.x * c - p.y * s;
                     double znew = - p.x * s - p.y * c;
                     p.y = p.z;
                     p.z = znew;
                     p.x = xnew;
                 }

                 if (pos == trav->begin()) {
                     // First point is move...
#ifdef DEBUG_CAD3D
                     printf("move to %9.2f %9.2f %9.2f\n",x,y,z);
#endif
                     fPendingMove = 1;
                 } else {
#ifdef DEBUG_CAD3D
                     printf("line to %9.2f %9.2f %9.2f\n", p.x, p.y, p.z);
#endif
                     filt->line(&p1, &p, false, fPendingMove);
                     fPendingMove = 0;
                 }
                 p1 = p;
             }
         }
      }
      if (surf_this_pass) {
          trav = mainfrm->surface_traverses_begin();
          tend = mainfrm->surface_traverses_end();
          for ( ; trav != tend; ++trav) {
             assert(trav->size() > 1);
             vector<PointInfo>::const_iterator pos = trav->begin();
             vector<PointInfo>::const_iterator end = trav->end();
             for ( ; pos != end; ++pos) {
                 p.x = pos->GetX();
                 p.y = pos->GetY();
                 p.z = pos->GetZ();

                 if (format == FMT_SKETCH) {
                     p.x -= min_x;
                     p.y -= min_y;
                     p.z -= min_z;
                 } else if (format == FMT_HPGL) {
                     p.x -= (min_x + max_x) / 2;
                     p.y -= (min_y + max_y) / 2;
                     p.z -= (min_z + max_z) / 2;
                 }

                 if (elevation) {
                     double xnew = p.x * c - p.y * s;
                     double znew = - p.x * s - p.y * c;
                     p.y = p.z;
                     p.z = znew;
                     p.x = xnew;
                 }

                 if (pos == trav->begin()) {
                     // First point is move...
#ifdef DEBUG_CAD3D
                     printf("surface move to %9.2f %9.2f %9.2f\n",x,y,z);
#endif
                     fPendingMove = 1;
                 } else {
#ifdef DEBUG_CAD3D
                     printf("surface line to %9.2f %9.2f %9.2f\n", p.x, p.y, p.z);
#endif
                     filt->line(&p1, &p, true, fPendingMove);
                     fPendingMove = 0;
                 }
                 p1 = p;
             }
         }
      }
      if (crosses_this_pass || labels_this_pass) {
          list<LabelInfo*>::const_iterator pos = mainfrm->GetLabels();
          list<LabelInfo*>::const_iterator end = mainfrm->GetLabelsEnd();
          for ( ; pos != end; ++pos) {
              p.x = (*pos)->GetX();
              p.y = (*pos)->GetY();
              p.z = (*pos)->GetZ();

              if (format == FMT_SKETCH) {
                  p.x -= min_x;
                  p.y -= min_y;
                  p.z -= min_z;
              } else if (format == FMT_HPGL) {
                  p.x -= (min_x + max_x) / 2;
                  p.y -= (min_y + max_y) / 2;
                  p.z -= (min_z + max_z) / 2;
              }

              if (elevation) {
                  double xnew = p.x * c - p.y * s;
                  double znew = - p.x * s - p.y * c;
                  p.y = p.z;
                  p.z = znew;
                  p.x = xnew;
              }
#ifdef DEBUG_CAD3D
              printf("label '%s' at %9.2f %9.2f %9.2f\n",(*pos)->GetText(),x,y,z);
#endif
              /* Use !UNDERGROUND as the criterion - we want stations where
               * a surface and underground survey meet to be in the
               * underground layer */
              if (labels_this_pass)
                  filt->label(&p, (*pos)->GetText().mb_str(), !(*pos)->IsUnderground());
              if (crosses_this_pass)
                  filt->cross(&p, !(*pos)->IsUnderground());
          }
      }
   }
   filt->footer();
   delete filt;
   osfree(htab);
   htab = NULL;
   return true;
}