La plupart des simulations interactives ont besoin d'un modèle de détection de collisions. Cette détection nécessite d'une part d'effectuer des requêtes de proximité entre les entités concernées et d'autre part de calculer un comportement à appliquer. Afin d'effectuer ces requêtes, les entités présentes dans une scène sont soit hiérarchisées dans un arbre ou dans un graphe de proximité, soit plongées dans une grille d'enregistrement. Nous présentons un nouveau modèle de détection de collisions s'appuyant sur deux piliers : une représentation de l'environnement par des cartes combinatoires multirésolutions et un suivi en temps réel de particules plongées dans ces cartes. Ce modèle nous permet de représenter des environnements complexes tout en suivant en temps réel les entités évoluant dans cet environnement. Nous présentons des outils d'enregistrement et de maintien de l'enregistrement de particules, d'arêtes et de surfaces dans des cartes combinatoires volumiques multirésolutions. / Most interactive simulations need a collision detection system. First, this system requires the querying of the proximity between the objects and then the computing of the behaviour to be applied. In order to perform these queries, the objects present in a scene are either classified in a tree, in a proximity graph, or embedded inside a registration grid.Our work present a new collision detection model based on two main concepts: representing the environment with a combinatorial multiresolution map, and tracking in real-time particles embedded inside this map. This model allows us to simulate complex environments while following in real-time the entities that are evolving within it.We present our framework used to register and update the registration of particles, edges and surfaces in volumetric combinatorial multiresolution maps. Results have been validated first in 2D with a crowd simulation application and then in 3D, in the medical field, with a percutaneous surgery simulation.

Nous étudions dans cette thèse la génération de maillages anisotropes basée sur la triangulation de Delaunay et le diagramme de Voronoi. Nous considérons tout d'abord les maillages anisotropes localement uniformes, développés par Boissonnat, Wormser et Yvinec. Bien que l'aspect théorique de cette approche soit connu, son utilité pratique n'a été que peu explorée. Une étude empirique exhaustive est présentée et révèle les avantages, mais aussi les inconvénients majeurs de cette méthode. Dans un second temps, nous étudions les diagrammes de Voronoi anisotropes définis par Labelle et Shewchuk. Nous donnons des conditions suffisantes sur un ensemble de points pour que le dual du diagramme soit une triangulation plongée en toute dimension ; un algorithme générant de tels ensembles est conçu. Ce diagramme est utilisé pour concevoir un algorithme qui génère efficacement un maillage anisotrope pour des domaines de dimension intrinsèque faible plongés dans des espaces de dimension large. Notre algorithme est prouvable, mais les résultats sont décevants. Enfin, nous présentons le diagramme de Voronoi Riemannien discret, qui utilise des avancées récentes dans l'estimation de distances géodésiques et dont le calcul est grandement accéléré par l'utilisation d'un graphe anisotrope. Nous donnons des conditions suffisantes pour que notre structure soit combinatoirement équivalente au diagramme de Voronoi Riemannien et que son dual utilisant des simplexes droits mais aussi courbes est une triangulation plongée en toute dimension. Nous obtenons de bien meilleurs résultats que pour nos autres techniques, mais dont l'utilité reste limitée / In this thesis, we study the generation of anisotropic meshes using the concepts of Delaunay triangulations and Voronoi diagrams. We first consider the framework of locally uniform anisotropic meshes introduced by Boissonnat, Wormser and Yvinec. Despite known theoretical guarantees, the practicality of this approach has only been hardly studied. An exhaustive empirical study is presented and reveals the strengths but also the overall impracticality of the method. In a second part, we investigate the anisotropic Voronoi diagram introduced by Labelle and Shewchuk and give conditions on a set of seeds such that the corresponding diagram has a dual that is an embedded triangulation in any dimension; an algorithm to generate such sets is devised. Using the same diagram, we propose an algorithm to generate efficiently anisotropic triangulations of low-dimensional manifolds embedded in high-dimensional spaces. Our algorithm is provable, but produces disappointing results. Finally, we study Riemannian Voronoi diagrams and introduce discrete Riemannian Voronoi diagrams, which employ recent developments in the numerical computation of geodesic distances and whose computation is accelerated through the use of an underlying anisotropic graph structure. We give conditions that guarantee that our discrete structure is combinatorially equivalent to the Riemannian Voronoi diagram and that its dual is an embedded triangulation, using both straight and curved simplices. We obtain significantly better results than with our other methods, but the overall utility of

