Embedded Files
Digital Animation of Powder-Snow Avalanches
F. Nascimento, F. S. Sousa, and A. Paiva
ACM Transactions on Graphics (SIGGRAPH 2025), 44 (4), 2025
Project Page ● DOI
Abstract: Powder-snow avalanches are natural phenomena that result from an instability in the snow cover on a mountain relief. It begins with a dense avalanche core moving fast down the mountain. During its evolution, the snow particles in the avalanche front mix with the air, forming a suspended turbulent cloud of snow dust surrounding the dense snow avalanche. This paper introduces a physically-based framework using the Finite Volume Method to simulate powder-snow avalanches under complex terrains. Specifically, the primary goal is to simulate the turbulent snow cloud dynamics within the avalanche in a visually realistic manner. Our approach relies on a multi-layer model that splits the avalanche into two main layers: dense and powder-snow. The dense-snow layer flow is simulated by solving a type of Shallow Water Equations suited for intricate basal surfaces, known as the Savage-Hutter model. The powder-snow layer flow is modeled as a two-phase mixture of miscible fluids and simulated using Navier-Stokes equations. Moreover, we propose a novel model for the transition layer, which is responsible for coupling the avalanche main layers, including the snow mass injected into the powder-snow cloud from the snow entrainment processes and its injection velocity. In brief, our framework comprehensively simulates powder-snow avalanches, allowing us to render convincing animations of one of the most complex gravity-driven flows.
ACM Transactions on Graphics (SIGGRAPH 2025), 44 (4), 2025
Project Page ● DOI
Abstract: Powder-snow avalanches are natural phenomena that result from an instability in the snow cover on a mountain relief. It begins with a dense avalanche core moving fast down the mountain. During its evolution, the snow particles in the avalanche front mix with the air, forming a suspended turbulent cloud of snow dust surrounding the dense snow avalanche. This paper introduces a physically-based framework using the Finite Volume Method to simulate powder-snow avalanches under complex terrains. Specifically, the primary goal is to simulate the turbulent snow cloud dynamics within the avalanche in a visually realistic manner. Our approach relies on a multi-layer model that splits the avalanche into two main layers: dense and powder-snow. The dense-snow layer flow is simulated by solving a type of Shallow Water Equations suited for intricate basal surfaces, known as the Savage-Hutter model. The powder-snow layer flow is modeled as a two-phase mixture of miscible fluids and simulated using Navier-Stokes equations. Moreover, we propose a novel model for the transition layer, which is responsible for coupling the avalanche main layers, including the snow mass injected into the powder-snow cloud from the snow entrainment processes and its injection velocity. In brief, our framework comprehensively simulates powder-snow avalanches, allowing us to render convincing animations of one of the most complex gravity-driven flows.
Parallel isogeometric boundary element analysis with T-splines on CUDA
M. A. Peres, G. Sanches, A. Paiva and. P. Pagliosa
Computer Methods in Applied Mechanics and Engineering, 432, Part A, 2024
Paper (PDF) ● Code ● BIB ● DOI
Abstract: We present a framework for parallel isogeometric boundary element analysis (BEA) of elastic solids on CUDA. To deal with traction discontinuities, we propose a BEA model that supports multiple nodes and discontinuous elements. The multiplicity of a node is defined by the number of regions containing any element influenced by the node. A region is a group of connected elements delimited by a closed crease curve. The default shape function of an element is determined by a linear operator applied to a set of basis functions. A BEA model is supposed to be generated from a watertight boundary representation of a solid. In this paper, we employ bicubic analysis-suitable T-splines. In this case, the shape of an element is defined by its Bezier extraction operator applied to the tensor product of Bernstein polynomials of degree 3. We describe the data structures of the BEA model and the main algorithms of the analysis pipeline on CUDA. In particular, we describe two strategies for parallel assembling of the linear system of equations. We also introduce a novel approach for inside integration based on the subdivision of the singular region in triangles with constant aspect ratio. In the T-splines context, we extend the B ́ezier extraction to handle unstructured T-meshes with crease edges. Moreover, we propose a scheme for embedding the influence of linked tangency handles on the shape of an element directly into the Bezier extraction operator. Such an embedding enables the removal of the corresponding nodes from the BEA model and the application of an alternative collocation method we discuss in the paper. We present several experiments to evaluate the accuracy and efficiency of the proposed framework. The results demonstrate that the GPU can be advantageously employed for parallelizing T-spline-based isogeometric analysis using boundary elements, achieving a speedup of up to 29x compared to the sequential code on a current laptop. We make the BEA code available in a prototype in MATLAB with a graphical interface that allows users to apply boundary conditions and visualize analysis results on the boundary.
Computer Methods in Applied Mechanics and Engineering, 432, Part A, 2024
Paper (PDF) ● Code ● BIB ● DOI
Abstract: We present a framework for parallel isogeometric boundary element analysis (BEA) of elastic solids on CUDA. To deal with traction discontinuities, we propose a BEA model that supports multiple nodes and discontinuous elements. The multiplicity of a node is defined by the number of regions containing any element influenced by the node. A region is a group of connected elements delimited by a closed crease curve. The default shape function of an element is determined by a linear operator applied to a set of basis functions. A BEA model is supposed to be generated from a watertight boundary representation of a solid. In this paper, we employ bicubic analysis-suitable T-splines. In this case, the shape of an element is defined by its Bezier extraction operator applied to the tensor product of Bernstein polynomials of degree 3. We describe the data structures of the BEA model and the main algorithms of the analysis pipeline on CUDA. In particular, we describe two strategies for parallel assembling of the linear system of equations. We also introduce a novel approach for inside integration based on the subdivision of the singular region in triangles with constant aspect ratio. In the T-splines context, we extend the B ́ezier extraction to handle unstructured T-meshes with crease edges. Moreover, we propose a scheme for embedding the influence of linked tangency handles on the shape of an element directly into the Bezier extraction operator. Such an embedding enables the removal of the corresponding nodes from the BEA model and the application of an alternative collocation method we discuss in the paper. We present several experiments to evaluate the accuracy and efficiency of the proposed framework. The results demonstrate that the GPU can be advantageously employed for parallelizing T-spline-based isogeometric analysis using boundary elements, achieving a speedup of up to 29x compared to the sequential code on a current laptop. We make the BEA code available in a prototype in MATLAB with a graphical interface that allows users to apply boundary conditions and visualize analysis results on the boundary.
Surface reconstruction method for particle-based fluids using discrete indicator functions
Filomen Incahuanaco and Afonso Paiva
Computers & Graphics (SIBGRAPI 2022), 114, pp. 26-35, 2023
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: Particle-based methods are practical for computing level-sets that represent liquid interfaces. However, these methods are computationally expensive to reconstruct the liquid surface when the number of particles increases considerably due to the massive amount of particle approximations. This paper introduces two simple and efficient surface reconstruction methods for particle-based fluids based on discrete indicator functions (DIFs). The first approach provides fast level-set approximation using a DIF defined by counting particles inside grid cells. The second approach generates a high-quality liquid surface using a DIF obtained by the particle distribution inside grid cells. Our DIF-based methods are fast, easy to code, and can be adapted straightforwardly in particle-based fluid solvers, even implemented in GPU. Moreover, we show the effectiveness of our approaches through experiments against prior surface reconstruction methods. Extended version of SIBGRAPI 2022 paper.
Computers & Graphics (SIBGRAPI 2022), 114, pp. 26-35, 2023
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: Particle-based methods are practical for computing level-sets that represent liquid interfaces. However, these methods are computationally expensive to reconstruct the liquid surface when the number of particles increases considerably due to the massive amount of particle approximations. This paper introduces two simple and efficient surface reconstruction methods for particle-based fluids based on discrete indicator functions (DIFs). The first approach provides fast level-set approximation using a DIF defined by counting particles inside grid cells. The second approach generates a high-quality liquid surface using a DIF obtained by the particle distribution inside grid cells. Our DIF-based methods are fast, easy to code, and can be adapted straightforwardly in particle-based fluid solvers, even implemented in GPU. Moreover, we show the effectiveness of our approaches through experiments against prior surface reconstruction methods. Extended version of SIBGRAPI 2022 paper.
Counting Particles: a simple and fast surface reconstruction method for particle-based fluids
Filomen Incahuanaco Quispe and Afonso Paiva
Proceedings of SIBGRAPI 2022, pp. 145-149, 2022
Honorable Mention in Computer Graphics/Visualization
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: We present a novel and efficient surface reconstruction for particle-based fluid methods. Although particle-based methods are practical for computing level-sets that represent liquid interfaces, these methods are computationally expensive when the number of particles increases considerably due to the intense usage of particle approximations. This paper introduces a simple level-set approximation using a discrete indicator function (DIF) defined by counting particles inside grid cells. Our method is fast, easy to code, and can be adapted straightforwardly in particle-based solvers, even implemented in GPU. Moreover, we show the effectiveness of our approach through a set of experiments against prior surface reconstruction methods.
Proceedings of SIBGRAPI 2022, pp. 145-149, 2022
Honorable Mention in Computer Graphics/Visualization
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: We present a novel and efficient surface reconstruction for particle-based fluid methods. Although particle-based methods are practical for computing level-sets that represent liquid interfaces, these methods are computationally expensive when the number of particles increases considerably due to the intense usage of particle approximations. This paper introduces a simple level-set approximation using a discrete indicator function (DIF) defined by counting particles inside grid cells. Our method is fast, easy to code, and can be adapted straightforwardly in particle-based solvers, even implemented in GPU. Moreover, we show the effectiveness of our approach through a set of experiments against prior surface reconstruction methods.
Region reconstruction with the sphere-of-influence diagram
Luiz Henrique de Figueiredo and Afonso Paiva
Computers & Graphics (SIBGRAPI 2022), 107, pp. 252-263, 2022
Paper (PDF) ● BIB ● DOI
Abstract: We describe a simple method for reconstructing regions in the plane from well-distributed point samples. For that, we introduce the sphere-of-influence diagram, a planar diagram extracted from the Delaunay triangulation and the sphere-of-influence graph. The sphere-of-influence diagram is simple to understand and to implement, and supports an intuitive scaling parameter to handle variations in the distribution of samples. We report experiments reconstructing regions from point clouds of varying shape and density. We compare our results with those obtained by α-shapes and by CT-shapes.
Computers & Graphics (SIBGRAPI 2022), 107, pp. 252-263, 2022
Paper (PDF) ● BIB ● DOI
Abstract: We describe a simple method for reconstructing regions in the plane from well-distributed point samples. For that, we introduce the sphere-of-influence diagram, a planar diagram extracted from the Delaunay triangulation and the sphere-of-influence graph. The sphere-of-influence diagram is simple to understand and to implement, and supports an intuitive scaling parameter to handle variations in the distribution of samples. We report experiments reconstructing regions from point clouds of varying shape and density. We compare our results with those obtained by α-shapes and by CT-shapes.
Narrow-band screen-space fluid rendering
Felipe Oliveira and Afonso Paiva
Computer Graphics Forum, 41 (6), pp. 82-93, 2022
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: This paper presents a novel and practical screen-space liquid rendering for particle-based fluids for real-time applications. Our rendering pipeline performs particle filtering only in a narrow-band around the boundary particles to provide a smooth liquid surface with volumetric rendering effects. We also introduce a novel boundary detection method allowing the user to select particle layers from the liquid interface. The proposed approach is simple, fast, memory-efficient, easy to code, and it can be adapted straightforwardly in the standard screen-space rendering methods, even in GPU architectures. We show through a set of experiments how the prior screen space techniques can be benefited and improved by our approach.
Computer Graphics Forum, 41 (6), pp. 82-93, 2022
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: This paper presents a novel and practical screen-space liquid rendering for particle-based fluids for real-time applications. Our rendering pipeline performs particle filtering only in a narrow-band around the boundary particles to provide a smooth liquid surface with volumetric rendering effects. We also introduce a novel boundary detection method allowing the user to select particle layers from the liquid interface. The proposed approach is simple, fast, memory-efficient, easy to code, and it can be adapted straightforwardly in the standard screen-space rendering methods, even in GPU architectures. We show through a set of experiments how the prior screen space techniques can be benefited and improved by our approach.
Foreword to the special section on SIBGRAPI 2021
Afonso Paiva and Gladimir Baranoski
Computers & Graphics, 103, pp. A10-A11, 2022
Paper (PDF) ● BIB ● DOI
Abstract: Special section of Computers & Graphcis and selection of best papers on Graphics and Visualization from SIBGRAPI 2021.
Computers & Graphics, 103, pp. A10-A11, 2022
Paper (PDF) ● BIB ● DOI
Abstract: Special section of Computers & Graphcis and selection of best papers on Graphics and Visualization from SIBGRAPI 2021.
STFT-LDA: An algorithm to facilitate the visual analysis of building seismic responses
Z. Zhao, D. Motta, M. Berger, J. A Levine, I. B. Kuzucu, R. B. Fleischman, A. Paiva and C. Scheidegger
Information Visualization, 20 (4), pp. 263-282, 2021
Paper (PDF) ● BIB ● DOI
Abstract: Civil engineers use numerical simulations of a building’s responses to seismic forces to understand the nature of building failures, the limitations of building codes, and how to determine the latter to prevent the former. Such simulations generate large ensembles of multivariate, multiattribute time series. Comprehensive understanding of this data requires techniques that support the multivariate nature of the time series and can compare behaviors that are both periodic and non-periodic across multiple time scales and multiple time series themselves. In this paper, we present a novel technique to extract such patterns from time series generated from simulations of seismic responses. The core of our approach is the use of topic modeling, where topics correspond to interpretable and discriminative features of the earthquakes. We transform the raw time series data into a time series of topics, and use this visual summary to compare temporal patterns in earthquakes, query earthquakes via the topics across arbitrary time scales, and enable details on demand by linking the topic visualization with the original earthquake data. We show, through a surrogate task and an expert study, that this technique allows analysts to more easily identify recurring patterns in such time series. By integrating this technique in a prototype system, we show how it enables novel forms of visual interaction.
Information Visualization, 20 (4), pp. 263-282, 2021
Paper (PDF) ● BIB ● DOI
Abstract: Civil engineers use numerical simulations of a building’s responses to seismic forces to understand the nature of building failures, the limitations of building codes, and how to determine the latter to prevent the former. Such simulations generate large ensembles of multivariate, multiattribute time series. Comprehensive understanding of this data requires techniques that support the multivariate nature of the time series and can compare behaviors that are both periodic and non-periodic across multiple time scales and multiple time series themselves. In this paper, we present a novel technique to extract such patterns from time series generated from simulations of seismic responses. The core of our approach is the use of topic modeling, where topics correspond to interpretable and discriminative features of the earthquakes. We transform the raw time series data into a time series of topics, and use this visual summary to compare temporal patterns in earthquakes, query earthquakes via the topics across arbitrary time scales, and enable details on demand by linking the topic visualization with the original earthquake data. We show, through a surrogate task and an expert study, that this technique allows analysts to more easily identify recurring patterns in such time series. By integrating this technique in a prototype system, we show how it enables novel forms of visual interaction.
On novelty detection for multi-class classification using non-linear metric learning
S. R. Silva, T. Vieira, D. Martínez and A. Paiva
Expert Systems with Applications, 167, pp. 114193, 2021
Paper (PDF) ● BIB ● DOI
Abstract: Novelty detection is a binary task aimed at identifying whether a test sample is novel or unusual compared to a previously observed training set. A typical approach is to consider distance as a criterion to detect such novelties. However, most previous work does not focus on finding an optimum distance for each particular problem. In this paper, we propose to detect novelties by exploiting non-linear distances learned from multi-class training data. For this purpose, we adopt a kernelization technique jointly with the Large Margin Nearest Neighbor (LMNN) metric learning algorithm. The optimum distance tries to keep each known class’s instances together while pushing instances from different known classes to remain reasonably distant. We propose a variant of the K-Nearest Neighbors (KNN) classifier that employs the learned distance to detect novelties. Besides, we use the learned distance to perform multi-class classification. We show quantitative and qualitative experiments conducted on synthetic and real data sets, revealing that the learned metrics are effective in improving novelty detection compared to other metrics. Our method also outperforms previous work regularly used for novelty detection.
Expert Systems with Applications, 167, pp. 114193, 2021
Paper (PDF) ● BIB ● DOI
Abstract: Novelty detection is a binary task aimed at identifying whether a test sample is novel or unusual compared to a previously observed training set. A typical approach is to consider distance as a criterion to detect such novelties. However, most previous work does not focus on finding an optimum distance for each particular problem. In this paper, we propose to detect novelties by exploiting non-linear distances learned from multi-class training data. For this purpose, we adopt a kernelization technique jointly with the Large Margin Nearest Neighbor (LMNN) metric learning algorithm. The optimum distance tries to keep each known class’s instances together while pushing instances from different known classes to remain reasonably distant. We propose a variant of the K-Nearest Neighbors (KNN) classifier that employs the learned distance to detect novelties. Besides, we use the learned distance to perform multi-class classification. We show quantitative and qualitative experiments conducted on synthetic and real data sets, revealing that the learned metrics are effective in improving novelty detection compared to other metrics. Our method also outperforms previous work regularly used for novelty detection.
CrimAnalyzer: understanding crime patterns in São Paulo City
G. Garcia, J. Silveira, J. Poco, A. Paiva, M. Nery, C. T. Silva, S. Adorno and L. G. Nonato
IEEE Transactions on Visualization and Computer Graphics (presented at VIS 2020), 27 (4), pp. 2313-2328, 2021
Paper (PDF) ● BIB ● DOI
Abstract: São Paulo is the largest city in South America, with criminality rates as large as the city. The number and type of crimes vary considerably around the city, assuming different patterns depending on urban and social characteristics of each particular location. Previous works have mostly focused on the analysis of crimes with the intent of uncovering patterns associated to social factors, time seasonality, and urban routine activities. Therefore, those studies and tools are more global in the sense that they are not designed to investigate specific regions of the city such as particular neighborhoods, avenues, or public areas. Enabling tools able to explore particular locations of the city is essential for domain experts to accomplish their analysis in a bottom-up fashion, more clearly revealing how local urban features related to mobility, passersby behavior, and presence of public infrastructures such as terminals of public transportation and schools, can influence the quantity and type of crimes. In this paper, we present CrimAnalyzer; a visualization assisted analytic tool that allows users to analyze the behavior of crimes in specific regions of a city.
IEEE Transactions on Visualization and Computer Graphics (presented at VIS 2020), 27 (4), pp. 2313-2328, 2021
Paper (PDF) ● BIB ● DOI
Abstract: São Paulo is the largest city in South America, with criminality rates as large as the city. The number and type of crimes vary considerably around the city, assuming different patterns depending on urban and social characteristics of each particular location. Previous works have mostly focused on the analysis of crimes with the intent of uncovering patterns associated to social factors, time seasonality, and urban routine activities. Therefore, those studies and tools are more global in the sense that they are not designed to investigate specific regions of the city such as particular neighborhoods, avenues, or public areas. Enabling tools able to explore particular locations of the city is essential for domain experts to accomplish their analysis in a bottom-up fashion, more clearly revealing how local urban features related to mobility, passersby behavior, and presence of public infrastructures such as terminals of public transportation and schools, can influence the quantity and type of crimes. In this paper, we present CrimAnalyzer; a visualization assisted analytic tool that allows users to analyze the behavior of crimes in specific regions of a city.
RBF Liquids: an adaptive PIC solver using RBF-FD
R. Nakanishi*, F. Nascimento*, R. Campos, P. Pagliosa and A. Paiva (* joint 1st authors)
ACM Transactions on Graphics (SIGGRAPH Asia 2020), 39 (6), 2020
Project Page ● DOI
Abstract: We introduce a novel liquid simulation approach that combines a spatially adaptive pressure projection solver with the Particle-in-Cell (PIC) method. The solver relies on a generalized version of the Finite Difference (FD) method to approximate the pressure field and its gradients in tree-based grid discretizations, possibly non-graded. In our approach, FD stencils are computed by using meshfree interpolations provided by a variant of Radial Basis Function (RBF), known as RBF-Finite-Difference (RBF-FD). This meshfree version of the FD produces differentiation weights on scattered nodes with high-order accuracy. Our method adapts a quadtree/octree dynamically in a narrow-band around the liquid interface, providing an adaptive particle sampling for the PIC advection step. Furthermore, RBF affords an accurate scheme for velocity transfer between the grid and particles, keeping the system’s stability and avoiding numerical dissipation. We also present a data structure that connects the spatial subdivision of a quadtree/octree with the topology of its corresponding dual-graph. Our data structure makes the setup of stencils straightforward, allowing its updating without the need to rebuild it from scratch at each time-step. We show the effectiveness and accuracy of our solver by simulating incompressible inviscid fluids and comparing results with regular PIC-based solvers available in the literature.
ACM Transactions on Graphics (SIGGRAPH Asia 2020), 39 (6), 2020
Project Page ● DOI
Abstract: We introduce a novel liquid simulation approach that combines a spatially adaptive pressure projection solver with the Particle-in-Cell (PIC) method. The solver relies on a generalized version of the Finite Difference (FD) method to approximate the pressure field and its gradients in tree-based grid discretizations, possibly non-graded. In our approach, FD stencils are computed by using meshfree interpolations provided by a variant of Radial Basis Function (RBF), known as RBF-Finite-Difference (RBF-FD). This meshfree version of the FD produces differentiation weights on scattered nodes with high-order accuracy. Our method adapts a quadtree/octree dynamically in a narrow-band around the liquid interface, providing an adaptive particle sampling for the PIC advection step. Furthermore, RBF affords an accurate scheme for velocity transfer between the grid and particles, keeping the system’s stability and avoiding numerical dissipation. We also present a data structure that connects the spatial subdivision of a quadtree/octree with the topology of its corresponding dual-graph. Our data structure makes the setup of stencils straightforward, allowing its updating without the need to rebuild it from scratch at each time-step. We show the effectiveness and accuracy of our solver by simulating incompressible inviscid fluids and comparing results with regular PIC-based solvers available in the literature.
Simple and reliable boundary detection for meshfree particle methods using interval analysis
M. Sandim, A. Paiva and L. H. de Figueiredo
Journal of Computational Physics, 420, pp. 109702, 2020
Paper (PDF) ● BIB ● DOI
Abstract: We present two novel algorithms for detecting boundary particles in 2D and 3D domains that are suitable for meshfree particle methods in Computational Fluid Dynamics. We combine a robust purely geometric sphere covering test based on interval analysis with an adaptive spatial subdivision of the sphere associated with a given particle. The methods are simple, fast, and easy to code. We report comparisons against state-of-the-art boundary detection methods in free-surface flow problems to demonstrate the effectiveness and accuracy of our approaches.
Journal of Computational Physics, 420, pp. 109702, 2020
Paper (PDF) ● BIB ● DOI
Abstract: We present two novel algorithms for detecting boundary particles in 2D and 3D domains that are suitable for meshfree particle methods in Computational Fluid Dynamics. We combine a robust purely geometric sphere covering test based on interval analysis with an adaptive spatial subdivision of the sphere associated with a given particle. The methods are simple, fast, and easy to code. We report comparisons against state-of-the-art boundary detection methods in free-surface flow problems to demonstrate the effectiveness and accuracy of our approaches.
Boundary particle resampling for surface reconstruction in liquid animation
M. Sandim*, N. Oe*, D. Cedrim, P. Pagliosa and A. Paiva (* joint 1st authors)
Computers & Graphics (SIBGRAPI 2019), 84, pp. 55-65, 2019
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: In this paper, we present a novel adaptive particle resampling method tailored for surface reconstruction of level-sets defined by the boundary particles from a particle-based liquid simulation. The proposed approach is simple and easy to implement, and only requires the positions of the particles to identify and refine regions with small and thin fluid features accurately. The method comprises four main stages: boundary detection, feature classification, particle refinement, and surface reconstruction. For each simulation frame, firstly the free-surface particles are captured through a boundary detection method. Then, the boundary particles are classified and labeled according to the deformation and the stretching of the free-surface computed from the Principal Component Analysis (PCA) of the particle positions. The particles placed at feature regions are then refined according to their feature classification. Finally, we extract the free-surface of the zero level-set defined by the resampled boundary particles and its normals. In order to render the free-surface, we demonstrate how the traditional methods of surface fitting in Computer Graphics and Computational Physics literature can benefit from the proposed resampling method. Furthermore, the results shown in the paper attest the effectiveness and robustness of our method when compared to state-of-the-art adaptive particle resampling techniques.
Computers & Graphics (SIBGRAPI 2019), 84, pp. 55-65, 2019
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: In this paper, we present a novel adaptive particle resampling method tailored for surface reconstruction of level-sets defined by the boundary particles from a particle-based liquid simulation. The proposed approach is simple and easy to implement, and only requires the positions of the particles to identify and refine regions with small and thin fluid features accurately. The method comprises four main stages: boundary detection, feature classification, particle refinement, and surface reconstruction. For each simulation frame, firstly the free-surface particles are captured through a boundary detection method. Then, the boundary particles are classified and labeled according to the deformation and the stretching of the free-surface computed from the Principal Component Analysis (PCA) of the particle positions. The particles placed at feature regions are then refined according to their feature classification. Finally, we extract the free-surface of the zero level-set defined by the resampled boundary particles and its normals. In order to render the free-surface, we demonstrate how the traditional methods of surface fitting in Computer Graphics and Computational Physics literature can benefit from the proposed resampling method. Furthermore, the results shown in the paper attest the effectiveness and robustness of our method when compared to state-of-the-art adaptive particle resampling techniques.
Vessel optimal transport for automated alignment of retinal fundus images
D. Motta, W. Casaca, and A. Paiva
IEEE Transactions on Image Processing, 28 (12), pp. 6154-6168, 2019
Paper (PDF) ● BIB ● DOI
Abstract: Optimal transport has emerged as a promising and useful tool for supporting modern image processing applications such as medical imaging and scientific visualization. Indeed, the optimal transport theory enables great flexibility in modeling problems related to image registration, as different optimization resources can be successfully used as well as the choice of suitable matching models to align the images. In this paper, we introduce an automated framework for fundus image registration which unifies optimal transport theory, image processing tools and graph matching schemes into a functional and concise methodology. Given two ocular fundus images, we construct representative graphs which embed in their structures spatial and topological information from the eye’s blood vessels. The graphs produced are then used as input by our optimal transport model in order to establish a correspondence between their sets of nodes. Finally, geometric transformations are performed between the images so as to accomplish the registration task properly.
IEEE Transactions on Image Processing, 28 (12), pp. 6154-6168, 2019
Paper (PDF) ● BIB ● DOI
Abstract: Optimal transport has emerged as a promising and useful tool for supporting modern image processing applications such as medical imaging and scientific visualization. Indeed, the optimal transport theory enables great flexibility in modeling problems related to image registration, as different optimization resources can be successfully used as well as the choice of suitable matching models to align the images. In this paper, we introduce an automated framework for fundus image registration which unifies optimal transport theory, image processing tools and graph matching schemes into a functional and concise methodology. Given two ocular fundus images, we construct representative graphs which embed in their structures spatial and topological information from the eye’s blood vessels. The graphs produced are then used as input by our optimal transport model in order to establish a correspondence between their sets of nodes. Finally, geometric transformations are performed between the images so as to accomplish the registration task properly.
Fundus image transformation revisited: towards determining more accurate registrations
D. Motta, W. Casaca, and A. Paiva
Proceedings of IEEE CBMS 2018, pp. 227-232, 2018
Paper (PDF) ● Code (Matlab) ● BIB ● DOI
Abstract: Image registration is an important pre-processing step in several computer vision applications, being crucial in medical imaging systems where patients are examined and diagnosed almost exclusively by images. Historically, geometric transformations derived from quadratic models have been widely used as a benchmark to perform registration on fundus images, but in this paper, we demonstrate that quadratic and other high-order mappings are not necessarily the best choices for this purpose, even for well-established state-of-the-art registration methods. From a novel overlapping metric designed to determine the best image transformation that maximizes the registration accuracy, we improve the assertiveness of several methods of the literature while still preserving the same computational burden initially reached by those methods.
Proceedings of IEEE CBMS 2018, pp. 227-232, 2018
Paper (PDF) ● Code (Matlab) ● BIB ● DOI
Abstract: Image registration is an important pre-processing step in several computer vision applications, being crucial in medical imaging systems where patients are examined and diagnosed almost exclusively by images. Historically, geometric transformations derived from quadratic models have been widely used as a benchmark to perform registration on fundus images, but in this paper, we demonstrate that quadratic and other high-order mappings are not necessarily the best choices for this purpose, even for well-established state-of-the-art registration methods. From a novel overlapping metric designed to determine the best image transformation that maximizes the registration accuracy, we improve the assertiveness of several methods of the literature while still preserving the same computational burden initially reached by those methods.
A user-friendly interactive framework for unsteady fluid flow segmentation and visualization
D. Motta, W. Casaca, P. Pagliosa and A. Paiva
Journal of Visualization, 21 (4), pp. 625-636, 2018
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: While vector fields are essential to simulate a great amount of natural phenomena, the difficulty to identify patterns and predict behaviors makes the visual segmentation in simulations an attractive and powerful tool. In this paper, we present a novel user-steered segmentation framework to cope with steady as well as unsteady vector fields on fluid flow simulations. Given a discrete vector field, our approach extracts multi-valued features from the field by exploiting its streamline structures so that these features are mapped to a visual space through a multidimensional projection technique. From an easy-to-handle interface, the user can interact with the projected data so as to partition and explore the most relevant vector features in a guidance frame of the simulation. Besides navigating and visually mining structures of interest, the interactivity with the projected data also allows us to progressively enhance the segmentation result according to the user insights. Finally, to successfully deal with unsteady simulations, the segments previously annotated by the user are used as a training set for a Support Vector Machine (SVM) approach that classifies the remaining frames in the flow. We attest the effectiveness and versatility of our methodology throughout a set of classical physical-inspired applications on fluid flow simulations as depicted in the experiment results section.
Journal of Visualization, 21 (4), pp. 625-636, 2018
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: While vector fields are essential to simulate a great amount of natural phenomena, the difficulty to identify patterns and predict behaviors makes the visual segmentation in simulations an attractive and powerful tool. In this paper, we present a novel user-steered segmentation framework to cope with steady as well as unsteady vector fields on fluid flow simulations. Given a discrete vector field, our approach extracts multi-valued features from the field by exploiting its streamline structures so that these features are mapped to a visual space through a multidimensional projection technique. From an easy-to-handle interface, the user can interact with the projected data so as to partition and explore the most relevant vector features in a guidance frame of the simulation. Besides navigating and visually mining structures of interest, the interactivity with the projected data also allows us to progressively enhance the segmentation result according to the user insights. Finally, to successfully deal with unsteady simulations, the segments previously annotated by the user are used as a training set for a Support Vector Machine (SVM) approach that classifies the remaining frames in the flow. We attest the effectiveness and versatility of our methodology throughout a set of classical physical-inspired applications on fluid flow simulations as depicted in the experiment results section.
Least-squares morphing of dynamic meshes
A. Medalha, L. Pagliosa, A. Paiva and P. Pagliosa
Proceedings of SIBGRAPI 2017, pp. 23-30, 2017
Best Paper Award in Computer Graphics/Visualization
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: Morphing is a technique that smoothly transforms a shape onto another. In this paper, we present a method for morphing of two dynamic meshes: mesh sequences representing the keyframes of animated shapes over time. The pipeline of the proposed method comprises two main stages: template-based cross-parameterization and dynamic mesh interpolation. In the cross-parameterization stage, we use a variation of least-squares (LS) meshes to provide a coarse approximation of the geometry of the source mesh onto the target mesh. In our method, the possible candidates for initial control points of the LS-mesh are detected using an approach based on the Heat Kernel Signature (HKS). Then, an iterative process of fine fitting adds new constraints in the LS-mesh processing. The cross-parameterization is performed just once for any two frames in order to establish a full correspondence between vertices of the source and target meshes. Next, we use such a correspondence in the dynamic mesh interpolation stage to produce the morphing results. The method is entirely mesh-based and does not demand the generation of skeletons, mesh segmentation or the use of any additional control structures.
Proceedings of SIBGRAPI 2017, pp. 23-30, 2017
Best Paper Award in Computer Graphics/Visualization
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: Morphing is a technique that smoothly transforms a shape onto another. In this paper, we present a method for morphing of two dynamic meshes: mesh sequences representing the keyframes of animated shapes over time. The pipeline of the proposed method comprises two main stages: template-based cross-parameterization and dynamic mesh interpolation. In the cross-parameterization stage, we use a variation of least-squares (LS) meshes to provide a coarse approximation of the geometry of the source mesh onto the target mesh. In our method, the possible candidates for initial control points of the LS-mesh are detected using an approach based on the Heat Kernel Signature (HKS). Then, an iterative process of fine fitting adds new constraints in the LS-mesh processing. The cross-parameterization is performed just once for any two frames in order to establish a full correspondence between vertices of the source and target meshes. Next, we use such a correspondence in the dynamic mesh interpolation stage to produce the morphing results. The method is entirely mesh-based and does not demand the generation of skeletons, mesh segmentation or the use of any additional control structures.
Boundary detection in particle-based fluids
M. Sandim, D. Cedrim, L. G. Nonato, P. Pagliosa and A. Paiva
Computer Graphics Forum (Eurographics 2016), 35 (2), pp. 215-224, 2016
Paper (PDF) ● Video (YouTube) ● Code (Python) ● BIB ● DOI
Abstract: This paper presents a novel method to detect free-surfaces on particle-based volume representation. In contrast to most particle-based free-surface detection methods, which perform the surface identification based on physical and geometrical properties derived from the underlying fluid flow simulation, the proposed approach only demands the spatial location of the particles to properly recognize surface particles, avoiding even the use of kernels. Boundary particles are identified through a Hidden Point Removal (HPR) operator used for visibility test. Our method is very simple, fast, easy to implement and robust to changes in the distribution of particles, even when facing large deformation of the free-surface. A set of comparisons against state-of-the-art boundary detection methods show the effectiveness of our approach. The good performance of our method is also attested in the context of fluid flow simulation involving free-surface, mainly when using level-sets for rendering purposes.
Computer Graphics Forum (Eurographics 2016), 35 (2), pp. 215-224, 2016
Paper (PDF) ● Video (YouTube) ● Code (Python) ● BIB ● DOI
Abstract: This paper presents a novel method to detect free-surfaces on particle-based volume representation. In contrast to most particle-based free-surface detection methods, which perform the surface identification based on physical and geometrical properties derived from the underlying fluid flow simulation, the proposed approach only demands the spatial location of the particles to properly recognize surface particles, avoiding even the use of kernels. Boundary particles are identified through a Hidden Point Removal (HPR) operator used for visibility test. Our method is very simple, fast, easy to implement and robust to changes in the distribution of particles, even when facing large deformation of the free-surface. A set of comparisons against state-of-the-art boundary detection methods show the effectiveness of our approach. The good performance of our method is also attested in the context of fluid flow simulation involving free-surface, mainly when using level-sets for rendering purposes.
Depth functions as a quality measure and for steering multidimensional projections
D. Cedrim, V. Vad, A. Paiva, E. Gröller, L. G. Nonato and A. Castelo
Computers & Graphics (SIBGRAPI 2016), 60, pp. 93-106, 2016
Paper (PDF) ● BIB ● DOI
Abstract: The analysis of multidimensional data has been for many years a topic of continuous research. This type of data can be found in several different areas of science. A common task analyzing such data is to investigate patterns by interacting with spatializations of the data in a visual domain. Understanding the relation between underlying dataset characteristics and the technique used to provide a visual representation of such dataset is of fundamental importance since it can provide a better intuition on what to expect from the spatialization. Motivated by this, in this paper we propose the usage of concepts from non-parametric statistics, namely depth functions, as quality measures for spatializations. We aim to evaluate the action of multidimensional projection (MP) techniques on such an estimate. In order to investigate this, we apply both qualitative and quantitative analyses on four different multidimensional techniques selected according to the properties they aim to preserve. We evaluate them with datasets of different characteristics: synthetic, real world, high dimensional and contaminated with outliers. As a straightforward application we propose to use depth information to guide MP techniques which rely on an intervention through the control point selection and positioning, showing that even for a technique which does no intend to preserve any centrality measure, interesting results can be achieved by separating regions possible contaminated with outliers.
Computers & Graphics (SIBGRAPI 2016), 60, pp. 93-106, 2016
Paper (PDF) ● BIB ● DOI
Abstract: The analysis of multidimensional data has been for many years a topic of continuous research. This type of data can be found in several different areas of science. A common task analyzing such data is to investigate patterns by interacting with spatializations of the data in a visual domain. Understanding the relation between underlying dataset characteristics and the technique used to provide a visual representation of such dataset is of fundamental importance since it can provide a better intuition on what to expect from the spatialization. Motivated by this, in this paper we propose the usage of concepts from non-parametric statistics, namely depth functions, as quality measures for spatializations. We aim to evaluate the action of multidimensional projection (MP) techniques on such an estimate. In order to investigate this, we apply both qualitative and quantitative analyses on four different multidimensional techniques selected according to the properties they aim to preserve. We evaluate them with datasets of different characteristics: synthetic, real world, high dimensional and contaminated with outliers. As a straightforward application we propose to use depth information to guide MP techniques which rely on an intervention through the control point selection and positioning, showing that even for a technique which does no intend to preserve any centrality measure, interesting results can be achieved by separating regions possible contaminated with outliers.
Visualizing and interacting with kernelized data
A. Barbosa, F. Paulovich, A. Paiva, S. Goldenstein, F. Petronetto and L. G. Nonato
IEEE Transactions on Visualization and Computer Graphics, 22 (3), pp. 1314-1325, 2016
Paper (PDF) ● BIB ● DOI
Abstract: Kernel-based methods have experienced a substantial progress in the last years, tuning out an essential mechanism for data classification, clustering and pattern recognition. The effectiveness of kernel-based techniques, though, depends largely on the capability of the underlying kernel to properly embed data in the feature space associated to the kernel. However, visualizing how a kernel embeds the data in a feature space is not so straightforward, as the embedding map and the feature space are implicitly defined by the kernel. In this work, we present a novel technique to visualize the action of a kernel, that is, how the kernel embeds data into a high-dimensional feature space. The proposed methodology relies on a solid mathematical formulation to map kernelized data onto a visual space. Our approach is faster and more accurate than most existing methods while still allowing interactive manipulation of the projection layout, a game-changing trait tha t other kernel-based projection techniques do not have.
IEEE Transactions on Visualization and Computer Graphics, 22 (3), pp. 1314-1325, 2016
Paper (PDF) ● BIB ● DOI
Abstract: Kernel-based methods have experienced a substantial progress in the last years, tuning out an essential mechanism for data classification, clustering and pattern recognition. The effectiveness of kernel-based techniques, though, depends largely on the capability of the underlying kernel to properly embed data in the feature space associated to the kernel. However, visualizing how a kernel embeds the data in a feature space is not so straightforward, as the embedding map and the feature space are implicitly defined by the kernel. In this work, we present a novel technique to visualize the action of a kernel, that is, how the kernel embeds data into a high-dimensional feature space. The proposed methodology relies on a solid mathematical formulation to map kernelized data onto a visual space. Our approach is faster and more accurate than most existing methods while still allowing interactive manipulation of the projection layout, a game-changing trait tha t other kernel-based projection techniques do not have.
Partial similarity of 3D shapes using cross recurrence plot
R. Nakanishi, J. P. Ono, P. Pagliosa, L. G. Nonato and A. Paiva
Proceedings of SIBGRAPI 2016, pp. 448-454, 2016
Paper (PDF) ● BIB ● DOI
Abstract: This paper presents a novel 3D partial shape retrieval algorithm based on time-series analysis. Given a piece of a 3D shape, the proposed method encodes the shape descriptor given by the Heat Kernel Signature (HKS) as a time-series, where the time is considered an ordered sequence of vertices provided by the Fiedler vector. Finally, a similarity metric is created using a well-known tool in time-series analysis called Cross Recurrence Plot (CRP). The good performance of our method is also attested in a large collection of shape models
Proceedings of SIBGRAPI 2016, pp. 448-454, 2016
Paper (PDF) ● BIB ● DOI
Abstract: This paper presents a novel 3D partial shape retrieval algorithm based on time-series analysis. Given a piece of a 3D shape, the proposed method encodes the shape descriptor given by the Heat Kernel Signature (HKS) as a time-series, where the time is considered an ordered sequence of vertices provided by the Fiedler vector. Finally, a similarity metric is created using a well-known tool in time-series analysis called Cross Recurrence Plot (CRP). The good performance of our method is also attested in a large collection of shape models
Particle-based fluids for viscous jet buckling
L.F.S. Andrade, M. Sandim, F. Petronetto, P. Pagliosa and A. Paiva
Computers & Graphics, 52, pp. 106-115, 2015
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: In this paper, we introduce a novel meshfree framework for animating free surface viscous liquids with jet buckling effects, such as coiling and folding. Our method is based on Smoothed Particle Hydrodynamics (SPH) fluids and allows more realistic and complex viscous behaviors than the previous SPH frameworks in computer animation literature. The viscous liquid is modeled by a non-Newtonian fluid flow and the variable viscosity under shear stress is achieved using a viscosity model known as Cross model. We demonstrate the efficiency and stability of our framework in a wide variety of animations, including scenarios with arbitrary geometries and high resolution of SPH particles. The interaction of the viscous liquid with complex solid obstacles is performed using boundary particles. Our framework is able to deal with different inlet velocity profiles and geometries of the injector, as well as moving inlet jet along trajectories given by cubic Hermite splines. Moreover, the simulation speed is significantly accelerated by using Computer Unified Device Architecture (CUDA) computing platform. Extended version of SIBGRAPI 2014 paper.
Computers & Graphics, 52, pp. 106-115, 2015
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: In this paper, we introduce a novel meshfree framework for animating free surface viscous liquids with jet buckling effects, such as coiling and folding. Our method is based on Smoothed Particle Hydrodynamics (SPH) fluids and allows more realistic and complex viscous behaviors than the previous SPH frameworks in computer animation literature. The viscous liquid is modeled by a non-Newtonian fluid flow and the variable viscosity under shear stress is achieved using a viscosity model known as Cross model. We demonstrate the efficiency and stability of our framework in a wide variety of animations, including scenarios with arbitrary geometries and high resolution of SPH particles. The interaction of the viscous liquid with complex solid obstacles is performed using boundary particles. Our framework is able to deal with different inlet velocity profiles and geometries of the injector, as well as moving inlet jet along trajectories given by cubic Hermite splines. Moreover, the simulation speed is significantly accelerated by using Computer Unified Device Architecture (CUDA) computing platform. Extended version of SIBGRAPI 2014 paper.
Concentric Radviz: visual exploration of multi-task classification
J. P. Ono, F. Sikansi, D. C. Correa, F. Paulovich, A. Paiva and L. G. Nonato
Proceedings of SIBGRAPI 2015, pp. 165-172, 2015
Honorable Mention in Computer Graphics/Visualization
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: The discovery of patterns in large data collections is a difficult task. Visualization and machine learning techniques have emerged as a way to facilitate data analysis, providing tools to uncover relevant patterns from the data. This paper presents Concentric Radviz, a general purpose class visualization system that takes into account multi-class, multi-label and multi-task classifiers. Concentric Radviz uses a force attenuation scheme, which minimizes cluttering and ambiguity in the visual layout. In addition, the user can add concentric circles to the layout in order to represent classification tasks. Our validation results and the application of Concentric Radviz for two real collections suggest that this tool can reveal important data patterns and relations. In our application, the user can interact with the visualization by selecting regions of interest according to specific criteria and changing projection parameters.
Proceedings of SIBGRAPI 2015, pp. 165-172, 2015
Honorable Mention in Computer Graphics/Visualization
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: The discovery of patterns in large data collections is a difficult task. Visualization and machine learning techniques have emerged as a way to facilitate data analysis, providing tools to uncover relevant patterns from the data. This paper presents Concentric Radviz, a general purpose class visualization system that takes into account multi-class, multi-label and multi-task classifiers. Concentric Radviz uses a force attenuation scheme, which minimizes cluttering and ambiguity in the visual layout. In addition, the user can add concentric circles to the layout in order to represent classification tasks. Our validation results and the application of Concentric Radviz for two real collections suggest that this tool can reveal important data patterns and relations. In our application, the user can interact with the visualization by selecting regions of interest according to specific criteria and changing projection parameters.
Exploratory segmentation of vector fields using multidimensional projection
D. Motta, M. L. Oliveira, P. Pagliosa, L. G. Nonato and A. Paiva
Proceedings of SIBGRAPI 2015, pp. 250-256, 2015
Paper (PDF) ● BIB ● DOI
Abstract: The difficulty to understand the complex behavior of vector fields makes its visual segmentation an area of constant interest in scientific visualization. In this paper, we present a novel interactive segmentation framework for discrete vector fields. In our method, the vector field domain is partitioned into multiple regions with same flow patterns. In order to accomplish this task, feature vectors are extracted from streamlines and mapped to a visual space using multidimensional projection. The interactivity with projected data in the visual space improving the results of the segmentation according to user’s knowledge. The provided results and comparisons show the flexibility and effectiveness of our framework.
Proceedings of SIBGRAPI 2015, pp. 250-256, 2015
Paper (PDF) ● BIB ● DOI
Abstract: The difficulty to understand the complex behavior of vector fields makes its visual segmentation an area of constant interest in scientific visualization. In this paper, we present a novel interactive segmentation framework for discrete vector fields. In our method, the vector field domain is partitioned into multiple regions with same flow patterns. In order to accomplish this task, feature vectors are extracted from streamlines and mapped to a visual space using multidimensional projection. The interactivity with projected data in the visual space improving the results of the segmentation according to user’s knowledge. The provided results and comparisons show the flexibility and effectiveness of our framework.
SHREC’15 Track: Non-rigid 3D Shape Retrieval
R. Nakanishi, A. Paiva, L. G. Nonato, et al.
Proceedings of EG Workshop on 3D Object Retrieval (3DOR), 2015
Paper (PDF) ● BIB ● DOI
Abstract: Non-rigid 3D shape retrieval has become a research hotpot in communities of computer graphics, computer vision, pattern recognition, etc. In this paper, we present the results of the SHREC’15 Track: Non-rigid 3D Shape Retrieval. The aim of this track is to provide a fair and effective platform to evaluate and compare the performance of current non-rigid 3D shape retrieval methods developed by different research groups around the world. The database utilized in this track consists of 1200 3D watertight triangle meshes which are equally classified into 50 categories. All models in the same category are generated from an original 3D mesh by implementing various pose transformations. The retrieval performance of a method is evaluated using 6 commonly-used measures (i.e., PR-plot, NN, FT, ST, E-measure and DCG.). Totally, there are 37 submissions and 11 groups taking part in this track. Evaluation results and comparison analyses described in this paper not only show the bright future in researches of non-rigid 3D shape retrieval but also point out several promising research directions in this topic.
Proceedings of EG Workshop on 3D Object Retrieval (3DOR), 2015
Paper (PDF) ● BIB ● DOI
Abstract: Non-rigid 3D shape retrieval has become a research hotpot in communities of computer graphics, computer vision, pattern recognition, etc. In this paper, we present the results of the SHREC’15 Track: Non-rigid 3D Shape Retrieval. The aim of this track is to provide a fair and effective platform to evaluate and compare the performance of current non-rigid 3D shape retrieval methods developed by different research groups around the world. The database utilized in this track consists of 1200 3D watertight triangle meshes which are equally classified into 50 categories. All models in the same category are generated from an original 3D mesh by implementing various pose transformations. The retrieval performance of a method is evaluated using 6 commonly-used measures (i.e., PR-plot, NN, FT, ST, E-measure and DCG.). Totally, there are 37 submissions and 11 groups taking part in this track. Evaluation results and comparison analyses described in this paper not only show the bright future in researches of non-rigid 3D shape retrieval but also point out several promising research directions in this topic.
All-in-focus imaging technique used to improve 3d retinal fundus image reconstruction
D. Motta, L. Matos, A. C. Souza, R. Marcato, A. Paiva and L. A. V. Carvalho
Proceedings of ACM SAC 2015, pp. 26-31, 2015.
Paper (PDF) ● BIB ● DOI
Abstract: In this paper we have applied the stacking technique on images from a medical device. There is an urgent need in ophthalmology for more cost effective instrumentation for the early diagnosis of glaucoma, one of the leading diseases in the cause of blindness worldwide. The current techniques involve expensive optical equipments generally called fundus camera, which most of the time capture a single high resolution frame for one eye at a time. In this research we have used stereoscopic videos of a state-of-the-art 3D retinal camera, which has simpler optics and electronics when compared to current monoscopic models. Nevertheless, the cross correlation algorithms for depth computation are very sensible to image noise and out of focus regions. We demonstrate the efficiency of our technique on experiments involving a sequence of images extracted from videos in simulations of optic nerves using artificial objects.
Proceedings of ACM SAC 2015, pp. 26-31, 2015.
Paper (PDF) ● BIB ● DOI
Abstract: In this paper we have applied the stacking technique on images from a medical device. There is an urgent need in ophthalmology for more cost effective instrumentation for the early diagnosis of glaucoma, one of the leading diseases in the cause of blindness worldwide. The current techniques involve expensive optical equipments generally called fundus camera, which most of the time capture a single high resolution frame for one eye at a time. In this research we have used stereoscopic videos of a state-of-the-art 3D retinal camera, which has simpler optics and electronics when compared to current monoscopic models. Nevertheless, the cross correlation algorithms for depth computation are very sensible to image noise and out of focus regions. We demonstrate the efficiency of our technique on experiments involving a sequence of images extracted from videos in simulations of optic nerves using artificial objects.
SPH fluids for viscous jet buckling
L.F.S. Andrade, M. Sandim, F. Petronetto, P. Pagliosa and A. Paiva
Proceedings of SIBGRAPI 2014, pp. 65-72, 2014
Best Paper Award in Computer Graphics
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: We present a novel meshfree technique for animating free surface viscous liquids with jet buckling effects, such as coiling and folding. Our technique is based on Smoothed Particle Hydrodynamics (SPH) fluids and allows more realistic and complex viscous behaviors than the preceding SPH frameworks in computer animation literature. The viscous liquid is modeled by a non-Newtonian fluid flow and the variable viscosity under shear stress is achieved using a viscosity model known as Cross model. The proposed technique is efficient and stable, and our framework can animate scenarios with high resolution of SPH particles in which the simulation speed is significantly accelerated by using Computer Unified Device Architecture (CUDA) computing platform. This work also includes several examples that demonstrate the ability of our technique.
Proceedings of SIBGRAPI 2014, pp. 65-72, 2014
Best Paper Award in Computer Graphics
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: We present a novel meshfree technique for animating free surface viscous liquids with jet buckling effects, such as coiling and folding. Our technique is based on Smoothed Particle Hydrodynamics (SPH) fluids and allows more realistic and complex viscous behaviors than the preceding SPH frameworks in computer animation literature. The viscous liquid is modeled by a non-Newtonian fluid flow and the variable viscosity under shear stress is achieved using a viscosity model known as Cross model. The proposed technique is efficient and stable, and our framework can animate scenarios with high resolution of SPH particles in which the simulation speed is significantly accelerated by using Computer Unified Device Architecture (CUDA) computing platform. This work also includes several examples that demonstrate the ability of our technique.
Approximating implicit curves on plane and surface triangulations with affine arithmetic
F. Nascimento, A. Paiva, L. H. de Figueiredo and J. Stolfi
Computers & Graphics, 40, pp. 36-48, 2014
Paper (PDF) ● BIB ● DOI
Abstract: We present a spatially and geometrically adaptive method for computing a robust polygonal approximation of an implicit curve defined on a planar region or on a triangulated surface. Our method uses affine arithmetic to identify regions where the curve lies inside a thin strip. Unlike other interval methods, even those based on affine arithmetic, our method works on both rectangular and triangular decompositions and can use any refinement scheme that the decomposition offers. Extended version of SIBGRAPI 2012 paper.
Computers & Graphics, 40, pp. 36-48, 2014
Paper (PDF) ● BIB ● DOI
Abstract: We present a spatially and geometrically adaptive method for computing a robust polygonal approximation of an implicit curve defined on a planar region or on a triangulated surface. Our method uses affine arithmetic to identify regions where the curve lies inside a thin strip. Unlike other interval methods, even those based on affine arithmetic, our method works on both rectangular and triangular decompositions and can use any refinement scheme that the decomposition offers. Extended version of SIBGRAPI 2012 paper.
MIST: multiscale information and summaries of texts
P. Pagliosa, R. M. Martins, D. Cedrim, A. Paiva, R. Minghim and L. G. Nonato
Proceedings of SIBGRAPI 2013, pp. 91-98, 2013
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: Combining distinct visual metaphors has been the mechanism adopted by several systems to enable the simultaneous visualization of multiple levels of information in a single layout. However, providing a meaningful layout while avoiding visual clutter is still a challenge. In this work we combine word clouds and a rigid-body simulation engine into an intuitive visualization tool that allows a user to visualize and interact with the content of document collections using a single overlap-free layout. The proposed force scheme ensures that neighboring documents are kept close to each other during and after layout change. Each group of neighboring documents formed on the layout generates a word cloud. A multi-seeded procedure guarantees a harmonious arrangement of distinct word clouds in visual space. The visual metaphor employs lenses to represent document instances where the size of each lens defines the importance of the document in the collection. To keep the visualization clean and intuitive, only the most relevant documents are depicted as lens while the remaining ones are either displayed as smaller glyphs to help convey density information or simply removed from the layout.
Proceedings of SIBGRAPI 2013, pp. 91-98, 2013
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: Combining distinct visual metaphors has been the mechanism adopted by several systems to enable the simultaneous visualization of multiple levels of information in a single layout. However, providing a meaningful layout while avoiding visual clutter is still a challenge. In this work we combine word clouds and a rigid-body simulation engine into an intuitive visualization tool that allows a user to visualize and interact with the content of document collections using a single overlap-free layout. The proposed force scheme ensures that neighboring documents are kept close to each other during and after layout change. Each group of neighboring documents formed on the layout generates a word cloud. A multi-seeded procedure guarantees a harmonious arrangement of distinct word clouds in visual space. The visual metaphor employs lenses to represent document instances where the size of each lens defines the importance of the document in the collection. To keep the visualization clean and intuitive, only the most relevant documents are depicted as lens while the remaining ones are either displayed as smaller glyphs to help convey density information or simply removed from the layout.
Normal correction towards smoothing point-based surfaces
P. Valdivia, D. Cedrim, F. Petronetto, A. Paiva and L. G. Nonato
Proceedings of SIBGRAPI 2013, pp. 187-194, 2013
Paper (PDF) ● BIB ● DOI
Abstract: Some of the main approaches for mesh denoising uses a two-step scheme, by filtering normals to later perform a point updating step to match the corrected normals. In this paper, we propose an adaptation of a two-step approach for point set surfaces, exploring three different weight schemes. Moreover, we explore three techniques for normal estimation, and how it improve the quality of filtering process. Towards a quantitative analysis, in addition of usual visual comparison, we evaluate the effectiveness of each combination by three measures and compare them against some state-of-art point-based denoising techniques.
Proceedings of SIBGRAPI 2013, pp. 187-194, 2013
Paper (PDF) ● BIB ● DOI
Abstract: Some of the main approaches for mesh denoising uses a two-step scheme, by filtering normals to later perform a point updating step to match the corrected normals. In this paper, we propose an adaptation of a two-step approach for point set surfaces, exploring three different weight schemes. Moreover, we explore three techniques for normal estimation, and how it improve the quality of filtering process. Towards a quantitative analysis, in addition of usual visual comparison, we evaluate the effectiveness of each combination by three measures and compare them against some state-of-art point-based denoising techniques.
Mesh-free discrete Laplace-Beltrami operator
F. Petronetto, A. Paiva, E. S. Helou, D. E. Stewart and L. G. Nonato
Computer Graphics Forum, 32 (6), pp. 214-226, 2013
Paper (PDF) ● BIB ● DOI
Abstract: In this work we propose a new discretization method for the Laplace-Beltrami operator defined on point-based surfaces. In contrast to the existing point-based discretization techniques, our approach does not rely on any triangle mesh structure, turning out truly meshfree. Based on a combination of Smoothed Particle Hydrodynamics and an optimization procedure to estimate area elements, our discretization method results in accurate solutions while still being robust when facing abrupt changes in the density of points. Moreover, the proposed scheme results in numerically stable discrete operators. The effectiveness of the proposed technique is brought to bear in many practical applications. In particular, we use the eigenstructure of the discrete operator for filtering and shape segmentation. Point-based surface deformation is another application that can be easily carried out from the proposed discretization method.
Computer Graphics Forum, 32 (6), pp. 214-226, 2013
Paper (PDF) ● BIB ● DOI
Abstract: In this work we propose a new discretization method for the Laplace-Beltrami operator defined on point-based surfaces. In contrast to the existing point-based discretization techniques, our approach does not rely on any triangle mesh structure, turning out truly meshfree. Based on a combination of Smoothed Particle Hydrodynamics and an optimization procedure to estimate area elements, our discretization method results in accurate solutions while still being robust when facing abrupt changes in the density of points. Moreover, the proposed scheme results in numerically stable discrete operators. The effectiveness of the proposed technique is brought to bear in many practical applications. In particular, we use the eigenstructure of the discrete operator for filtering and shape segmentation. Point-based surface deformation is another application that can be easily carried out from the proposed discretization method.
Spectral image segmentation using image decomposition and inner product-based metric
W. Casaca, A. Paiva, E. G. Nieto, P. Joia and L. G. Nonato
Journal of Mathematical Imaging and Vision, 45 (3), pp. 227-238, 2013
Paper (PDF) ● BIB ● DOI
Abstract: Image segmentation is an indispensable tool in computer vision applications, such as recognition, detection, and tracking. In this work, we introduce a novel user-assisted image segmentation technique which combines image decomposition, inner product-based similarity metric, and spectral graph theory into a concise and unified framework.First, we perform an image decomposition to split the image into texture and cartoon components. Then, an affinity graph is generated and the weights are assigned to its edges according to a gradient-based inner-product function. From the eigenstructure of the affinity graph, the image is partitioned through the spectral cut of the underlying graph. The computational effort of our framework is alleviated by an image coarsening process, which reduces the graph size considerably. Finally, a coarse-to-fine interpolation is applied in order to assemble the partition back onto the original image. The efficiency of the proposed methodology is attested by comparisons with state-of-the-art spectral segmentation methods through a qualitative and quantitative analysis of the results. Extended version of SIBGRAPI 2011 paper.
Journal of Mathematical Imaging and Vision, 45 (3), pp. 227-238, 2013
Paper (PDF) ● BIB ● DOI
Abstract: Image segmentation is an indispensable tool in computer vision applications, such as recognition, detection, and tracking. In this work, we introduce a novel user-assisted image segmentation technique which combines image decomposition, inner product-based similarity metric, and spectral graph theory into a concise and unified framework.First, we perform an image decomposition to split the image into texture and cartoon components. Then, an affinity graph is generated and the weights are assigned to its edges according to a gradient-based inner-product function. From the eigenstructure of the affinity graph, the image is partitioned through the spectral cut of the underlying graph. The computational effort of our framework is alleviated by an image coarsening process, which reduces the graph size considerably. Finally, a coarse-to-fine interpolation is applied in order to assemble the partition back onto the original image. The efficiency of the proposed methodology is attested by comparisons with state-of-the-art spectral segmentation methods through a qualitative and quantitative analysis of the results. Extended version of SIBGRAPI 2011 paper.
Approximating implicit curves on triangulations with affine arithmetic
F. Nascimento, A. Paiva, L. H. de Figueiredo and J. Stolfi
Proceedings of SIBGRAPI 2012, pp. 94-101, 2012
Paper (PDF) ● BIB ● DOI
Abstract: We present an adaptive method for computing a robust polygonal approximation of an implicit curve in the plane that uses affine arithmetic to identify regions where the curve lies inside a thin strip. Unlike other interval methods, even those based on affine arithmetic, our method works on triangulations, not only on rectangular quadtrees.
Proceedings of SIBGRAPI 2012, pp. 94-101, 2012
Paper (PDF) ● BIB ● DOI
Abstract: We present an adaptive method for computing a robust polygonal approximation of an implicit curve in the plane that uses affine arithmetic to identify regions where the curve lies inside a thin strip. Unlike other interval methods, even those based on affine arithmetic, our method works on triangulations, not only on rectangular quadtrees.
Colorization by multidimensional projection
W. Casaca, E. G. Nieto, C. Ferreira, G. Tavares, P. Pagliosa, F. Paulovich L. G. Nonato and A. Paiva
Proceedings of SIBGRAPI 2012, pp. 32-38, 2012.
Paper (PDF) ● BIB ● DOI
Abstract: Most image colorization techniques assign colors to grayscale images by embedding image pixels into a highdimensional feature space and applying a color pattern to each cluster of high-dimensional data. A main drawback of such an approach is that, depending on texture patterns and image complexity, clusters of similar pixels can hardly be defined automatically, rendering existing methods prone to fail. In this work we present a novel approach to colorize grayscale images that allows for user intervention. Our methodology makes use of multidimensional projection to map high-dimensional data to a visual space. User can manipulate projected data in the visual space so as to further improve clusters and thus the colorization result. Different from other methods, our interactive tool is ease of use while still being flexible enough to enable local color modification. We show the effectiveness of our approach through a set of examples and comparisons against existing colorization methods.
Proceedings of SIBGRAPI 2012, pp. 32-38, 2012.
Paper (PDF) ● BIB ● DOI
Abstract: Most image colorization techniques assign colors to grayscale images by embedding image pixels into a highdimensional feature space and applying a color pattern to each cluster of high-dimensional data. A main drawback of such an approach is that, depending on texture patterns and image complexity, clusters of similar pixels can hardly be defined automatically, rendering existing methods prone to fail. In this work we present a novel approach to colorize grayscale images that allows for user intervention. Our methodology makes use of multidimensional projection to map high-dimensional data to a visual space. User can manipulate projected data in the visual space so as to further improve clusters and thus the colorization result. Different from other methods, our interactive tool is ease of use while still being flexible enough to enable local color modification. We show the effectiveness of our approach through a set of examples and comparisons against existing colorization methods.
Class-specific metrics for multidimensional data projection applied to CBIR
P. Joia, E. G. Nieto, J. B. Neto, W. Casaca, G. Botelho, A. Paiva and L. G. Nonato
The Visual Computer, 28 (10), pp. 1027-1037, 2012
Paper (PDF) ● BIB ● DOI
Abstract: Content-based image retrieval is still a challenging issue due to the inherent complexity of images and choice of the most discriminant descriptors. Recent developments in the field have introduced multidimensional projections to burst accuracy in the retrieval process, but many issues such as introduction of pattern recognition tasks and deeper user intervention to assist the process of choosing the most discriminant features still remain unaddressed. In this paper we present a novel framework to CBIR that combines pattern recognition tasks, class-specific metrics and multidimensional projection to devise an effective and interactive image retrieval system. User interaction plays an essential role in the computation of the final multidimensional projection from which image retrieval will be attained. Results have shown that the proposed approach outperforms existing methods, turning out to be a very attractive alternativefor managing image data sets. Extended version of SIBGRAPI 2011 paper.
The Visual Computer, 28 (10), pp. 1027-1037, 2012
Paper (PDF) ● BIB ● DOI
Abstract: Content-based image retrieval is still a challenging issue due to the inherent complexity of images and choice of the most discriminant descriptors. Recent developments in the field have introduced multidimensional projections to burst accuracy in the retrieval process, but many issues such as introduction of pattern recognition tasks and deeper user intervention to assist the process of choosing the most discriminant features still remain unaddressed. In this paper we present a novel framework to CBIR that combines pattern recognition tasks, class-specific metrics and multidimensional projection to devise an effective and interactive image retrieval system. User interaction plays an essential role in the computation of the final multidimensional projection from which image retrieval will be attained. Results have shown that the proposed approach outperforms existing methods, turning out to be a very attractive alternativefor managing image data sets. Extended version of SIBGRAPI 2011 paper.
Sketch-based adaptive mesh augmentation using stellar operators
A. Paiva, R. Amorim, L. Velho and M. C. Sousa
Proceedings of SIBGRAPI 2011, pp. 78-85, 2011.
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: In this paper we present a new method for modeling and editing surface detail using free-form curves and a natural interface. It combines in a original way an adaptive multiresolution mesh structure with a simple, intuitive sketch-based interface. One of the novel contributions of this work is the curve sensitive mesh resolution control, which allows the definition of a rich set of operators that locally modify the surface geometry. Furthermore, the present framework provides the basic functionality to build a complete feature based modeling system.
Proceedings of SIBGRAPI 2011, pp. 78-85, 2011.
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: In this paper we present a new method for modeling and editing surface detail using free-form curves and a natural interface. It combines in a original way an adaptive multiresolution mesh structure with a simple, intuitive sketch-based interface. One of the novel contributions of this work is the curve sensitive mesh resolution control, which allows the definition of a rich set of operators that locally modify the surface geometry. Furthermore, the present framework provides the basic functionality to build a complete feature based modeling system.
Spectral segmentation using cartoon-texture decomposition and inner product-based metric
W. Casaca, A. Paiva and L. G. Nonato
Proceedings of SIBGRAPI 2011, pp. 266-273, 2011
Best Paper Award in Image Processing
Paper (PDF) ● BIB ● DOI
Abstract: This paper presents a user-assisted image partition technique that combines Cartoon-Texture Decomposition (CTD), inner product-based similarity metric, and spectral cut into a unified framework. The CTD is used to first split the image into textured and texture-free components, the latter being used to define a gradient-based inner-product function. An affinity graph is then derived and weights are assigned to its edges according to the inner product-based metric. Spectral cut is computed on the affinity graph so as to partition the image. The computational burden of the spectral cut is mitigated by a fine-to-coarse image representation process, which enables moderate size graphs that can be handled more efficiently. The partitioning can be steered by interactively by changing the weights of the graph through user strokes. Weights are updated by combining the texture component computed in the first stage of our pipeline and a recent harmonic analysis technique that captures waving patterns. Finally, a coarse-to-fine interpolation is applied in order to project the partition back onto the original image. The suitable performance of the proposed methodology is attested by comparisons against state-of-art spectral segmentation methods.
Proceedings of SIBGRAPI 2011, pp. 266-273, 2011
Best Paper Award in Image Processing
Paper (PDF) ● BIB ● DOI
Abstract: This paper presents a user-assisted image partition technique that combines Cartoon-Texture Decomposition (CTD), inner product-based similarity metric, and spectral cut into a unified framework. The CTD is used to first split the image into textured and texture-free components, the latter being used to define a gradient-based inner-product function. An affinity graph is then derived and weights are assigned to its edges according to the inner product-based metric. Spectral cut is computed on the affinity graph so as to partition the image. The computational burden of the spectral cut is mitigated by a fine-to-coarse image representation process, which enables moderate size graphs that can be handled more efficiently. The partitioning can be steered by interactively by changing the weights of the graph through user strokes. Weights are updated by combining the texture component computed in the first stage of our pipeline and a recent harmonic analysis technique that captures waving patterns. Finally, a coarse-to-fine interpolation is applied in order to project the partition back onto the original image. The suitable performance of the proposed methodology is attested by comparisons against state-of-art spectral segmentation methods.
Projection-based image retrieval using class-specific metrics
P. Joia, E. G. Nieto, G. Botelho, J. B. Neto, A. Paiva and L. G. Nonato
Proceedings of SIBGRAPI 2011, pp. 125-132, 2011
Best Paper Award in Computer Graphics/Visualization
Paper (PDF) ● BIB ● DOI
Abstract: Content-based image classification/retrieval based on image descriptors has become an essential component in most database systems. However, most existing systems do not provide mechanisms that enable interactive multi-objective queries, hampering the user experience. In this paper we present a novel methodology capable of accomplishing multi-objective searches while still being interactive. Our approach relies on a combination of class-specific metrics and multidimensional projection to devise an effective and interactive image retrieval system. Besides allowing visual exploration of image data sets, the provided results and comparisons show that the proposed approach outperforms existing methods, turning out to be a veryattractive alternative for managing image data sets.
Proceedings of SIBGRAPI 2011, pp. 125-132, 2011
Best Paper Award in Computer Graphics/Visualization
Paper (PDF) ● BIB ● DOI
Abstract: Content-based image classification/retrieval based on image descriptors has become an essential component in most database systems. However, most existing systems do not provide mechanisms that enable interactive multi-objective queries, hampering the user experience. In this paper we present a novel methodology capable of accomplishing multi-objective searches while still being interactive. Our approach relies on a combination of class-specific metrics and multidimensional projection to devise an effective and interactive image retrieval system. Besides allowing visual exploration of image data sets, the provided results and comparisons show that the proposed approach outperforms existing methods, turning out to be a veryattractive alternative for managing image data sets.
Meshless Helmholtz-Hodge decomposition
F. Petronetto, A. Paiva, M. Lage, G. Tavares, H. Lopes and T. Lewiner
IEEE Transactions on Visualization and Computer Graphics, 16 (2), pp. 338-349, 2010
Paper (PDF) ● BIB ● DOI
Abstract: Vector fields analysis traditionally distinguishes conservative (curl-free) from mass preserving (divergence-free) components. The Helmholtz-Hodge decomposition allows separating any vector field into the sum of three uniquely defined components: curl-free, divergence-free and harmonic. This decomposition is usually achieved by using mesh-based methods such as finite differences or finite elements. This work presents a new meshless approach to the Helmholtz-Hodge decomposition for the analysis of 2D discrete vector fields. It embeds into the SPH particle-based framework. The proposed method is efficient and can be applied to extract features from a 2D discrete vector field and to multiphase fluid flow simulation to ensure incompressibility.
IEEE Transactions on Visualization and Computer Graphics, 16 (2), pp. 338-349, 2010
Paper (PDF) ● BIB ● DOI
Abstract: Vector fields analysis traditionally distinguishes conservative (curl-free) from mass preserving (divergence-free) components. The Helmholtz-Hodge decomposition allows separating any vector field into the sum of three uniquely defined components: curl-free, divergence-free and harmonic. This decomposition is usually achieved by using mesh-based methods such as finite differences or finite elements. This work presents a new meshless approach to the Helmholtz-Hodge decomposition for the analysis of 2D discrete vector fields. It embeds into the SPH particle-based framework. The proposed method is efficient and can be applied to extract features from a 2D discrete vector field and to multiphase fluid flow simulation to ensure incompressibility.
Fluid-based hatching for tone mapping in line illustrations
A. Paiva, E. V. Brazil, F. Petronetto and M. C. Sousa
The Visual Computer (CGI 2009), 25 (5-7), pp. 519-527, 2009
Paper (PDF) ● BIB ● DOI
Abstract: This paper presents a novel physically-based, meshless framework for line art rendering of surfaces with complex geometry and arbitrary topology. We apply an inviscid fluid flow simulation using Smoothed Particles Hydrodynamics to compute the global velocity and cross fields over the surface model. These fields guide the automatic placement of strokes while marinating the geometric and topological coherence with the model. Target tones are matched by tonal value maps allowing different hatching and cross-hatching effects. We demonstrate the simplicity and effectiveness of our method with sample renderings obtained for a variety of models.
The Visual Computer (CGI 2009), 25 (5-7), pp. 519-527, 2009
Paper (PDF) ● BIB ● DOI
Abstract: This paper presents a novel physically-based, meshless framework for line art rendering of surfaces with complex geometry and arbitrary topology. We apply an inviscid fluid flow simulation using Smoothed Particles Hydrodynamics to compute the global velocity and cross fields over the surface model. These fields guide the automatic placement of strokes while marinating the geometric and topological coherence with the model. Target tones are matched by tonal value maps allowing different hatching and cross-hatching effects. We demonstrate the simplicity and effectiveness of our method with sample renderings obtained for a variety of models.
Particle-based viscoplastic fluid/solid simulation
A. Paiva, F. Petronetto, T. Lewiner and G. Tavares
Computer-Aided Design, 41(4), pp. 306-314, 2009
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: Simulations of viscoplastic materials are traditionally governed by continuum mechanics. The viscous behavior is typically modeled as an internal force, defined by diverse quantities. This work introduces a fluid model to simulate the viscoplastic effect of solid materials, such as plastic, wax, clay and polymer. Our method consists in modeling a solid object through a non-Newtonian fluid with high viscosity. This fluid simulation uses the Smoothed Particle Hydrodynamics and the viscosity is formulated by using the General Newtonian Fluid model. This model concentrates the viscoplasticity in a single parameter. Our results show clear effects of creep, melting, hardening and flowing. Extended version of SIBGRAPI 2006 paper.
Computer-Aided Design, 41(4), pp. 306-314, 2009
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: Simulations of viscoplastic materials are traditionally governed by continuum mechanics. The viscous behavior is typically modeled as an internal force, defined by diverse quantities. This work introduces a fluid model to simulate the viscoplastic effect of solid materials, such as plastic, wax, clay and polymer. Our method consists in modeling a solid object through a non-Newtonian fluid with high viscosity. This fluid simulation uses the Smoothed Particle Hydrodynamics and the viscosity is formulated by using the General Newtonian Fluid model. This model concentrates the viscoplasticity in a single parameter. Our results show clear effects of creep, melting, hardening and flowing. Extended version of SIBGRAPI 2006 paper.
Robust adaptive meshes for implicit surfaces
A. Paiva, H. Lopes, T. Lewiner and L.H. de Figueiredo
Proceedings of SIBGRAPI 2006, pp. 205-212, 2006
Paper (PDF) ● BIB ● DOI
Abstract: This work introduces a robust algorithm for computing good polygonal approximations of implicit surfaces, where robustness entails recovering the exact topology of the implicit surface. Furthermore, the approximate triangle mesh adapts to the geometry and to the topology of the real implicit surface. This method generates an octree subdivided according to the interval evaluation of the implicit function in order to guarantee the robustness, and to the interval automatic differentiation in order to adapt the octree to the geometry of the implicit surface. The triangle mesh is then generated from that octree through an enhanced dual marching.
Proceedings of SIBGRAPI 2006, pp. 205-212, 2006
Paper (PDF) ● BIB ● DOI
Abstract: This work introduces a robust algorithm for computing good polygonal approximations of implicit surfaces, where robustness entails recovering the exact topology of the implicit surface. Furthermore, the approximate triangle mesh adapts to the geometry and to the topology of the real implicit surface. This method generates an octree subdivided according to the interval evaluation of the implicit function in order to guarantee the robustness, and to the interval automatic differentiation in order to adapt the octree to the geometry of the implicit surface. The triangle mesh is then generated from that octree through an enhanced dual marching.
Vector field reconstruction from sparse samples with applications
M. Lage, F. Petronetto, A. Paiva, H. Lopes, T. Lewiner and G. Tavares
Proceedings of SIBGRAPI 2006, pp. 297-304, 2006.
Paper (PDF) ● BIB ● DOI
Abstract: This work presents a novel algorithm for 2D vector field reconstruction from sparse set of points-vectors pairs. Our approach subdivides the domain adaptively in order to make local piecewise polynomial approximations for the field. It uses partition of unity to blend those local approximations together, generating a global approximation for the field. The flexibility of this scheme allows handling data from very different sources. In particular, this work presents important applications of the proposed method to velocity and acceleration fields analysis, in particular for fluid dynamics visualization. Front cover of the Proceedings of SIBGRAPI 2006.
Proceedings of SIBGRAPI 2006, pp. 297-304, 2006.
Paper (PDF) ● BIB ● DOI
Abstract: This work presents a novel algorithm for 2D vector field reconstruction from sparse set of points-vectors pairs. Our approach subdivides the domain adaptively in order to make local piecewise polynomial approximations for the field. It uses partition of unity to blend those local approximations together, generating a global approximation for the field. The flexibility of this scheme allows handling data from very different sources. In particular, this work presents important applications of the proposed method to velocity and acceleration fields analysis, in particular for fluid dynamics visualization. Front cover of the Proceedings of SIBGRAPI 2006.
Particle-based non-Newtonian fluid animation for melting objects
A. Paiva, F. Petronetto, T. Lewiner and G. Tavares
Proceedings of SIBGRAPI 2006, pp. 78-85, 2006
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: This paper presents a new visually realistic animation technique for objects that melt and flow. It simulates viscoplastic properties of materials such as metal, plastic, wax, polymer and lava. The technique consists in modeling the object by the transition of a non-Newtonian fluid with high viscosity to a liquid of low viscosity. During the melting, the viscosity is formulated using the General Newtonian fluids model, whose properties depend on the local temperature. The phase transition is then driven by the heat equation. The fluid simulation framework uses a variation of the Lagrangian method called Smoothed Particle Hydrodynamics. This paper also includes several schemes that improve the efficiency and the numerical stability of the equations.
Proceedings of SIBGRAPI 2006, pp. 78-85, 2006
Paper (PDF) ● Video (YouTube) ● BIB ● DOI
Abstract: This paper presents a new visually realistic animation technique for objects that melt and flow. It simulates viscoplastic properties of materials such as metal, plastic, wax, polymer and lava. The technique consists in modeling the object by the transition of a non-Newtonian fluid with high viscosity to a liquid of low viscosity. During the melting, the viscosity is formulated using the General Newtonian fluids model, whose properties depend on the local temperature. The phase transition is then driven by the heat equation. The fluid simulation framework uses a variation of the Lagrangian method called Smoothed Particle Hydrodynamics. This paper also includes several schemes that improve the efficiency and the numerical stability of the equations.
Robust visualization of strange attractors using affine arithmetic
A. Paiva, L.H. de Figueiredo and J. Stolfi
Computers & Graphics, 30 (6), pp. 1020-1026, 2006
Paper (PDF) ● BIB ● DOI
Abstract: We propose the use of affine arithmetic in cell-mapping methods for the robust visualization of strange attractors and show that the resulting cellular approximations converge faster than those produced by cell-mapping methods based on classical interval arithmetic.
Computers & Graphics, 30 (6), pp. 1020-1026, 2006
Paper (PDF) ● BIB ● DOI
Abstract: We propose the use of affine arithmetic in cell-mapping methods for the robust visualization of strange attractors and show that the resulting cellular approximations converge faster than those produced by cell-mapping methods based on classical interval arithmetic.
Page updated
Report abuse