]> git.mxchange.org Git - simgear.git/blobdiff - simgear/io/sg_binobj.cxx
Introduce SGBinaryFile
[simgear.git] / simgear / io / sg_binobj.cxx
index c6cfc63686eaa23a2a9181bb210a3a82f4cacd7b..1c604a28b39fb7e2555324050b3197c026071cae 100644 (file)
@@ -2,7 +2,7 @@
 //
 // Written by Curtis Olson, started January 2000.
 //
-// Copyright (C) 2000  Curtis L. Olson  - curt@flightgear.org
+// Copyright (C) 2000  Curtis L. Olson  - http://www.flightgear.org/~curt
 //
 // 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
 //
 // 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., 675 Mass Ave, Cambridge, MA 02139, USA.
+// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
 //
 // $Id$
 //
 
 
+#ifdef HAVE_CONFIG_H
+#  include <simgear_config.h>
+#endif
+
 #include <simgear/compiler.h>
+#include <simgear/debug/logstream.hxx>
 
 #include <stdio.h>
 #include <time.h>
+#include <cstring>
+#include <cstdlib> // for system()
+#include <cassert>
 
 #include <vector>
-#include STL_STRING
+#include <string>
+#include <iostream>
+#include <bitset>
 
 #include <simgear/bucket/newbucket.hxx>
+#include <simgear/misc/sg_path.hxx>
+#include <simgear/math/SGGeometry.hxx>
+#include <simgear/structure/exception.hxx>
 
 #include "lowlevel.hxx"
 #include "sg_binobj.hxx"
 
 
-SG_USING_STD( string );
-SG_USING_STD( vector );
-
-#if !defined (SG_HAVE_NATIVE_SGI_COMPILERS)
-SG_USING_STD( cout );
-SG_USING_STD( endl );
-#endif
-
+using std::string;
+using std::vector;
+using std::cout;
+using std::endl;
 
-enum {
+enum sgObjectTypes {
     SG_BOUNDING_SPHERE = 0,
 
     SG_VERTEX_LIST = 1,
     SG_NORMAL_LIST = 2,
     SG_TEXCOORD_LIST = 3,
+    SG_COLOR_LIST = 4,
+    SG_VA_FLOAT_LIST = 5,
+    SG_VA_INTEGER_LIST = 6,
+
+    SG_POINTS = 9,
 
     SG_TRIANGLE_FACES = 10,
     SG_TRIANGLE_STRIPS = 11,
     SG_TRIANGLE_FANS = 12
-} tgObjectTypes;
+};
+
+enum sgIndexTypes {
+    SG_IDX_VERTICES =    0x01,
+    SG_IDX_NORMALS =     0x02,
+    SG_IDX_COLORS =      0x04,
+    SG_IDX_TEXCOORDS_0 = 0x08,
+    SG_IDX_TEXCOORDS_1 = 0x10,
+    SG_IDX_TEXCOORDS_2 = 0x20,
+    SG_IDX_TEXCOORDS_3 = 0x40,
+};
 
-enum {
-    SG_MATERIAL = 0
-} tgPropertyTypes;
+enum sgVertexAttributeTypes {
+    // vertex attributes
+    SG_VA_INTEGER_0 = 0x00000001,
+    SG_VA_INTEGER_1 = 0x00000002,
+    SG_VA_INTEGER_2 = 0x00000004,
+    SG_VA_INTEGER_3 = 0x00000008,
+
+    SG_VA_FLOAT_0 =   0x00000100,
+    SG_VA_FLOAT_1 =   0x00000200,
+    SG_VA_FLOAT_2 =   0x00000400,
+    SG_VA_FLOAT_3 =   0x00000800,
+};
+
+enum sgPropertyTypes {
+    SG_MATERIAL = 0,
+    SG_INDEX_TYPES = 1,
+    SG_VERT_ATTRIBS = 2
+};
 
 
 class sgSimpleBuffer {
@@ -68,120 +107,442 @@ private:
 
     char *ptr;
     unsigned int size;
-
+    size_t offset;
 public:
 
-    inline sgSimpleBuffer( unsigned int s )
+    sgSimpleBuffer( unsigned int s = 0) :
+        ptr(NULL),
+        size(0),
+        offset(0)
     {
-       size = 1;
-       while ( size < s ) {
-           size *= 2;
-       }
-       cout << "Creating a new buffer of size = " << size << endl;
-       ptr = new char[size];
-    }
-
-    inline ~sgSimpleBuffer() {
-       delete [] ptr;
-    }
-
-    inline unsigned int get_size() const { return size; }
-    inline char *get_ptr() const { return ptr; }
-    inline void resize( unsigned int s ) {
-       if ( s > size ) {
-           if ( ptr != NULL ) {
-               delete [] ptr;
-           }
-           while ( size < s ) {
-               size *= 2;
-           }
-           cout << "resizing buffer to size = " << size << endl;
-           ptr = new char[size];
-       }
+        resize(s);
     }
-};
 
+    ~sgSimpleBuffer()
+    {
+        delete [] ptr;
+    }
 
-// calculate the center of a list of points, by taking the halfway
-// point between the min and max points.
-static Point3D calc_center( point_list& wgs84_nodes ) {
-    Point3D p, min, max;
+    unsigned int get_size() const { return size; }
+    char *get_ptr() const { return ptr; }
+    
+    void reset()
+    {
+        offset = 0;
+    }
+    
+    void resize( unsigned int s )
+    {
+        if ( s > size ) {
+            if ( ptr != NULL ) {
+                delete [] ptr;
+            }
+          
+            if ( size == 0) {
+                size = 16;
+            }
+          
+            while ( size < s ) {
+              size = size << 1;
+            }
+            ptr = new char[size];
+        }
+    }
+    
+    float readInt()
+    {
+        unsigned int* p = reinterpret_cast<unsigned int*>(ptr + offset);
+        if ( sgIsBigEndian() ) {            
+            sgEndianSwap((uint32_t *) p);
+        }
+        
+        offset += sizeof(unsigned int);
+        return *p;
+    }
+    
+    SGVec3d readVec3d()
+    {
+        double* p = reinterpret_cast<double*>(ptr + offset);
+        
+        if ( sgIsBigEndian() ) {            
+            sgEndianSwap((uint64_t *) p + 0);
+            sgEndianSwap((uint64_t *) p + 1);
+            sgEndianSwap((uint64_t *) p + 2);
+        }
+        
+        offset += 3 * sizeof(double);
+        return SGVec3d(p);
+    }
+    
+    float readFloat()
+    {
+        float* p = reinterpret_cast<float*>(ptr + offset);
+        if ( sgIsBigEndian() ) {            
+            sgEndianSwap((uint32_t *) p);
+        }
+        
+        offset += sizeof(float);
+        return *p;
+    }
+    
+    SGVec2f readVec2f()
+    {
+        float* p = reinterpret_cast<float*>(ptr + offset);
+        
+        if ( sgIsBigEndian() ) {            
+            sgEndianSwap((uint32_t *) p + 0);
+            sgEndianSwap((uint32_t *) p + 1);
+        }
+        
+        offset += 2 * sizeof(float);
+        return SGVec2f(p);
+    }
+    
+    SGVec3f readVec3f()
+    {
+        float* p = reinterpret_cast<float*>(ptr + offset);
+        
+        if ( sgIsBigEndian() ) {            
+            sgEndianSwap((uint32_t *) p + 0);
+            sgEndianSwap((uint32_t *) p + 1);
+            sgEndianSwap((uint32_t *) p + 2);
+        }
+        
+        offset += 3 * sizeof(float);
+        return SGVec3f(p);
+    }
+    
+    SGVec4f readVec4f()
+    {
+        float* p = reinterpret_cast<float*>(ptr + offset);
+        
+        if ( sgIsBigEndian() ) {            
+            sgEndianSwap((uint32_t *) p + 0);
+            sgEndianSwap((uint32_t *) p + 1);
+            sgEndianSwap((uint32_t *) p + 2);
+            sgEndianSwap((uint32_t *) p + 3);
+        }
+        
+        offset += 4 * sizeof(float);
+        return SGVec4f(p);
+    }
+};
 
