]> git.mxchange.org Git - simgear.git/blobdiff - simgear/io/sg_binobj.cxx
Remove deprecated vector classes - finally!
[simgear.git] / simgear / io / sg_binobj.cxx
index 0e7082d08ad585beb7ee82445b13c90883aabff3..faed1e68f067ce9bb8e6e5087c11548208c75e0b 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 <config.h>
+#  include <simgear_config.h>
 #endif
 
 #include <simgear/compiler.h>
+#include <simgear/debug/logstream.hxx>
 
 #include <stdio.h>
 #include <time.h>
 
 #include <vector>
-#include STL_STRING
+#include <string>
+#include <iostream>
 
 #include <simgear/bucket/newbucket.hxx>
+#include <simgear/misc/sg_path.hxx>
+#include <simgear/math/SGGeometry.hxx>
 
 #include "lowlevel.hxx"
 #include "sg_binobj.hxx"
 
 
-FG_USING_STD( string );
-FG_USING_STD( vector );
-FG_USING_STD( cout );
-FG_USING_STD( endl );
+using std::string;
+using std::vector;
+using std::cout;
+using std::endl;
 
-
-enum {
+enum sgObjectTypes {
     SG_BOUNDING_SPHERE = 0,
 
     SG_VERTEX_LIST = 1,
+    SG_COLOR_LIST = 4,
     SG_NORMAL_LIST = 2,
     SG_TEXCOORD_LIST = 3,
 
+    SG_POINTS = 9,
+
     SG_TRIANGLE_FACES = 10,
     SG_TRIANGLE_STRIPS = 11,
     SG_TRIANGLE_FANS = 12
-} tgObjectTypes;
+};
 
-enum {
-    SG_MATERIAL = 0
-} tgPropertyTypes;
+enum sgIndexTypes {
+    SG_IDX_VERTICES =  0x01,
+    SG_IDX_NORMALS =   0x02,
+    SG_IDX_COLORS =    0x04,
+    SG_IDX_TEXCOORDS = 0x08
+};
+
+enum sgPropertyTypes {
+    SG_MATERIAL = 0,
+    SG_INDEX_TYPES = 1
+};
 
 
 class sgSimpleBuffer {
@@ -78,7 +92,7 @@ public:
        while ( size < s ) {
            size *= 2;
        }
-       cout << "Creating a new buffer of size = " << size << endl;
+        SG_LOG(SG_EVENT, SG_DEBUG, "Creating a new buffer of size = " << size);
        ptr = new char[size];
     }
 
@@ -96,90 +110,178 @@ public:
            while ( size < s ) {
                size *= 2;
            }
-           cout << "resizing buffer to size = " << size << endl;
+            SG_LOG(SG_EVENT, SG_DEBUG, "resizing buffer to size = " << size);
            ptr = new char[size];
        }
     }
 };
 
+// read object properties
+static void read_object( gzFile fp,
+                        int obj_type,
+                        int nproperties,
+                        int nelements,
+                        group_list *vertices,
+                        group_list *normals,
+                        group_list *colors,
+                        group_list *texcoords,
+                        string_list *materials )
+{
+    unsigned int nbytes;
+    unsigned char idx_mask;
+    int idx_size;
+    bool do_vertices, do_normals, do_colors, do_texcoords;
+    int j, k, idx;
+    sgSimpleBuffer buf( 32768 );  // 32 Kb
+    char material[256];
 
-// 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;
-
-    if ( wgs84_nodes.size() ) {
-       min = max = wgs84_nodes[0];
+    // default values
+    if ( obj_type == SG_POINTS ) {
+       idx_size = 1;
+       idx_mask = SG_IDX_VERTICES;
+       do_vertices = true;
+       do_normals = false;
+       do_colors = false;
+       do_texcoords = false;
     } else {
-       min = max = Point3D( 0 );
+       idx_size = 2;
+       idx_mask = (char)(SG_IDX_VERTICES | SG_IDX_TEXCOORDS);
+       do_vertices = true;
+       do_normals = false;
+       do_colors = false;
+       do_texcoords = true;
     }
 
-    for ( int i = 0; i < (int)wgs84_nodes.size(); ++i ) {
-       p = wgs84_nodes[i];
-
-       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() ); }
-
-       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() ); }
+    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;
+       } else if ( prop_type == SG_INDEX_TYPES ) {
+           idx_mask = ptr[0];
+           // cout << "idx_mask = " << (int)idx_mask << endl;
+           idx_size = 0;
+           do_vertices = false;
+           do_normals = false;
+           do_colors = false;
+           do_texcoords = false;
+           if ( idx_mask & SG_IDX_VERTICES ) {
+               do_vertices = true;
+               ++idx_size;
+           }
+           if ( idx_mask & SG_IDX_NORMALS ) {
+               do_normals = true;
+               ++idx_size;
+           }
+           if ( idx_mask & SG_IDX_COLORS ) {
+               do_colors = true;
+               ++idx_size;
+           }
+           if ( idx_mask & SG_IDX_TEXCOORDS ) {
+               do_texcoords = true;
+               ++idx_size;
+           }
+       }
     }
 
