Compact Representations for Fast Nonrigid Registration of Medical Images

We develop efficient techniques for the non-rigid registration of medical images by using representations that adapt to the anatomy found in such images. Images of anatomical structures typically have uniform intensity interiors and smooth boundaries. We create methods to represent such regions compactly using tetrahedra. Unlike voxel-based representations, tetrahedra can accurately describe the expected smooth surfaces of medical objects. Furthermore, the interior of such objects can be represented using a small number of tetrahedra. Rather than describing a medical object using tens of thousands of voxels, our representations generally contain only a few thousand elements. Tetrahedra facilitate the creation of efficient non-rigid registration algorithms based on finite element methods (FEM). We create a fast, FEM-based method to non-rigidly register segmented anatomical structures from two subjects. Using our compact tetrahedral representations, this method generally requires less than one minute of processing time on a desktop PC. We also create a novel method for the non-rigid registration of gray scale images. To facilitate a fast method, we create a tetrahedral representation of a displacement field that automatically adapts to both the anatomy in an image and to the displacement field. The resulting algorithm has a computational cost that is dominated by the number of nodes in the mesh (about 10,000), rather than the number of voxels in an image (nearly 10,000,000). For many non-rigid registration problems, we can find a transformation from one image to another in five minutes. This speed is important as it allows use of the algorithm during surgery. We apply our algorithms to find correlations between the shape of anatomical structures and the presence of schizophrenia. We show that a study based on our representations outperforms studies based on other representations. We also use the results of our non-rigid registration algorithm as the basis of a segmentation algorithm. That algorithm also outperforms other methods in our tests, producing smoother segmentations and more accurately reproducing manual segmentations.

Image Segmentation Using Information Bottleneck Method

In image processing, segmentation algorithms constitute one of the main focuses of research. In this paper, new image segmentation algorithms based on a hard version of the information bottleneck method are presented. The objective of this method is to extract a compact representation of a variable, considered the input, with minimal loss of mutual information with respect to another variable, considered the output. First, we introduce a split-and-merge algorithm based on the definition of an information channel between a set of regions (input) of the image and the intensity histogram bins (output). From this channel, the maximization of the mutual information gain is used to optimize the image partitioning. Then, the merging process of the regions obtained in the previous phase is carried out by minimizing the loss of mutual information. From the inversion of the above channel, we also present a new histogram clustering algorithm based on the minimization of the mutual information loss, where now the input variable represents the histogram bins and the output is given by the set of regions obtained from the above split-and-merge algorithm. Finally, we introduce two new clustering algorithms which show how the information bottleneck method can be applied to the registration channel obtained when two multimodal images are correctly aligned. Different experiments on 2-D and 3-D images show the behavior of the proposed algorithms

A fast two-grid and finite section method for a class of integral equations on the real line with application to an acoustic scattering problem in the half-plane

We consider the numerical treatment of second kind integral equations on the real line of the form ∅(s) = ∫_(-∞)^(+∞)▒〖κ(s-t)z(t)ϕ(t)dt,s=R〗 (abbreviated ϕ= ψ+K_z ϕ) in which K ϵ L_1 (R), z ϵ L_∞ (R) and ψ ϵ BC(R), the space of bounded continuous functions on R, are assumed known and ϕ ϵ BC(R) is to be determined. We first derive sharp error estimates for the finite section approximation (reducing the range of integration to [-A, A]) via bounds on (1-K_z )^(-1)as an operator on spaces of weighted continuous functions. Numerical solution by a simple discrete collocation method on a uniform grid on R is then analysed: in the case when z is compactly supported this leads to a coefficient matrix which allows a rapid matrix-vector multiply via the FFT. To utilise this possibility we propose a modified two-grid iteration, a feature of which is that the coarse grid matrix is approximated by a banded matrix, and analyse convergence and computational cost. In cases where z is not compactly supported a combined finite section and two-grid algorithm can be applied and we extend the analysis to this case. As an application we consider acoustic scattering in the half-plane with a Robin or impedance boundary condition which we formulate as a boundary integral equation of the class studied. Our final result is that if z (related to the boundary impedance in the application) takes values in an appropriate compact subset Q of the complex plane, then the difference between ϕ(s)and its finite section approximation computed numerically using the iterative scheme proposed is ≤C_1 [kh log⁡〖(1⁄kh)+(1-Θ)^((-1)⁄2) (kA)^((-1)⁄2) 〗 ] in the interval [-ΘA,ΘA](Θ<1) for kh sufficiently small, where k is the wavenumber and h the grid spacing. Moreover this numerical approximation can be computed in ≤C_2 N log⁡N operations, where N = 2A/h is the number of degrees of freedom. The values of the constants C1 and C2 depend only on the set Q and not on the wavenumber k or the support of z.

