public double getValeurPente()
{
if (p42_valeurPente == null)
{
p42_valeurPente = FLabServices.createCalculLow().GetPenteFromPoints3D(p11_pointDbt.p10_coord, p12_pointFin.p10_coord);
}
return((double)p42_valeurPente);
}
public double getPente()
{
if (p33_penteValeur == null)
{
p33_penteValeur = FLabServices.createCalculLow().GetPente(getVecteurPenteMaxi());
}
return((double)p33_penteValeur);
}
public double[] getVecteurPenteMaxi()
{
if (p32_penteVecteur == null)
{
bool v_normaliser_vf = true;
p32_penteVecteur = FLabServices.createCalculLow().GetVecteurPenteMaxi(getNormaleDuPlan(), v_normaliser_vf);
}
return(p32_penteVecteur);
}
public double[] getNormaleDuPlan()
{
if (p31_normaleDuPlan == null)
{
bool p_exceptionSiNotAPlanElseVecteur000_vf = false;
p31_normaleDuPlan = FLabServices.createCalculLow().GetNormaleDuPlan(p01_pointsDeFacette[0].p10_coord, p01_pointsDeFacette[2].p10_coord, p01_pointsDeFacette[1].p10_coord, p_exceptionSiNotAPlanElseVecteur000_vf);
}
return(p31_normaleDuPlan);
}
public enum_qualificationMorpho_arc GetQualificationMorphoDeLArc(BeanArc_internal p_arc)
{
enum_qualificationMorpho_arc v_result;
try
{
//L'arc sépare 2 facettes (sauf en frontière)
//On va exprimer les 2 points n'appartenant pas à l'arc dans le plan de plus forte pente
//Si les 2 points sont 'au-dessus' alors on considère que l'arc est 'talweg'
//Si les 2 points sont 'au-dessous' alors on considère que l'arc est 'ligne de crète'
//Les autres cas correspondent, au plus, à des ruptures de pente
//-On détermine le plan.
//Ce plan doit être tel que:
//-sa pente est celle de l'arc
//-si on translate le vecteur correspondant selon un vecteur de même élévation (=selon la courbe de niveau),
//alors ce vecteur doit être inclus dans ce plan.
//Du coup, pour déterminer ce plan, on détermine on utilise le vecteur de l'arc et un vecteur normal en xy à ce vecteur
double[] v_vecteurArc = FLabServices.createCalculLow().GetVectorBrutFromTwoPoints(p_arc.p11_pointDbt.p10_coord, p_arc.p12_pointFin.p10_coord);
double[] v_vecteurNormalXy = new double[3] {
-1 * v_vecteurArc[1], v_vecteurArc[0], 0
};
double[] v_pointNormal = new double[3] {
v_vecteurNormalXy[0] + p_arc.p11_pointDbt.p10_coord[0], v_vecteurNormalXy[1] + p_arc.p11_pointDbt.p10_coord[1], p_arc.p11_pointDbt.p10_coord[2]
};
double[] v_vecteurNormalAuPlanDePenteMaxi;
bool p_exceptionSiNotAPlanElseVecteur000_vf = false;
v_vecteurNormalAuPlanDePenteMaxi = FLabServices.createCalculLow().GetNormaleDuPlan(p_arc.p11_pointDbt.p10_coord, p_arc.p12_pointFin.p10_coord, v_pointNormal, p_exceptionSiNotAPlanElseVecteur000_vf);
//=>On calcule la matrice inverse...
v_vecteurArc = FLabServices.createCalculLow().GetNormalisationVecteurXYZ(v_vecteurArc);
v_vecteurNormalXy = FLabServices.createCalculLow().GetNormalisationVecteurXYZ(v_vecteurNormalXy);
double[,] v_matriceDeRotation = FLabServices.createCalculLow().GetMatriceInverse3x3(v_vecteurArc, v_vecteurNormalXy, v_vecteurNormalAuPlanDePenteMaxi);
//...on l'applique sur les 2 points
List <double[]> v_coordonnees = new List <double[]>();
double[] v_coordDansLeRepereDuPlan;
BeanPoint_internal v_pointATester;
double[] v_vectorATester;
if (p_arc.p21_facetteGauche != null)
{
v_pointATester = p_arc.p21_facetteGauche.p01_pointsDeFacette.Where(c => c.p00_id != p_arc.p11_pointDbt.p00_id && c.p00_id != p_arc.p12_pointFin.p00_id).First();
v_vectorATester = FLabServices.createCalculLow().GetVectorBrutFromTwoPoints(p_arc.p11_pointDbt.p10_coord, v_pointATester.p10_coord);
v_coordDansLeRepereDuPlan = FLabServices.createCalculLow().GetProduitMatriceParVector(v_matriceDeRotation, v_vectorATester);
v_coordonnees.Add(v_coordDansLeRepereDuPlan);
}
if (p_arc.p22_facetteDroite != null)
{
v_pointATester = p_arc.p22_facetteDroite.p01_pointsDeFacette.Where(c => c.p00_id != p_arc.p11_pointDbt.p00_id && c.p00_id != p_arc.p12_pointFin.p00_id).First();
v_vectorATester = FLabServices.createCalculLow().GetVectorBrutFromTwoPoints(p_arc.p11_pointDbt.p10_coord, v_pointATester.p10_coord);
v_coordDansLeRepereDuPlan = FLabServices.createCalculLow().GetProduitMatriceParVector(v_matriceDeRotation, v_vectorATester);
v_coordonnees.Add(v_coordDansLeRepereDuPlan);
}
//On exploite l'info:
//(pour l'instant, on ne traite pas les arcs frontières)
if (v_coordonnees.Where(c => c[2] > 0).Count() > 1)
{
v_result = enum_qualificationMorpho_arc.talweg;
return(v_result);
}
if (v_coordonnees.Where(c => c[2] < 0).Count() > 1)
{
v_result = enum_qualificationMorpho_arc.crete;
return(v_result);
}
return(enum_qualificationMorpho_arc.autre);;
}
catch (Exception)
{
throw;
}
}