-    return ( min + max ) / 2.0;
-}
-
-// 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;
-    
-    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 < 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 / (idx_size * sizeof(unsigned short));
+       unsigned short *sptr = (unsigned short *)ptr;
+       int_list vs; vs.clear();
+       int_list ns; ns.clear();
+       int_list cs; cs.clear();
+       int_list tcs; tcs.clear();
+       for ( k = 0; k < count; ++k ) {
+           if ( sgIsBigEndian() ) {
+               for ( idx = 0; idx < idx_size; ++idx ) {
+                   sgEndianSwap( (uint16_t *)&(sptr[idx]) );
+               }
+           }
+           idx = 0;
+           if ( do_vertices ) {
+               vs.push_back( sptr[idx++] );
+           }
+           if ( do_normals ) {
+               ns.push_back( sptr[idx++] );
+                   }
+           if ( do_colors ) {
+               cs.push_back( sptr[idx++] );
+           }
+           if ( do_texcoords ) {
+               tcs.push_back( sptr[idx++] );
+           }
+           // cout << sptr[0] << " ";
+           sptr += idx_size;
+       }
+       // cout << endl;
+       vertices->push_back( vs );
+       normals->push_back( ns );
+       colors->push_back( cs );
+       texcoords->push_back( tcs );
+       materials->push_back( material );
     }
-
-    return sqrt(radius_squared);
 }
 
 
 // read a binary file and populate the provided structures.
 bool SGBinObject::read_bin( const string& file ) {
-    Point3D p;
+    SGVec3d p;
     int i, j, k;
-    char material[256];
     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_tc.clear();
+    pt_materials.clear();
+
     tris_v.clear();
+    tris_n.clear();
+    tris_c.clear();
     tris_tc.clear();
     tri_materials.clear();
 
     strips_v.clear();
+    strips_n.clear();
+    strips_c.clear();
     strips_tc.clear();
     strip_materials.clear();
 
     fans_v.clear();
+    fans_n.clear();
+    fans_c.clear();
     fans_tc.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;
+            SG_LOG( SG_EVENT, SG_ALERT,
+               "ERROR: opening " << file << " or " << filegz << " for reading!");
+
            return false;
        }
     }
