Home
Events
Education
Research
Committees
Teaching
Advising
Software
Publications
Music
     

Tao Ju

Assistant Professor
Department of Computer Science and Engineering
Washington University in St. Louis
Office: Lopata Hall 517
Phone: 314-935-6648
Email: taoju at cse.wustl.edu
Events
Education
Research
I am part of the Media and Machines Lab. My research interests lie in the field of computer graphics and its applications in bio-medical research. I am particularly interested in geometric modeling, mesh processing, visualization, and bio-medical modeling.

Publications by year or by topic:

  Volumetric methods in mesh processing
  Barycentric coordinates and mesh deformation
  Protein structure identification using skeletons
  Deformable imaging using anatomical atlases
Program Committees
  • Eurographics (2009).
  • Siggraph Asia (2008).
  • International Symposium on Visual Computing (2008).
  • Pacific Graphics (2008).
  • Computer Animation and Social Agents (2008).
  • ACM Symposium on Solid and Physical Modeling (2008).
  • Geometric Modeling and Processing (2008).
  • Pacific Graphics (Co-chair, 2007).
  • International Symposium on Visual Computing (Co-chair, 2007).
  • ACM Symposium on Solid and Physical Modeling (2007).
  • Eurographics (2007).
  • Pacific Graphics (2006).
  • International Symposium on Visual Computing (2006).
  • Geometric Modeling and Processing (2006).
  • The 24th Computer Graphics International Conference (2006),
Teaching
Advising
Graphics is fun...and cool! If you are an undergraduate interested in fun and rewarding research, or a master student looking for project ideas, send me an email or come by my office.
Software
  • PolyMender: A robust, efficient program for removing cracks, holes, T-juctions and self-intersections from arbitrary polygonal models.
  • StackAligner: A graphical toolkit for reconstructing smooth 3D volumes from a stack of serial tissue sections containing both global and local distortions.
  • GeneAtlas.org: A web-based spatial database of gene expression patterns in the mouse brain featuring fully-customized graphical query interface.
Publications
2008
Interactive Separation of Segmented Bones in CT Volumes Using Graph Cut
Lecture Notes in Computer Science (Proceedings of MICCAI 2008), 5241:296-304 (Paper, Poster)
L. Liu, D. Raber, D. Nopachai, P. Commean, D. Sinacore, F. Prior, R. Pless, and T. Ju

We present a fast, interactive method for separating bones that have been collectively segmented from a CT volume. Given userprovided seed points, the method computes the separation as a multiway cut on a weighted graph constructed from the binary, segmented volume. By properly designing and weighting the graph, we show that the resulting cut can accurately be placed at bone-interfaces using only a small number of seed points even when the data is noisy. The method has been implemented with an interactive graphical interface, and used to separate the 12 human foot bones in 10 CT volumes. The interactive tool produced compatible result with a ground-truth separation, generated by a completely manual labelling procedure, while reducing the human interaction time from a mean of 2.4 hours per volume in manual labelling down to approximately 18 minutes.
Segmentation-free skeletonization of grayscale volumes for shape understanding.
IEEE International Conference on Shape Modeling and Applications 2008, pp. 63-71 (Paper)
S. Abeysinghe, M. Baker, W. Chiu, T. Ju

Medical imaging has produced a large number of volumetric images capturing biological structures in 3D. Computer-based understanding of these structures can often benefit from the knowledge of shape components, particularly rod-like and plate-like parts, in such volumes. Previously, skeletons have been a common tool for identifying these shape components in a solid object. However, obtaining skeletons of a grayscale volume poses new challenges due to the lack of a clear boundary between object and background. In this paper, we present a new skeletonization algorithm on grayscale volumes typical to medical imaging (e.g., MRI, CT and EM scans), for the purpose of identifying shape components. Our algorithm does not require an explicit segmentation of the volume into object and background, and is capable of producing skeletal curves and surfaces that lie centered at rod-shaped and plate-shaped parts in the grayscale volume. Our method is demonstrated on both synthetic and medical data.
Surface Reconstruction From Non-parallel Curve Networks.
Computer Graphics Forum (Proceedings of Eurographics 2008), 27(2):155-163 (Paper)
L. Liu, C. Bajaj, J.O. Deasy, D.A. Low, T. Ju

