I am having diffulty designing the logic for partitioning an octree for viewing frustum culling and could use some brainstorming ideas …

A little background on my project …

I have several massive databases of terrain contour data. The databases represent elevation data using 2 paradigms: spot elevations, and line segments. The quantity of spot elevations can reach into the hundreds of thousands and the quantity of line segments can reach into the millions.

I defined a set of classes to define the geometry like in the following pseudocode:

// ----------------------------------------------------------------------------

// Define a vector/vertexstruct MgVector3D

{

float x,y,z;

};// ----------------------------------------------------------------------------

// Define a simple bounding boxclass MgBoundingBox

{

MgVector3D minCoord;

MgVector3D maxCoord;

};// ----------------------------------------------------------------------------

// Define a list (a simple array) of vertices

// This is just an encapsulation of MgVector3D* pArrayclass VertexList : public TMgArray

{

public:

const MgBoundingBox& boundingBox() const;

…

}// ----------------------------------------------------------------------------

// Define a simple spot elevation (actually the indices of the spot’s

// vertex in a master vertex listtypedef size_t SpotElevation;

// ----------------------------------------------------------------------------

// Define a list (a simple array) of spot elevations

// This is just an encapsulation of SpotElevation* pArrayclass SpotElevationList : public TMgArray

{

…

};// ----------------------------------------------------------------------------

// Define a simple line segmentstruct LineSegment

{

/// Index of starting vertex in the external vertex array.

unsigned int startIndex;

/// Index of ending vertex in the external vertex array.

unsigned int endIndex;

/// Precomputed length of line

T length;

/// Construct given the indices and line length

LineSegment(unsigned int s=0, unsigned int e=0, T len=0)

: startIndex(s), endIndex(e), length(len)

{

}

};// ----------------------------------------------------------------------------

// Define a list (a simple array) of line segments

// This is just an encapsulation of LineSegment* pArrayclass LineSegmentList : public TMgArray

{

public:

const MgBoundingBox& boundingBox() const;

…

};// ----------------------------------------------------------------------------

// Define the contour model that holds everythingclass ContourModel

{

const MgBoundingBox& boundingBox() const { return m_Vertices.boundingBox(); }

private:

/// Master pool of vertices.

VertexList m_Vertices;

/// Collection of line segments.

LineSegmentList m_LineSegments;

/// Collection of spot elevations.

SpotElevationList m_SpotElevations;

…

};

I am experimenting with an Octree for partitioning the geometry and have run across a problem that I do not quite understand how to solve …

When building the Octree, I can easily place the spot elevations based on their position in space. However, deciding which octants to include the line segments in is a lot more difficult. The lines range in length from 0.1 meters to 980 meters (databases range up to 6000x4000 meters).

I thought about doing a line/box intersection test with the line segments and each octant but the lines could easily span many neighboring octants.

I thought about computing the midpoint of each line segment and using that as the basis for assigning them to octants however, the longer segments would frequently be culled even when they should be displayed.

I am new to Octrees and could use any help … Any ideas?

[This message has been edited by pleopard (edited 01-24-2002).]