@@ -188,7 +290,6 @@ bool SGBinObject::read_bin( const string& file ) {
 
     // read headers
     unsigned int header;
-    unsigned short version;
     sgReadUInt( fp, &header );
     if ( ((header & 0xFF000000) >> 24) == 'S' &&
         ((header & 0x00FF0000) >> 16) == 'G' ) {
@@ -197,31 +298,56 @@ bool SGBinObject::read_bin( const string& file ) {
        version = (header & 0x0000FFFF);
        // cout << "File version = " << version << endl;
     } else {
+       // close the file before we return
+       gzclose(fp);
+
        return false;
     }
 
     // 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
+    // disable it again once the debugging task is finished.
     struct tm *local_tm;
     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 );
+    unsigned short nobjects;
+    if ( version >= 7 ) {
+        sgReadUShort( fp, &nobjects );
+    } else {
+        short tmp;
+        sgReadShort( fp, &tmp );
+        nobjects = tmp;
+    }
     // cout << "Total objects to read = " << nobjects << endl;
 
     // read in objects
     for ( i = 0; i < nobjects; ++i ) {
        // read object header
        char obj_type;
-       short nproperties, nelements;
+       unsigned short nproperties, nelements;
        sgReadChar( fp, &obj_type );
-       sgReadShort( fp, &nproperties );
-       sgReadShort( fp, &nelements );
+        if ( version >= 7 ) {
+            sgReadUShort( fp, &nproperties );
+            sgReadUShort( fp, &nelements );
+        } else {
+            short tmp;
+            sgReadShort( fp, &tmp );
+            nproperties = tmp;
+            sgReadShort( fp, &tmp );
+            nelements = tmp;
+        }
 
        // cout << "object " << i << " = " << (int)obj_type << " props = "
        //      << nproperties << " elements = " << nelements << endl;
@@ -248,11 +374,19 @@ bool SGBinObject::read_bin( const string& file ) {
                sgReadBytes( fp, nbytes, ptr );
 
                double *dptr = (double *)ptr;
-               gbs_center = Point3D( dptr[0], dptr[1], dptr[2] );
+               if ( sgIsBigEndian() ) {
+                   sgEndianSwap( (uint64_t *)&(dptr[0]) );
+                   sgEndianSwap( (uint64_t *)&(dptr[1]) );
+                   sgEndianSwap( (uint64_t *)&(dptr[2]) );
+               }
+               gbs_center = SGVec3d( dptr[0], dptr[1], dptr[2] );
                // cout << "Center = " << gbs_center << endl;
                ptr += sizeof(double) * 3;
                
                float *fptr = (float *)ptr;
+               if ( sgIsBigEndian() ) {
+                   sgEndianSwap( (uint32_t *)fptr );
+               }
                gbs_radius = fptr[0];
                // cout << "Bounding radius = " << gbs_radius << endl;
            }
@@ -278,15 +412,19 @@ bool SGBinObject::read_bin( const string& file ) {
                sgReadBytes( fp, nbytes, ptr );
                int count = nbytes / (sizeof(float) * 3);
                float *fptr = (float *)ptr;
+               wgs84_nodes.reserve( count );
                for ( k = 0; k < count; ++k ) {
-                   p = Point3D( fptr[0], fptr[1], fptr[2] );
-                   // cout << "node = " << p << endl;
-                   wgs84_nodes.push_back( p );
+                   if ( sgIsBigEndian() ) {
+                       sgEndianSwap( (uint32_t *)&(fptr[0]) );
+                       sgEndianSwap( (uint32_t *)&(fptr[1]) );
+                       sgEndianSwap( (uint32_t *)&(fptr[2]) );
+                   }
+                   wgs84_nodes.push_back( SGVec3d(fptr[0], fptr[1], fptr[2]) );
                    fptr += 3;
                }
            }
-       } else if ( obj_type == SG_NORMAL_LIST ) {
-           // read normal list properties
+       } else if ( obj_type == SG_COLOR_LIST ) {
+           // read color list properties
            for ( j = 0; j < nproperties; ++j ) {
                char prop_type;
                sgReadChar( fp, &prop_type );
@@ -298,54 +436,30 @@ bool SGBinObject::read_bin( const string& file ) {
                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 ) {
-                   p = Point3D( ptr[0] / 128.0 - 1.0,
-                                ptr[1] / 128.0 - 1.0,
-                                ptr[2] / 128.0 - 1.0 );
-                   // 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
+           // read color 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);
+               int count = nbytes / (sizeof(float) * 4);
                float *fptr = (float *)ptr;
+               colors.reserve(count);
                for ( k = 0; k < count; ++k ) {
-                   p = Point3D( fptr[0], fptr[1], 0 );
-                   // cout << "texcoord = " << p << endl;
-                   texcoords.push_back( p );
-                   fptr += 2;
+                   if ( sgIsBigEndian() ) {
+                       sgEndianSwap( (uint32_t *)&(fptr[0]) );
+                       sgEndianSwap( (uint32_t *)&(fptr[1]) );
+                       sgEndianSwap( (uint32_t *)&(fptr[2]) );
+                       sgEndianSwap( (uint32_t *)&(fptr[3]) );
+                   }
+                   SGVec4f color( fptr[0], fptr[1], fptr[2], fptr[3] );
+                   colors.push_back( color );
+                   fptr += 4;
                }
            }
-       } else if ( obj_type == SG_TRIANGLE_FACES ) {
-           // read triangle face properties
+       } else if ( obj_type == SG_NORMAL_LIST ) {
+           // read normal list properties
            for ( j = 0; j < nproperties; ++j ) {
                char prop_type;
                sgReadChar( fp, &prop_type );
@@ -355,37 +469,28 @@ bool SGBinObject::read_bin( const string& file ) {
                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
+           // 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 ); }
-               char *ptr = buf.get_ptr();
+               unsigned char *ptr = (unsigned char *)(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();
+               int count = nbytes / 3;
+               normals.reserve( count );
                for ( k = 0; k < count; ++k ) {
-                   vs.push_back( sptr[0] );
-                   tcs.push_back( sptr[1] );
-                   // cout << sptr[0] << "/" << sptr[1] << " ";
-                   sptr += 2;
+                    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;
                }
-               // 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
+       } else if ( obj_type == SG_TEXCOORD_LIST ) {
+           // read texcoord list properties
            for ( j = 0; j < nproperties; ++j ) {
                char prop_type;
                sgReadChar( fp, &prop_type );
@@ -395,75 +500,44 @@ bool SGBinObject::read_bin( const string& file ) {
                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
+           // 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(short) * 2);
-               short *sptr = (short *)ptr;
-               int_list vs, tcs;
-               vs.clear(); tcs.clear();
+               int count = nbytes / (sizeof(float) * 2);
+               float *fptr = (float *)ptr;
+               texcoords.reserve(count);
                for ( k = 0; k < count; ++k ) {
-                   vs.push_back( sptr[0] );
-                   tcs.push_back( sptr[1] );
-                   // cout << sptr[0] << "/" << sptr[1] << " ";
-                   sptr += 2;
+                   if ( sgIsBigEndian() ) {
+                       sgEndianSwap( (uint32_t *)&(fptr[0]) );
+                       sgEndianSwap( (uint32_t *)&(fptr[1]) );
+                   }
+                   texcoords.push_back( SGVec2f( fptr[0], fptr[1] ) );
+                   fptr += 2;
                }
-               // cout << endl;
-               strips_v.push_back( vs );
-               strips_tc.push_back( tcs );
-               strip_materials.push_back( material );
            }
+       } else if ( obj_type == SG_POINTS ) {
+           // read point elements
+           read_object( fp, SG_POINTS, nproperties, nelements,
+                        &pts_v, &pts_n, &pts_c, &pts_tc, &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_tc, &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_tc,
+                        &strip_materials );
        } 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 ) {
-                   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 );
-           }
+           read_object( fp, SG_TRIANGLE_FANS, nproperties, nelements,
+                        &fans_v, &fans_n, &fans_c, &fans_tc, &fan_materials );
        } else {
            // unknown object type, just skip
 
@@ -508,30 +582,24 @@ bool SGBinObject::read_bin( const string& file ) {
 bool SGBinObject::write_bin( const string& base, const string& name,
                             const SGBucket& b )
 {
-    Point3D p;
-    sgVec2 t;
-    sgVec3 pt;
     int i, j;
+    unsigned char idx_mask;
+    int idx_size;
 
-    string dir = base + "/" + b.gen_base_path();
-    string command = "mkdir -p " + dir;
-#ifdef _MSC_VER
-    system( (string("mkdir ") + dir).c_str() );
-#else
-    system(command.c_str());
-#endif
-
-    string file = dir + "/" + name + ".gz";
-    cout << "Output file = " << file << endl;
+    SGPath file = base + "/" + b.gen_base_path() + "/" + name + ".gz";
+    file.create_dir( 0755 );
+    cout << "Output file = " << file.str() << endl;
 
     gzFile fp;
     if ( (fp = gzopen( file.c_str(), "wb9" )) == NULL ) {
-       cout << "ERROR: opening " << file << " for writing!" << endl;
+       cout << "ERROR: opening " << file.str() << " for writing!" << endl;
        return false;
     }
 
     sgClearWriteError();
 
+    cout << "points size = " << pts_v.size() << "  pt_materials = " 
+        << pt_materials.size() << endl;
     cout << "triangles size = " << tris_v.size() << "  tri_materials = " 
         << tri_materials.size() << endl;
     cout << "strips size = " << strips_v.size() << "  strip_materials = " 
@@ -539,26 +607,43 @@ bool SGBinObject::write_bin( const string& base, const string& name,
     cout << "fans size = " << fans_v.size() << "  fan_materials = " 
         << fan_materials.size() << endl;
 
-    cout << "points = " << wgs84_nodes.size() << endl;
+    cout << "nodes = " << wgs84_nodes.size() << endl;
+    cout << "colors = " << colors.size() << endl;
+    cout << "normals = " << normals.size() << endl;
     cout << "tex coords = " << texcoords.size() << endl;
 
     // write header magic
     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;
+    unsigned short nobjects = 0;
     nobjects++;                        // for gbs
     nobjects++;                        // for vertices
+    nobjects++;                        // for colors
     nobjects++;                        // for normals
     nobjects++;                        // for texcoords
 
+    // points
+    unsigned short npts = 0;
+    start = 0; end = 1;
+    while ( start < (int)pt_materials.size() ) {
+       material = pt_materials[start];
+       while ( (end < (int)pt_materials.size()) &&
+               (material == pt_materials[end]) ) {
+           end++;
+       }
+       npts++;
+       start = end; end = start + 1;
+    }
+    nobjects += npts;
+
     // tris
-    short ntris = 0;
+    unsigned short ntris = 0;
     start = 0; end = 1;
     while ( start < (int)tri_materials.size() ) {
        material = tri_materials[start];
@@ -572,7 +657,7 @@ bool SGBinObject::write_bin( const string& base, const string& name,
     nobjects += ntris;
 
     // strips
-    short nstrips = 0;
+    unsigned short nstrips = 0;
     start = 0; end = 1;
     while ( start < (int)strip_materials.size() ) {
        material = strip_materials[start];
@@ -586,7 +671,7 @@ bool SGBinObject::write_bin( const string& base, const string& name,
     nobjects += nstrips;
 
     // fans
-    short nfans = 0;
+    unsigned short nfans = 0;
     start = 0; end = 1;
     while ( start < (int)fan_materials.size() ) {
        material = fan_materials[start];
@@ -600,12 +685,12 @@ bool SGBinObject::write_bin( const string& base, const string& name,
     nobjects += nfans;
 
     cout << "total top level objects = " << nobjects << endl;
-    sgWriteShort( fp, nobjects );
+    sgWriteUShort( fp, nobjects );
 
     // write bounding sphere
     sgWriteChar( fp, (char)SG_BOUNDING_SPHERE );        // type
-    sgWriteShort( fp, 0 );                             // nproperties
-    sgWriteShort( fp, 1 );                             // nelements
+    sgWriteUShort( fp, 0 );                            // nproperties
+    sgWriteUShort( fp, 1 );                            // nelements
 
     sgWriteUInt( fp, sizeof(double) * 3 + sizeof(float) ); // nbytes
     sgdVec3 center;
@@ -615,38 +700,105 @@ bool SGBinObject::write_bin( const string& base, const string& name,
 
     // dump vertex list
     sgWriteChar( fp, (char)SG_VERTEX_LIST );             // type
-    sgWriteShort( fp, 0 );                              // nproperties
-    sgWriteShort( fp, 1 );                              // nelements
+    sgWriteUShort( fp, 0 );                             // nproperties
+    sgWriteUShort( fp, 1 );                             // nelements
     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 );
+        SGVec3f p = toVec3f(wgs84_nodes[i] - gbs_center);
+       sgWriteVec3( fp, p.data() );
+    }
+
+    // dump vertex color list
+    sgWriteChar( fp, (char)SG_COLOR_LIST );             // type
+    sgWriteUShort( fp, 0 );                             // nproperties
+    sgWriteUShort( fp, 1 );                             // nelements
+    sgWriteUInt( fp, colors.size() * sizeof(float) * 4 ); // nbytes
+    for ( i = 0; i < (int)colors.size(); ++i ) {
+       sgWriteVec4( fp, colors[i].data() );
     }
 
     // dump vertex normal list
     sgWriteChar( fp, (char)SG_NORMAL_LIST );            // type
-    sgWriteShort( fp, 0 );                             // nproperties
-    sgWriteShort( fp, 1 );                             // nelements
+    sgWriteUShort( fp, 0 );                            // nproperties
+    sgWriteUShort( fp, 1 );                            // nelements
     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) * 128);
-       normal[1] = (unsigned char)((p.y() + 1.0) * 128);
-       normal[2] = (unsigned char)((p.z() + 1.0) * 128);
+        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
+    sgWriteUShort( fp, 0 );                            // nproperties
+    sgWriteUShort( fp, 1 );                            // nelements
     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 );
+       sgWriteVec2( fp, texcoords[i].data() );
+    }
+
+    // dump point groups if they exist
+    if ( pts_v.size() > 0 ) {
+       int start = 0;
+       int end = 1;
+       string material;
+       while ( start < (int)pt_materials.size() ) {
+           // find next group
+           material = pt_materials[start];
+           while ( (end < (int)pt_materials.size()) && 
+                   (material == pt_materials[end]) )
+               {
+                   // cout << "end = " << end << endl;
+                   end++;
+               }
+           // cout << "group = " << start << " to " << end - 1 << endl;
+
+           // write group headers
+           sgWriteChar( fp, (char)SG_POINTS );          // type
+           sgWriteUShort( fp, 2 );                      // nproperties
+           sgWriteUShort( fp, end - start );             // nelements
+
+           sgWriteChar( fp, (char)SG_MATERIAL );        // property
+           sgWriteUInt( fp, material.length() );        // nbytes
+           sgWriteBytes( fp, material.length(), material.c_str() );
+
+           idx_mask = 0;
+           idx_size = 0;
+           if ( pts_v.size() ) { idx_mask |= SG_IDX_VERTICES; ++idx_size; }
+           if ( pts_n.size() ) { idx_mask |= SG_IDX_NORMALS; ++idx_size; }
+           if ( pts_c.size() ) { idx_mask |= SG_IDX_COLORS; ++idx_size; }
+           if ( pts_tc.size() ) { idx_mask |= SG_IDX_TEXCOORDS; ++idx_size; }
+           sgWriteChar( fp, (char)SG_INDEX_TYPES );     // property
+           sgWriteUInt( fp, 1 );                        // nbytes
+           sgWriteChar( fp, idx_mask );
+
+           // write strips
+           for ( i = start; i < end; ++i ) {
+               // nbytes
+               sgWriteUInt( fp, pts_v[i].size() * idx_size
+                                 * sizeof(unsigned short) );
+               for ( j = 0; j < (int)pts_v[i].size(); ++j ) {
+                   if ( pts_v.size() ) { 
+                       sgWriteUShort( fp, (unsigned short)pts_v[i][j] );
+                   }
+                   if ( pts_n.size() ) { 
+                       sgWriteUShort( fp, (unsigned short)pts_n[i][j] );
+                   }
+                   if ( pts_c.size() ) { 
+                       sgWriteUShort( fp, (unsigned short)pts_c[i][j] );
+                   }
+                   if ( pts_tc.size() ) { 
+                       sgWriteUShort( fp, (unsigned short)pts_tc[i][j] );
+                   }
+               }
+           }
+           
+           start = end;
+           end = start + 1;
+       }
     }
 
     // dump individual triangles if they exist
@@ -658,7 +810,8 @@ bool SGBinObject::write_bin( const string& base, const string& name,
            // find next group
            material = tri_materials[start];
            while ( (end < (int)tri_materials.size()) && 
-                   (material == tri_materials[end]) )
+                   (material == tri_materials[end]) &&
+                   3*(end-start) < 32760 )
            {
                // cout << "end = " << end << endl;
                end++;
@@ -667,20 +820,42 @@ bool SGBinObject::write_bin( const string& base, const string& name,
 
            // write group headers
            sgWriteChar( fp, (char)SG_TRIANGLE_FACES ); // type
-           sgWriteShort( fp, 1 );                      // nproperties
-           sgWriteShort( fp, 1 );                      // nelements
+           sgWriteUShort( fp, 2 );                     // nproperties
+           sgWriteUShort( 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
+           idx_mask = 0;
+           idx_size = 0;
+           if ( tris_v.size() ) { idx_mask |= SG_IDX_VERTICES; ++idx_size; }
+           if ( tris_n.size() ) { idx_mask |= SG_IDX_NORMALS; ++idx_size; }
+           if ( tris_c.size() ) { idx_mask |= SG_IDX_COLORS; ++idx_size; }
+           if ( tris_tc.size() ) { idx_mask |= SG_IDX_TEXCOORDS; ++idx_size; }
+           sgWriteChar( fp, (char)SG_INDEX_TYPES );     // property
+           sgWriteUInt( fp, 1 );                        // nbytes
+           sgWriteChar( fp, idx_mask );
+
+           // nbytes
+           sgWriteUInt( fp, (end - start) * 3 * idx_size
+                             * sizeof(unsigned short) );
 
            // 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] );
+                   if ( tris_v.size() ) {
+                       sgWriteUShort( fp, (unsigned short)tris_v[i][j] );
+                   }
+                   if ( tris_n.size() ) {
+                       sgWriteUShort( fp, (unsigned short)tris_n[i][j] );
+                   }
+                   if ( tris_c.size() ) {
+                       sgWriteUShort( fp, (unsigned short)tris_c[i][j] );
+                   }
+                   if ( tris_tc.size() ) {
+                       sgWriteUShort( fp, (unsigned short)tris_tc[i][j] );
+                   }
                }
            }
 
@@ -707,20 +882,41 @@ bool SGBinObject::write_bin( const string& base, const string& name,
 
            // write group headers
            sgWriteChar( fp, (char)SG_TRIANGLE_STRIPS ); // type
-           sgWriteShort( fp, 1 );                       // nproperties
-           sgWriteShort( fp, end - start );             // nelements
+           sgWriteUShort( fp, 2 );                      // nproperties
+           sgWriteUShort( fp, end - start );             // nelements
 
            sgWriteChar( fp, (char)SG_MATERIAL );        // property
            sgWriteUInt( fp, material.length() );        // nbytes
            sgWriteBytes( fp, material.length(), material.c_str() );
 
+           idx_mask = 0;
+           idx_size = 0;
+           if ( strips_v.size() ) { idx_mask |= SG_IDX_VERTICES; ++idx_size; }
+           if ( strips_n.size() ) { idx_mask |= SG_IDX_NORMALS; ++idx_size; }
+           if ( strips_c.size() ) { idx_mask |= SG_IDX_COLORS; ++idx_size; }
+           if ( strips_tc.size() ) { idx_mask |= SG_IDX_TEXCOORDS; ++idx_size;}
+           sgWriteChar( fp, (char)SG_INDEX_TYPES );     // property
+           sgWriteUInt( fp, 1 );                        // nbytes
+           sgWriteChar( fp, idx_mask );
+
            // write strips
            for ( i = start; i < end; ++i ) {
                // nbytes
-               sgWriteUInt( fp, strips_v[i].size() * 2 * sizeof(short) );
+               sgWriteUInt( fp, strips_v[i].size() * idx_size
+                                 * sizeof(unsigned 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] );
+                   if ( strips_v.size() ) { 
+                       sgWriteUShort( fp, (unsigned short)strips_v[i][j] );
+                   }
+                   if ( strips_n.size() ) { 
+                       sgWriteUShort( fp, (unsigned short)strips_n[i][j] );
+                   }
+                   if ( strips_c.size() ) { 
+                       sgWriteUShort( fp, (unsigned short)strips_c[i][j] );
+                   }
+                   if ( strips_tc.size() ) { 
+                       sgWriteUShort( fp, (unsigned short)strips_tc[i][j] );
+                   }
                }
            }
            
@@ -747,20 +943,41 @@ bool SGBinObject::write_bin( const string& base, const string& name,
 
            // write group headers
            sgWriteChar( fp, (char)SG_TRIANGLE_FANS );   // type
-           sgWriteShort( fp, 1 );                       // nproperties
-           sgWriteShort( fp, end - start );             // nelements
+           sgWriteUShort( fp, 2 );                      // nproperties
+           sgWriteUShort( fp, end - start );             // nelements
 
            sgWriteChar( fp, (char)SG_MATERIAL );       // property
            sgWriteUInt( fp, material.length() );        // nbytes
            sgWriteBytes( fp, material.length(), material.c_str() );
 
+           idx_mask = 0;
+           idx_size = 0;
+           if ( fans_v.size() ) { idx_mask |= SG_IDX_VERTICES; ++idx_size; }
+           if ( fans_n.size() ) { idx_mask |= SG_IDX_NORMALS; ++idx_size; }
+           if ( fans_c.size() ) { idx_mask |= SG_IDX_COLORS; ++idx_size; }
+           if ( fans_tc.size() ) { idx_mask |= SG_IDX_TEXCOORDS; ++idx_size; }
+           sgWriteChar( fp, (char)SG_INDEX_TYPES );     // property
+           sgWriteUInt( fp, 1 );                        // nbytes
+           sgWriteChar( fp, idx_mask );
+
            // write fans
            for ( i = start; i < end; ++i ) {
                // nbytes
-               sgWriteUInt( fp, fans_v[i].size() * 2 * sizeof(short) );
+               sgWriteUInt( fp, fans_v[i].size() * idx_size
+                                 * sizeof(unsigned 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] );
+                   if ( fans_v.size() ) {
+                       sgWriteUShort( fp, (unsigned short)fans_v[i][j] );
+                   }
+                   if ( fans_n.size() ) {
+                       sgWriteUShort( fp, (unsigned short)fans_n[i][j] );
+                   }
+                   if ( fans_c.size() ) {
+                       sgWriteUShort( fp, (unsigned short)fans_c[i][j] );
+                   }
+                   if ( fans_tc.size() ) {
+                       sgWriteUShort( fp, (unsigned short)fans_tc[i][j] );
+                   }
                }
            }
            
@@ -787,24 +1004,15 @@ bool SGBinObject::write_bin( const string& base, const string& name,
 bool SGBinObject::write_ascii( const string& base, const string& name,
                               const SGBucket& b )
 {
-    Point3D p;
     int i, j;
 
-    string dir = base + "/" + b.gen_base_path();
-    string command = "mkdir -p " + dir;
-#ifdef _MSC_VER
-    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;
+       cout << "ERROR: opening " << file.str() << " for writing!" << endl;
        return false;
     }
 
@@ -837,7 +1045,7 @@ bool SGBinObject::write_ascii( const string& base, const string& name,
     // dump vertex list
     fprintf(fp, "# vertex list\n");
     for ( i = 0; i < (int)wgs84_nodes.size(); ++i ) {
-       p = wgs84_nodes[i] - gbs_center;
+        SGVec3d p = wgs84_nodes[i] - gbs_center;
        
        fprintf(fp,  "v %.5f %.5f %.5f\n", p.x(), p.y(), p.z() );
     }
@@ -845,7 +1053,7 @@ bool SGBinObject::write_ascii( const string& base, const string& name,
 
     fprintf(fp, "# vertex normal list\n");
     for ( i = 0; i < (int)normals.size(); ++i ) {
-       p = normals[i];
+        SGVec3f p = normals[i];
        fprintf(fp,  "vn %.5f %.5f %.5f\n", p.x(), p.y(), p.z() );
     }
     fprintf(fp, "\n");
@@ -853,7 +1061,7 @@ bool SGBinObject::write_ascii( const string& base, const string& name,
     // dump texture coordinates
     fprintf(fp, "# texture coordinate list\n");
     for ( i = 0; i < (int)texcoords.size(); ++i ) {
-       p = texcoords[i];
+        SGVec2f p = texcoords[i];
        fprintf(fp,  "vt %.5f %.5f\n", p.x(), p.y() );
     }
     fprintf(fp, "\n");
@@ -876,19 +1084,16 @@ bool SGBinObject::write_ascii( const string& base, const string& name,
            }
            // 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 );
-               }
-           }
-
+      SGSphered d;
+      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());
@@ -927,18 +1132,16 @@ bool SGBinObject::write_ascii( const string& base, const string& name,
                }
            // 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 );
-               }
-           }
+
+      SGSphered d;
+      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");
@@ -963,7 +1166,7 @@ bool SGBinObject::write_ascii( const string& base, const string& name,
     // close the file
     fclose(fp);
 
-    command = "gzip --force --best " + file;
+    string command = "gzip --force --best " + file.str();
     system(command.c_str());
 
     return true;