Building surfaces from cross-section curves has wide applications including bio-medical modeling. Previous work in this area has mostly focused on connecting simple closed curves on parallel cross-sections. Here we consider the more general problem where input data may lie on non-parallel cross-sections and consist of curve networks that represent the segmentation of the underlying object by different material or tissue types (e.g., skin, muscle, bone, etc.) on each cross-section. The desired output is a surface network that models both the exterior surface and the internal partitioning of the object. We introduce an algorithm that is capable of handling curve networks of arbitrary shape and topology on cross-section planes with arbitrary orientations. Our algorithm is simple to implement and is guaranteed to produce a closed surface network that interpolates the curve network on each cross-section. Our method is demonstrated on both synthetic and bio-medical examples.
Geometric interpretation of the Gamma dose distribution comparison technique: interpolation-free calculation
Medical Physics, 35(3):879-887 (Paper)
T. Ju, T. Simpson, J.O. Deasy, D.A. Low

The Gamma dose distribution comparison tool has been used by numerous investigators to quantitatively compare multidimensional dose distributions. The tool simultaneously evaluates the dose difference and distance-to-agreement of two dose distributions. One of the weaknesses of the tool is that the comparison requires one of the dose distributions to have a relatively high spatial resolution, due to the fact that the Gamma tool measures the closest pixel in one of the dose distributions with individual pixels of another, and this closest distance can not be accurately measured unless the pixels are finely spaced, which requires time-consuming interpolation. We provide a reinterpretation of the Gamma distance as the geometric distance between two 3D or 4D meshes representing the two 2D or 3D dose distributions. The geometric approach avoids the drastic growth of calculation time incurred by interpolation and makes the Gamma tool more practical and more accurate.
Shape modeling and matching in identifying 3D protein structures.
Computer Aided-Design, 40:708-720 (Paper)
S. Abeysinghe, T. Ju, M. L. Baker, W. Chiu

In this paper, we describe a novel geometric approach in the process of recovering 3D protein structures from scalar volumes. The input to our method is a sequence of alpha-helices that make up a protein, and a low-resolution protein density volume where possible locations of alpha-helices have been detected. Our task is to identify the correspondence between the two sets of helices, which will shed light on how the protein folds in space. The central theme of our approach is to cast the correspondence problem as that of shape matching between the 3D volume and the 1D sequence. We model both shapes as attributed relational graphs, and formulate a constrained inexact graph matching problem. To compute the matching, we developed an optimal algorithm based on the A*-search with several choices of heuristic functions. As demonstrated in a suite of synthetic and authentic inputs, the shape-modeling approach is capable of identifying helix correspondences in noise-abundant volumes at high accuracy with minimal or no user intervention.
Tarsal and Metatarsal Bone Mineral Density Measurement using Volumetric Quantitative Computed Tomography.
Journal of Digital Imaging, To Appear
P. K. Commean, T. Ju, L. Liu, D. R. Sinacore, M. K. Hastings, M. J. Mueller

A new method for measuring bone mineral density (BMD) of the tarsal and metatarsals is described using volumetric quantitative computed tomography (VQCT) in conjunction with geometric subdivision in subjects with diabetes mellitus and peripheral neuropathy. In addition to whole-bone segmentation and measurement, we performed atlas-based partitioning of sub-regions within the second metatarsal for all subjects, from which the volumes and BMDs were obtained for each sub-region. The sub-region measurement BMD errors (root mean square coefficient of variation) within the shaft, proximal end and distal end were shown to vary by approximately 1% between the two scans of each subject. These methods can provide an important outcome measure for clinical research trials investigating the effects of interventions, aging or disease progression on bone loss or gain in individual foot bones.
2007
Editing The Topology of 3D Models by Sketching.
ACM Transactions on Graphics (Proceedings of SIGGRAPH 2007), 26(3): 42 (Paper, Video (33MB))
T. Ju, Q-Y Zhou, S-M Hu