De nombreux scanners laser permettent d'obtenir la surface 3D a partir d'un objet. Néanmoins, la surface reconstruite est souvent lisse, ce qui est du au débruitage interne du scanner et aux décalages entre les scans. Cette these utilise des scans haute precision et choisit de ne pas perdre ni alterer les echantillons initiaux au cours du traitement afin de les visualiser. C'est en effet la seule façon de decouvrir les imperfections (trous, decalages de scans). De plus, comme les donnees haute precision capturent meme le plus leger detail, tout debruitage ou sous-echantillonnage peut amener a perdre ces details.La these s'attache a prouver que l'on peut trianguler le nuage de point initial en ne perdant presque aucun echantillon. Le probleme de la visualisation exacte sur des donnees de plus de 35 millions de points et de 300 scans differents est ainsi resolu. Deux problemes majeurs sont traites: le premier est l'orientation du nuage de point brut complet et la creation d'un maillage. Le second est la correction des petits decalages entre les scans qui peuvent creer un tres fort aliasing et compromettre la visualisation de la surface. Le second developpement de la these est une decomposition des nuages de points en hautes/basses frequences. Ainsi, des methodes classiques pour l'analyse d'image, l'arbre des ensembles de niveau et la representation MSER, sont etendues aux maillages, ce qui donne une methode intrinseque de segmentation de maillages. Une analyse mathematiques d'operateurs differentiels discrets, proposes dans la litterature et operant sur des nuages de points est realisee. En considerant les developpements asymptotiques de ces operateurs sur une surface reguliere, ces operateurs peuvent etre classifies. Cette analyse amene au developpement d'un operateur discret consistant avec Ie mouvement par courbure moyenne (l'equation de la chaleur intrinseque) definissant ainsi un espace-echelle numerique simple et remarquablement robuste. Cet espace-echelle permet de resoudre de maniere unifiee tous les problemes mentionnes auparavant (orientation et triangulation du nuage de points, fusion de scans, segmentation de maillages) qui sont ordinairement traites avec des techniques distinctes. / Many laser devices acquire directly 3D objects and reconstruct their surface. Nevertheless, the final reconstructed surface is usually smoothed out as a result of the scanner internal de-noising process and the offsets between different scans. This thesis, working on results from high precision scans, adopts the somewhat extreme conservative position, not to loose or alter any raw sample throughout the whole processing pipeline, and to attempt to visualize them. Indeed, it is the only way to discover all surface imperfections (holes, offsets). Furthermore, since high precision data can capture the slightest surface variation, any smoothing and any sub-sampling can incur in the loss of textural detail.The thesis attempts to prove that one can triangulate the raw point cloud with almost no sample loss. It solves the exact visualization problem on large data sets of up to 35 million points made of 300 different scan sweeps and more. Two major problems are addressed. The first one is the orientation of the complete raw point set, an the building of a high precision mesh. The second one is the correction of the tiny scan misalignments which can cause strong high frequency aliasing and hamper completely a direct visualization.The second development of the thesis is a general low-high frequency decomposition algorithm for any point cloud. Thus classic image analysis tools, the level set tree and the MSER representations, are extended to meshes, yielding an intrinsic mesh segmentation method.The underlying mathematical development focuses on an analysis of a half dozen discrete differential operators acting on raw point clouds which have been proposed in the literature. By considering the asymptotic behavior of these operators on a smooth surface, a classification by their underlying curvature operators is obtained.This analysis leads to the development of a discrete operator consistent with the mean curvature motion (the intrinsic heat equation) defining a remarkably simple and robust numerical scale space. By this scale space all of the above mentioned problems (point set orientation, raw point set triangulation, scan merging, segmentation), usually addressed by separated techniques, are solved in a unified framework.

L'acquisition de modèles 3D des scènes réelles trouve son utilité dans de nombreuses applications pratiques, comme l'archivage numérique, les jeux vidéo, l'ingénierie, la publicité. Il existe principalement deux méthodes pour acquérir un modèle 3D: la reconstruction avec un scanner laser (méthode active) et la reconstruction à partir de plusieurs photographies d'une même scène prise dans des points de vues différentes (méthode passive). Si la méthode active permet d'acquérir des modèles avec une grande précision, il est cependant coûteux et difficile à mettre en place pour de grandes scènes extérieures. La méthode passive, ou la stéréo multi-vues est en revanche plus flexible, facile à mettre en oeuvre et surtout moins coûteuse que la méthode active. Cette thèse s'attaque au problème de la reconstruction de stéréo multi-vues à grande échelle et précise pour les scènes extérieures. Nous améliorons des méthodes précédentes et les assemblons pour créer une chaîne de stéréo multi-vues efficace tirant parti de l'accélération de cartes graphiques. La chaîne produit des maillages de qualité à partir d'images de haute résolution, ce qui permet d'atteindre les meilleurs scores dans de nombreuses évaluations. Aux plus grandes échelles, nous développons d'une part des techniques de type diviser-pour-régner pour reconstruire des morceaux partiaux de la scène. D'autre part, pour combiner ces résultats séparés, nous créons une nouvelle méthode qui fusionne rapidement des centaines de maillages. Nous réussissons à reconstruire de beaux maillages urbains et des monuments historiques précis à partir de grandes collections d'images (environ 1600 images de 5M Pixel) / Acquisition of 3D model of real objects and scenes is indispensable and useful in many practical applications, such as digital archives, game and entertainment industries, engineering, advertisement. There are 2 main methods for 3D acquisition : laser-based reconstruction (active method) and image-based reconstruction from multiple images of the scene in different points of view (passive method). While laser-based reconstruction achieves high accuracy, it is complex, expensive and difficult to set up for large-scale outdoor reconstruction. Image-based, or multi-view stereo methods are more versatile, easier, faster and cheaper. By the time we begin this thesis, most multi-view methods could handle only low resolution images under controlled environment. This thesis targets multi-view stereo both both in large scale and high accuracy issues. We significantly improve some previous methods and combine them into a remarkably effective multi-view pipeline with GPU acceleration. From high-resolution images, we produce highly complete and accurate meshes that achieve best scores in many international recognized benchmarks. Aiming even larger scale, on one hand, we develop Divide and Conquer approaches in order to reconstruct many small parts of a big scene. On the other hand, to combine separate partial results, we create a new merging method, which can merge automatically and quickly hundreds of meshes. With all these components, we are successful to reconstruct highly accurate water-tight meshes for cities and historical monuments from large collections of high-resolution images (around 1600 images of 5 M Pixel images)

Dans ce travail, nous étudions le cadre de la morphologie mathématique sur les complexes simpliciaux. Complexes simpliciaux sont une structure versatile et largement utilisée pour représenter des données multidimensionnelles, telles que des maillages, qui sont des complexes tridimensionnels, ou des graphes, qui peuvent être interprétées comme des complexes bidimensionnels. La morphologie mathématique est l'un des cadres les plus puissants pour le traitement de l'image, y compris le traitement des structures numériques, et est largement utilisé pour de nombreuses applications. Toutefois, les opérateurs de morphologie mathématique sur des espaces complexes simpliciaux n'est pas un concept entièrement développé dans la littérature. Dans ce travail, nous passons en revue certains opérateurs classiques des complexes simpliciaux sous la lumière de la morphologie mathématique, de montrer qu'ils sont des opérateurs de morphologie. Nous définissons certains treillis de base et les opérateurs agissant sur ces treillis: dilatations, érosions, ouvertures, fermetures et filtres alternés séquentiels, et aussi leur extension à simplexes pondérés. Cependant, les principales contributions de ce travail sont ce que nous appelions les opérateurs dimensionnels, petites et polyvalents opérateurs qui peuvent être utilisés pour définir de nouveaux opérateurs sur les complexes simpliciaux, qui garde les propriétés de la morphologie mathématique. Ces opérateurs peuvent également être utilisés pour exprimer pratiquement n'importe quel opérateur dans la littérature. Nous illustrons les opérateurs définis et nous comparons les filtres alternés séquentiels contre filtres définis dans la littérature, où nos filtres présentent de meilleurs résultats pour l'enlèvement du petit, intense bruit des images binaires / In this work we study the framework of mathematical morphology on simplicial complex spaces. Simplicial complexes are a versatile and widely used structure to represent multidimensional data, such as meshes, that are tridimensional complexes, or graphs, that can be interpreted as bidimensional complexes. Mathematical morphology is one of the most powerful frameworks for image processing, including the processing of digital structures, and is heavily used for many applications. However, mathematical morphology operators on simplicial complex spaces is not a concept fully developped in the literature. In this work, we review some classical operators from simplicial complexes under the light of mathematical morphology, to show that they are morphology operators. We define some basic lattices and operators acting on these lattices: dilations, erosions, openings, closings and alternating sequential filters, including their extension to weighted simplexes. However, the main contributions of this work are what we called dimensional operators, small, versatile operators that can be used to define new operators on simplicial complexes, while mantaining properties from mathematical morphology. These operators can also be used to express virtually any operator from the literature. We illustrate all the defined operators and compare the alternating sequential filters against filters defined in the literature, where our filters show better results for removal of small, intense, noise from binary images

Ces dernières décennies, de nouvelles méthodes numériques connues sous le nom de « méthodes sans maillage » ont été développées. Contrairement à la MEF, ces méthodes n'utilisent qu'un ensemble de noeuds répartis dans le domaine sans demander un maillage de celui-ci. Jusqu'à présent, aucune de ces méthodes n'est parvenue à satisfaire les utilisateurs de la MEF. Dans cette thèse, nous proposons une méthode sans maillage, utilisant les approximations de Taylor. Cette méthode a l'avantage de n'utiliser que des points sur la frontière. En effet, l'EDP est résolue sous sa forme forte dans le domaine et les conditions aux limites sont appliquées par la méthode des moindres carrés. Cette méthode a été introduite, il y a 3 ans par S. Zeze dans sa thèse. Les tests académiques effectués en linéaire ont montré que cette méthode est très précise et que la convergence est améliorée en augmentant le degré, comme dans la p-version des EF. Nos travaux de thèse sont une suite des travaux de S. Zeze et ils visent à rendre plus robuste la méthode et aussi à élargir son champ d'application. Dans un premier temps, nous faisons une analyse mathématique de la méthode. Cette analyse passe par l'analyse des séries calculées. Le but de cette analyse est d'évaluer le domaine de convergence de la solution. Les résultats obtenus ont montré que pour certains problèmes, il faut subdiviser le domaine en quelques sous domaines et faire une résolution par sous domaine. La suite de nos travaux a donc été d'établir une technique de raccordement qui permettra d'assurer les conditions de transmission aux interfaces, dans le cas d'une résolution par sous domaine. En dernière partie, nous étendons l'application de la méthode aux problèmes non linéaires, en la couplant à une méthode de linéarisation / In these last decades, new numerical methods known as « meshless methods » have been developped. Contrary to the FEM, these methods uses only a set of nodes in the domain, without need of any mesh. Until now, any of these methods has convinced users of FEM. In this paper, we present a new meshless method using Taylor series expansion. In this method, the PDE is solved quasi exactly in the domain and the boundary conditions are applied by using a least square method. Then only the boundary discretisation is needed so the proposed method is a « true boundary meshless method ». This technique has been proposed for the first time by S. Zeze in his PhD thesis. The study of some linear problems has shown that this technique leads to a very good accuracy and that the convergence can be improved by increasing approximation degree. Our work is a continuation of S. Zeze work, and it consists to make the proposed method more robust and to extend its range of application. For that, we first make an analysis of the series computed by the method. The aim of this analysis was to evaluate the domain of validity of these series. This analysis showed that, for some problems, an accuracy cannot be obtained without splitting the domain in subdomains and making a resolution by subdomains. Therefore the second part of our work was to define a technique which will ensure the continuity at the interface between subdomains, in the case of a resolution by subdomains. The last part of our work was dedicated to non-linear problems. We establish an algorithm to show how the proposed method can deal with nonlinear-problems

Cette thèse s'intéresse à la conception de méthodes de discrétisation non-conforme pour un modèle de poromécanique. Le but de ce travail est de simplifier les couplages liant la géomécanique d'un milieu poreux à l'écoulement polyphasique compositionnel ayant cours en son sein tels qu'ils sont réalisés actuellement dans l'industrie pétrolière, en discrétisant sur un même maillage, typiquement non-conforme car à l'image de la lithologie, la mécanique et l'écoulement. La nouveauté consiste donc à traiter la mécanique par une méthode d'approximation non-conforme sur maillages généraux. Dans cette thèse, nous nous concentrons sur un modèle d'élasticité linéaire. Les difficultés inhérentes à son approximation non-conforme sont son manque de coercivité (se traduisant par la nécessité de satisfaire une inégalité de Korn sur un espace discret discontinu), ainsi que le phénomène de verrouillage numérique lorsque le matériau tend à devenir incompressible. Dans une première partie, nous construisons un espace d'approximation sur maillages généraux, s'apparentant à une extension de l'espace de Crouzeix-Raviart. Nous explicitons ses propriétés d'approximation et de conformité, et montrons que ce dernier est adapté à une discrétisation primale coercive et robuste au locking du modèle d'élasticité sur maillages généraux. La méthode proposée est moins coûteuse que son équivalent éléments finis (en termes de propriétés) P2. Nous nous intéressons dans une deuxième partie à l'approximation non-conforme d'un modèle couplé de poroélasticité. Nous étudions la convergence d'une famille de schémas numériques dont la discrétisation en espace utilise le formalisme des schémas Gradient, auquel appartient la méthode développée pour la mécanique. Nous prouvons la convergence de telles approximations vers la solution de régularité minimale du problème continu, indépendamment des paramètres physiques du système / This manuscript focuses on the conception of nonconforming discretization methods for a poromechanical model. The aim of this work is to ease the coupling between the geomechanics and the multiphase compositional Darcy flow in porous media by discretizing mechanics and flow on the same mesh, typically nonconforming as it represents the lithology. Hence, the novelty hinges on a nonconforming treatment of mechanics on general meshes. In this work, we focus on a linear elasticity model. The nonconforming approximation of such a model is not straightforward owing to its lack of coercivity (meaning that a discrete Korn's inequality must hold on a discontinuous discrete space) and to the numerical locking phenomenon occurring as the material becomes incompressible. In a first part, we design an approximation space on general meshes, which can be viewed as an extension of the so-called Crouzeix-Raviart space. We study its approximation and conformity properties, and prove that this latter is well-adapted to the design of a primal, coercive, and locking-free discretization of the elasticity model on general meshes. The proposed method is less costly than its finite element equivalent (in terms of properties) P2. In a second part, we tackle the nonconforming approximation of a coupled poroelasticity model. We study the convergence of a family of numerical schemes whose space discretization relies on the Gradient schemes framework, to which belongs the method developed for mechanics. We prove the convergence of such approximations toward the minimal regularity solution of the continuous problem, and independently of the choice of physical parameters