-    if ( wgs84_nodes.size() ) {
-       min = max = wgs84_nodes[0];
-    } else {
-       min = max = Point3D( 0 );
+template <class T>
+static void read_indices(char* buffer, 
+                         size_t bytes,
+                         int indexMask,
+                         int vaMask,
+                         int_list& vertices, 
+                         int_list& normals,
+                         int_list& colors,
+                         tci_list& texCoords,
+                         vai_list& vas
+                        )
+{
+    const int indexSize = sizeof(T) * std::bitset<32>((int)indexMask).count();
+    const int vaSize = sizeof(T) * std::bitset<32>((int)vaMask).count();
+    const int count = bytes / (indexSize + vaSize);
+    
+    // fix endian-ness of the whole lot, if required
+    if (sgIsBigEndian()) {
+        int indices = bytes / sizeof(T);
+        T* src = reinterpret_cast<T*>(buffer);
+        for (int i=0; i<indices; ++i) {
+            sgEndianSwap(src++);
+        }
+    }
+    
+    T* src = reinterpret_cast<T*>(buffer);
+    for (int i=0; i<count; ++i) {
+        if (indexMask & SG_IDX_VERTICES) vertices.push_back(*src++);
+        if (indexMask & SG_IDX_NORMALS) normals.push_back(*src++);
+        if (indexMask & SG_IDX_COLORS) colors.push_back(*src++);
+        if (indexMask & SG_IDX_TEXCOORDS_0) texCoords[0].push_back(*src++);
+        if (indexMask & SG_IDX_TEXCOORDS_1) texCoords[1].push_back(*src++);
+        if (indexMask & SG_IDX_TEXCOORDS_2) texCoords[2].push_back(*src++);
+        if (indexMask & SG_IDX_TEXCOORDS_3) texCoords[3].push_back(*src++);
+        
+        if ( vaMask ) {
+            if (vaMask & SG_VA_INTEGER_0) vas[0].push_back(*src++);
+            if (vaMask & SG_VA_INTEGER_1) vas[1].push_back(*src++);
+            if (vaMask & SG_VA_INTEGER_2) vas[2].push_back(*src++);
+            if (vaMask & SG_VA_INTEGER_3) vas[3].push_back(*src++);
+            if (vaMask & SG_VA_FLOAT_0) vas[4].push_back(*src++);
+            if (vaMask & SG_VA_FLOAT_1) vas[5].push_back(*src++);
+            if (vaMask & SG_VA_FLOAT_2) vas[6].push_back(*src++);
+            if (vaMask & SG_VA_FLOAT_3) vas[7].push_back(*src++);
+        }
+    } // of elements in the index
+    
+    // WS2.0 fix : toss zero area triangles
+    if ( ( count == 3 ) && (indexMask & SG_IDX_VERTICES) ) {
+        if ( (vertices[0] == vertices[1]) || 
+             (vertices[1] == vertices[2]) || 
+             (vertices[2] == vertices[0]) ) {
+            vertices.clear();
+        }
     }
+}
+
+template <class T>
+void write_indice(gzFile fp, T value)
+{
+    sgWriteBytes(fp, sizeof(T), &value);
+}
 
-    for ( int i = 0; i < (int)wgs84_nodes.size(); ++i ) {
-       p = wgs84_nodes[i];
+// specialize template to call endian-aware conversion methods
+template <>
+void write_indice(gzFile fp, uint16_t value)
+{
+    sgWriteUShort(fp, value);
+}
 
-       if ( p.x() < min.x() ) { min.setx( p.x() ); }
-       if ( p.y() < min.y() ) { min.sety( p.y() ); }
-       if ( p.z() < min.z() ) { min.setz( p.z() ); }
+template <>
+void write_indice(gzFile fp, uint32_t value)
+{
+    sgWriteUInt(fp, value);
+}
 
-       if ( p.x() > max.x() ) { max.setx( p.x() ); }
-       if ( p.y() > max.y() ) { max.sety( p.y() ); }
-       if ( p.z() > max.z() ) { max.setz( p.z() ); }
-    }
 
-    return ( min + max ) / 2.0;
+template <class T>
+void write_indices(gzFile fp, 
+    unsigned char indexMask, 
+    unsigned int vaMask,
+    const int_list& vertices, 
+    const int_list& normals, 
+    const int_list& colors,
+    const tci_list& texCoords,
+    const vai_list& vas )
+{
+    unsigned int count = vertices.size();
+    const int indexSize = sizeof(T) * std::bitset<32>((int)indexMask).count();
+    const int vaSize = sizeof(T) * std::bitset<32>((int)vaMask).count();
+    sgWriteUInt(fp, (indexSize + vaSize) * count);
+            
+    for (unsigned int i=0; i < count; ++i) {
+        write_indice(fp, static_cast<T>(vertices[i]));
+        
+        if (indexMask & SG_IDX_NORMALS) {
+            write_indice(fp, static_cast<T>(normals[i]));
+        }
+        if (indexMask & SG_IDX_COLORS) {
+            write_indice(fp, static_cast<T>(colors[i]));
+        }
+        if (indexMask & SG_IDX_TEXCOORDS_0) {
+            write_indice(fp, static_cast<T>(texCoords[0][i]));
+        }
+        if (indexMask & SG_IDX_TEXCOORDS_1) {
+            write_indice(fp, static_cast<T>(texCoords[1][i]));
+        }
+        if (indexMask & SG_IDX_TEXCOORDS_2) {
+            write_indice(fp, static_cast<T>(texCoords[2][i]));
+        }
+        if (indexMask & SG_IDX_TEXCOORDS_3) {
+            write_indice(fp, static_cast<T>(texCoords[3][i]));
+        }
+        
+        if (vaMask) {
+            if (vaMask & SG_VA_INTEGER_0) {
+                write_indice(fp, static_cast<T>(vas[0][i]));
+            }
+            if (vaMask & SG_VA_INTEGER_1) {
+                write_indice(fp, static_cast<T>(vas[1][i]));
+            }
+            if (vaMask & SG_VA_INTEGER_2) {
+                write_indice(fp, static_cast<T>(vas[2][i]));
+            }
+            if (vaMask & SG_VA_INTEGER_3) {
+                write_indice(fp, static_cast<T>(vas[3][i]));
+            }
+
+            if (vaMask & SG_VA_FLOAT_0) {
+                write_indice(fp, static_cast<T>(vas[4][i]));
+            }
+            if (vaMask & SG_VA_FLOAT_1) {
+                write_indice(fp, static_cast<T>(vas[5][i]));
+            }
+            if (vaMask & SG_VA_FLOAT_2) {
+                write_indice(fp, static_cast<T>(vas[6][i]));
+            }
+            if (vaMask & SG_VA_FLOAT_3) {
+                write_indice(fp, static_cast<T>(vas[7][i]));
+            }
+        }
+    }
 }
 
-// calculate the bounding sphere.  Center is the center of the
-// tile and zero elevation
-double sgCalcBoundingRadius( Point3D center, point_list& wgs84_nodes ) {
-    double dist_squared;
-    double radius_squared = 0;
+
+// read object properties
+void SGBinObject::read_object( gzFile fp,
+                         int obj_type,
+                         int nproperties,
+                         int nelements,
+                         group_list& vertices, 
+                         group_list& normals,
+                         group_list& colors,
+                         group_tci_list& texCoords,
+                         group_vai_list& vertexAttribs,
+                         string_list& materials)
+{
+    unsigned int  nbytes;
+    unsigned char idx_mask;
+    unsigned int  vertex_attrib_mask;
+    int j;
+    sgSimpleBuffer buf( 32768 );  // 32 Kb
+    char material[256];
+
+    // default values
+    if ( obj_type == SG_POINTS ) {
+        idx_mask = SG_IDX_VERTICES;
+    } else {
+        idx_mask = (char)(SG_IDX_VERTICES | SG_IDX_TEXCOORDS_0);
+    }
+    vertex_attrib_mask = 0;
     
-    for ( int i = 0; i < (int)wgs84_nodes.size(); ++i ) {
-        dist_squared = center.distance3Dsquared( wgs84_nodes[i] );
-       if ( dist_squared > radius_squared ) {
-            radius_squared = dist_squared;
+    for ( j = 0; j < nproperties; ++j ) {
+        char prop_type;
+        sgReadChar( fp, &prop_type );
+        sgReadUInt( fp, &nbytes );
+        
+        buf.resize(nbytes);
+        char *ptr = buf.get_ptr();
+        
+        switch( prop_type )
+        {
+            case SG_MATERIAL:
+                sgReadBytes( fp, nbytes, ptr );
+                if (nbytes > 255) {
+                    nbytes = 255;
+                }
+                strncpy( material, ptr, nbytes );
+                material[nbytes] = '\0';
+                break;
+                
+            case SG_INDEX_TYPES:
+                if (nbytes == 1) {
+                    sgReadChar( fp, (char *)&idx_mask );                    
+                } else {
+                    sgReadBytes( fp, nbytes, ptr );                    
+                }
+                break;
+
+            case SG_VERT_ATTRIBS:
+                if (nbytes == 4) {
+                    sgReadUInt( fp, &vertex_attrib_mask );
+                } else {
+                    sgReadBytes( fp, nbytes, ptr );
+                }
+                break;
+                
+            default:
+                sgReadBytes( fp, nbytes, ptr );
+                SG_LOG(SG_IO, SG_ALERT, "Found UNKNOWN property type with nbytes == " << nbytes << " mask is " << (int)idx_mask );
+                break;
         }
     }
 
-    return sqrt(radius_squared);
+    if ( sgReadError() ) {
+        throw sg_exception("Error reading object properties");
+    }
+    
+    size_t indexCount = std::bitset<32>((int)idx_mask).count();
+    if (indexCount == 0) {
+        throw sg_exception("object index mask has no bits set");
+    }
+    
+    for ( j = 0; j < nelements; ++j ) {
+        sgReadUInt( fp, &nbytes );
+        if ( sgReadError() ) {
+            throw sg_exception("Error reading element size");
+        }
+        
+        buf.resize( nbytes );
+        char *ptr = buf.get_ptr();
+        sgReadBytes( fp, nbytes, ptr );
+        
+        if ( sgReadError() ) {
+            throw sg_exception("Error reading element bytes");
+        }
+                
+        int_list vs;
+        int_list ns;
+        int_list cs;
+        tci_list tcs;
+        vai_list vas;
+
+        if (version >= 10) {
+            read_indices<uint32_t>(ptr, nbytes, idx_mask, vertex_attrib_mask, vs, ns, cs, tcs, vas );
+        } else {
+            read_indices<uint16_t>(ptr, nbytes, idx_mask, vertex_attrib_mask, vs, ns, cs, tcs, vas );
+        }
+
+        // Fix for WS2.0 - ignore zero area triangles
+        if ( !vs.empty() ) {
+            vertices.push_back( vs );
+            normals.push_back( ns );
+            colors.push_back( cs );
+            texCoords.push_back( tcs );
+            vertexAttribs.push_back( vas );
+            materials.push_back( material );
+        }
+    } // of element iteration
 }
 
 
 // read a binary file and populate the provided structures.
 bool SGBinObject::read_bin( const string& file ) {
-    Point3D p;
-    int i, j, k;
-    char material[256];
+    SGVec3d p;
+    int i, k;
+    size_t j;
     unsigned int nbytes;
-    static sgSimpleBuffer buf( 32768 );  // 32 Kb
+    sgSimpleBuffer buf( 32768 );  // 32 Kb
 
     // zero out structures
-    gbs_center = Point3D( 0 );
+    gbs_center = SGVec3d(0, 0, 0);
     gbs_radius = 0.0;
 
     wgs84_nodes.clear();
     normals.clear();
     texcoords.clear();
 
+    pts_v.clear();
+    pts_n.clear();
+    pts_c.clear();
+    pts_tcs.clear();
+    pts_vas.clear();
+    pt_materials.clear();
+
     tris_v.clear();
-    tris_tc.clear();
+    tris_n.clear();
+    tris_c.clear();
+    tris_tcs.clear();
+    tris_vas.clear();
     tri_materials.clear();
 
     strips_v.clear();
-    strips_tc.clear();
+    strips_n.clear();
+    strips_c.clear();
+    strips_tcs.clear();
+    strips_vas.clear();
     strip_materials.clear();
 
     fans_v.clear();
-    fans_tc.clear();
+    fans_n.clear();
+    fans_c.clear();
+    fans_tcs.clear();
+    fans_vas.clear();
     fan_materials.clear();
-   
+
     gzFile fp;
     if ( (fp = gzopen( file.c_str(), "rb" )) == NULL ) {
-       string filegz = file + ".gz";
-       if ( (fp = gzopen( filegz.c_str(), "rb" )) == NULL ) {
-           // cout << "ERROR: opening " << file << " or " << filegz
-           //      << "for reading!" << endl;
+        string filegz = file + ".gz";
+        if ( (fp = gzopen( filegz.c_str(), "rb" )) == NULL ) {
+            SG_LOG( SG_EVENT, SG_ALERT,
+               "ERROR: opening " << file << " or " << filegz << " for reading!");
 
-           return false;
-       }
+            throw sg_io_exception("Error opening for reading (and .gz)", sg_location(file));
+        }
     }
 
     sgClearReadError();
@@ -190,23 +551,22 @@ bool SGBinObject::read_bin( const string& file ) {
     unsigned int header;
     sgReadUInt( fp, &header );
     if ( ((header & 0xFF000000) >> 24) == 'S' &&
-        ((header & 0x00FF0000) >> 16) == 'G' ) {
-       // cout << "Good header" << endl;
-       // read file version
-       version = (header & 0x0000FFFF);
-       // cout << "File version = " << version << endl;
+         ((header & 0x00FF0000) >> 16) == 'G' ) {
+    
+        // read file version
+        version = (header & 0x0000FFFF);
     } else {
-       // close the file before we return
-       gzclose(fp);
-
-       return false;
+        // close the file before we return
+        gzclose(fp);
+        throw sg_io_exception("Bad BTG magic/version", sg_location(file));
     }
-
+    
     // read creation time
-    time_t calendar_time;
-    sgReadLong( fp, &calendar_time );
+    unsigned int foo_calendar_time;
+    sgReadUInt( fp, &foo_calendar_time );
 
 #if 0
+    time_t calendar_time = foo_calendar_time;
     // The following code has a global effect on the host application
     // and can screws up the time elsewhere.  It should be avoided
     // unless you need this for debugging in which case you should
@@ -215,609 +575,497 @@ bool SGBinObject::read_bin( const string& file ) {
     local_tm = localtime( &calendar_time );
     char time_str[256];
     strftime( time_str, 256, "%a %b %d %H:%M:%S %Z %Y", local_tm);
-    cout << "File created on " << time_str << endl;
+    SG_LOG( SG_EVENT, SG_DEBUG, "File created on " << time_str);
 #endif
 
     // read number of top level objects
-    short nobjects;
-    sgReadShort( fp, &nobjects );
-    // cout << "Total objects to read = " << nobjects << endl;
+    int nobjects;
+    if ( version >= 10) { // version 10 extends everything to be 32-bit
+        sgReadInt( fp, &nobjects );
+    } else if ( version >= 7 ) {
+        uint16_t v;
+        sgReadUShort( fp, &v );
+        nobjects = v;
+    } else {
+        int16_t v;
+        sgReadShort( fp, &v );
+        nobjects = v;
+    }
+     
+    SG_LOG(SG_IO, SG_DEBUG, "SGBinObject::read_bin Total objects to read = " << nobjects);
 
+    if ( sgReadError() ) {
+        throw sg_io_exception("Error reading BTG file header", sg_location(file));
+    }
+    
     // read in objects
     for ( i = 0; i < nobjects; ++i ) {
-       // read object header
-       char obj_type;
-       short nproperties, nelements;
-       sgReadChar( fp, &obj_type );
-       sgReadShort( fp, &nproperties );
-       sgReadShort( fp, &nelements );
-
-       // cout << "object " << i << " = " << (int)obj_type << " props = "
-       //      << nproperties << " elements = " << nelements << endl;
-           
-       if ( obj_type == SG_BOUNDING_SPHERE ) {
-           // read bounding sphere properties
-           for ( j = 0; j < nproperties; ++j ) {
-               char prop_type;
-               sgReadChar( fp, &prop_type );
-
-               sgReadUInt( fp, &nbytes );
-               // cout << "property size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-           }
-
-           // read bounding sphere elements
-           for ( j = 0; j < nelements; ++j ) {
-               sgReadUInt( fp, &nbytes );
-               // cout << "element size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-
-               double *dptr = (double *)ptr;
-               if ( sgIsBigEndian() ) {
-                   sgEndianSwap( (uint64 *)&(dptr[0]) );
-                   sgEndianSwap( (uint64 *)&(dptr[1]) );
-                   sgEndianSwap( (uint64 *)&(dptr[2]) );
-               }
-               gbs_center = Point3D( dptr[0], dptr[1], dptr[2] );
-               // cout << "Center = " << gbs_center << endl;
-               ptr += sizeof(double) * 3;
-               
-               float *fptr = (float *)ptr;
-               if ( sgIsBigEndian() ) {
-                   sgEndianSwap( (unsigned int *)fptr );
-               }
-               gbs_radius = fptr[0];
-               // cout << "Bounding radius = " << gbs_radius << endl;
-           }
-       } else if ( obj_type == SG_VERTEX_LIST ) {
-           // read vertex list properties
-           for ( j = 0; j < nproperties; ++j ) {
-               char prop_type;
-               sgReadChar( fp, &prop_type );
-
-               sgReadUInt( fp, &nbytes );
-               // cout << "property size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-           }
-
-           // read vertex list elements
-           for ( j = 0; j < nelements; ++j ) {
-               sgReadUInt( fp, &nbytes );
-               // cout << "element size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-               int count = nbytes / (sizeof(float) * 3);
-               float *fptr = (float *)ptr;
-               for ( k = 0; k < count; ++k ) {
-                   if ( sgIsBigEndian() ) {
-                       sgEndianSwap( (unsigned int *)&(fptr[0]) );
-                       sgEndianSwap( (unsigned int *)&(fptr[1]) );
-                       sgEndianSwap( (unsigned int *)&(fptr[2]) );
-                   }
-                   p = Point3D( fptr[0], fptr[1], fptr[2] );
-                   // cout << "node = " << p << endl;
-                   wgs84_nodes.push_back( p );
-                   fptr += 3;
-               }
-           }
-       } else if ( obj_type == SG_NORMAL_LIST ) {
-           // read normal list properties
-           for ( j = 0; j < nproperties; ++j ) {
-               char prop_type;
-               sgReadChar( fp, &prop_type );
-
-               sgReadUInt( fp, &nbytes );
-               // cout << "property size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-           }
-
-           // read normal list elements
-           for ( j = 0; j < nelements; ++j ) {
-               sgReadUInt( fp, &nbytes );
-               // cout << "element size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               unsigned char *ptr = (unsigned char *)(buf.get_ptr());
-               sgReadBytes( fp, nbytes, ptr );
-               int count = nbytes / 3;
-               for ( k = 0; k < count; ++k ) {
-                    sgdVec3 normal;
-                    sgdSetVec3( normal,
-                               (ptr[0]) / 127.5 - 1.0,
-                               (ptr[1]) / 127.5 - 1.0,
-                               (ptr[2]) / 127.5 - 1.0 );
-                    sgdNormalizeVec3( normal );
-
-                   p = Point3D( normal[0], normal[1], normal[2] );
-                   // cout << "normal = " << p << endl;
-                   normals.push_back( p );
-                   ptr += 3;
-               }
-           }
-       } else if ( obj_type == SG_TEXCOORD_LIST ) {
-           // read texcoord list properties
-           for ( j = 0; j < nproperties; ++j ) {
-               char prop_type;
-               sgReadChar( fp, &prop_type );
-
-               sgReadUInt( fp, &nbytes );
-               // cout << "property size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-           }
-
-           // read texcoord list elements
-           for ( j = 0; j < nelements; ++j ) {
-               sgReadUInt( fp, &nbytes );
-               // cout << "element size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-               int count = nbytes / (sizeof(float) * 2);
-               float *fptr = (float *)ptr;
-               for ( k = 0; k < count; ++k ) {
-                   if ( sgIsBigEndian() ) {
-                       sgEndianSwap( (unsigned int *)&(fptr[0]) );
-                       sgEndianSwap( (unsigned int *)&(fptr[1]) );
-                   }
-                   p = Point3D( fptr[0], fptr[1], 0 );
-                   // cout << "texcoord = " << p << endl;
-                   texcoords.push_back( p );
-                   fptr += 2;
-               }
-           }
-       } else if ( obj_type == SG_TRIANGLE_FACES ) {
-           // read triangle face properties
-           for ( j = 0; j < nproperties; ++j ) {
-               char prop_type;
-               sgReadChar( fp, &prop_type );
-
-               sgReadUInt( fp, &nbytes );
-               // cout << "property size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-               if ( prop_type == SG_MATERIAL ) {
-                   strncpy( material, ptr, nbytes );
-                   material[nbytes] = '\0';
-                   // cout << "material type = " << material << endl;
-               }
-           }
-
-           // read triangle face elements
-           for ( j = 0; j < nelements; ++j ) {
-               sgReadUInt( fp, &nbytes );
-               // cout << "element size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-               int count = nbytes / (sizeof(short) * 2);
-               short *sptr = (short *)ptr;
-               int_list vs, tcs;
-               vs.clear(); tcs.clear();
-               for ( k = 0; k < count; ++k ) {
-                   if ( sgIsBigEndian() ) {
-                       sgEndianSwap( (unsigned short *)&(sptr[0]) );
-                       sgEndianSwap( (unsigned short *)&(sptr[1]) );
-                   }
-                   vs.push_back( sptr[0] );
-                   tcs.push_back( sptr[1] );
-                   // cout << sptr[0] << "/" << sptr[1] << " ";
-                   sptr += 2;
-               }
-               // cout << endl;
-               tris_v.push_back( vs );
-               tris_tc.push_back( tcs );
-               tri_materials.push_back( material );
-           }
-       } else if ( obj_type == SG_TRIANGLE_STRIPS ) {
-           // read triangle strip properties
-           for ( j = 0; j < nproperties; ++j ) {
-               char prop_type;
-               sgReadChar( fp, &prop_type );
-
-               sgReadUInt( fp, &nbytes );
-               // cout << "property size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-               if ( prop_type == SG_MATERIAL ) {
-                   strncpy( material, ptr, nbytes );
-                   material[nbytes] = '\0';
-                   // cout << "material type = " << material << endl;
-               }
-           }
-
-           // read triangle strip elements
-           for ( j = 0; j < nelements; ++j ) {
-               sgReadUInt( fp, &nbytes );
-               // cout << "element size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-               int count = nbytes / (sizeof(short) * 2);
-               short *sptr = (short *)ptr;
-               int_list vs, tcs;
-               vs.clear(); tcs.clear();
-               for ( k = 0; k < count; ++k ) {
-                   if ( sgIsBigEndian() ) {
-                       sgEndianSwap( (unsigned short *)&(sptr[0]) );
-                       sgEndianSwap( (unsigned short *)&(sptr[1]) );
-                   }
-                   vs.push_back( sptr[0] );
-                   tcs.push_back( sptr[1] );
-                   // cout << sptr[0] << "/" << sptr[1] << " ";
-                   sptr += 2;
-               }
-               // cout << endl;
-               strips_v.push_back( vs );
-               strips_tc.push_back( tcs );
-               strip_materials.push_back( material );
-           }
-       } else if ( obj_type == SG_TRIANGLE_FANS ) {
-           // read triangle fan properties
-           for ( j = 0; j < nproperties; ++j ) {
-               char prop_type;
-               sgReadChar( fp, &prop_type );
-
-               sgReadUInt( fp, &nbytes );
-               // cout << "property size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-               if ( prop_type == SG_MATERIAL ) {
-                   strncpy( material, ptr, nbytes );
-                   material[nbytes] = '\0';
-                   // cout << "material type = " << material << endl;
-               }
-           }
-
-           // read triangle fan elements
-           for ( j = 0; j < nelements; ++j ) {
-               sgReadUInt( fp, &nbytes );
-               // cout << "element size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-               int count = nbytes / (sizeof(short) * 2);
-               short *sptr = (short *)ptr;
-               int_list vs, tcs;
-               vs.clear(); tcs.clear();
-               for ( k = 0; k < count; ++k ) {
-                   if ( sgIsBigEndian() ) {
-                       sgEndianSwap( (unsigned short *)&(sptr[0]) );
-                       sgEndianSwap( (unsigned short *)&(sptr[1]) );
-                   }
-                   vs.push_back( sptr[0] );
-                   tcs.push_back( sptr[1] );
-                   // cout << sptr[0] << "/" << sptr[1] << " ";
-                   sptr += 2;
-               }
-               // cout << endl;
-               fans_v.push_back( vs );
-               fans_tc.push_back( tcs );
-               fan_materials.push_back( material );
-           }
-       } else {
-           // unknown object type, just skip
-
-           // read properties
-           for ( j = 0; j < nproperties; ++j ) {
-               char prop_type;
-               sgReadChar( fp, &prop_type );
-
-               sgReadUInt( fp, &nbytes );
-               // cout << "property size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-           }
-
-           // read elements
-           for ( j = 0; j < nelements; ++j ) {
-               sgReadUInt( fp, &nbytes );
-               // cout << "element size = " << nbytes << endl;
-               if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
-               char *ptr = buf.get_ptr();
-               sgReadBytes( fp, nbytes, ptr );
-           }
-       }
+        // read object header
+        char obj_type;
+        uint32_t nproperties, nelements;
+        sgReadChar( fp, &obj_type );
+        if ( version >= 10 ) {
+            sgReadUInt( fp, &nproperties );
+            sgReadUInt( fp, &nelements );
+        } else if ( version >= 7 ) {
+            uint16_t v;
+            sgReadUShort( fp, &v );
+            nproperties = v;
+            sgReadUShort( fp, &v );
+            nelements = v;
+        } else {
+            int16_t v;
+            sgReadShort( fp, &v );
+            nproperties = v;
+            sgReadShort( fp, &v );
+            nelements = v;
+        }
+
+        SG_LOG(SG_IO, SG_DEBUG, "SGBinObject::read_bin object " << i << 
+                " = " << (int)obj_type << " props = " << nproperties << 
+                " elements = " << nelements);
+         
+        if ( obj_type == SG_BOUNDING_SPHERE ) {
+            // read bounding sphere properties
+            read_properties( fp, nproperties );
+            
+            // read bounding sphere elements
+            for ( j = 0; j < nelements; ++j ) {
+                sgReadUInt( fp, &nbytes );
+                buf.resize( nbytes );
+                buf.reset();
+                char *ptr = buf.get_ptr();
+                sgReadBytes( fp, nbytes, ptr );
+                gbs_center = buf.readVec3d();
+                gbs_radius = buf.readFloat();
+            }
+        } else if ( obj_type == SG_VERTEX_LIST ) {
+            // read vertex list properties
+            read_properties( fp, nproperties );
+
+            // read vertex list elements
+            for ( j = 0; j < nelements; ++j ) {
+                sgReadUInt( fp, &nbytes );
+                buf.resize( nbytes );
+                buf.reset();
+                char *ptr = buf.get_ptr();
+                sgReadBytes( fp, nbytes, ptr );
+                int count = nbytes / (sizeof(float) * 3);
+                wgs84_nodes.reserve( count );
+                for ( k = 0; k < count; ++k ) {
+                    SGVec3f v = buf.readVec3f();
+                    // extend from float to double, hmmm
+                    wgs84_nodes.push_back( SGVec3d(v[0], v[1], v[2]) );
+                }
+            }
+        } else if ( obj_type == SG_COLOR_LIST ) {
+            // read color list properties
+            read_properties( fp, nproperties );
+
+            // read color list elements
+            for ( j = 0; j < nelements; ++j ) {
+                sgReadUInt( fp, &nbytes );
+                buf.resize( nbytes );
+                buf.reset();
+                char *ptr = buf.get_ptr();
+                sgReadBytes( fp, nbytes, ptr );
+                int count = nbytes / (sizeof(float) * 4);
+                colors.reserve(count);
+                for ( k = 0; k < count; ++k ) {
+                    colors.push_back( buf.readVec4f() );
+                }
+            }
+        } else if ( obj_type == SG_NORMAL_LIST ) {
+            // read normal list properties
+            read_properties( fp, nproperties );
+
+            // read normal list elements
+            for ( j = 0; j < nelements; ++j ) {
+                sgReadUInt( fp, &nbytes );
+                buf.resize( nbytes );
+                buf.reset();
+                unsigned char *ptr = (unsigned char *)(buf.get_ptr());
+                sgReadBytes( fp, nbytes, ptr );
+                int count = nbytes / 3;
+                normals.reserve( count );
+                for ( k = 0; k < count; ++k ) {
+                    SGVec3f normal( (ptr[0]) / 127.5 - 1.0,
+                                    (ptr[1]) / 127.5 - 1.0, 
+                                    (ptr[2]) / 127.5 - 1.0);
+                    normals.push_back(normalize(normal));
+                    ptr += 3;
+                }
+            }
+        } else if ( obj_type == SG_TEXCOORD_LIST ) {
+            // read texcoord list properties
+            read_properties( fp, nproperties );
+
+            // read texcoord list elements
+            for ( j = 0; j < nelements; ++j ) {
+                sgReadUInt( fp, &nbytes );
+                buf.resize( nbytes );
+                buf.reset();
+                char *ptr = buf.get_ptr();
+                sgReadBytes( fp, nbytes, ptr );
+                int count = nbytes / (sizeof(float) * 2);
+                texcoords.reserve(count);
+                for ( k = 0; k < count; ++k ) {
+                    texcoords.push_back( buf.readVec2f() );
+                }
+            }
+        } else if ( obj_type == SG_VA_FLOAT_LIST ) {
+            // read vertex attribute (float) properties
+            read_properties( fp, nproperties );
+            
+            // read vertex attribute list elements
+            for ( j = 0; j < nelements; ++j ) {
+                sgReadUInt( fp, &nbytes );
+                buf.resize( nbytes );
+                buf.reset();
+                char *ptr = buf.get_ptr();
+                sgReadBytes( fp, nbytes, ptr );
+                int count = nbytes / (sizeof(float));
+                va_flt.reserve(count);
+                for ( k = 0; k < count; ++k ) {
+                    va_flt.push_back( buf.readFloat() );
+                }
+            }
+        } else if ( obj_type == SG_VA_INTEGER_LIST ) {
+            // read vertex attribute (integer) properties
+            read_properties( fp, nproperties );
+            
+            // read vertex attribute list elements
+            for ( j = 0; j < nelements; ++j ) {
+                sgReadUInt( fp, &nbytes );
+                buf.resize( nbytes );
+                buf.reset();
+                char *ptr = buf.get_ptr();
+                sgReadBytes( fp, nbytes, ptr );
+                int count = nbytes / (sizeof(unsigned int));
+                va_int.reserve(count);
+                for ( k = 0; k < count; ++k ) {
+                    va_int.push_back( buf.readInt() );
+                }
+            }
+        } else if ( obj_type == SG_POINTS ) {
+            // read point elements
+            read_object( fp, SG_POINTS, nproperties, nelements,
+                         pts_v, pts_n, pts_c, pts_tcs,
+                         pts_vas, pt_materials );
+        } else if ( obj_type == SG_TRIANGLE_FACES ) {
+            // read triangle face properties
+            read_object( fp, SG_TRIANGLE_FACES, nproperties, nelements,
+                         tris_v, tris_n, tris_c, tris_tcs,
+                         tris_vas, tri_materials );
+        } else if ( obj_type == SG_TRIANGLE_STRIPS ) {
+            // read triangle strip properties
+            read_object( fp, SG_TRIANGLE_STRIPS, nproperties, nelements,
+                         strips_v, strips_n, strips_c, strips_tcs,
+                         strips_vas, strip_materials );
+        } else if ( obj_type == SG_TRIANGLE_FANS ) {
+            // read triangle fan properties
+            read_object( fp, SG_TRIANGLE_FANS, nproperties, nelements,
+                         fans_v, fans_n, fans_c, fans_tcs,
+                         fans_vas, fan_materials );
+        } else {
+            // unknown object type, just skip
+            read_properties( fp, nproperties );
+
+            // read elements
+            for ( j = 0; j < nelements; ++j ) {
+                sgReadUInt( fp, &nbytes );
+                // cout << "element size = " << nbytes << endl;
+                if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
+                char *ptr = buf.get_ptr();
+                sgReadBytes( fp, nbytes, ptr );
+            }
+        }
+        
+        if ( sgReadError() ) {
+            throw sg_io_exception("Error while reading object", sg_location(file, i));
+        }
     }
 
     // close the file
     gzclose(fp);
 
-    if ( sgReadError() ) {
-       cout << "We detected an error while reading the file." << endl;
-       return false;
+    return true;
+}
+
+void SGBinObject::write_header(gzFile fp, int type, int nProps, int nElements)
+{
+    sgWriteChar(fp, (unsigned char) type);
+    if (version == 7) {
+        sgWriteUShort(fp, nProps);
+        sgWriteUShort(fp, nElements);
+    } else {
+        sgWriteUInt(fp, nProps);
+        sgWriteUInt(fp, nElements);
     }
+}
 
-    return true;
+unsigned int SGBinObject::count_objects(const string_list& materials)
+{
+    unsigned int result = 0;
+    unsigned int start = 0, end = 1;
+    unsigned int count = materials.size();
+    string m;
+    
+    while ( start < count ) {
+        m = materials[start];
+        for (end = start+1; (end < count) && (m == materials[end]); ++end) { }     
+        ++result;
+        start = end; 
+    }
+    
+    return result;
 }
 
+void SGBinObject::write_objects(gzFile fp, int type, 
+                                const group_list& verts,
+                                const group_list& normals, 
+                                const group_list& colors, 
+                                const group_tci_list& texCoords,
+                                const group_vai_list& vertexAttribs, 
+                                const string_list& materials)
+{
+    if (verts.empty()) {
+        return;
+    }
+        
+    unsigned int start = 0, end = 1;
+    string m;
+    int_list emptyList;
+    
+    while (start < materials.size()) {
+        m = materials[start];
+        // find range of objects with identical material, write out as a single object
+        for (end = start+1; (end < materials.size()) && (m == materials[end]); ++end) {}
+    
+        // calc the number of elements
+        const int count = end - start;
+
+        // calc the number of properties
+        unsigned int va_mask = 0;
+        unsigned int va_count = vertexAttribs.size();
+        for ( unsigned int va=0; va<va_count; va++ ) {
+            if ( !vertexAttribs[va].empty() && !vertexAttribs[va].front().empty() ) {
+                va_mask |= ( 1 << va );
+            }
+        }
+
+        if ( va_mask ) {
+            write_header(fp, type, 3, count);
+        } else {
+            write_header(fp, type, 2, count);
+        }
+        
+    // properties
+        // material property
+        sgWriteChar( fp, (char)SG_MATERIAL );        // property
+        sgWriteUInt( fp, m.length() );        // nbytes
+        sgWriteBytes( fp, m.length(), m.c_str() );
+    
+        // index mask property
+        unsigned char idx_mask = 0;
+        if ( !verts.empty() && !verts[start].empty()) idx_mask |= SG_IDX_VERTICES;
+        if ( !normals.empty() && !normals[start].empty()) idx_mask |= SG_IDX_NORMALS;
+        if ( !colors.empty() && !colors[start].empty()) idx_mask |= SG_IDX_COLORS;
+        if ( !texCoords.empty() && !texCoords[start][0].empty()) idx_mask |= SG_IDX_TEXCOORDS_0;
+        if ( !texCoords.empty() && !texCoords[start][1].empty()) idx_mask |= SG_IDX_TEXCOORDS_1;
+        if ( !texCoords.empty() && !texCoords[start][2].empty()) idx_mask |= SG_IDX_TEXCOORDS_2;
+        if ( !texCoords.empty() && !texCoords[start][3].empty()) idx_mask |= SG_IDX_TEXCOORDS_3;
+        
+        if (idx_mask == 0) {
+            SG_LOG(SG_IO, SG_ALERT, "SGBinObject::write_objects: object with material:"
+                << m << "has no indices set");
+        }
+
+        sgWriteChar( fp, (char)SG_INDEX_TYPES );     // property
+        sgWriteUInt( fp, 1 );                        // nbytes
+        sgWriteChar( fp, idx_mask );
+
+        // vertex attribute property
+        if (va_mask != 0) {
+            sgWriteChar( fp, (char)SG_VERT_ATTRIBS );    // property
+            sgWriteUInt( fp, 4 );                        // nbytes
+            sgWriteChar( fp, va_mask );
+        }
+        
+    // elements
+        for (unsigned int i=start; i < end; ++i) {
+            const int_list& va(verts[i]);
+            const int_list& na((idx_mask & SG_IDX_NORMALS) ? normals[i] : emptyList);
+            const int_list& ca((idx_mask & SG_IDX_COLORS) ? colors[i] : emptyList);
+
+            // pass the whole texcoord array - we'll figure out which indicies to write 
+            // in write_indices
+            const tci_list& tca( texCoords[i] );
+
+            // pass the whole vertex array - we'll figure out which indicies to write 
+            // in write_indices
+            const vai_list& vaa( vertexAttribs[i] );
+            
+            if (version == 7) {
+                write_indices<uint16_t>(fp, idx_mask, va_mask, va, na, ca, tca, vaa);
+            } else {
+                write_indices<uint32_t>(fp, idx_mask, va_mask, va, na, ca, tca, vaa);
+            }
+        }
+    
+        start = end;
+    } // of materials iteration
+}
 
 // write out the structures to a binary file.  We assume that the
 // groups come to us sorted by material property.  If not, things
 // don't break, but the result won't be as optimal.
 bool SGBinObject::write_bin( const string& base, const string& name,
-                            const SGBucket& b )
+                             const SGBucket& b )
 {
-    Point3D p;
-    sgVec2 t;
-    sgVec3 pt;
-    int i, j;
 
-    string dir = base + "/" + b.gen_base_path();
-    string command = "mkdir -p " + dir;
-#if defined(_MSC_VER) || defined(__MINGW32__)
-    system( (string("mkdir ") + dir).c_str() );
-#else
-    system(command.c_str());
-#endif
+    SGPath file = base + "/" + b.gen_base_path() + "/" + name + ".gz";
+    return write_bin_file(file);
+}
+
+static unsigned int max_object_size( const string_list& materials )
+{
+    unsigned int max_size = 0;
+
+    for (unsigned int start=0; start < materials.size();) {
+        string m = materials[start];
+        unsigned int end = start + 1;
+        // find range of objects with identical material, calc its size
+        for (; (end < materials.size()) && (m == materials[end]); ++end) {}
+
+        unsigned int cur_size = end - start;
+        max_size = std::max(max_size, cur_size);
+        start = end;
+    }
+
+    return max_size;
+}
 
-    string file = dir + "/" + name + ".gz";
-    cout << "Output file = " << file << endl;
+const unsigned int VERSION_7_MATERIAL_LIMIT = 0x7fff;
+
+bool SGBinObject::write_bin_file(const SGPath& file)
+{
+    int i;
+    
+    SGPath file2(file);
+    file2.create_dir( 0755 );
 
     gzFile fp;
     if ( (fp = gzopen( file.c_str(), "wb9" )) == NULL ) {
-       cout << "ERROR: opening " << file << " for writing!" << endl;
-       return false;
+        cout << "ERROR: opening " << file.str() << " for writing!" << endl;
+        return false;
     }
 
     sgClearWriteError();
 
-    cout << "triangles size = " << tris_v.size() << "  tri_materials = " 
-        << tri_materials.size() << endl;
-    cout << "strips size = " << strips_v.size() << "  strip_materials = " 
-        << strip_materials.size() << endl;
-    cout << "fans size = " << fans_v.size() << "  fan_materials = " 
-        << fan_materials.size() << endl;
-
-    cout << "points = " << wgs84_nodes.size() << endl;
-    cout << "tex coords = " << texcoords.size() << endl;
+    SG_LOG(SG_IO, SG_DEBUG, "points size = " << pts_v.size() 
+         << "  pt_materials = " << pt_materials.size() );
+    SG_LOG(SG_IO, SG_DEBUG, "triangles size = " << tris_v.size() 
+         << "  tri_materials = " << tri_materials.size() );
+    SG_LOG(SG_IO, SG_DEBUG, "strips size = " << strips_v.size() 
+         << "  strip_materials = " << strip_materials.size() );
+    SG_LOG(SG_IO, SG_DEBUG, "fans size = " << fans_v.size() 
+         << "  fan_materials = " << fan_materials.size() );
+
+    SG_LOG(SG_IO, SG_DEBUG, "nodes = " << wgs84_nodes.size() );
+    SG_LOG(SG_IO, SG_DEBUG, "colors = " << colors.size() );
+    SG_LOG(SG_IO, SG_DEBUG, "normals = " << normals.size() );
+    SG_LOG(SG_IO, SG_DEBUG, "tex coords = " << texcoords.size() );
+
+    version = 10;
+    bool shortMaterialsRanges = 
+        (max_object_size(pt_materials) < VERSION_7_MATERIAL_LIMIT) &&
+        (max_object_size(fan_materials) < VERSION_7_MATERIAL_LIMIT) &&
+        (max_object_size(strip_materials) < VERSION_7_MATERIAL_LIMIT) &&
+        (max_object_size(tri_materials) < VERSION_7_MATERIAL_LIMIT);
+                    
+    if ((wgs84_nodes.size() < 0xffff) &&
+        (normals.size() < 0xffff) &&
+        (texcoords.size() < 0xffff) &&
+        shortMaterialsRanges) {
+        version = 7; // use smaller indices if possible
+    }
 
     // write header magic
+    
+    /** Magic Number for our file format */
+    #define SG_FILE_MAGIC_NUMBER  ( ('S'<<24) + ('G'<<16) + version )
+    
     sgWriteUInt( fp, SG_FILE_MAGIC_NUMBER );
     time_t calendar_time = time(NULL);
-    sgWriteLong( fp, (long int)calendar_time );
+    sgWriteLong( fp, (int32_t)calendar_time );
 
     // calculate and write number of top level objects
-    string material;
-    int start;
-    int end;
-    short nobjects = 0;
-    nobjects++;                        // for gbs
-    nobjects++;                        // for vertices
-    nobjects++;                        // for normals
-    nobjects++;                        // for texcoords
-
-    // tris
-    short ntris = 0;
-    start = 0; end = 1;
-    while ( start < (int)tri_materials.size() ) {
-       material = tri_materials[start];
-       while ( (end < (int)tri_materials.size()) &&
-               (material == tri_materials[end]) ) {
-           end++;
-       }
-       ntris++;
-       start = end; end = start + 1;
-    }
-    nobjects += ntris;
-
-    // strips
-    short nstrips = 0;
-    start = 0; end = 1;
-    while ( start < (int)strip_materials.size() ) {
-       material = strip_materials[start];
-       while ( (end < (int)strip_materials.size()) &&
-               (material == strip_materials[end]) ) {
-           end++;
-       }
-       nstrips++;
-       start = end; end = start + 1;
-    }
-    nobjects += nstrips;
-
-    // fans
-    short nfans = 0;
-    start = 0; end = 1;
-    while ( start < (int)fan_materials.size() ) {
-       material = fan_materials[start];
-       while ( (end < (int)fan_materials.size()) &&
-               (material == fan_materials[end]) ) {
-           end++;
-       }
-       nfans++;
-       start = end; end = start + 1;
-    }
-    nobjects += nfans;
-
-    cout << "total top level objects = " << nobjects << endl;
-    sgWriteShort( fp, nobjects );
+    int nobjects = 5; // gbs, vertices, colors, normals, texcoords
+    nobjects += count_objects(pt_materials);
+    nobjects += count_objects(tri_materials);
+    nobjects += count_objects(strip_materials);
+    nobjects += count_objects(fan_materials);
 
-    // write bounding sphere
-    sgWriteChar( fp, (char)SG_BOUNDING_SPHERE );        // type
-    sgWriteShort( fp, 0 );                             // nproperties
-    sgWriteShort( fp, 1 );                             // nelements
+    SG_LOG(SG_IO, SG_DEBUG, "total top level objects = " << nobjects);
 
+    if (version == 7) {
+        sgWriteUShort( fp, (uint16_t) nobjects );
+    } else {
+        sgWriteInt( fp, nobjects );
+    } 
+    
+    // write bounding sphere
+    write_header( fp, SG_BOUNDING_SPHERE, 0, 1);
     sgWriteUInt( fp, sizeof(double) * 3 + sizeof(float) ); // nbytes
-    sgdVec3 center;
-    sgdSetVec3( center, gbs_center.x(), gbs_center.y(), gbs_center.z() );
-    sgWritedVec3( fp, center );
+    sgWritedVec3( fp, gbs_center );
     sgWriteFloat( fp, gbs_radius );
 
     // dump vertex list
-    sgWriteChar( fp, (char)SG_VERTEX_LIST );             // type
-    sgWriteShort( fp, 0 );                              // nproperties
-    sgWriteShort( fp, 1 );                              // nelements
+    write_header( fp, SG_VERTEX_LIST, 0, 1);
     sgWriteUInt( fp, wgs84_nodes.size() * sizeof(float) * 3 ); // nbytes
     for ( i = 0; i < (int)wgs84_nodes.size(); ++i ) {
-       p = wgs84_nodes[i] - gbs_center;
-       sgSetVec3( pt, p.x(), p.y(), p.z() );
-       sgWriteVec3( fp, pt );
+        sgWriteVec3( fp, toVec3f(wgs84_nodes[i] - gbs_center));
+    }
+
+    // dump vertex color list
+    write_header( fp, SG_COLOR_LIST, 0, 1);
+    sgWriteUInt( fp, colors.size() * sizeof(float) * 4 ); // nbytes
+    for ( i = 0; i < (int)colors.size(); ++i ) {
+      sgWriteVec4( fp, colors[i]);
     }
 
     // dump vertex normal list
-    sgWriteChar( fp, (char)SG_NORMAL_LIST );            // type
-    sgWriteShort( fp, 0 );                             // nproperties
-    sgWriteShort( fp, 1 );                             // nelements
+    write_header( fp, SG_NORMAL_LIST, 0, 1);
     sgWriteUInt( fp, normals.size() * 3 );              // nbytes
     char normal[3];
     for ( i = 0; i < (int)normals.size(); ++i ) {
-       p = normals[i];
-       normal[0] = (unsigned char)((p.x() + 1.0) * 127.5);
-       normal[1] = (unsigned char)((p.y() + 1.0) * 127.5);
-       normal[2] = (unsigned char)((p.z() + 1.0) * 127.5);
-       sgWriteBytes( fp, 3, normal );
+        SGVec3f p = normals[i];
+        normal[0] = (unsigned char)((p.x() + 1.0) * 127.5);
+        normal[1] = (unsigned char)((p.y() + 1.0) * 127.5);
+        normal[2] = (unsigned char)((p.z() + 1.0) * 127.5);
+        sgWriteBytes( fp, 3, normal );
     }
 
     // dump texture coordinates
-    sgWriteChar( fp, (char)SG_TEXCOORD_LIST );          // type
-    sgWriteShort( fp, 0 );                             // nproperties
-    sgWriteShort( fp, 1 );                             // nelements
+    write_header( fp, SG_TEXCOORD_LIST, 0, 1);
     sgWriteUInt( fp, texcoords.size() * sizeof(float) * 2 ); // nbytes
     for ( i = 0; i < (int)texcoords.size(); ++i ) {
-       p = texcoords[i];
-       sgSetVec2( t, p.x(), p.y() );
-       sgWriteVec2( fp, t );
-    }
-
-    // dump individual triangles if they exist
-    if ( tris_v.size() > 0 ) {
-       int start = 0;
-       int end = 1;
-       string material;
-       while ( start < (int)tri_materials.size() ) {
-           // find next group
-           material = tri_materials[start];
-           while ( (end < (int)tri_materials.size()) && 
-                   (material == tri_materials[end]) )
-           {
-               // cout << "end = " << end << endl;
-               end++;
-           }
-           // cout << "group = " << start << " to " << end - 1 << endl;
-
-           // write group headers
-           sgWriteChar( fp, (char)SG_TRIANGLE_FACES ); // type
-           sgWriteShort( fp, 1 );                      // nproperties
-           sgWriteShort( fp, 1 );                      // nelements
-
-           sgWriteChar( fp, (char)SG_MATERIAL );       // property
-           sgWriteUInt( fp, material.length() );        // nbytes
-           sgWriteBytes( fp, material.length(), material.c_str() );
-
-           sgWriteUInt( fp, (end - start) * 3 * 2 * sizeof(short) ); // nbytes
-
-           // write group
-           for ( i = start; i < end; ++i ) {
-               for ( j = 0; j < 3; ++j ) {
-                   sgWriteShort( fp, (short)tris_v[i][j] );
-                   sgWriteShort( fp, (short)tris_tc[i][j] );
-               }
-           }
-
-           start = end;
-           end = start + 1;
-       }
-    }
-
-    // dump triangle strips
-    if ( strips_v.size() > 0 ) {
-       int start = 0;
-       int end = 1;
-       string material;
-       while ( start < (int)strip_materials.size() ) {
-           // find next group
-           material = strip_materials[start];
-           while ( (end < (int)strip_materials.size()) && 
-                   (material == strip_materials[end]) )
-               {
-                   // cout << "end = " << end << endl;
-                   end++;
-               }
-           // cout << "group = " << start << " to " << end - 1 << endl;
-
-           // write group headers
-           sgWriteChar( fp, (char)SG_TRIANGLE_STRIPS ); // type
-           sgWriteShort( fp, 1 );                       // nproperties
-           sgWriteShort( fp, end - start );             // nelements
-
-           sgWriteChar( fp, (char)SG_MATERIAL );        // property
-           sgWriteUInt( fp, material.length() );        // nbytes
-           sgWriteBytes( fp, material.length(), material.c_str() );
-
-           // write strips
-           for ( i = start; i < end; ++i ) {
-               // nbytes
-               sgWriteUInt( fp, strips_v[i].size() * 2 * sizeof(short) );
-               for ( j = 0; j < (int)strips_v[i].size(); ++j ) {
-                   sgWriteShort( fp, (short)strips_v[i][j] );
-                   sgWriteShort( fp, (short)strips_tc[i][j] );
-               }
-           }
-           
-           start = end;
-           end = start + 1;
-       }
-    }
-
-    // dump triangle fans
-    if ( fans_v.size() > 0 ) {
-       int start = 0;
-       int end = 1;
-       string material;
-       while ( start < (int)fan_materials.size() ) {
-           // find next group
-           material = fan_materials[start];
-           while ( (end < (int)fan_materials.size()) && 
-                   (material == fan_materials[end]) )
-               {
-                   // cout << "end = " << end << endl;
-                   end++;
-               }
-           // cout << "group = " << start << " to " << end - 1 << endl;
-
-           // write group headers
-           sgWriteChar( fp, (char)SG_TRIANGLE_FANS );   // type
-           sgWriteShort( fp, 1 );                       // nproperties
-           sgWriteShort( fp, end - start );             // nelements
-
-           sgWriteChar( fp, (char)SG_MATERIAL );       // property
-           sgWriteUInt( fp, material.length() );        // nbytes
-           sgWriteBytes( fp, material.length(), material.c_str() );
-
-           // write fans
-           for ( i = start; i < end; ++i ) {
-               // nbytes
-               sgWriteUInt( fp, fans_v[i].size() * 2 * sizeof(short) );
-               for ( j = 0; j < (int)fans_v[i].size(); ++j ) {
-                   sgWriteShort( fp, (short)fans_v[i][j] );
-                   sgWriteShort( fp, (short)fans_tc[i][j] );
-               }
-           }
-           
-           start = end;
-           end = start + 1;
-       }
+      sgWriteVec2( fp, texcoords[i]);
     }
 
+    write_objects(fp, SG_POINTS, pts_v, pts_n, pts_c, pts_tcs, pts_vas, pt_materials);
+    write_objects(fp, SG_TRIANGLE_FACES, tris_v, tris_n, tris_c, tris_tcs, tris_vas, tri_materials);
+    write_objects(fp, SG_TRIANGLE_STRIPS, strips_v, strips_n, strips_c, strips_tcs, strips_vas, strip_materials);
+    write_objects(fp, SG_TRIANGLE_FANS, fans_v, fans_n, fans_c, fans_tcs, fans_vas, fan_materials);
+    
     // close the file
     gzclose(fp);
 
     if ( sgWriteError() ) {
-       cout << "We detected an error while writing the file." << endl;
-       return false;
+        cout << "Error while writing file " << file.str() << endl;
+        return false;
     }
 
     return true;
@@ -828,35 +1076,26 @@ bool SGBinObject::write_bin( const string& base, const string& name,
 // groups come to us sorted by material property.  If not, things
 // don't break, but the result won't be as optimal.
 bool SGBinObject::write_ascii( const string& base, const string& name,
-                              const SGBucket& b )
+                               const SGBucket& b )
 {
-    Point3D p;
     int i, j;
 
-    string dir = base + "/" + b.gen_base_path();
-    string command = "mkdir -p " + dir;
-#if defined(_MSC_VER) || defined(__MINGW32__)
-    system( (string("mkdir ") + dir).c_str() );
-#else
-    system(command.c_str());
-#endif
-
-    // string file = dir + "/" + b.gen_index_str();
-    string file = dir + "/" + name;
-    cout << "Output file = " << file << endl;
+    SGPath file = base + "/" + b.gen_base_path() + "/" + name;
+    file.create_dir( 0755 );
+    cout << "Output file = " << file.str() << endl;
 
     FILE *fp;
     if ( (fp = fopen( file.c_str(), "w" )) == NULL ) {
-       cout << "ERROR: opening " << file << " for writing!" << endl;
-       return false;
+        cout << "ERROR: opening " << file.str() << " for writing!" << endl;
+        return false;
     }
 
     cout << "triangles size = " << tris_v.size() << "  tri_materials = " 
-        << tri_materials.size() << endl;
+         << tri_materials.size() << endl;
     cout << "strips size = " << strips_v.size() << "  strip_materials = " 
-        << strip_materials.size() << endl;
+         << strip_materials.size() << endl;
     cout << "fans size = " << fans_v.size() << "  fan_materials = " 
-        << fan_materials.size() << endl;
+         << fan_materials.size() << endl;
 
     cout << "points = " << wgs84_nodes.size() << endl;
     cout << "tex coords = " << texcoords.size() << endl;
@@ -874,140 +1113,181 @@ bool SGBinObject::write_ascii( const string& base, const string& name,
 
     // write bounding sphere
     fprintf(fp, "# gbs %.5f %.5f %.5f %.2f\n",
-           gbs_center.x(), gbs_center.y(), gbs_center.z(), gbs_radius);
+            gbs_center.x(), gbs_center.y(), gbs_center.z(), gbs_radius);
     fprintf(fp, "\n");
 
     // dump vertex list
     fprintf(fp, "# vertex list\n");
     for ( i = 0; i < (int)wgs84_nodes.size(); ++i ) {
-       p = wgs84_nodes[i] - gbs_center;
-       
-       fprintf(fp,  "v %.5f %.5f %.5f\n", p.x(), p.y(), p.z() );
+        SGVec3d p = wgs84_nodes[i] - gbs_center;
+        
+        fprintf(fp,  "v %.5f %.5f %.5f\n", p.x(), p.y(), p.z() );
     }
     fprintf(fp, "\n");
 
     fprintf(fp, "# vertex normal list\n");
     for ( i = 0; i < (int)normals.size(); ++i ) {
-       p = normals[i];
-       fprintf(fp,  "vn %.5f %.5f %.5f\n", p.x(), p.y(), p.z() );
+        SGVec3f p = normals[i];
+        fprintf(fp,  "vn %.5f %.5f %.5f\n", p.x(), p.y(), p.z() );
     }
     fprintf(fp, "\n");
 
     // dump texture coordinates
     fprintf(fp, "# texture coordinate list\n");
     for ( i = 0; i < (int)texcoords.size(); ++i ) {
-       p = texcoords[i];
-       fprintf(fp,  "vt %.5f %.5f\n", p.x(), p.y() );
+        SGVec2f p = texcoords[i];
+        fprintf(fp,  "vt %.5f %.5f\n", p.x(), p.y() );
     }
     fprintf(fp, "\n");
 
     // dump individual triangles if they exist
-    if ( tris_v.size() > 0 ) {
-       fprintf(fp, "# triangle groups\n");
-
-       int start = 0;
-       int end = 1;
-       string material;
-       while ( start < (int)tri_materials.size() ) {
-           // find next group
-           material = tri_materials[start];
-           while ( (end < (int)tri_materials.size()) && 
-                   (material == tri_materials[end]) )
-           {
-               // cout << "end = " << end << endl;
-               end++;
-           }
-           // cout << "group = " << start << " to " << end - 1 << endl;
-
-           // make a list of points for the group
-           point_list group_nodes;
-           group_nodes.clear();
-           Point3D bs_center;
-           double bs_radius = 0;
-           for ( i = start; i < end; ++i ) {
-               for ( j = 0; j < (int)tris_v[i].size(); ++j ) {
-                   group_nodes.push_back( wgs84_nodes[ tris_v[i][j] ] );
-                   bs_center = calc_center( group_nodes );
-                   bs_radius = sgCalcBoundingRadius( bs_center, group_nodes );
-               }
-           }
-
-           // write group headers
-           fprintf(fp, "\n");
-           fprintf(fp, "# usemtl %s\n", material.c_str());
-           fprintf(fp, "# bs %.4f %.4f %.4f %.2f\n",
-                   bs_center.x(), bs_center.y(), bs_center.z(), bs_radius);
-
-           // write groups
-           for ( i = start; i < end; ++i ) {
-               fprintf(fp, "f");
-               for ( j = 0; j < (int)tris_v[i].size(); ++j ) {
-                   fprintf(fp, " %d/%d", tris_v[i][j], tris_tc[i][j] );
-               }
-               fprintf(fp, "\n");
-           }
-
-           start = end;
-           end = start + 1;
-       }
+    if ( ! tris_v.empty() ) {
+        fprintf(fp, "# triangle groups\n");
+
+        int start = 0;
+        int end = 1;
+        string material;
+        while ( start < (int)tri_materials.size() ) {
+            // find next group
+            material = tri_materials[start];
+            while ( (end < (int)tri_materials.size()) && 
+                    (material == tri_materials[end]) )
+            {
+                // cout << "end = " << end << endl;
+                end++;
+            }
+            // cout << "group = " << start << " to " << end - 1 << endl;
+
+      SGSphered d( SGVec3d(0.0, 0.0, 0.0), -1.0 );
+      for ( i = start; i < end; ++i ) {
+        for ( j = 0; j < (int)tris_v[i].size(); ++j ) {
+          d.expandBy(wgs84_nodes[ tris_v[i][j] ]);
+        }
+      }
+      
+      SGVec3d bs_center = d.getCenter();
+      double bs_radius = d.getRadius();
+            
+            // write group headers
+            fprintf(fp, "\n");
+            fprintf(fp, "# usemtl %s\n", material.c_str());
+            fprintf(fp, "# bs %.4f %.4f %.4f %.2f\n",
+                    bs_center.x(), bs_center.y(), bs_center.z(), bs_radius);
+
+            // write groups
+            for ( i = start; i < end; ++i ) {
+                fprintf(fp, "f");
+                for ( j = 0; j < (int)tris_v[i].size(); ++j ) {
+                    fprintf(fp, " %d/%d", tris_v[i][j], tris_tcs[i][0][j] );
+                }
+                fprintf(fp, "\n");
+            }
+
+            start = end;
+            end = start + 1;
+        }
     }
 
     // dump triangle groups
-    if ( strips_v.size() > 0 ) {
-       fprintf(fp, "# triangle strips\n");
-
-       int start = 0;
-       int end = 1;
-       string material;
-       while ( start < (int)strip_materials.size() ) {
-           // find next group
-           material = strip_materials[start];
-           while ( (end < (int)strip_materials.size()) && 
-                   (material == strip_materials[end]) )
-               {
-                   // cout << "end = " << end << endl;
-                   end++;
-               }
-           // cout << "group = " << start << " to " << end - 1 << endl;
-
-           // make a list of points for the group
-           point_list group_nodes;
-           group_nodes.clear();
-           Point3D bs_center;
-           double bs_radius = 0;
-           for ( i = start; i < end; ++i ) {
-               for ( j = 0; j < (int)strips_v[i].size(); ++j ) {
-                   group_nodes.push_back( wgs84_nodes[ strips_v[i][j] ] );
-                   bs_center = calc_center( group_nodes );
-                   bs_radius = sgCalcBoundingRadius( bs_center, group_nodes );
-               }
-           }
-
-           // write group headers
-           fprintf(fp, "\n");
-           fprintf(fp, "# usemtl %s\n", material.c_str());
-           fprintf(fp, "# bs %.4f %.4f %.4f %.2f\n",
-                   bs_center.x(), bs_center.y(), bs_center.z(), bs_radius);
-
-           // write groups
-           for ( i = start; i < end; ++i ) {
-               fprintf(fp, "ts");
-               for ( j = 0; j < (int)strips_v[i].size(); ++j ) {
-                   fprintf(fp, " %d/%d", strips_v[i][j], strips_tc[i][j] );
-               }
-               fprintf(fp, "\n");
-           }
-           
-           start = end;
-           end = start + 1;
-       }
+    if ( ! strips_v.empty() ) {
+        fprintf(fp, "# triangle strips\n");
+
+        int start = 0;
+        int end = 1;
+        string material;
+        while ( start < (int)strip_materials.size() ) {
+            // find next group
+            material = strip_materials[start];
+            while ( (end < (int)strip_materials.size()) && 
+                    (material == strip_materials[end]) )
+                {
+                    // cout << "end = " << end << endl;
+                    end++;
+                }
+            // cout << "group = " << start << " to " << end - 1 << endl;
+
+
+      SGSphered d( SGVec3d(0.0, 0.0, 0.0), -1.0 );
+      for ( i = start; i < end; ++i ) {
+        for ( j = 0; j < (int)tris_v[i].size(); ++j ) {
+          d.expandBy(wgs84_nodes[ tris_v[i][j] ]);
+        }
+      }
+      
+      SGVec3d bs_center = d.getCenter();
+      double bs_radius = d.getRadius();
+
+            // write group headers
+            fprintf(fp, "\n");
+            fprintf(fp, "# usemtl %s\n", material.c_str());
+            fprintf(fp, "# bs %.4f %.4f %.4f %.2f\n",
+                    bs_center.x(), bs_center.y(), bs_center.z(), bs_radius);
+
+            // write groups
+            for ( i = start; i < end; ++i ) {
+                fprintf(fp, "ts");
+                for ( j = 0; j < (int)strips_v[i].size(); ++j ) {
+                    fprintf(fp, " %d/%d", strips_v[i][j], strips_tcs[i][0][j] );
+                }
+                fprintf(fp, "\n");
+            }
+            
+            start = end;
+            end = start + 1;
+        }
     }
 
     // close the file
     fclose(fp);
 
-    command = "gzip --force --best " + file;
-    system(command.c_str());
+    string command = "gzip --force --best " + file.str();
+    int err = system(command.c_str());
+    if (err)
+    {
+        cout << "ERROR: gzip " << file.str() << " failed!" << endl;
+    }
+
+    return (err == 0);
+}
 
+void SGBinObject::read_properties(gzFile fp, int nproperties)
+{
+    sgSimpleBuffer buf;
+    uint32_t nbytes;
+    
+    // read properties
+    for ( int j = 0; j < nproperties; ++j ) {
+        char prop_type;
+        sgReadChar( fp, &prop_type );
+        sgReadUInt( fp, &nbytes );
+        // cout << "property size = " << nbytes << endl;
+        if ( nbytes > buf.get_size() ) { buf.resize( nbytes ); }
+        char *ptr = buf.get_ptr();
+        sgReadBytes( fp, nbytes, ptr );
+    }
+}
+
+bool SGBinObject::add_point( const SGBinObjectPoint& pt )
+{
+    // add the point info
+    pt_materials.push_back( pt.material );
+    
+    pts_v.push_back( pt.v_list );
+    pts_n.push_back( pt.n_list );    
+    pts_c.push_back( pt.c_list );
+    
     return true;
 }
+
+bool SGBinObject::add_triangle( const SGBinObjectTriangle& tri )
+{
+    // add the triangle info and keep lists aligned
+    tri_materials.push_back( tri.material );
+    tris_v.push_back( tri.v_list );
+    tris_n.push_back( tri.n_list );
+    tris_c.push_back( tri.c_list );
+    tris_tcs.push_back( tri.tc_list );
+    tris_vas.push_back( tri.va_list );
+    
+    return true;
+}
\ No newline at end of file