We present a method for modifying the topology of a 3D model with user control. The heart of our method is a guided topology editing algorithm. Given a source model and a user-provided target shape, the algorithm modifies the source so that the resulting model is topologically consistent with the target. Our algorithm permits removing or adding various topological features (e.g., handles, cavities and islands) in a common framework and ensures that each topological change is made by minimal modification to the source model. To create the target shape, we have also designed a convenient 2D sketching interface for drawing 3D line skeletons. As demonstrated in a suite of examples, the use of sketching allows more accurate removal of topological artifacts than previous methods, and enables creative designs with specific topological goals.
Real-time homogenous translucent material editing.
Computer Graphics Forum (Proceedings of Eurographics 2007), 26(3):545-552 (Paper)
K. Xu, Y. Gao, Y. Li, T. Ju, S.-M. Hu

This paper presents a novel method for real-time homogenous translucent material editing under fixed illumination. We consider the complete analytic BSSRDF model proposed by Jensen et al. [JMLH01], including both multiple scattering and single scattering. Our method allows the user to adjust the analytic parameters of BSSRDF and provides high-quality, real-time rendering feedback. Inspired by recently developed Precomputed Radiance Transfer (PRT) techniques, we approximate both the multiple scattering diffuse reflectance function and the single scattering exponential attenuation function in the analytic model using basis functions, so that re-computing the outgoing radiance at each vertex as parameters change reduces to simple dot products. In addition, using a non-uniform piecewise polynomial basis, we are able to achieve smaller approximation error than using bases adopted in previous PRT-based works, such as spherical harmonics and wavelets. Using hardware acceleration, we demonstrate that our system generates images comparable to [JMLH01] at real-time frame-rates.
Shape modeling and matching in identifying protein structure from low-resolution images.
ACM Symposium on Solid and Physical Modeling 2007 (Paper,Talk)
S. Abeysinghe, T. Ju, M. Baker, W. Chiu

In this paper, we describe a novel, shape-modeling approach to recovering 3D protein structures from volumetric images. The input to our method is a sequence of a-helices that make up a protein, and a low-resolution volumetric image of the protein where possible locations of a-helices have been detected. Our task is to identify the correspondence between the two sets of helices, which will shed light on how the protein folds in space. The central theme of our approach is to cast the correspondence problem as that of shape matching between the 3D volume and the 1D sequence. We model both the shapes as attributed relational graphs, and formulate a constrained inexact graph matching problem. To compute the matching, we developed an optimal algorithm based on the A*-search with several choices of heuristic functions. As demonstrated in a suite of real protein data, the shape-modeling approach is capable of correctly identifying helix correspondences in noise-abundant volumes with minimal or no user intervention.
Computing a family of skeletons of volumetric models for shape description
Computer-Aided Design, 39(5):352-360 (Paper)
T. Ju, M. Baker, W. Chiu

Skeletons are important shape descriptors in object representation and recognition. Typically, skeletons of volumetric models are computed using iterative thinning. However, traditional thinning methods often generate skeletons with complex structures that are unsuitable for shape description, and appropriate pruning methods are lacking. In this paper, we present a new method for computing skeletons of volumetric models by alternating thinning and a novel skeleton pruning routine. Our method creates a family of skeletons parameterized by two user-specified numbers that determine respectively the size of curve and surface features on the skeleton. As demonstrated on both real-world models and protein images in bio-medical research, our method generates skeletons with simple and meaningful structures that are particularly suitable for describing cylindrical and plate-like shapes.
Learning-based Segmentation Framework for Tissue Images Containing Gene Expression Data
IEEE Transactions on Medical Imaging, 26(5):728-744 (Paper)
M. Bello, T. Ju, J. Carson, J. Warren, W. Chiu, I.A. Kakadiaris