L'objectif de cette thèse était de proposer, puis de valider expérimentalement, un ensemble de méthodes algorithmiques permettant le remaillage parallèle de maillages distribués, en s'appuyant sur une méthode séquentielle de remaillage préexistante. Cet objectif a été atteint par étapes : définition de structures de données et de schémas de communication adaptés aux maillages distribués, permettant le déplacement à moindre coût des interfaces entre sous-domaines sur les processeurs d'une architecture à mémoire distribuée ; utilisation d'algorithmes de répartition dynamique de la charge adaptés aux techniques parallèles de remaillage ; conception d'algorithmes parallèles permettant de scinder le problème global de remaillage parallèle en plusieurs sous-tâches séquentielles, susceptibles de s'exécuter concurremment sur les processeurs de la machine parallèle. Ces contributions ont été mises en oeuvre au sein de la bibliothèque parallèle PaMPA, en s'appuyant sur les briques logicielles MMG3D (remaillage séquentiel de maillages tétraédriques) et PT-Scotch (repartitionnement parallèle de graphes). La bibliothèque PaMPA offre ainsi les fonctionnalités suivantes : communication transparente entre processeurs voisins des valeurs portées par les noeuds, les éléments, etc. ;remaillage, selon des critères fournis par l'utilisateur, de portions du maillage distribué, en offrant une qualité constante, que les éléments à remailler soient portés par un unique processeur ou bien répartis sur plusieurs d'entre eux ; répartition et redistribution de la charge des maillages pour préserver l'efficacité des simulations après remaillage. / The purpose of this thesis was to propose and to validate experimentally a set of algorithmic methods for the parallel remeshing of distributed meshes, based on a preexisting sequential remeshing method. This goal has been achieved through several steps : definition of data structures and of communication schemes suitable for distributed meshes, allowing for cheap migration of subdomain interfaces across the processors of a distributed-memory architecture ; use of dynamic load balancing algorithms suitable for parallel remeshing techniques ; design of parallel algorithms for splitting the global remeshing problem into several independent sequential tasks, susceptible to be executed concurrently across the processors of the parallel machine. These contributions have been implemented into the PaMPA parallel library, taking advantage of the MMG3D (sequential anisotropic tetrahedral remesher) PT-Scotch (parallel graph repartitioning) software. The PaMPA library consequently provides the following features : transparent communication across neighboring processors of data borne by nodes, elements, etc.; remeshing, according to used-defined criteria, of portions of the distributed mesh, that yields constant quality, irrespective of whether elements to be remeshed are located on a single processor or distributed across several of them ; balancing and redistribution of the workload of the mesh, to preserve the efficiency of simulations after the remeshing phase.

