00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024 #if defined(_MSC_VER)
00025 #pragma once
00026 #endif
00027
00028 #ifndef PBRT_VOLUMES_VOLUMEGRID_H
00029 #define PBRT_VOLUMES_VOLUMEGRID_H
00030
00031
00032 #include "volume.h"
00033
00034
00035 class VolumeGridDensity : public DensityRegion {
00036 public:
00037
00038 VolumeGridDensity(const Spectrum &sa, const Spectrum &ss, float gg,
00039 const Spectrum &emit, const BBox &e, const Transform &v2w,
00040 int x, int y, int z, const float *d)
00041 : DensityRegion(sa, ss, gg, emit, v2w), nx(x), ny(y), nz(z), extent(e) {
00042 density = new float[nx*ny*nz];
00043 memcpy(density, d, nx*ny*nz*sizeof(float));
00044 }
00045 ~VolumeGridDensity() { delete[] density; }
00046 BBox WorldBound() const { return Inverse(WorldToVolume)(extent); }
00047 bool IntersectP(const Ray &r, float *t0, float *t1) const {
00048 Ray ray = WorldToVolume(r);
00049 return extent.IntersectP(ray, t0, t1);
00050 }
00051 float Density(const Point &Pobj) const;
00052 float D(int x, int y, int z) const {
00053 x = Clamp(x, 0, nx-1);
00054 y = Clamp(y, 0, ny-1);
00055 z = Clamp(z, 0, nz-1);
00056 return density[z*nx*ny + y*nx + x];
00057 }
00058 private:
00059
00060 float *density;
00061 const int nx, ny, nz;
00062 const BBox extent;
00063 };
00064
00065
00066 VolumeGridDensity *CreateGridVolumeRegion(const Transform &volume2world,
00067 const ParamSet ¶ms);
00068
00069 #endif // PBRT_VOLUMES_VOLUMEGRID_H