Associating specific gene activity with functional locations in the brain results in a greater understanding of the role of the gene. To perform such an association for the over 20,000 genes in the mammalian genome, reliable automated methods that characterize the distribution of gene expression in relation to a standard anatomical model are required. In this paper, we propose a new automatic method that results in the segmentation of gene expression images into distinct anatomical regions in which the expression can be quantified and compared with other images. Our contribution is a novel hybrid atlas that utilizes a statistical shape model based on a subdivision mesh, texture differentiation at region boundaries, and features of anatomical landmarks to delineate boundaries of anatomical regions in gene expression images. This atlas, which provides a common coordinate system for internal brain data, is being used to create a searchable database of gene expression patterns in the adult mouse brain. Our framework annotates the images about four times faster and has achieved a median spatial overlap of up to 0.92 compared with expert segmentation in 64 images tested. This tool is intended to help scientists interpret large-scale gene expression patterns more efficiently.
Topology Repair of Solid Models Using Skeletons
IEEE Transactions on Visualization and Computer Graphics, 13(4):675-685 (Paper)
Q-Y Zhou, T. Ju, S-M. Hu

We present a method for repairing topological errors on solid models in the form of small surface handles, which often arise from surface reconstruction algorithms. We utilize a skeleton representation that offers a new mechanism for identifying and measuring handles. Our method presents two unique advantages over previous approaches. First, handle removal is guaranteed not to introduce invalid geometry or additional handles. Second, by using an adaptive grid structure, our method is capable of processing huge models efficiently at high resolutions.
A general geometric construction of coordinates in a convex simplicial polytope
Computer Aided Geometric Design, 24(3): 161-178 (Paper, Talk)
T. Ju, P. Liepa, J. Warren

Barycentric coordinates are a fundamental concept in computer graphics and geometric modeling. We extend the geometric construction of Floater's mean value coordinates to a general form that is capable of constructing a family of coordinates in a convex 2D polygon, 3D triangular polyhedron, or a higher-dimensional simplicial polytope. This family unifies previously known coordinates, including Wachspress coordinates, mean value coordinates and discrete harmonic coordinates, in a simple geometric framework. Using the construction, we are able to create a new set of coordinates in 3D and higher dimensions and study its relation with known coordinates. We show that our general construction is complete, that is, the resulting family includes all possible coordinates in any convex simplicial polytope.
Manifold Dual Contouring
IEEE Transactions on Visualization and Computer Graphics, 13(3):610-619 (Paper)
S. Schaefer, T. Ju, J. Warren

Dual Contouring is a feature-preserving iso-surfacing method that extracts crack-free surfaces from both uniform and adaptive octree grids. We present an extension of Dual Contouring that further guarantees that the mesh generated is a manifold even under adaptive simplification. Our main contribution is an octree-based, topology-preserving vertex clustering algorithm for adaptive contouring. The contoured surface generated by our method contains only manifold vertices and edges, preserves sharp features, and possesses much better adaptivity than those generated by other iso-surfacing methods under topologically safe simplification.
Identification of Secondary Structure Elements in Intermediate Resolution Density Maps
Structure, 15(1):7-19, 2007 (Paper)
M. L. Baker, T. Ju, W. Chiu

An increasing number of structural studies of large macromolecular complexes, both in X-ray crystallography and electron cryomicroscopy, have resulted in intermediate resolution (5-10 A) structures. Despite being limited in resolution, significant structural and functional information may be extractable from these maps. To aid in the analysis and annotation of these complexes, we have developed SSEhunter, a tool for the quantitative detection of alpha-helices and beta-sheets. Based on density skeletonization, local geometry calculations and a template-based search, SSEhunter has been tested and validated on a variety of simulated and authentic subnanometer resolution density maps. The result is a robust, user-friendly approach that allows users to quickly visualize, assess and annotate intermediate resolution density maps. Beyond secondary structure element identification, the skeletonization algorithm in SSEhunter provides secondary structure topology, potentially useful in leading to structural models of individual molecular components directly from the density.
2006
A Unified, Integral Construction For Coordinates Over Closed Curves
Computer-Aided Geometric Design, 24(8-9):481-493 (Paper)
S. Schaefer, T. Ju and J. Warren