La simulation numérique est largement utilisée pour évaluer les performances aérodynamiques des aéronefs ainsi qu'en optimisation de forme. Ainsi l'objectif de ces simulations est souvent le calcul de fonctions aérodynamiques. L'objet de cette thèse est d'étudier des méthodes d'adaptation de maillages basées sur la dérivée totale de ces fonctions par rapport aux coordonnées du maillage (notée dJ/dX). Celle-ci pouvant être calculée par la méthode adjointe discrète. La première partie de cette étude concerne l'application de méthodes d'adaptation de maillages appliquées à des écoulements de fluides parfaits. Le senseur qui détecte les zones de maillage à raffiner s'appuie sur la norme de cette dérivée pour adapter des maillages pour le calcul d'une fonction J. La seconde partie du travail est la construction et l'étude de critères plus fiables basés sur dJ/dX pour d'une part adapter des maillages et d'autre part estimer si un maillage est bien adapté ou non pour le calcul de la fonction J. De plus une méthode de remaillage plus efficace basée sur une EDP elliptique est aussi présentée. Cette nouvelle méthode est appliquée pour des écoulements bidimensionnels de fluides parfaits ainsi que pour un écoulement décrit par les équations RANS. La dernière partie de l'étude est consacrée à l'application de la méthode proposée à des cas tridimensionnels d'écoulement RANS sur des géométries d'intérêt industriel. / Numerical simulation is widely used for the assessment of aircraft aerodynamic performances and shape optimizations. Hence the objective of these simulations is often to compute aerodynamic outputs. The purpose of this thesis is to study mesh adaptation methods based on the total derivative of the outputs with respect to mesh coordinates (denoted dJ/dX). This derivative can be computed using the discrete adjoint method. The first part of this study is about the application of mesh adaptation methods applied for Eulerian flows. The mesh locations to refine are detected using a sensor based on the norm of the derivative dJ/dX. This study confirmed that this derivative is relevant in order to adapt a mesh for the computation of the output J. The second part of this work is the construction and the study of reliable criteria based on dJ/dX for both mesh adaptation and the quality assessment of a given mesh for the computation of the output J. Moreover a more efficient remeshing method based on an elliptic PDE is presented too. This new method is applied for both two-dimensional Eulerian flows and a flow described by the RANS equations. The last part of the study is devoted to the application of the proposed method to three-dimensional RANS flows on geometries of industrial interest.

