Esempio n. 1
0
        /// <summary>
        /// Carica una mesh poligonale da un Object File Format su una struttura DCEL.
        /// </summary>
        /// <param name="filename">Il percorso completo del file da caricare.</param>
        /// <returns>Una mesh DCEL se caricata correttamente, null altrimenti.</returns>
        public static DCELMesh LoadFromOFF(string filename)
        {
            Dictionary <int, List <int> > leavingEdges = new Dictionary <int, List <int> >();
            Hashtable    vIndexes = new Hashtable();
            DCELMesh     mesh     = new DCELMesh();
            StreamReader reader;

            char[] separators = new char[] { ' ', '\t' };
            int    vertex_count = 0, face_count = 0;

            try
            {
                reader = new StreamReader(filename);
            }
            catch (FileNotFoundException)
            {
                return(null);
            }
            string line = reader.ReadLine();

            if (line == "OFF" || line == "NOFF")
            {
                while ((line = reader.ReadLine()) != null)
                {
                    if (CheckLine(line))
                    {
                        string[] values = line.Split(separators, StringSplitOptions.RemoveEmptyEntries);

                        vertex_count = Convert.ToInt32(values[0]);
                        face_count   = Convert.ToInt32(values[1]);

                        //MESH INITIALIZATION

                        #region Vertices Creation
                        //
                        for (int i = 0; i < vertex_count; i++)
                        {
                            line = reader.ReadLine();

                            if (CheckLine(line))
                            {
                                values = line.Split(separators, StringSplitOptions.RemoveEmptyEntries);

                                // x, y, z coordinates
                                mesh.AddVertex(new DCELVertex(new Point3D(
                                                                  Double.Parse(values[0].Replace(',', '.')),
                                                                  Double.Parse(values[1].Replace(',', '.')),
                                                                  Double.Parse(values[2].Replace(',', '.'))), null));

                                leavingEdges.Add(i, new List <int>());

                                try
                                {
                                    vIndexes.Add(mesh.VertexList[i], i);
                                }
                                catch (Exception)
                                {
                                    //non possono esistere due vertici con le stesse coordinate
                                    return(null);
                                }
                            }
                            else
                            {
                                i--;
                            }
                        }
                        //
                        #endregion

                        #region Faces Creation
                        //
                        int he_counter = 0;

                        for (int i = 0; i < face_count; i++)
                        {
                            line = reader.ReadLine();

                            if (CheckLine(line))
                            {
                                values = line.Split(separators, StringSplitOptions.RemoveEmptyEntries);

                                //k è il numero di vertici della faccia corrente
                                int k = Convert.ToInt32(values[0]);

                                int  starting_index = he_counter;
                                bool check_first    = true;
                                mesh.AddFace(new DCELFace());

                                //for each vertex
                                for (int j = 1; j <= k; j++)
                                {
                                    //@vertex_index è l'indice del vertice all'interno della collezione VertexList
                                    int vertex_index = Convert.ToInt32(values[j]);

                                    #region Half Edges Calculation

                                    if (check_first)
                                    {
                                        check_first = false;
                                        //creo un nuovo halfedge e lo aggiungo alla lista
                                        mesh.AddHalfEdge(new DCELHalfEdge());
                                        //imposto come vertice di origine il point3d corrispondente all'indice vertex_index in VertexList
                                        mesh.HalfEdgeList[he_counter].Origin = mesh.VertexList[vertex_index];
                                        //imposto il primo halfedge come halfedge di riferimento della faccia corrente
                                        mesh.FaceList[i].Edge = mesh.HalfEdgeList[he_counter];
                                        //l'halfedge corrente si collega alla faccia corrente
                                        mesh.HalfEdgeList[he_counter].Face = mesh.FaceList[i];
                                        //
                                        leavingEdges[vertex_index].Add(he_counter);
                                        //incremento il contatore
                                        he_counter += 1;
                                    }
                                    else
                                    {
                                        //creo un nuovo halfedge e lo aggiungo alla lista
                                        mesh.AddHalfEdge(new DCELHalfEdge());
                                        //l'halfedge successivo al precedente è quello corrente. Incremento face_count e assegno
                                        mesh.HalfEdgeList[he_counter - 1].Next = mesh.HalfEdgeList[he_counter];
                                        //assegno il punto d'origine corrispondente all'indice vertex_index in Mesh.VertexList
                                        mesh.HalfEdgeList[he_counter].Origin = mesh.VertexList[vertex_index];
                                        //collego l'halfedge alla faccia corrente
                                        mesh.HalfEdgeList[he_counter].Face = mesh.FaceList[i];
                                        //
                                        leavingEdges[vertex_index].Add(he_counter);
                                        //incremento il contatore
                                        he_counter += 1;
                                    }
                                    //finito di iterare sui vertici, l'ultimo halfedge ha come next l'halfedge di partenza
                                    mesh.HalfEdgeList[he_counter - 1].Next = mesh.HalfEdgeList[starting_index];

                                    #endregion
                                }
                                //face normal
                                MathUtils.CalculateNormal(mesh.FaceList[i]);
                            }
                            else
                            {
                                i--;
                            }
                        }
                        //
                        #endregion

                        #region Leaving Edges and Vertex Normals Calculation

                        for (int i = 0; i < vertex_count; i++)
                        {
                            if (leavingEdges[i].Count > 0)
                            {
                                mesh.VertexList[i].Leaving = mesh.HalfEdgeList[leavingEdges[i][0]];
                                //face normal
                                mesh.VertexList[i].Normal = new Vector3D();
                                //for each leaving edge from the current vertex
                                foreach (var lEdgeIndex in leavingEdges[i])
                                {
                                    //sum the face normals for every triangle containing the vertex
                                    mesh.VertexList[i].Normal += mesh.HalfEdgeList[lEdgeIndex].Face.Normal;
                                }

                                //get average normal
                                mesh.VertexList[i].Normal /= leavingEdges[i].Count;
                                mesh.VertexList[i].Normal  = MathUtils.Normalize(mesh.VertexList[i].Normal);
                            }
                        }

                        #endregion

                        #region Twin Edges Calculation

                        foreach (var edge in mesh.HalfEdgeList)
                        {
                            //index è la posizione in VertexList del vertice dell'edge successivo a quello corrente
                            int index = (int)vIndexes[edge.Next.Origin];

                            //edge_index è la posizione in EdgeList dell'edge con origine nel vertice dell'indice di cui sopra!
                            foreach (var edge_index in leavingEdges[index])
                            {
                                if (edge.Origin == mesh.HalfEdgeList[edge_index].Next.Origin)
                                {
                                    edge.Twin = mesh.HalfEdgeList[edge_index];
                                    break;
                                }
                            }
                        }
                        //
                        #endregion

                        reader.Close();
                        //restituisco la mesh caricata
                        return(mesh);
                    }
                }
            }
            //if null is returned, something is gone wrong!
            return(null);
        }