We propose a simple generalization of Shephard's interpolation to piecewise smooth, convex closed curves that yields a family of boundary interpolants with linear precision. Two instances of this family reduce to previously known interpolants: one based on a generalization of Wachspress coordinates to smooth curves and the other an integral version of mean value coordinates for smooth curves. A third instance of this family yields a previously unknown generalization of discrete harmonic coordinates to smooth curves. For closed, piecewise linear curves, we prove that our interpolant reproduces a general family of barycentric coordinates considered by Floater, Hormann and Kos that includes Wachspress coordinates, mean value coordinates and discrete harmonic coordinates.
Probing 3'-ssDNA Loop Formation in E. coli RecBCD/RecBC-DNA Complexes using Non-natural DNA: A Model for "Chi" Recognition Complexes
Journal of Molecular Biology, 362(1):26-43 (Paper)
C. Jason Wong, Rachel L. Rice, Nathan A. Baker, Tao Ju and Timothy M. Lohman

The equilibrium binding of E. coli RecBC and RecBCD helicases to duplex DNA ends containing varying lengths of polyethylene glycol (PEG) spacers within pre-formed 3'-single-stranded (ss) DNA ((dT)n) tails were studied. These studies were designed to test a previous proposal that the 3'-(dT)n tail can be looped out upon binding RecBC and RecBCD for 3'-ssDNA tails with n>=6 nucleotides. Equilibrium binding of protein to unlabeled DNA substrates with ends containing PEG-substituted 3'-ssDNA tails was examined by competition with a Cy3-labeled reference DNA which undergoes a Cy3 fluorescence enhancement upon protein binding. We find that the binding affinities of both RecBC and RecBCD for a DNA end are unaffected upon substituting PEG for the ssDNA between the sixth and the final two nucleotides of the 3'-(dT)n tail. However, placing PEG at the end of the 3'-(dT)n tail increases the binding affinities to their maximum values (i.e. the same as binding constants for RecBC or RecBCD to a DNA end with only a 3'-(dT)6 tail). Equilibrium binding studies of a RecBC mutant containing a nuclease domain deletion, RecBC suggest that looping of the 3'-tail (when n>=6 nucleotides) occurs even in the absence of the RecB nuclease domain, the nuclease domain stabilizes such loop formation. Computer modeling of the RecBCD-DNA complexes suggests that the loop in the 3'-ssDNA tail may form at the RecB/RecC interface. Based on these results we suggest a model for how a loop in the 3'-ssDNA tail might form upon encounter of a "Chi" recognition sequence during unwinding of DNA by the RecBCD helicase.
Intersection-free Contouring on An Octree Grid
Proceedings of Pacific Graphics, 2006 (Paper)
T. Ju and T. Udeshi

A method for extracting intersection-free iso-surfaces from volumetric data with an octree structure is presented. Unlike contouring techniques designed for uniform grids (such as Marching Cubes), adaptive contouring methods (such as Dual Contouring) can and do often generate intersecting polygons. Our main contribution is a polygon generation algorithm that produces triangles enclosed in nonoverlapping volumes, which guarantees an intersection-free mesh. Like other adaptive contouring methods, this new method generates crack-free and feature-preserving surfaces on both uniform and octree grids. We demonstrate the method on both scanned objects and industrial models.
Computing a family of skeletons of volumetric models for shape description
Proceedings of Geometric Modeling and Processing 2006, pp. 235 - 247 (Paper,Talk)
T. Ju, M. Baker, and W. Chiu

Skeletons are important shape descriptors in object representation and recognition. Typically, skeletons of volumetric models are computed via an iterative thinning process. However, traditional thinning methods often generate skeletons with complex structures that are unsuitable for shape description, and appropriate pruning methods are lacking. In this paper, we present a new method for computing skeletons on volumes by alternating thinning and a novel skeleton pruning routine. Our method creates a family of skeletons parameterized by two user-specified numbers that determine respectively the size of curve and surface features on the skeleton. As demonstrated on both real-world models and medical images, our method generates skeletons with simple and meaningful structures that are particularly suitable for describing cylindrical and plate-like shapes.
3D Volume Reconstruction of a Mouse Brain from Histological Sections using Warp Filtering
Journal of Neuroscience Methods, 156(1-2):84-100 (Paper,Program)
T. Ju, J. Warren, J. Carson, M. Bello, I. Kakadiaris, W. Chiu, C. Thaller and G. Eichele

Sectioning tissues for optical microscopy often introduces upon the resulting sec- tions distortions that make 3d reconstruction di”Ącult. Here we present an automatic method for producing a smooth 3D volume from distorted 2D sections in the ab- sence of any undistorted references. The method is based on pairwise elastic image warps between successive tissue sections, which can be computed by 2D image reg- istration. Using a Gaussian ”„lter, an average warp is computed for each section from the pairwise warps in a group of its neighboring sections. The average warps deform each section to match its neighboring sections, thus creating a smooth vol- ume where corresponding features on successive sections lie close to each other. The proposed method can be used with any existing 2D image registration method for 3D reconstruction. In particular, we present a novel image warping algorithm based on dynamic programming that extends Dynamic Time Warping in 1D speech recog- nition to compute pairwise warps between high-resolution 2D images. The warping algorithm e”Ąciently computes a restricted class of 2D local deformations that are characteristic between successive tissue sections. Finally, a validation framework is proposed and applied to evaluate the quality of reconstruction using both real sections and a synthetic volume.
2005
Building 3D surface networks from 2D curve networks with application to anatomical modeling
Proceedings of Pacific Graphics 2005, 21(8-10):764-773
The Visual Computer, (Paper)
T. Ju, J. Warren, J. Carson, G. Eichele, C. Thaller, W. Chiu, M. Bello and I. Kakadiaris

We present a novel method that automatically constructs a surface network from curve networks with arbitrary topology and partitioning an arbitrary number of materials. The surface network exactly interpolates the curve network on each plane and is guaranteed to be free of gaps or self-intersections. In addition, our method provides a flexible framework for user interaction so that the surface topology can be modified conveniently when necessary. As an application, we applied the method to build a high-resolution 3D model of the mouse brain from 2D anatomical boundaries defined on 350 tissue sections.
A Digital Atlas to Characterize the Mouse Brain Transcriptome
PLoS Computational Biology, 1(4): e41, 2005 (Paper)
J. Carson, T. Ju, H. Lu, C. Thaller, M. Xu, S. Pallas, M. C. Crair, J. Warren, W. Chiu and G. Eichele

Here we have developed a computational method for annotating gene expression patterns in the context of a digital atlas to facilitate custom user-queries and comparisons of this type of data. This procedure has been applied to 200 genes in the postnatal mouse brain. As an illustration of utility, we identify candidate genes that may be related to Parkinson's disease by using the expression of a dopamine transporter in the substantia nigra as a search query pattern. In addition, we discover that transcription factor Rorb is down-regulated in the barrelless mutant relative to control mice by quantitative comparison of expression patterns in layer IV somatosensory cortex.
Hybrid Segmentation Framework for Tissue Images Containing Gene Expression Data
Proceedings of the International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI), p254-261, 2005 (Paper)
M. Bello, T. Ju, J. Warren, J. Carson, W. Chiu, C. Thaller, G. Eichele and I. Kakadiaris

