Go to the documentation of this file.00001
00002 #ifndef vgl_point_2d_h
00003 #define vgl_point_2d_h
00004 #ifdef VCL_NEEDS_PRAGMA_INTERFACE
00005 #pragma interface
00006 #endif
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019 #include <vcl_iosfwd.h>
00020 #include <vgl/vgl_fwd.h>
00021 #include <vgl/vgl_vector_2d.h>
00022 #include <vcl_cassert.h>
00023 #include <vcl_vector.h>
00024
00025
00026 template <class Type>
00027 class vgl_point_2d
00028 {
00029
00030 Type x_;
00031 Type y_;
00032
00033 public:
00034
00035
00036
00037
00038 inline vgl_point_2d () {}
00039
00040
00041 inline vgl_point_2d (Type px, Type py) : x_(px), y_(py) {}
00042
00043
00044 inline vgl_point_2d (Type const v[2]) : x_(v[0]), y_(v[1]) {}
00045
00046
00047 vgl_point_2d (vgl_homg_point_2d<Type> const& p);
00048
00049
00050 vgl_point_2d (vgl_line_2d<Type> const& l1,
00051 vgl_line_2d<Type> const& l2);
00052
00053 #if 0 // The compiler defaults for these are doing what they should do:
00054
00055 inline vgl_point_2d(vgl_point_2d<Type> const& p) : x_(p.x()), y_(p.y()) {}
00056
00057 inline ~vgl_point_2d () {}
00058
00059 inline vgl_point_2d<Type>& operator=(const vgl_point_2d<Type>& p)
00060 { x_ = p.x(); y_ = p.y(); return *this; }
00061 #endif
00062
00063
00064 inline bool operator==(vgl_point_2d<Type>const& p) const
00065 { return this==&p || (x_==p.x() && y_==p.y()); }
00066 inline bool operator!=(vgl_point_2d<Type>const& p)const
00067 { return !operator==(p); }
00068
00069
00070
00071 inline Type &x() {return x_;}
00072 inline Type &y() {return y_;}
00073
00074 inline Type x() const {return x_;}
00075 inline Type y() const {return y_;}
00076
00077
00078
00079 inline void set (Type px, Type py){ x_ = px; y_ = py; }
00080
00081
00082 inline void set (Type const p[2]) { x_ = p[0]; y_ = p[1]; }
00083
00084
00085
00086 inline bool ideal(Type = (Type)0) const { return false; }
00087
00088
00089
00090
00091
00092 vcl_istream& read(vcl_istream& is);
00093 };
00094
00095
00096
00097
00098
00099 template <class Type>
00100 vcl_ostream& operator<<(vcl_ostream& s, vgl_point_2d<Type> const& p);
00101
00102
00103
00104
00105
00106
00107 template <class Type>
00108 vcl_istream& operator>>(vcl_istream& s, vgl_point_2d<Type>& p);
00109
00110
00111
00112
00113
00114 template <class Type> inline
00115 bool is_ideal(vgl_point_2d<Type> const&, Type = 0) { return false; }
00116
00117
00118
00119 template <class Type> inline
00120 vgl_vector_2d<Type> operator-(vgl_point_2d<Type> const& p1,
00121 vgl_point_2d<Type> const& p2)
00122 { return vgl_vector_2d<Type>(p1.x()-p2.x(), p1.y()-p2.y()); }
00123
00124
00125
00126
00127 template <class Type> inline
00128 vgl_point_2d<Type> operator+(vgl_point_2d<Type> const& p,
00129 vgl_vector_2d<Type> const& v)
00130 { return vgl_point_2d<Type>(p.x()+v.x(), p.y()+v.y()); }
00131
00132
00133
00134 template <class Type> inline
00135 vgl_point_2d<Type>& operator+=(vgl_point_2d<Type>& p,
00136 vgl_vector_2d<Type> const& v)
00137 { p.set(p.x()+v.x(), p.y()+v.y()); return p; }
00138
00139
00140
00141 template <class Type> inline
00142 vgl_point_2d<Type> operator-(vgl_point_2d<Type> const& p,
00143 vgl_vector_2d<Type> const& v)
00144 { return p + (-v); }
00145
00146
00147
00148 template <class Type> inline
00149 vgl_point_2d<Type>& operator-=(vgl_point_2d<Type>& p,
00150 vgl_vector_2d<Type> const& v)
00151 { return p += (-v); }
00152
00153
00154
00155
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173 template <class T>
00174 double cross_ratio(vgl_point_2d<T>const& p1, vgl_point_2d<T>const& p2,
00175 vgl_point_2d<T>const& p3, vgl_point_2d<T>const& p4);
00176
00177
00178
00179 template <class Type> inline
00180 bool collinear(vgl_point_2d<Type> const& p1,
00181 vgl_point_2d<Type> const& p2,
00182 vgl_point_2d<Type> const& p3)
00183 { return parallel(p1-p2, p1-p3); }
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193 template <class Type> inline
00194 double ratio(vgl_point_2d<Type> const& p1,
00195 vgl_point_2d<Type> const& p2,
00196 vgl_point_2d<Type> const& p3)
00197 { return (p3-p1)/(p2-p1); }
00198
00199
00200
00201
00202
00203
00204
00205 template <class Type> inline
00206 vgl_point_2d<Type> midpoint(vgl_point_2d<Type> const& p1,
00207 vgl_point_2d<Type> const& p2,
00208 Type f = (Type)0.5)
00209 {
00210 return vgl_point_2d<Type>((Type)((1-f)*p1.x() + f*p2.x()),
00211 (Type)((1-f)*p1.y() + f*p2.y()));
00212 }
00213
00214
00215
00216
00217
00218 template <class Type> inline
00219 vgl_point_2d<Type> centre(vgl_point_2d<Type> const& p1,
00220 vgl_point_2d<Type> const& p2)
00221 {
00222 return vgl_point_2d<Type>((p1.x() + p2.x())/2 ,
00223 (p1.y() + p2.y())/2 );
00224 }
00225
00226
00227
00228 template <class Type> inline
00229 vgl_point_2d<Type> centre(vgl_point_2d<Type> const& p1,
00230 vgl_point_2d<Type> const& p2,
00231 vgl_point_2d<Type> const& p3)
00232 {
00233 return vgl_point_2d<Type>((p1.x() + p2.x() + p3.x())/3 ,
00234 (p1.y() + p2.y() + p3.y())/3 );
00235 }
00236
00237
00238
00239 template <class Type> inline
00240 vgl_point_2d<Type> centre(vgl_point_2d<Type> const& p1,
00241 vgl_point_2d<Type> const& p2,
00242 vgl_point_2d<Type> const& p3,
00243 vgl_point_2d<Type> const& p4)
00244 {
00245 return vgl_point_2d<Type>((p1.x() + p2.x() + p3.x() + p4.x())/4 ,
00246 (p1.y() + p2.y() + p3.y() + p4.y())/4 );
00247 }
00248
00249
00250
00251
00252 template <class Type> inline
00253 vgl_point_2d<Type> centre(vcl_vector<vgl_point_2d<Type> > const& v)
00254 {
00255 int n=v.size();
00256 assert(n>0);
00257 Type x = 0, y = 0;
00258 for (int i=0; i<n; ++i) x+=v[i].x(), y+=v[i].y();
00259 return vgl_point_2d<Type>(x/n,y/n);
00260 }
00261
00262 #define VGL_POINT_2D_INSTANTIATE(T) extern "please include vgl/vgl_point_2d.txx first"
00263
00264 #endif // vgl_point_2d_h