Changeset 2517 in Sophya for trunk/SophyaPI/PIext/pintuple.cc

Timestamp:

Mar 17, 2004, 7:04:57 PM (22 years ago)

Author:

cmv

Message:

GetDistanceToPoint cmv 17/03/04

File:

• trunk/SophyaPI/PIext/pintuple.cc (modified) (5 diffs)

trunk/SophyaPI/PIext/pintuple.cc

@@ -4 +4 @@
 
 #include <stdio.h>
+#include <stdlib.h>
+#include <iostream>
+#include <math.h>
 #include "pintuple.h"
 
@@ -55 +58 @@
   SelectLabel(NULL);
   SetName("NTupleDrw");
+  NptDraw = 0;
+
 }
 
@@ -189 +194 @@
 xp = yp = xl = yl = 0;
 PIGrCoord xpolyg[NMXMULTP_LOCAL], ypolyg[NMXMULTP_LOCAL];
-npolyg = 0;  
+npolyg = 0;
+NptDraw = 0;
 for (long i=0; i<(long)mNT->NbLines(); i++) {
   xl = xp;  yl = yp;
@@ -198 +204 @@
   if ( (xp < xmin) || (xp > xmax) || (yp < ymin) || (yp > ymax) )  continue;
 
+  NptDraw++;
 // Taille - couleur de marker en fonction du poids
   if (wK >= 0) wp = mNT->GetCell(i, wK);
@@ -391 +398 @@
 
 
-
+/* --Methode-- */
+double PINTuple::GetDistanceToPoint(double x, double y)
+{
+ if(!mNT) return 1.e+9;
+ if( xK<0 || yK<0 ) return 1.e+9;
+
+ const int nessai = 100;
+ long inc = (NptDraw>nessai) ? (long)(NptDraw/nessai)+1 : 1;
+
+ double dist = -1.e+18;
+ long n = 0;
+ for(long i=0; i<(long)mNT->NbLines(); i++) {
+   double xp=mNT->GetCell(i,xK);
+   if(xp<XMin() || xp>XMax()) continue;
+   double yp=mNT->GetCell(i,yK);
+   if(yp<YMin() || yp>YMax()) continue;
+   if(n%inc==0) {
+     xp = (xp-x)/(XMax()-XMin())/0.5;
+     yp = (yp-y)/(YMax()-YMin())/0.5;
+     xp = xp*xp+yp*yp;
+     if(dist<0. || xp<dist) dist = xp;
+   }
+   n++;
+ }
+ dist=sqrt(fabs(dist));
+ //cout<<"PINTuple: xlim="<<XMin()<<","<<XMax()<<" ylim="<<YMin()<<","<<YMax()
+ //    <<" NbLines="<<mNT->NbLines()<<" inc="<<inc<<endl;
+ //cout<<"....d="<<dist<<" x="<<x<<" y="<<y<<" NptDraw="<<NptDraw<<endl;
+
+ return dist;
+}