In this work, we propose a new automatic method that results in the segmentation of gene expression images into distinct anatomical regions in which the expression can be quantified and compared with other images. Our method utilizes models of shape of training images, texture differentiation at region boundaries, and features of anatomical landmarks to deform a subdivision mesh-based atlas to fit gene expression images.
Geometric Construction of Coordinates for Convex Polyhedra using Polar Duals
Proceedings of Eurographics Symposium on Geometry Processing, p181-186, 2005 (Paper)
T. Ju, S. Schaefer, J. Warren, M.Desbrun

A fundamental problem in geometry processing is that of expressing a point inside a convex polyhedron as a combination of the vertices of the polyhedron. In this paper, we present a uned geometric construction for building these weighted combinations using the notion of polar duals. We show that our method yields a simple geometric construction for Wachspress's barycentric coordinates, as well as for constructing Colin de Verdire matrices from convex polyhedra, critical step in Lovasz's method with applications to parameterizations.
Mean value coordinates for closed triangular meshes
Proceedings of ACM SIGGRAPH, 2005
ACM Transactions on Graphics, 24(3):561-566 (Paper)
T. Ju, S. Schaefer, J. Warren

In this paper, we generalize mean value coordinates from closed 2D polygons to closed triangular meshes. Given such a mesh P, we show that these coordinates are continuous everywhere and smooth on the interior of P. The coordinates are linear on the triangles of P and can reproduce linear functions on the interior of P. To illustrate their usefulness, we conclude by considering several interesting applications including constructing volumetric textures and surface deformation.


