OpenVDB  8.1.1
Transform.h
Go to the documentation of this file.
1 // Copyright Contributors to the OpenVDB Project
2 // SPDX-License-Identifier: MPL-2.0
3 
4 #ifndef OPENVDB_MATH_TRANSFORM_HAS_BEEN_INCLUDED
5 #define OPENVDB_MATH_TRANSFORM_HAS_BEEN_INCLUDED
6 
7 #include "Maps.h"
8 #include <openvdb/Types.h>
9 #include <iosfwd>
10 
11 namespace openvdb {
13 namespace OPENVDB_VERSION_NAME {
14 namespace math {
15 
16 // Forward declaration
17 class Transform;
18 
19 
20 // Utility methods
21 
25 OPENVDB_API void
26 calculateBounds(const Transform& t, const Vec3d& minWS, const Vec3d& maxWS,
27  Vec3d& minIS, Vec3d& maxIS);
28 
31 //void calculateBounds(const Transform& t, const Vec3d& center, const Real radius,
32 // Vec3d& minIS, Vec3d& maxIS);
33 
34 
36 
37 
40 {
41 public:
44 
45  Transform(): mMap(MapBase::Ptr(new ScaleMap())) {}
46  Transform(const MapBase::Ptr&);
47  Transform(const Transform&);
49 
50  Ptr copy() const { return Ptr(new Transform(mMap->copy())); }
51 
53  static Transform::Ptr createLinearTransform(double voxelSize = 1.0);
55  static Transform::Ptr createLinearTransform(const Mat4R&);
56  static Transform::Ptr createFrustumTransform(const BBoxd&, double taper,
57  double depth, double voxelSize = 1.0);
59 
61  bool isLinear() const { return mMap->isLinear(); }
62 
64  bool isIdentity() const ;
66  Name mapType() const { return mMap->type(); }
67 
68 
70  void preRotate(double radians, const Axis axis = X_AXIS);
76  void preTranslate(const Vec3d&);
77  void preScale(const Vec3d&);
78  void preScale(double);
79  void preShear(double shear, Axis axis0, Axis axis1);
80  void preMult(const Mat4d&);
81  void preMult(const Mat3d&);
82 
83  void postRotate(double radians, const Axis axis = X_AXIS);
84  void postTranslate(const Vec3d&);
85  void postScale(const Vec3d&);
86  void postScale(double);
87  void postShear(double shear, Axis axis0, Axis axis1);
88  void postMult(const Mat4d&);
89  void postMult(const Mat3d&);
91 
93  Vec3d voxelSize() const { return mMap->voxelSize(); }
97  Vec3d voxelSize(const Vec3d& xyz) const { return mMap->voxelSize(xyz); }
98 
100  double voxelVolume() const { return mMap->determinant(); }
102  double voxelVolume(const Vec3d& xyz) const { return mMap->determinant(xyz); }
104  bool hasUniformScale() const { return mMap->hasUniformScale(); }
105 
107  Vec3d indexToWorld(const Vec3d& xyz) const { return mMap->applyMap(xyz); }
109  Vec3d indexToWorld(const Coord& ijk) const { return mMap->applyMap(ijk.asVec3d()); }
110  Vec3d worldToIndex(const Vec3d& xyz) const { return mMap->applyInverseMap(xyz); }
111  Coord worldToIndexCellCentered(const Vec3d& xyz) const {return Coord::round(worldToIndex(xyz));}
112  Coord worldToIndexNodeCentered(const Vec3d& xyz) const {return Coord::floor(worldToIndex(xyz));}
114 
116  BBoxd indexToWorld(const CoordBBox&) const;
119  BBoxd indexToWorld(const BBoxd&) const;
121 
122  BBoxd worldToIndex(const BBoxd&) const;
125  CoordBBox worldToIndexCellCentered(const BBoxd&) const;
126  CoordBBox worldToIndexNodeCentered(const BBoxd&) const;
128 
130  MapBase::ConstPtr baseMap() const { return mMap; }
132  MapBase::Ptr baseMap() { return mMap; }
134 
136  template<typename MapType> typename MapType::Ptr map();
139  template<typename MapType> typename MapType::ConstPtr map() const;
140  template<typename MapType> typename MapType::ConstPtr constMap() const;
142 
144  void read(std::istream&);
146  void write(std::ostream&) const;
147 
151  void print(std::ostream& os = std::cout, const std::string& indent = "") const;
152 
153  bool operator==(const Transform& other) const;
154  inline bool operator!=(const Transform& other) const { return !(*this == other); }
155 
156 private:
157  MapBase::Ptr mMap;
158 }; // class Transform
159 
160 
161 OPENVDB_API std::ostream& operator<<(std::ostream&, const Transform&);
162 
163 
165 
166 
167 template<typename MapType>
168 inline typename MapType::Ptr
169 Transform::map()
170 {
171  if (mMap->type() == MapType::mapType()) {
172  return StaticPtrCast<MapType>(mMap);
173  }
174  return typename MapType::Ptr();
175 }
176 
177 
178 template<typename MapType>
179 inline typename MapType::ConstPtr
180 Transform::map() const
181 {
182  return ConstPtrCast<const MapType>(
183  const_cast<Transform*>(this)->map<MapType>());
184 }
185 
186 
187 template<typename MapType>
188 inline typename MapType::ConstPtr
189 Transform::constMap() const
190 {
191  return map<MapType>();
192 }
193 
194 
196 
197 
199 template<typename ResolvedMapType, typename OpType>
200 inline void
201 doProcessTypedMap(Transform& transform, OpType& op)
202 {
203  ResolvedMapType& resolvedMap = *transform.map<ResolvedMapType>();
204  op.template operator()<ResolvedMapType>(resolvedMap);
205 }
206 
208 template<typename ResolvedMapType, typename OpType>
209 inline void
210 doProcessTypedMap(const Transform& transform, OpType& op)
211 {
212  const ResolvedMapType& resolvedMap = *transform.map<ResolvedMapType>();
213  op.template operator()<ResolvedMapType>(resolvedMap);
214 }
215 
216 
231 template<typename TransformType, typename OpType>
232 bool
233 processTypedMap(TransformType& transform, OpType& op)
234 {
235  using namespace openvdb;
236 
237  const Name mapType = transform.mapType();
238  if (mapType == UniformScaleMap::mapType()) {
239  doProcessTypedMap<UniformScaleMap, OpType>(transform, op);
240 
241  } else if (mapType == UniformScaleTranslateMap::mapType()) {
242  doProcessTypedMap<UniformScaleTranslateMap, OpType>(transform, op);
243 
244  } else if (mapType == ScaleMap::mapType()) {
245  doProcessTypedMap<ScaleMap, OpType>(transform, op);
246 
247  } else if (mapType == ScaleTranslateMap::mapType()) {
248  doProcessTypedMap<ScaleTranslateMap, OpType>(transform, op);
249 
250  } else if (mapType == UnitaryMap::mapType()) {
251  doProcessTypedMap<UnitaryMap, OpType>(transform, op);
252 
253  } else if (mapType == AffineMap::mapType()) {
254  doProcessTypedMap<AffineMap, OpType>(transform, op);
255 
256  } else if (mapType == TranslationMap::mapType()) {
257  doProcessTypedMap<TranslationMap, OpType>(transform, op);
258 
259  } else if (mapType == NonlinearFrustumMap::mapType()) {
260  doProcessTypedMap<NonlinearFrustumMap, OpType>(transform, op);
261  } else {
262  return false;
263  }
264  return true;
265 }
266 
267 } // namespace math
268 } // namespace OPENVDB_VERSION_NAME
269 } // namespace openvdb
270 
271 #endif // OPENVDB_MATH_TRANSFORM_HAS_BEEN_INCLUDED
Coord worldToIndexCellCentered(const Vec3d &xyz) const
Apply this transformation to the given coordinates.
Definition: Transform.h:111
Axis-aligned bounding box of signed integer coordinates.
Definition: Coord.h:248
Vec3d voxelSize(const Vec3d &xyz) const
Return the size of a voxel at position (x, y, z).
Definition: Transform.h:97
#define OPENVDB_API
Helper macros for defining library symbol visibility.
Definition: Platform.h:240
bool isLinear() const
Return true if the transformation map is exclusively linear/affine.
Definition: Transform.h:61
Transform()
Definition: Transform.h:45
SharedPtr< Transform > Ptr
Definition: Transform.h:42
double voxelVolume() const
Return the voxel volume of the linear component of the map.
Definition: Transform.h:100
Axis
Definition: Math.h:904
OPENVDB_API void calculateBounds(const Transform &t, const Vec3d &minWS, const Vec3d &maxWS, Vec3d &minIS, Vec3d &maxIS)
Calculate an axis-aligned bounding box in index space from an axis-aligned bounding box in world spac...
std::string Name
Definition: Name.h:17
MapBase::Ptr baseMap()
Return a base pointer to the transformation map.
Definition: Transform.h:132
SharedPtr< const Transform > ConstPtr
Definition: Transform.h:43
double voxelVolume(const Vec3d &xyz) const
Return the voxel volume at position (x, y, z).
Definition: Transform.h:102
void doProcessTypedMap(const Transform &transform, OpType &op)
Helper function used internally by processTypedMap()
Definition: Transform.h:210
Name mapType() const
Return the transformation map&#39;s type-name.
Definition: Transform.h:66
Vec3d voxelSize() const
Return the size of a voxel using the linear component of the map.
Definition: Transform.h:93
Vec3d asVec3d() const
Definition: Coord.h:144
bool isIdentity(const MatType &m)
Determine if a matrix is an identity matrix.
Definition: Mat.h:882
std::shared_ptr< T > SharedPtr
Definition: openvdb/Types.h:110
MapType::Ptr map()
Return the result of downcasting the base map pointer to a MapType pointer, or return a null pointer ...
Definition: Transform.h:169
bool hasUniformScale() const
Return true if the voxels in world space are uniformly sized cubes.
Definition: Transform.h:104
Vec3< double > Vec3d
Definition: Vec3.h:668
SharedPtr< MapBase > Ptr
Definition: Maps.h:137
bool operator!=(const Transform &other) const
Definition: Transform.h:154
bool processTypedMap(TransformType &transform, OpType &op)
Utility function that, given a generic map pointer, calls a functor on the fully-resoved map...
Definition: Transform.h:233
Definition: openvdb/Exceptions.h:13
MatType shear(Axis axis0, Axis axis1, typename MatType::value_type shear)
Set the matrix to a shear along axis0 by a fraction of axis1.
Definition: Mat.h:710
Vec3d worldToIndex(const Vec3d &xyz) const
Apply this transformation to the given coordinates.
Definition: Transform.h:110
Coord worldToIndexNodeCentered(const Vec3d &xyz) const
Apply this transformation to the given coordinates.
Definition: Transform.h:112
Definition: Math.h:905
Abstract base class for maps.
Definition: Maps.h:134
Vec3d indexToWorld(const Coord &ijk) const
Apply this transformation to the given coordinates.
Definition: Transform.h:109
bool operator==(const Vec3< T0 > &v0, const Vec3< T1 > &v1)
Equality operator, does exact floating point comparisons.
Definition: Vec3.h:477
Signed (x, y, z) 32-bit integer coordinates.
Definition: Coord.h:25
Ptr copy() const
Definition: Transform.h:50
OPENVDB_API std::ostream & operator<<(std::ostream &, const Transform &)
~Transform()
Definition: Transform.h:48
SharedPtr< const MapBase > ConstPtr
Definition: Maps.h:138
#define OPENVDB_VERSION_NAME
The version namespace name for this library version.
Definition: version.h.in:116
void print(const ast::Node &node, const bool numberStatements=true, std::ostream &os=std::cout, const char *indent=" ")
Writes a descriptive printout of a Node hierarchy into a target stream.
A specialized Affine transform that scales along the principal axis the scaling need not be uniform i...
Definition: Maps.h:669
#define OPENVDB_USE_VERSION_NAMESPACE
Definition: version.h.in:178
Definition: Transform.h:39