From 9f188bce5ec41135563a1787b428179d75218f66 Mon Sep 17 00:00:00 2001
From: Boris Sedji <sedji.boris@hotmail.com>
Date: Fri, 18 Dec 2009 14:40:00 +0000
Subject: [PATCH]
---
contrib/arc/ImportLSImage.cpp | 92 +++++++++++++++++++++++++++++++++++
1 file changed, 92 insertions(+)
create mode 100644 contrib/arc/ImportLSImage.cpp
diff --git a/contrib/arc/ImportLSImage.cpp b/contrib/arc/ImportLSImage.cpp
new file mode 100644
index 0000000000..715a80504d
--- /dev/null
+++ b/contrib/arc/ImportLSImage.cpp
@@ -0,0 +1,92 @@
+//
+// Description:
+//
+//
+// Author: <Boris Sedji>, 12/2009
+//
+// Copyright: See COPYING file that comes with this distribution
+//
+//
+
+
+#ifndef _ITK_H_
+
+#include "itkImage.h"
+#include "itkImageFileReader.h"
+#include "itkImageFileWriter.h"
+#define _ITK_H_
+
+#endif
+
+#include "Octree.h"
+#include <stdio.h>
+#include "Gmsh.h"
+#include "GModel.h"
+#include "MElement.h"
+#include "gmshLevelset.h"
+#include "OctreeLSImage/OctreeLSImage.cpp"
+
+
+
+int main( int argc, char *argv[] )
+{
+ if( argc < 4 )
+ {
+ std::cerr << "Missing Parameters " << std::endl;
+ std::cerr << "Usage: " << argv[0];
+ std::cerr << " inputLevelSetImage MeshSizeMax MeshSizeMin"<< std::endl;
+ return 1;
+ }
+
+ typedef float PixelTypeFloat;
+ const unsigned int Dimension = 3;
+ typedef itk::Image< PixelTypeFloat, Dimension > ImageTypeFloat;
+ typedef itk::ImageFileReader< ImageTypeFloat > ImageReaderTypeFloat;
+ ImageReaderTypeFloat::Pointer reader = ImageReaderTypeFloat::New();
+
+ reader->SetFileName(argv[1]);
+
+ ImageTypeFloat::Pointer image = reader->GetOutput();
+ image->Update();
+
+ ImageTypeFloat::RegionType region;
+ region = image->GetLargestPossibleRegion ();
+
+ std::cout<<"\nImage dimensions : " << region.GetSize(0) << " x " << region.GetSize(1) << " x " << region.GetSize(2);
+
+ OctreeLSImage octree(image);
+
+ int sizemax = atoi(argv[2]);
+ int sizemin = atoi(argv[3]);
+
+ // Mesh with conditions on mesh size
+ octree.Mesh(sizemax,sizemin);
+
+ // Smoothing mesh to avoid too much generation between adjacent elements
+ bool statut=false;
+ int k = 0;
+ std::cout<<"\nLeaf Number : "<< (octree.GetLeafNumber())<<"\n";
+ while((!statut) & (k < 20)){
+ statut = octree.Smooth();
+ std::cout<<"\nk : "<< k;
+ std::cout<<"\nsmoothed : " << (int)statut;
+ std::cout<<"\nLeaf Number : "<< (octree.GetLeafNumber())<<"\n";
+ k++;
+ }
+
+ // Create GModel with the octree mesh
+ GModel *m = octree.CreateGModel();
+
+ // Write a .msh file with the mesh
+ std::string ModelName = "OctreeMesh.msh" ;
+ m->writeMSH(ModelName,2.1,false,false);
+
+ // Write a .msh file with the level set values as postview data
+ PView *pv = octree.CreateLSPView(m);
+ bool useadapt = true;
+ pv->getData(useadapt)->writeMSH("LSPView.msh", false);
+
+ std::cout<<"\n";
+
+ return 0;
+}
--
GitLab