Isocontour.h Source File

00001 /***************************************************************************
00002  *cr
00003  *cr            (C) Copyright 1995-2019 The Board of Trustees of the
00004  *cr                        University of Illinois
00005  *cr                         All Rights Reserved
00006  *cr
00007  ***************************************************************************/
00008 
00009 #ifndef ISOCONTOUR_H
00010 #define ISOCONTOUR_H
00011 
00012 #include "BaseMolecule.h"
00013 #include "ResizeArray.h"
00014 
00016 typedef struct {
00017    float x;
00018    float y;
00019 } XY;
00020 
00022 typedef struct {
00023    XY p[2];     
00024 } LINE;
00025 
00027 typedef struct {
00028    XY p[4];     
00029    float val[4]; 
00030    int cubeindex;
00031 } SQUARECELL;
00032 
00034 class IsoContour {
00035  public:
00036    int numtriangles;              
00037    ResizeArray<float>  v;         
00038    ResizeArray<int>    f;         
00039    const VolumetricData * vol;    
00040    float xax[3], yax[3], zax[3];  
00041 
00042  public:
00043    IsoContour(); 
00044 
00046    int compute(const VolumetricData *, float isovalue, int step); 
00047 
00048    void clear(); 
00049 
00050  /* private: */
00051  public:
00052    int DoCell(int x, int y, int z, float isovalue, int step);
00053    int Polygonise(const SQUARECELL grid, const float isolevel, LINE *lines);
00054    void VertexInterp(float isolevel, const SQUARECELL grid, int ind1, int ind2, XY *vert); 
00055 };
00056 
00057 
00058 #if defined(ISOCONTOUR_INTERNAL)
00059 
00060 static int edgeTable[16]={
00061 0x0  , 0x109, 0x203, 0x30a, 0x406, 0x50f, 0x605, 0x70c,
00062 0x70c, 0x605, 0x50f, 0x406, 0x30a, 0x203, 0x109, 0x0   };
00063 
00064 int lineTable[16][16] =
00065 {{-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00066  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00067  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00068  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00069  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00070  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00071  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00072  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00073  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00074  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00075  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00076  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00077  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00078  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00079  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1},
00080  {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1}};
00081 
00082 #endif
00083 
00084 #endif