Salome HOME
rnc: Import Picture and Feature Detection dialogs small modifications
[modules/geom.git] / src / GEOM_I / GEOM_IMeasureOperations_i.cc
1 // Copyright (C) 2007-2011  CEA/DEN, EDF R&D, OPEN CASCADE
2 //
3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 //
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License.
10 //
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 // Lesser General Public License for more details.
15 //
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 //
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 //
22
23 #include <Standard_Stream.hxx>
24
25 #include "GEOM_IMeasureOperations_i.hh"
26
27 #include "utilities.h"
28 #include "OpUtil.hxx"
29
30 #include "GEOM_Engine.hxx"
31 #include "GEOM_Object.hxx"
32
33 //=============================================================================
34 /*!
35  *   constructor:
36  */
37 //=============================================================================
38 GEOM_IMeasureOperations_i::GEOM_IMeasureOperations_i (PortableServer::POA_ptr thePOA,
39                                                       GEOM::GEOM_Gen_ptr theEngine,
40                                                       ::GEOMImpl_IMeasureOperations* theImpl)
41 :GEOM_IOperations_i(thePOA, theEngine, theImpl)
42 {
43   MESSAGE("GEOM_IMeasureOperations_i::GEOM_IMeasureOperations_i");
44 }
45
46 //=============================================================================
47 /*!
48  *  destructor
49  */
50 //=============================================================================
51 GEOM_IMeasureOperations_i::~GEOM_IMeasureOperations_i()
52 {
53   MESSAGE("GEOM_IMeasureOperations_i::~GEOM_IMeasureOperations_i");
54 }
55
56 //=============================================================================
57 /*!
58  *  KindOfShape
59  */
60 //=============================================================================
61 GEOM::GEOM_IKindOfShape::shape_kind GEOM_IMeasureOperations_i::KindOfShape
62                                    (GEOM::GEOM_Object_ptr  theShape,
63                                     GEOM::ListOfLong_out   theIntegers,
64                                     GEOM::ListOfDouble_out theDoubles)
65 {
66   GEOMImpl_IMeasureOperations::ShapeKind aKind = GEOMImpl_IMeasureOperations::SK_NO_SHAPE;
67
68   // allocate the CORBA arrays
69   GEOM::ListOfLong_var anIntegersArray = new GEOM::ListOfLong();
70   GEOM::ListOfDouble_var aDoublesArray = new GEOM::ListOfDouble();
71
72   //Get the reference shape
73   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
74
75   if (!aShape.IsNull()) {
76     Handle(TColStd_HSequenceOfInteger) anIntegers = new TColStd_HSequenceOfInteger;
77     Handle(TColStd_HSequenceOfReal)    aDoubles   = new TColStd_HSequenceOfReal;
78
79     // Detect kind of shape and parameters
80     aKind = GetOperations()->KindOfShape(aShape, anIntegers, aDoubles);
81
82     int nbInts = anIntegers->Length();
83     int nbDbls = aDoubles->Length();
84
85     anIntegersArray->length(nbInts);
86     aDoublesArray->length(nbDbls);
87
88     for (int ii = 0; ii < nbInts; ii++) {
89       anIntegersArray[ii] = anIntegers->Value(ii + 1);
90     }
91     for (int id = 0; id < nbDbls; id++) {
92       aDoublesArray[id] = aDoubles->Value(id + 1);
93     }
94   }
95
96   // initialize out-parameters with local arrays
97   theIntegers = anIntegersArray._retn();
98   theDoubles  = aDoublesArray._retn();
99   return (GEOM::GEOM_IKindOfShape::shape_kind)aKind;
100 }
101
102 //=============================================================================
103 /*!
104  *  GetPosition
105  */
106 //=============================================================================
107 void GEOM_IMeasureOperations_i::GetPosition
108                  (GEOM::GEOM_Object_ptr theShape,
109                   CORBA::Double& Ox, CORBA::Double& Oy, CORBA::Double& Oz,
110                   CORBA::Double& Zx, CORBA::Double& Zy, CORBA::Double& Zz,
111                   CORBA::Double& Xx, CORBA::Double& Xy, CORBA::Double& Xz)
112 {
113   //Set a not done flag
114   GetOperations()->SetNotDone();
115
116   //Set default values: global CS
117   Ox = Oy = Oz = Zx = Zy = Xy = Xz = 0.;
118   Zz = Xx = 1.;
119
120   //Get the reference shape
121   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
122   if (aShape.IsNull()) return;
123
124   // Get shape parameters
125   GetOperations()->GetPosition(aShape, Ox,Oy,Oz, Zx,Zy,Zz, Xx,Xy,Xz);
126 }
127
128 //=============================================================================
129 /*!
130  *  GetCentreOfMass
131  */
132 //=============================================================================
133 GEOM::GEOM_Object_ptr GEOM_IMeasureOperations_i::GetCentreOfMass
134                                               (GEOM::GEOM_Object_ptr theShape)
135 {
136   GEOM::GEOM_Object_var aGEOMObject;
137
138   //Set a not done flag
139   GetOperations()->SetNotDone();
140
141   //Get the reference shape
142   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
143   if (aShape.IsNull()) return aGEOMObject._retn();
144
145   // Make Point - centre of mass of theShape
146   Handle(GEOM_Object) anObject = GetOperations()->GetCentreOfMass(aShape);
147   if (!GetOperations()->IsDone() || anObject.IsNull())
148     return aGEOMObject._retn();
149
150   return GetObject(anObject);
151 }
152
153 //=============================================================================
154 /*!
155  *  GetVertexByIndex
156  */
157 //=============================================================================
158 GEOM::GEOM_Object_ptr GEOM_IMeasureOperations_i::GetVertexByIndex
159   (GEOM::GEOM_Object_ptr theShape, CORBA::Long theIndex)
160 {
161   GEOM::GEOM_Object_var aGEOMObject;
162
163   //Set a not done flag
164   GetOperations()->SetNotDone();
165
166   //Get the reference shape
167   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
168   if ( aShape.IsNull() ) return aGEOMObject._retn();
169
170   // Get vertex by index
171   Handle(GEOM_Object) anObject = GetOperations()->GetVertexByIndex(aShape, theIndex);
172   if (!GetOperations()->IsDone() || anObject.IsNull())
173     return aGEOMObject._retn();
174
175   return GetObject(anObject);
176 }
177
178 //=============================================================================
179 /*!
180  *  GetNormal
181  */
182 //=============================================================================
183 GEOM::GEOM_Object_ptr GEOM_IMeasureOperations_i::GetNormal
184                                        (GEOM::GEOM_Object_ptr theFace,
185                                         GEOM::GEOM_Object_ptr theOptionalPoint)
186 {
187   GEOM::GEOM_Object_var aGEOMObject;
188
189   //Set a not done flag
190   GetOperations()->SetNotDone();
191
192   //Get the reference shape
193   Handle(GEOM_Object) aFace = GetObjectImpl(theFace);
194   if (aFace.IsNull()) return aGEOMObject._retn();
195
196   // Get the OptionalPoint (can be not defined)
197   Handle(GEOM_Object) anOptionalPoint = GetObjectImpl(theOptionalPoint);
198
199   // Make Vector - normal to theFace (in point theOptionalPoint if the face is not planar)
200   Handle(GEOM_Object) anObject = GetOperations()->GetNormal(aFace, anOptionalPoint);
201   if (!GetOperations()->IsDone() || anObject.IsNull())
202     return aGEOMObject._retn();
203
204   return GetObject(anObject);
205 }
206
207 //=============================================================================
208 /*!
209  *  GetBasicProperties
210  */
211 //=============================================================================
212 void GEOM_IMeasureOperations_i::GetBasicProperties (GEOM::GEOM_Object_ptr theShape,
213                                                     CORBA::Double& theLength,
214                                                     CORBA::Double& theSurfArea,
215                                                     CORBA::Double& theVolume)
216 {
217   //Set a not done flag
218   GetOperations()->SetNotDone();
219
220   //Get the reference shape
221   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
222   if (aShape.IsNull()) return;
223
224   // Get shape parameters
225   GetOperations()->GetBasicProperties(aShape, theLength, theSurfArea, theVolume);
226 }
227
228 //=============================================================================
229 /*!
230  *  GetInertia
231  */
232 //=============================================================================
233 void GEOM_IMeasureOperations_i::GetInertia
234   (GEOM::GEOM_Object_ptr theShape,
235    CORBA::Double& I11, CORBA::Double& I12, CORBA::Double& I13,
236    CORBA::Double& I21, CORBA::Double& I22, CORBA::Double& I23,
237    CORBA::Double& I31, CORBA::Double& I32, CORBA::Double& I33,
238    CORBA::Double& Ix , CORBA::Double& Iy , CORBA::Double& Iz)
239 {
240   //Set a not done flag
241   GetOperations()->SetNotDone();
242
243   //Get the reference shape
244   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
245   if (aShape.IsNull()) return;
246
247   // Get shape parameters
248   GetOperations()->GetInertia(aShape,
249                               I11, I12, I13,
250                               I21, I22, I23,
251                               I31, I32, I33,
252                               Ix , Iy , Iz);
253 }
254
255 //=============================================================================
256 /*!
257  *  GetBoundingBox
258  */
259 //=============================================================================
260 void GEOM_IMeasureOperations_i::GetBoundingBox (GEOM::GEOM_Object_ptr theShape,
261                                                 CORBA::Double& Xmin, CORBA::Double& Xmax,
262                                                 CORBA::Double& Ymin, CORBA::Double& Ymax,
263                                                 CORBA::Double& Zmin, CORBA::Double& Zmax)
264 {
265   //Set a not done flag
266   GetOperations()->SetNotDone();
267
268   //Get the reference shape
269   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
270   if (aShape.IsNull()) return;
271
272   // Get shape parameters
273   GetOperations()->GetBoundingBox(aShape, Xmin, Xmax, Ymin, Ymax, Zmin, Zmax);
274 }
275
276 //=============================================================================
277 /*!
278  *  GetTolerance
279  */
280 //=============================================================================
281 void GEOM_IMeasureOperations_i::GetTolerance
282                                 (GEOM::GEOM_Object_ptr theShape,
283                                  CORBA::Double& FaceMin, CORBA::Double& FaceMax,
284                                  CORBA::Double& EdgeMin, CORBA::Double& EdgeMax,
285                                  CORBA::Double& VertMin, CORBA::Double& VertMax)
286 {
287   //Set a not done flag
288   GetOperations()->SetNotDone();
289
290   //Get the reference shape
291   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
292   if (aShape.IsNull()) return;
293
294   // Get shape parameters
295   GetOperations()->GetTolerance(aShape,
296                                 FaceMin, FaceMax,
297                                 EdgeMin, EdgeMax,
298                                 VertMin, VertMax);
299 }
300
301 //=============================================================================
302 /*!
303  *  CheckShape
304  */
305 //=============================================================================
306 CORBA::Boolean GEOM_IMeasureOperations_i::CheckShape (GEOM::GEOM_Object_ptr theShape,
307                                                       CORBA::String_out     theDescription)
308 {
309   //Set a not done flag
310   GetOperations()->SetNotDone();
311
312   if (CORBA::is_nil(theShape))
313   {
314     theDescription = CORBA::string_dup("null");
315     return 0;
316   }
317
318   //Get the reference shape
319   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
320
321   if (aShape.IsNull())
322   {
323     theDescription = CORBA::string_dup("null2");
324     return 0;
325   }
326
327   // Get shape parameters
328   TCollection_AsciiString aDump;
329   if (GetOperations()->CheckShape(aShape, /*check_geom = */false, aDump))
330   {
331     theDescription = CORBA::string_dup("OK");
332     return 1;
333   }
334   theDescription = CORBA::string_dup(aDump.ToCString());
335   return 0;
336 }
337
338 CORBA::Boolean GEOM_IMeasureOperations_i::CheckShapeWithGeometry (GEOM::GEOM_Object_ptr theShape,
339                                                                   CORBA::String_out     theDescription)
340 {
341   //Set a not done flag
342   GetOperations()->SetNotDone();
343
344   if (CORBA::is_nil(theShape))
345   {
346     theDescription = CORBA::string_dup("null");
347     return 0;
348   }
349
350   //Get the reference shape
351   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
352
353   if (aShape.IsNull())
354   {
355     theDescription = CORBA::string_dup("null2");
356     return 0;
357   }
358
359   // Get shape parameters
360   TCollection_AsciiString aDump;
361   if (GetOperations()->CheckShape(aShape, /*check_geom = */true, aDump))
362   {
363     theDescription = CORBA::string_dup("OK");
364     return 1;
365   }
366   theDescription = CORBA::string_dup(aDump.ToCString());
367   return 0;
368 }
369
370 //=============================================================================
371 /*!
372  *  IsGoodForSolid
373  */
374 //=============================================================================
375 char* GEOM_IMeasureOperations_i::IsGoodForSolid (GEOM::GEOM_Object_ptr theShape)
376 {
377   //Set a not done flag
378   GetOperations()->SetNotDone();
379
380   //Get the reference shape
381   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
382   if (aShape.IsNull()) return CORBA::string_dup("WRN_NULL_OBJECT_OR_SHAPE");
383
384   // Get shape parameters
385   TCollection_AsciiString aDescription = GetOperations()->IsGoodForSolid(aShape);
386   return CORBA::string_dup(aDescription.ToCString());
387 }
388
389 //=============================================================================
390 /*!
391  *  WhatIs
392  */
393 //=============================================================================
394 char* GEOM_IMeasureOperations_i::WhatIs (GEOM::GEOM_Object_ptr theShape)
395 {
396   //Set a not done flag
397   GetOperations()->SetNotDone();
398
399   //Get the reference shape
400   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
401   if (aShape.IsNull()) return NULL;
402
403   // Get shape parameters
404   TCollection_AsciiString aDescription = GetOperations()->WhatIs(aShape);
405   return CORBA::string_dup(aDescription.ToCString());
406 }
407
408 //=============================================================================
409 /*!
410  *  AreCoordsInside
411  */
412 //=============================================================================
413 GEOM::ListOfBool* GEOM_IMeasureOperations_i::AreCoordsInside (GEOM::GEOM_Object_ptr theShape,
414                                                               const GEOM::ListOfDouble& theCoords,
415                                                               CORBA::Double tolerance)
416 {
417   //Set a not done flag
418   GetOperations()->SetNotDone();
419   
420   unsigned int nb_points = theCoords.length()/3;
421
422   GEOM::ListOfBool_var aResults = new GEOM::ListOfBool;
423   aResults->length(nb_points);
424   
425   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
426   
427   std::vector<double> tmp(3*nb_points);
428   for (int i = 0; i < 3*nb_points; i++)
429     tmp[i] = theCoords[i];
430   std::vector<bool> res = GetOperations()->AreCoordsInside(aShape, tmp, tolerance);
431   for (int i = 0; i < nb_points; i++)
432     aResults[i] = i < res.size() ? res[i] : false;
433   return aResults._retn();
434 }
435
436 //=============================================================================
437 /*!
438  *  GetMinDistance
439  */
440 //=============================================================================
441 CORBA::Double GEOM_IMeasureOperations_i::GetMinDistance
442   (GEOM::GEOM_Object_ptr theShape1, GEOM::GEOM_Object_ptr theShape2,
443    CORBA::Double& X1, CORBA::Double& Y1, CORBA::Double& Z1,
444    CORBA::Double& X2, CORBA::Double& Y2, CORBA::Double& Z2)
445 {
446   //Set a not done flag
447   GetOperations()->SetNotDone();
448
449   //Get the reference shape
450   Handle(GEOM_Object) aShape1 = GetObjectImpl(theShape1);
451   Handle(GEOM_Object) aShape2 = GetObjectImpl(theShape2);
452   if (aShape1.IsNull() || aShape2.IsNull()) return -1.0;
453
454   // Get shape parameters
455   return GetOperations()->GetMinDistance(aShape1, aShape2, X1, Y1, Z1, X2, Y2, Z2);
456 }
457
458 //=============================================================================
459 /*!
460  *  PointCoordinates
461  */
462 //=============================================================================
463 void GEOM_IMeasureOperations_i::PointCoordinates (GEOM::GEOM_Object_ptr theShape,
464                                                   CORBA::Double& X, CORBA::Double& Y, CORBA::Double& Z)
465
466 {
467   //Set a not done flag
468   GetOperations()->SetNotDone();
469
470   //Get the reference shape
471   Handle(GEOM_Object) aShape = GetObjectImpl(theShape);
472   if (aShape.IsNull())
473     return;
474
475   // Get shape parameters
476   GetOperations()->PointCoordinates( aShape, X, Y, Z );
477 }
478
479 //=============================================================================
480 /*!
481  *  GetAngle
482  */
483 //=============================================================================
484 CORBA::Double GEOM_IMeasureOperations_i::GetAngle (GEOM::GEOM_Object_ptr theShape1,
485                                                    GEOM::GEOM_Object_ptr theShape2)
486 {
487   //Set a not done flag
488   GetOperations()->SetNotDone();
489
490   //Get the reference shapes
491   Handle(GEOM_Object) aShape1 = GetObjectImpl(theShape1);
492   Handle(GEOM_Object) aShape2 = GetObjectImpl(theShape2);
493   if (aShape1.IsNull() || aShape2.IsNull()) return -1.0;
494
495   // Get the angle
496   return GetOperations()->GetAngle(aShape1, aShape2);
497 }
498
499 //=============================================================================
500 /*!
501  *  GetAngle
502  */
503 //=============================================================================
504 CORBA::Double GEOM_IMeasureOperations_i::GetAngleBtwVectors (GEOM::GEOM_Object_ptr theShape1,
505                                                              GEOM::GEOM_Object_ptr theShape2)
506 {
507   //Set a not done flag
508   GetOperations()->SetNotDone();
509
510   //Get the reference shapes
511   Handle(GEOM_Object) aShape1 = GetObjectImpl(theShape1);
512   Handle(GEOM_Object) aShape2 = GetObjectImpl(theShape2);
513   if (aShape1.IsNull() || aShape2.IsNull()) return -1.0;
514
515   // Get the angle
516   return GetOperations()->GetAngleBtwVectors(aShape1, aShape2);
517 }
518
519
520 //=============================================================================
521 /*!
522  *  CurveCurvatureByParam
523  */
524 //=============================================================================
525 CORBA::Double GEOM_IMeasureOperations_i::CurveCurvatureByParam
526                        (GEOM::GEOM_Object_ptr theCurve, CORBA::Double theParam)
527 {
528   //Set a not done flag
529   GetOperations()->SetNotDone();
530
531   //Get the reference shape
532   Handle(GEOM_Object) aShape = GetObjectImpl(theCurve);
533   if(aShape.IsNull()) return -1.0;
534
535   return GetOperations()->CurveCurvatureByParam(aShape,theParam);
536 }
537
538 //=============================================================================
539 /*!
540  *  CurveCurvatureByPoint
541  */
542 //=============================================================================
543 CORBA::Double GEOM_IMeasureOperations_i::CurveCurvatureByPoint
544                (GEOM::GEOM_Object_ptr theCurve, GEOM::GEOM_Object_ptr thePoint)
545 {
546   //Set a not done flag
547   GetOperations()->SetNotDone();
548
549   //Get the reference shape
550   Handle(GEOM_Object) aShape = GetObjectImpl(theCurve);
551   Handle(GEOM_Object) aPoint = GetObjectImpl(thePoint);
552   if( aShape.IsNull() || aPoint.IsNull() ) return -1.0;
553
554   return GetOperations()->CurveCurvatureByPoint(aShape,aPoint);
555 }
556
557
558 //=============================================================================
559 /*!
560  *  MaxSurfaceCurvatureByParam
561  */
562 //=============================================================================
563 CORBA::Double GEOM_IMeasureOperations_i::MaxSurfaceCurvatureByParam
564                                                 (GEOM::GEOM_Object_ptr theSurf,
565                                                  CORBA::Double theUParam,
566                                                  CORBA::Double theVParam)
567 {
568   //Set a not done flag
569   GetOperations()->SetNotDone();
570
571   //Get the reference shape
572   Handle(GEOM_Object) aShape = GetObjectImpl(theSurf);
573   if(aShape.IsNull()) return -1.0;
574
575   return GetOperations()->MaxSurfaceCurvatureByParam(aShape,theUParam,theVParam);
576 }
577
578 //=============================================================================
579 /*!
580  *  MaxSurfaceCurvatureByPoint
581  */
582 //=============================================================================
583 CORBA::Double GEOM_IMeasureOperations_i::MaxSurfaceCurvatureByPoint
584                 (GEOM::GEOM_Object_ptr theSurf, GEOM::GEOM_Object_ptr thePoint)
585 {
586   //Set a not done flag
587   GetOperations()->SetNotDone();
588
589   //Get the reference shape
590   Handle(GEOM_Object) aShape = GetObjectImpl(theSurf);
591   Handle(GEOM_Object) aPoint = GetObjectImpl(thePoint);
592   if( aShape.IsNull() || aPoint.IsNull() ) return -1.0;
593
594   return GetOperations()->MaxSurfaceCurvatureByPoint(aShape,aPoint);
595 }
596
597 //=============================================================================
598 /*!
599  *  MinSurfaceCurvatureByParam
600  */
601 //=============================================================================
602 CORBA::Double GEOM_IMeasureOperations_i::MinSurfaceCurvatureByParam
603                                                 (GEOM::GEOM_Object_ptr theSurf,
604                                                  CORBA::Double theUParam,
605                                                  CORBA::Double theVParam)
606 {
607   //Set a not done flag
608   GetOperations()->SetNotDone();
609
610   //Get the reference shape
611   Handle(GEOM_Object) aShape = GetObjectImpl(theSurf);
612   if (aShape.IsNull()) return -1.0;
613
614   return GetOperations()->MinSurfaceCurvatureByParam(aShape,theUParam,theVParam);
615 }
616
617 //=============================================================================
618 /*!
619  *  MinSurfaceCurvatureByPoint
620  */
621 //=============================================================================
622 CORBA::Double GEOM_IMeasureOperations_i::MinSurfaceCurvatureByPoint
623                 (GEOM::GEOM_Object_ptr theSurf, GEOM::GEOM_Object_ptr thePoint)
624 {
625   //Set a not done flag
626   GetOperations()->SetNotDone();
627
628   //Get the reference shape
629   Handle(GEOM_Object) aShape = GetObjectImpl(theSurf);
630   Handle(GEOM_Object) aPoint = GetObjectImpl(thePoint);
631   if (aShape.IsNull() || aPoint.IsNull()) return -1.0;
632
633   return GetOperations()->MinSurfaceCurvatureByPoint(aShape,aPoint);
634 }