ENH: add dataNode methods using directly the OGR api -> need to use the OGRAdapters (wip)

Code/Common/otbDataNode.h

@@ -266,6 +266,17 @@ public:
   double EuclideanDistance(const DataNode* node);
 
   double EuclideanDistance(const PointType point);
+
+   /**
+    * \return true if the two datanodes intersects
+    */
+  bool Intersects(const DataNode* node);
+
+   /**
+    * \return true if the two datanodes intersects
+    */
+  bool Within(const DataNode* node);
+
 /**
  * Clear all fields.
  */

Code/Common/otbDataNode.txx

@@ -636,6 +636,32 @@ DataNode<TPrecision, VDimension, TValuePrecision>
   return currentGeometry->Distance(&ogrPointSrc);
 }
 
+template <class TPrecision, unsigned int VDimension, class TValuePrecision>
+bool
+DataNode<TPrecision, VDimension, TValuePrecision>
+::Intersects(const DataNode* node)
+{
+  // Convert the nodes to OGRGeometries
+  OGRGeometry * dstGeomtery     = this->ConvertDataNodeToOGRGeometry(node);
+  OGRGeometry * currentGeometry = this->ConvertDataNodeToOGRGeometry(this);
+  
+  // 
+  return currentGeometry->Intersects(dstGeomtery);
+}
+
+template <class TPrecision, unsigned int VDimension, class TValuePrecision>
+bool
+DataNode<TPrecision, VDimension, TValuePrecision>
+::Within(const DataNode* node)
+{
+  // Convert the nodes to OGRGeometries
+  OGRGeometry * dstGeomtery     = this->ConvertDataNodeToOGRGeometry(node);
+  OGRGeometry * currentGeometry = this->ConvertDataNodeToOGRGeometry(this);
+  
+  // 
+  return currentGeometry->Within(dstGeomtery);
+}
+
 
 } // end namespace otb