Construire une surface lisse d'un objet 3D est un problème important dans de nombreuses applications graphiques. En particulier, les méthodes de subdivision permettent de passer facilement d'un maillage discret à une surface continue. Un problème général résultant de la subdivision de deux maillages initialement connectés le long d'un bord est l'apparition de fissures ou de trous entre eux. Ces fissures produisent non seulement des formes indésirables, mais induisent aussi des difficultés pour les traitements ultérieurs. Il faut donc réparer ces défauts de sorte que la surface obtenue soit lisse et puisse être subdivisée ou modifiée. Nous proposons de nouvelles méthodes pour relier deux maillages avec des résolutions différentes en utilisant une transformée en ondelettes B-splines et une approximation locale ou une interpolation locale à l'aide de fonctions de base radiales (RBF). Ces procédés génèrent un maillage de connexion où la continuité est contrôlée. La résolution du maillage est ajustable pour respecter le changement de résolution entre les zones grossières et fines. En outre, nous présentons des méthodes pour combler les trous à n-côtés, et le raffinement des maillages grâce à un schéma de subdivision adaptative. Nous avons conçu, implémenté et testé les algorithmes en MatLab pour illustrer nos méthodes et montrer des résultats expérimentaux. Ces algorithmes sont mis en oeuvre sur de nombreux modèles d'objets 3D avec des formes complexes. En outre, nous avons fourni des approches différentes pour chaque problème. Ainsi, les résultats des différentes approches sont comparés et évalués afin d'exploiter les avantages et les inconvénients de ces approches. / Constructing a smooth surface of a 3D object is an important problem in many graphical applications. In particular, subdivision methods permit to pass easily from a discrete mesh to a continuous surface. A generic problem arising from subdividing two meshes initially connected along a common boundary is the occurrence of cracks or holes between them. These cracks not only produce undesired shapes, but also bring serious trouble for further mesh processing. They must be removed or filled so that the produced surface is smooth and can be further subdivided or edited. In order to remove cracks, we propose new methods for joining two meshes with different resolutions using a Lifted B-spline wavelet transform and a local approximation or radial basis function (RBF) local interpolation. These methods generate a connecting mesh where continuity is controlled from one boundary to the other and the connecting mesh can change gradually in resolution between coarse and fine areas. Additionally, we introduce methods for filling n-sided holes, and refining meshes with an adaptive subdivision scheme. We have designed, implemented, and tested the algorithms in MatLab to illustrate our proposed methods and show experimental results. These algorithms are implemented on many 3D object models with complex shapes. Additionally, we have provided some different approaches for each problem. Thus, results from the different approaches are compared and evaluated to exploit the advantages and disadvantages of these approaches.