Search

Whole brain neuroanatomy using tera-voxel light-microscopic data sets is of much current interest. A fundamental problem in this field is the mapping of individual brain data sets to a reference space. Previous work has not rigorously quantified the distortions in brain geometry from in-vivo to ex-vivo brains due to the tissue processing, which will be important when computing properties such as local cell and process densities at the voxel level in creating reference brain maps. Further, existing approaches focus on registering uni-modal volumetric data; however, given the increasing interest in the marmoset model for neuroscience research, it is necessary to cross-register multi-modal data sets including MRIs and multiple histological series that can help address individual variations in brain architecture. Here we present a computational approach for same-subject multimodal MRI guided reconstruction of a histological series, jointly with diffeomorphic mapping to a reference atlas. We quantify the scale change during the different stages of histological processing of the brains using the Jacobian determinant of the diffeomorphic transformations involved. There are two major steps in the histology process with associated scale distortions (a) brain perfusion (b) histological sectioning and reassembly. By mapping the final image stacks to the ex-vivo post fixation MRI, we show that tape-transfer histology can be reassembled accurately into 3D volumes with a local scale change of 2.0 $\pm$ 0.4% per axis dimension. In contrast, the perfusion step, as assessed by mapping the in-vivo MRIs to the ex-vivo post fixation MRIs, shows a larger local scale change of 6.9 $\pm$ 2.1% per axis dimension. This is the first systematic quantification of the local metric distortions associated with whole-brain histological processing, and we expect that the results will generalize to other species.

This paper presents a variational framework for dense diffeomorphic atlas-mapping onto high-throughput histology stacks at the 20 um meso-scale. The observed sections are modelled as Gaussian random fields conditioned on a sequence of unknown section by section rigid motions and unknown diffeomorphic transformation of a three-dimensional atlas. To regularize over the high-dimensionality of our parameter space (which is a product space of the rigid motion dimensions and the diffeomorphism dimensions), the histology stacks are modelled as arising from a first order Sobolev space smoothness prior. We show that the joint maximum a-posteriori, penalized-likelihood estimator of our high dimensional parameter space emerges as a joint optimization interleaving rigid motion estimation for histology restacking and large deformation diffeomorphic metric mapping to atlas coordinates. We show that joint optimization in this parameter space solves the classical curvature non-identifiability of the histology stacking problem. The algorithms are demonstrated on a collection of whole-brain histological image stacks from the Mouse Brain Architecture Project.

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects., Comment: 17 pages, 13 figures, accepted by ApJ

We identify galaxy groups and clusters in volume-limited samples of the SDSS redshift survey, using a redshift-space friends-of-friends algorithm. We optimize the friends-of-friends linking lengths to recover galaxy systems that occupy the same dark matter halos, using a set of mock catalogs created by populating halos of N-body simulations with galaxies. Extensive tests with these mock catalogs show that no combination of perpendicular and line-of-sight linking lengths is able to yield groups and clusters that simultaneously recover the true halo multiplicity function, projected size distribution, and velocity dispersion. We adopt a linking length combination that yields, for galaxy groups with ten or more members: a group multiplicity function that is unbiased with respect to the true halo multiplicity function; an unbiased median relation between the multiplicities of groups and their associated halos; a spurious group fraction of less than ~1%; a halo completeness of more than ~97%; the correct projected size distribution as a function of multiplicity; and a velocity dispersion distribution that is ~20% too low at all multiplicities. These results hold over a range of mock catalogs that use different input recipes of populating halos with galaxies. We apply our group-finding algorithm to the SDSS data and obtain three group and cluster catalogs for three volume-limited samples that cover 3495.1 square degrees on the sky. We correct for incompleteness caused by fiber collisions and survey edges, and obtain measurements of the group multiplicity function, with errors calculated from realistic mock catalogs. These multiplicity function measurements provide a key constraint on the relation between galaxy populations and dark matter halos., Comment: 26 emulateapj pages including 19 figures. Replaced with final ApJ version

We present a new catalog of isolated galaxies obtained through an automated systematic search. These 2980 isolated galaxies were found in approximately 2099 sq deg of sky in the Sloan Digital Sky Survey Data Release 1 (SDSS DR1) photometry. The selection algorithm, implementing a variation on the criteria developed by Karachentseva in 1973, proved to be very efficient and fast. This catalog will be useful for studies of the general galaxy characteristics. Here we report on our results., Comment: 67 pages, which includes 14 figures. Accepted for publication by AJ

We identify close pairs of galaxies from 278 deg^2 of Sloan Digital Sky Survey commissioning imaging data. The pairs are drawn from a sample of 330,041 galaxies with 18 < r^* < 20. We determine the angular correlation function of galaxy pairs, and find it to be stronger than the correlation function of single galaxies by a factor of 2.9 +/- 0.4. The two correlation functions have the same logarithmic slope of 0.77. We invert Limber's equation to estimate the three-dimensional correlation functions; we find clustering lengths of r_0= 4.2 +/- 0.4 h^{-1} Mpc for galaxies and 7.8 +/- 0.7 h^{-1} Mpc for galaxy pairs. These results agree well with the global richness dependence of the correlation functions of galaxy systems., Comment: 12 pages. ApJ, in press

The discovery of an optical counterpart to GRB010222 (detected by BeppoSAX; Piro 2001) was announced 4.4 hrs after the burst by Henden (2001a). The Sloan Digital Sky Survey's 0.5m photometric telescope (PT) and 2.5m survey telescope were used to observe the afterglow of GRB010222 starting 4.8 hours after the GRB. The 0.5m PT observed the afterglow in five, 300 sec g' band exposures over the course of half an hour, measuring a temporal decay rate in this short period of F_nu \propto t^{-1.0+/-0.5}. The 2.5m camera imaged the counterpart nearly simultaneously in five filters (u' g' r' i' z'), with r' = 18.74+/-0.02 at 12:10 UT. These multicolor observations, corrected for reddening and the afterglow's temporal decay, are well fit by the power-law F_nu \propto nu^{-0.90+/-0.03} with the exception of the u' band UV flux which is 20% below this slope. We examine possible interpretations of this spectral shape, including source extinction in a star forming region., Comment: 8 pages, 4 figures, accepted for publication in ApJ. Two figures added, minor changes to text in this draft. Related material can be found at: http://sdss.fnal.gov:8000/grb/