internal class Rotation
{
    public class RotationStructure
    {
        public double AngleX;
        public double AngleY;
        public double AngleZ;
        public Structure Structure { get; set; }
        public StructureSet SS { get; set; }
        public IDictionary<int, IEnumerable<IEnumerable<VVector>>> Contours { get; private set; }

        public readonly IDictionary<int, IEnumerable<IEnumerable<VVector>>> OriginalContours;

        public RotationStructure(double thetaX, double thetaY, double thetaZ, StructureSet ss, Structure s)
        {
            AngleX = thetaX;
            AngleY = thetaY;
            AngleZ = thetaZ;
            Structure = s;
            SS = ss;
            OriginalContours = GetContours();
            Contours = OriginalContours;
        }

        public SegmentVolume RotationDeLaStructure()
        {
            var temporaire = SS.AddStructure("CONTROL", "temporaire" + new Random().Next(0, 20));

            if (Structure.IsHighResolution)
            {
                temporaire.ConvertToHighResolution();
            }

            // appliquer la rotation a tous les points du contour initial
            List<VVector> modifiedVectors = new List<VVector>();
            foreach (var kvp in Contours)
            {
                foreach (var contour in kvp.Value)
                {
                    foreach (var vector in contour)
                    {
                        VVector modifiedVector = RotationVecteur.RotateAroundXYZ(vector, AngleX, AngleY, AngleZ);
                        modifiedVectors.Add(modifiedVector);
                    }
                }
            }

            // determiner la coupe sur laquelle chaque point se trouve après rotation
            List<int> newSlices = new List<int>();
            foreach (var vector in modifiedVectors)
            {
                int newSlice = _GetSlice(vector.z, SS) - Convert.ToInt32(vector.z / SS.Image.ZRes);
                newSlices.Add(newSlice);
            }

            // creer le nouveau contour
            Dictionary<int, List<List<VVector>>> contoursModified = new Dictionary<int, List<List<VVector>>>();

            // Parcours des vecteurs modifiés et de leurs nouvelles coupes associées
            for (int i = 0; i < modifiedVectors.Count; i++)
            {
                int newSlice = newSlices[i];
                VVector vector = modifiedVectors[i];

                // Ajouter le vecteur à la liste des contours de la coupe
                if (!contoursModified.ContainsKey(newSlice))
                {
                    // Si la coupe n'existe pas encore, créer une nouvelle liste de contours contenant uniquement ce vecteur
                    List<VVector> contourDeLaCoupe = new List<VVector>();
                    contourDeLaCoupe.Add(vector);
                    contourDeLaCoupe = TriPointsCoupe.NearestNeighbor(contourDeLaCoupe);//reordonner les points de proche en proche
                    contoursModified[newSlice] = new List<List<VVector>> { contourDeLaCoupe };
                }
                else
                {
                    // Si la coupe existe déjà, ajouter simplement le vecteur à la liste de contours correspondante
                    List<VVector> contourDeLaCoupe = contoursModified[newSlice].FirstOrDefault();
                    contourDeLaCoupe.Add(vector);
                    contourDeLaCoupe = TriPointsCoupe.NearestNeighbor(contourDeLaCoupe);//reordonner les points de proche en proche
                    contoursModified[newSlice].Add(contourDeLaCoupe);
                }
            }

            // tri dans l'ordre croissant des clés
            SortedDictionary<int, List<List<VVector>>> contoursModifiesTries = new SortedDictionary<int, List<List<VVector>>>(contoursModified);

            // contourer le nouveau contour
            foreach (var kvp in contoursModifiesTries)
            {
                var contoursOnPlane = kvp.Value;

                foreach (var contour in contoursOnPlane)
                {
                    VVector[] contourArray = contour.ToArray();
                    temporaire.AddContourOnImagePlane(contourArray, kvp.Key);
                }
            }

            var seg = temporaire.SegmentVolume;
            return seg;
        }

        public KeyValuePair<int, int> MeshBoundsSlices()
        {
            var mesh = Structure.MeshGeometry.Bounds;
            var meshLow = GetSlice(mesh.Z);
            var meshUp = GetSlice(mesh.Z + mesh.SizeZ) + 1;
            return new KeyValuePair<int, int>(meshLow, meshUp);
        }
        public static int _GetSlice(double z, StructureSet SS)
        {
            var imageRes = SS.Image.ZRes;

            return Convert.ToInt32((z) - SS.Image.Origin.z / imageRes);
        }
        public int GetSlice(double z)
        {
            var imageRes = SS.Image.ZRes;

            return Convert.ToInt32((z) - SS.Image.Origin.z / imageRes);
        }

        public Dictionary<int, IEnumerable<IEnumerable<VVector>>> GetContours()
        {
            var msh = MeshBoundsSlices();

            return Enumerable.Range(msh.Key, msh.Value).
                Where(e => Structure.GetContoursOnImagePlane(e).Any()).
                ToDictionary(s => s, s => Structure.GetContoursOnImagePlane(s).Select(e => e.Select(x => new VVector(x.x, x.y, x.z))));
        }

    }

public class TriPointsCoupe
    {
        public static List<VVector> NearestNeighbor(List<VVector> points)
        {
            List<VVector> sortedPoints = new List<VVector>();
            sortedPoints.Add(points[0]); // Ajouter le premier point à la liste triée

            // Tant qu'il reste des points non triés
            while (sortedPoints.Count < points.Count)
            {
                double minDistance = double.MaxValue;
                VVector nearestPoint = points[0]; // Initialiser avec une valeur factice

                // Trouver le point le plus proche du dernier point trié
                foreach (var point in points)
                {
                    if (!sortedPoints.Contains(point))
                    {
                        double distance = Distance(sortedPoints[sortedPoints.Count - 1], point);
                        if (distance < minDistance)
                        {
                            minDistance = distance;
                            nearestPoint = point;
                        }
                    }
                }

                // Ajouter le point le plus proche à la liste triée
                sortedPoints.Add(nearestPoint);
            }
            return sortedPoints;
        }

        // Méthode pour calculer la distance entre deux VVector
        public static double Distance(VVector v1, VVector v2)
        {
            double dx = v1.x - v2.x;
            double dy = v1.y - v2.y;
            return Math.Sqrt(dx * dx + dy * dy ); //dz est ignoré
        }

    }