On asymptotic behavior at infinity and the finite section method for integral equations on the half-line

e consider integral equations on the half-line of the form and the finite section approximation to x obtained by replacing the infinite limit of integration by the finite limit β. We establish conditions under which, if the finite section method is stable for the original integral equation (i.e. exists and is uniformly bounded in the space of bounded continuous functions for all sufficiently large β), then it is stable also for a perturbed equation in which the kernel k is replaced by k + h. The class of perturbations allowed includes all compact and some non-compact perturbations of the integral operator. Using this result we study the stability and convergence of the finite section method in the space of continuous functions x for which ()()()=−∫∞dttxt,sk)s(x0()syβxβx()sxsp+1 is bounded. With the additional assumption that ()(tskt,sk−≤ where ()()(),qsomefor,sassOskandRLkq11>+∞→=∈− we show that the finite-section method is stable in the weighted space for ,qp≤≤0 provided it is stable on the space of bounded continuous functions. With these results we establish error bounds in weighted spaces for x - xβ and precise information on the asymptotic behaviour at infinity of x. We consider in particular the case when the integral operator is a perturbation of a Wiener-Hopf operator and illustrate this case with a Wiener-Hopf integral equation arising in acoustics.

A new method for the characterisation and verification of local spatial predictability for convective scale ensembles

The use of kilometre-scale ensembles in operational forecasting provides new challenges for forecast interpretation and evaluation to account for uncertainty on the convective scale. A new neighbourhood based method is presented for evaluating and characterising the local predictability variations from convective scale ensembles. Spatial scales over which ensemble forecasts agree (agreement scales, S^A) are calculated at each grid point ij, providing a map of the spatial agreement between forecasts. By comparing the average agreement scale obtained from ensemble member pairs (S^A(mm)_ij), with that between members and radar observations (S^A(mo)_ij), this approach allows the location-dependent spatial spread-skill relationship of the ensemble to be assessed. The properties of the agreement scales are demonstrated using an idealised experiment. To demonstrate the methods in an operational context the S^A(mm)_ij and S^A(mo)_ij are calculated for six convective cases run with the Met Office UK Ensemble Prediction System. The S^A(mm)_ij highlight predictability differences between cases, which can be linked to physical processes. Maps of S^A(mm)_ij are found to summarise the spatial predictability in a compact and physically meaningful manner that is useful for forecasting and for model interpretation. Comparison of S^A(mm)_ij and S^A(mo)_ij demonstrates the case-by-case and temporal variability of the spatial spread-skill, which can again be linked to physical processes.

Role of Polyelectrolyte for Layer-by-Layer Compact TiO(2) Films in Efficiency Enhanced Dye-Sensitized Solar Cells

Efficient compact TiO(2) films using different polyeleetrolytes are prepared by the layer-by-layer technique (LbL) and applied as an effective contact and blocking film in dye-sensitized solar cells (DSCs). The polyanion thermal stability plays a major role on the compact layers, which decreases back electron transfer processes and current losses at the FTO/TiO(2) interface. FESEM images show that polyelectrolytes such is sodium sullonated polystyrene (PSS) and sulfonated lignin (SE), in comparison to poly(acrylic acid) (FAA), ensure an adequate morphology for the LbL TiO(2) layer deposited before the mesoporous film, even triter the sintering step at 450 degrees C. The so treated photoanode in DSCs leads to a 30% improvement On the overall conversion efficiency. Electrochemical impedance spectroscopy (EIS) is employed to ascertain the role of die compact films with such polyelectrolytes. The significant increase in V(oc) of the solar cells with adequate polyelectrolytes in the LbL TiO(2) films shows their pivotal role in decreasing the electron recombination at the FTO surface and enhancing the electrical contact of FTO with the mesoporous TiO(2) layer.

A grid-based index method for time warping distance

Recently DTW (dynamic time warping) has been recognized as the most robust distance function to measure the similarity between two time series, and this fact has spawned a flurry of research on this topic. Most indexing methods proposed for DTW are based on the R-tree structure. Because of high dimensionality and loose lower bounds for time warping distance, the pruning power of these tree structures are quite weak, resulting in inefficient search. In this paper, we propose a dimensionality reduction method motivated by observations about the inherent character of each time series. A very compact index file is constructed. By scanning the index file, we can get a very small candidate set, so that the number of page access is dramatically reduced. We demonstrate the effectiveness of our approach on real and synthetic datasets.

DDR: an index method for large time-series datasets

The tree index structure is a traditional method for searching similar data in large datasets. It is based on the presupposition that most sub-trees are pruned in the searching process. As a result, the number of page accesses is reduced. However, time-series datasets generally have a very high dimensionality. Because of the so-called dimensionality curse, the pruning effectiveness is reduced in high dimensionality. Consequently, the tree index structure is not a suitable method for time-series datasets. In this paper, we propose a two-phase (filtering and refinement) method for searching time-series datasets. In the filtering step, a quantizing time-series is used to construct a compact file which is scanned for filtering out irrelevant. A small set of candidates is translated to the second step for refinement. In this step, we introduce an effective index compression method named grid-based datawise dimensionality reduction (DRR) which attempts to preserve the characteristics of the time-series. An experimental comparison with existing techniques demonstrates the utility of our approach.

Compact stacked planar inverted-F antenna for passive deep brain stimulation implants

A compact meandered three-layer stacked circular planar inverted-F antenna is designed and simulated at the UHF band (902.75 – 927.25 MHz) for passive deep brain stimulation implants. The UHF band is used because it offers small antenna size, and high data rate. The top and middle radiating layers are meandered, and low cost substrate and superstrate materials are used to limit the radius and height of the antenna to 5 mm and 1.64 mm, respectively. A dielectric substrate of FR-4 of εr= 4.7 and δ= 0.018, and a biocompatible superstrate of silicone of er= 3.7 and d= 0.003 with thickness of 0.2 mm are used in the design. The resonance frequency of the proposed antenna is 918 MHz with a bandwidth of 24 MHz at return loss of −10 dB in free space. The antenna parameter such as 3D gain pattern of the designed antenna within a skin-tissue model is evaluated by using the finite element method. The compactness, wide bandwidth, round shape, and stable characteristics in skin make this antenna suitable for DBS. The feasibility of the wireless power transmission to the implant in the human head is also examined.

Investigation of flat spectrum radio sources by the interplanetary scintillation method at 111 MHz

Abstract. Interplanetary scintillation observations of 48 of the 55 Augusto et al. (1998) flat spectrum radio sources were carried out at 111 MHz using the interplanetary scintillation method on the Large Phased Array (LPA) in Russia. Due to the large size of the LPA beam (1◦ × 0.5◦) a careful inspection of all possible confusion sources was made using extant large radio surveys: 37 of the 48 sources are not confused. We were able to estimate the scintillating flux densities of 13 sources, getting upper limits for the remaining 35. Gathering more or improving extant VLBI data on these sources might significantly improve our results. This proof-of-concept project tells us that compact (<1 ) flat spectrum radio sources show strong enough scintillations at 111 MHz to establish/constrain their spectra (low-frequency end). Key words. galaxies: general – galaxies: active – galaxies: quasars: general

The influence of preparation method of OC1OC6-PPV films on the photo-oxidation process

Poly(p-phenylene vinylene) (PPV) derivatives are well known for their applications in polymer light emitting diodes (PLEDs). PPV derivatives are highly susceptible to photo-oxidation though, which is mainly caused by the scission of the vinyl double bond on the polymer backbone. In this work, we show that Langmuir-Blodgett (LB) films are less degraded than cast films of a PPV derivative (OC1OC6-PPV). Both films had similar thickness (similar to 50 nm) to allow for a more realistic comparison. Degradation was monitored with UV-vis and FTIR spectroscopies. The results indicated that cast films were completely degraded in ca. 400 min, while LB took longer time, i.e. about four times the values for the cast films. The differences can be attributed to the more compact morphology in the LB than in the cast films. With a compact morphology the diffusion of oxygen in the LB film is hampered and this causes a delay in the degradation process. (c) 2006 Elsevier Ltd. All rights reserved.

A case study of compact cogeneration using various fuels

A method for studying the technical and economic feasibility of absorption refrigeration systems in compact cogenerators is presented. The system studied consists of an internal combustion engine, an electric generator and a heat exchanger to recover residual heat from the refrigeration water and exhaust gases. As an application, a computer program simulates the cogeneration system in a building which already has 75 kW of installed electric power. The maximum electric and refrigeration demands are 45 kW and 76 kW respectively. This study simulates the system performance, utilizing diesel oil, sugar cane alcohol and natural gas as possible fuels. (C) 1997 Elsevier B.V. Ltd.

Effective field theory methods to model compact binaries

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

AÇÃO REGULATÓRIA DO GnRH NO DESENVOLVIMENTO EMBRIONÁRIO PRECOCE E VITRIFICAÇÃO DE EMBRIÕES SUÍNOS PELO MÉTODO DE MICROGOTA (ROLE OF GnRH DURING EARLY EMBRYONIC DEVELOPMENT AND DEVELOPMENT OF SWINE EMBRYOS FOLLOWING VITRIFICATION BY MICRODROPLET METHOD)

The objective of this study was to investigate the role of GnRH on the preimplantation development of mouse embryos in vitro. GnRH-I, GnRH-II, and GnRH agonists: Des-Gly, Des-Trp and histrelin did not improve embryo development. However, treatment with the specific GnRH antagonist SB-75 blocked embryo development at morula stage. The inhibition of embryo development by SB-75 could be rescued by the addition of histrelin. To determine which intracellular signaling cascade is involved following binding of GnRH to the GnRHR, embryos were cultured in the presence of specific PKC (GFX) or PKA (SQ22536) inhibitors. The PKC inhibitor blocked embryo development at a similar stage as SB-75, whereas SQ22536 had an inhibitory effect, diminishing blastocyst formation and hatched rates. There are evidences that GnRH has an essential autocrine effect on mouse embryonic development via GnRHR, probably by activating PKC signaling cascade while the inhibition of the GnRH signaling does not activate apoptotic mechanisms involving caspase-3. In another experiment, development in vitro of embryos from Chinese Meishan (M) and occidental white crossbred (WC) females were investigated after improving the vitrification protocol for pig embryos. Efficient cryopreservation of zona pellucida-intact porcine embryos and studies of the difference among breeds could greatly impact the swine industry. The percentage of embryos surviving 24 h after cryopreservation without lysis or degeneration was higher for M (72%) than WC (44%). However, in vitro development of embryos that survived cryopreservation was not different between M and WC at the expanded (64%) or hatched (22%) blastocyst stages. Developmental rates were significantly higher for control embryos than frozen embryos from both breeds at expanded blastocyst stage, but not at hatched blastocyst stage. Rates of expanded blastocyst formation did not differ between M and WC control embryos (98 and 95%, respectively). With a new procedure to warm vitrified pig embryos, the survival rates may be improved. The optimal stages to vitrify pig embryos using the microdroplet method ranges from late compact morula to early expanded blastocyst. The results suggest that M embryos have a higher capacity to survive the vitrification process than WC embryos. O objetivo do presente estudo foi investigar a importância do GnRH no desenvolvimento embrionário precoce em camundongos. GnRH-I, GnRH-II e os GnRH agonistas: Des-Gly, Des-Trp e histrelina não incrementaram o desenvolvimento embrionário. Entretanto, o tratamento com SB-75, um antagonista específico do GnRH, bloqueou o desenvolvimento embrionário no estádio de mórula. A inibição do desenvolvimento embrionário pelo SB-75 pôde ser revertida com a adição de histrelina. Para determinar a cascata do sinal intracelular desencadeada pela ligação do GnRH com o seu receptor, embriões foram cultivados na presença de inibidores específicos da PKC (GFX) e da PKA (SQ22536). O inibidor da PKC bloqueou o desenvolvimento embrionário em estádio similar ao bloqueio mediado pelo SB- 75, enquanto o SQ22536 teve efeito inibitório diminuindo a formação de blastocisto e taxas de eclosão. Os resultados sugerem que o GnRH tem um efeito autócrino essencial no desenvolvimento embrionário através do GnRHR, provavelmente, ativando a cascata da PKC. Por outro lado, a inibição do sinal do GnRH não ativa mecanismos apoptóticos que involvam caspase-3. Em outro experimento, foi investigado o desenvolvimento in vitro de embriões da raça Meishan (M) e branco cruzado (WC) após vitrificação pelo método microgota. O desenvolvimento de protocolos eficientes para criopreservação de embriões suínos com a zona pelúcida intacta e a avaliação das diferenças entre raças pode ter um significativo impacto na suinocultura. A percentagem de embriões que sobreviveram à criopreservação depois de 24 h foi maior na M (72%) do que na WC (44%). No entanto, o desenvolvimento in vitro dos embriões que sobreviveram à criopreservação não foi diferente entre M e WC nos estádios de blastocisto expandido (64%) ou eclodido (22%). Os índices de desenvolvimento foram significativamente mais altos para os embriões controle do que para os embriões vitrificados nas duas raças no estádio de blastocisto expandido, porém não foram diferentes para o estádio de blastocisto eclodido. A formação de blastocisto expandido não diferiu entre os embriões controle M e WC (98 e 95%, respectivamente). Com o novo procedimento (“hot warm”) para descongelar embriões vitrificados pelo método de microgota, pode-se aumentar dos índices de sobrevivência. Os melhores estádios embrionários para a vitrificação de embriões suínos variam de mórula compacta tardia até blastocisto expandido inicial. Os resultados sugerem que embriões M têm mais capacidade de sobreviver ao processo de vitrificação do que embriões WC.

MODEL PARTITIONS AND COMPACT TEST CASE SUITES

We present a generalized test case generation method, called the G method. Although inspired by the W method, the G method, in contrast, allows for test case suite generation even in the absence of characterization sets for the specification models. Instead, the G method relies on knowledge about the index of certain equivalences induced at the implementation models. We show that the W method can be derived from the G method as a particular case. Moreover, we discuss some naturally occurring infinite classes of FSM models over which the G method generates test suites that are exponentially more compact than those produced by the W method.