2004
Robust Repair of Polygonal Models
Proceedings of ACM SIGGRAPH, 2004
ACM Transactions on Graphics, 23(3):888-895 (Paper, Slides, Program)
Tao Ju

We present a robust method for repairing arbitrary polygon models. The method is guaranteed to produce a closed surface that partitions the space into disjoint internal and external volumes. Our novel algorithm can efficiently process large models containing millions of polygons and is capable of reproducing sharp features in the original geometry.
Automated Characterization of Gene Expression Patterns with an Atlas of the Mouse Brain
Proceedings of IEEE International Conference of the Engineering in Medicine and Biology Society (EMBS), p2917-2920, 2004 (Paper)
J. P. Carson, T. Ju, C. Thaller, J. Warren, M. Bello, I. Kakadiaris, W. Chiu, G. Eichele

A spatio-temporal map of gene activity in the brain would be an important contribution to the understanding of brain development, disease, and function. Such a resource is now possible using high-throughput in situ hybridization, a method for transcriptome-wide acquisition of cellular resolution gene expression patterns in serial tissue sections. However, querying an enormous quantity of image data requires computational methods for describing and organizing gene expression patterns in a consistent manner. In addressing this, we have developed procedures for automated annotation of gene expression patterns in the postnatal mouse brain.
Landmark-driven, Atlas-based Segmentation of Mouse Brain Tissue Images Containing Gene Expression Data
Proceedings of the International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI), p192-199, 2004 (Paper)
Ioannis A. Kakadiaris, Musodiq Bello, Shiva Arunachalam, Wei Kang, Tao Ju, Joe Warren, James Carson, Wah Chiu, Christina Thaller, and Gregor Eichele

In this paper, we present an anatomical landmark detection method that has been incorporated into an atlasbased segmentation. The addition of this technique significantly increases the accuracy of automated atlas-deformation. The resulting large-scale annotation will help scientists interpret gene expression patterns more rapidly and accurately.
Turtle Geometry in Computer Graphics and Computer Aided Design
Computer-Aided Design, 36(14): 1471-1482, 2004 (Paper)
Ron Goldman, Scott Schaefer, and Tao Ju

LOGO is a programming language incorporating turtle graphics, originally devised for teaching computing to young children in elementary and middle schools. Here we advocate the use of LOGO to help introduce some of the basic concepts of computer graphics and computer aided design to undergraduate and graduate students in colleges and universities. We shall show how to motivate affine coordinates and affine transformations, fractal curves and iterated function systems, relaxation methods and subdivision schemes from elementary notions in turtle geometry and turtle programming.
Recursive Turtle Programs and Iterated Affine Transformations
Computer and Graphics, 28(6): 991-1004, 2004. (Paper)
Tao Ju, Scott Schaefer and Ron Goldman.

Recursive turtle programs (RTP) and iterated affine transformations (IAT) are two popular methods for generating fractals. We show that these two models are equivalent in their expressive power. Conversion algorithms in both directions are presented explicitly from the structure of the RTP and the affine transformations in the IAT.
2003
A geometry database for gene expression data
Proceedings of Eurographics Symposium on Geomtry Processing, 2003. (Paper, Slides)
Tao Ju, Joe Warren, Gregor Eichele, Christina Thaller, Wah Chiu and James Carson

