wzvoxmap.cxx

Delaunay三角形的网格剖分程序

//#include <iostream.h>
//#include "wzoutput.hxx"
//#include <stdio.h>
//#include <stdlib.h>
#include "wzvoxmap.hxx"       

wzGrayScaleTable::wzGrayScaleTable()
  :ranges(3*sizeof(wzGrayValue))
  ,min(ranges.base)
  ,max(ranges.base,sizeof(wzGrayValue))
  ,value(ranges.base,2*sizeof(wzGrayValue))
  ,mvalue(0)
{;}

wzIndex wzGrayScaleTable::addRange(wzByte V, wzGrayValue Min, wzGrayValue Max)
{
  wzIndex range = ranges.create();
  min[range] = Min;
  max[range] = Max;
  value[range] = V;
  if(V > mvalue) mvalue = V;
  return range;
}

wzByte wzGrayScaleTable::operator()(wzGrayValue v) const
{
  wzIndex range;
  wzRangeLoop(ranges, range){
    if(v<min[range]) continue;
    if(v>max[range]) continue;
    return value[range];
  }
  return 0;
}

wzGrayScaleVoxmap::wzGrayScaleVoxmap(int lx, int ly, int lz)
:Lx(lx),Ly(ly),Lz(lz)
{g(lx*ly*lz);}

wzByteVoxmap::wzByteVoxmap(int lx, int ly, int lz)
:Lx(lx),Ly(ly),Lz(lz)
{map(lx*ly*lz);}

wzByteVoxmap::wzByteVoxmap(const wzGrayScaleVoxmap& tmap,
 		 const wzGrayScaleTable& table)
:Lx(tmap.lx()),Ly(tmap.ly()),Lz(tmap.lz())
{
  wzIndex ix,iy,iz;
  map(Lx*Ly*Lz);
  for(ix=0;ix<Lx;ix++)
    for(iy=0;iy<Ly;iy++)
      for(iz=0;iz<Lz;iz++){
	//	wzGrayValue v = tmap(ix,iy,iz);
	//	  wzByte b = (wzByte) table((wzGrayValue) v);
	operator()(ix,iy,iz) = table(tmap(ix,iy,iz));
      }
}