In this paper, we describe the structure of a geometric database for the mouse brain that allows biologists to organize and search gene expression data in the mouse brain. The central component of this database is a standard atlas, represented as a subdivision mesh, that explicitly partitions the mouse brain into key anatomical subregions. Due to this partitioning, user queries comparing expression data between various genes can be restricted to anatomical subregions without difficulty while the multi-resolution structure of the subdivision mesh allows these queries to be processed efficiently. The database and searching tools are available at www.geneatlas.org.
Morphing of rational b-spline curves and surfaces using mass distributions
Proceedings of Eurographics, short papers, 2003. (Paper, Slides)
Tao Ju and Ron Goldman.

A rational B-spline curve or surface is a collection of points associated with a mass (weight) distribution. These mass distributions can be used to exert local control over the morph between two rational B-spline curves or surfaces. Here we propose a technique for designing customized morphs by attaching appropriate mass distributions to target B-spline curves and surfaces. We also develop a user interface for this morphing method that is easy to use and requires no knowledge of B-splines on the part of the user.
Convex Contouring on Volumetric Data
The Visual Computer, 19: 513-525, 2003. (Paper, Slides)
Tao Ju, Scott Schaefer, Joe Warren.

In this paper we present a fast, table-driven isosurface extraction technique on volumetric data. Unlike Marching Cubes or other cell-based algorithms, the proposed polygonization generates convex negative space inside individual cells, enabling fast collision detection on the triangulated isosurface. In our implementation, we are able to perform over 2 million point classifications per second. The algorithm is driven by an automatically constructed look-up table that stores compact decision trees by sign configurations. Using the same technique, we can perform fast, crack-free multi-resolution contouring on nested grids of volumetric data.
2002
Dual Contouring on Hermite Data
Proceedings of ACM SIGGRAPH, 2002. (Paper)
Tao Ju, Frank Losasso, Scott Schaefer and Joe Warren

This paper describes a new method for contouring a signed grid whose edges are tagged by Hermite data (i.e; exact intersection points and normals). We extend this contouring method to the case of multi-signed functions and demonstrate how to model textured contours using multi-signed functions. Using a new, numerically stable representation for quadratic error functions, we develop an octreebased method for simplifying these contours and their textured regions. We next extend our contouring method to these simplified octrees.
2001
Modifying the shape of NURBS surfaces with geometric constraints
Computer Aided Design, 33(12): 903-912, 2001. (Paper)
Shi-Min Hu, You-Fu Li, Tao Ju and Xiang Zhu

NURBS surfaces are among the most commonly used parametric surfaces in CAGD and Computer Graphics. This paper investigates shape modification of NURBS surfaces with geometric constraints, such as point, normal vector, and curve constraints. Two new methods are presented by constrained optimization and energy minimization. The former is based on minimizing changes in control net of surfaces, whereas the latter is based on strain energy minimization. By these two methods, we change control points and weights of an original surface, such that the modified surface satisfies the given constraints. Comparison results and practical examples are also given.
Approximate merging of a pair of Bezier curves
Computer Aided Design, 33(2): 125-136, 2001. (Paper)
Hu Shi-Min, Tong Ruo-Feng, Ju Tao and Sun Jia-Guang

This paper deals with the merging problem, i.e. to approximate two adjacent Bezier curves by a single Bezier curve. A novel approach for approximate merging is introduced in the paper by using the constrained optimization method. The basic idea of this method is to find conditions for the precise merging of Bezier curves first, and then compute the constrained optimization solution by moving the control points. "Discrete" coefficient norm in L2 sense and "squared difference integral" norm are used in our method. Continuity at the endpoints of curves are considered in the merging process, and approximate merging with points constraints are also discussed. Further, it is shown that the degree elevation of original Bezier curves will reduce the merging error.

Comments or suggestions: taoju at cs.wustl.edu