A novel cell nuclei segmentation method for 3D C. elegans embryonic time-lapse images
© Chen et al.; licensee BioMed Central Ltd. 2013
Received: 19 June 2013
Accepted: 8 November 2013
Published: 19 November 2013
Recently a series of algorithms have been developed, providing automatic tools for tracing C. elegans embryonic cell lineage. In these algorithms, 3D images collected from a confocal laser scanning microscope were processed, the output of which is cell lineage with cell division history and cell positions with time. However, current image segmentation algorithms suffer from high error rate especially after 350-cell stage because of low signal-noise ratio as well as low resolution along the Z axis (0.5-1 microns). As a result, correction of the errors becomes a huge burden. These errors are mainly produced in the segmentation of nuclei. Thus development of a more accurate image segmentation algorithm will alleviate the hurdle for automated analysis of cell lineage.
This paper presents a new type of nuclei segmentation method embracing an bi-directional prediction procedure, which can greatly reduce the number of false negative errors, the most common errors in the previous segmentation. In this method, we first use a 2D region growing technique together with the level-set method to generate accurate 2D slices. Then a modified gradient method instead of the existing 3D local maximum method is adopted to detect all the 2D slices located in the nuclei center, each of which corresponds to one nucleus. Finally, the bi-directional pred- iction method based on the images before and after the current time point is introduced into the system to predict the nuclei in low quality parts of the images. The result of our method shows a notable improvement in the accuracy rate. For each nucleus, its precise location, volume and gene expression value (gray value) is also obtained, all of which will be useful in further downstream analyses.
The result of this research demonstrates the advantages of the bi-directional prediction method in the nuclei segmentation over that of StarryNite/MatLab StarryNite. Several other modifications adopted in our nuclei segmentation system are also discussed.
The development of the live-cell imaging microscopy and fluorescent tagging provides us unprecedented opportunity to observe gene expression, nuclei movement and nuclei division processes during embryogenesis at the single cell level [1-3]. Customized algorithms were developed using high temporal resolution of image stacks to automatically trace cell divisions during animal development.
One system named StarryNite has been developed for automated cell lineage tracing and gene expression profiling during C.elegans embryogenesis [4, 5]. A lineage tree can be generated with StarryNite automatically based on the 3D time-lapse image of a developing C. elegans embryo. To obtain cell nuclei with high signal/noise ratio, the C. elegans embryos are first tagged with a fluorescent protein. Then the authors took images every 60 or 90 seconds with a confocal microscope during C. elegans embryogenesis. At every time point there are 41 image planes and the resolution in the Z-axis is 0.71 microns. A complete 3D time-lapse image series contains 180-240 time points when the cell number of the embryo is over 500. So there are over 6000 images in a single image series, which record approximately 600 cell divisions and 50000 nuclei at different time points. With these images as input to StarryNite system, information on reporter expression, nucleus division and movement can be efficiently extracted from a huge mass of image data. [Additional file 1: Table S1].
However, due to the low resolution in the Z-axis and high noise in the 3D time-lapse image of C. elegans, StarryNite/MatLab StarryNite output suffers from high error rate especially during late embryogenesis. It takes around 2-5 human hours for manually annotating StarryNite’s output up to 350 cell stage, but for most cells, at least one division process cannot be correctly traced in this kind of lineage tree [6-8].
StarryNite/MatLab StarryNite can be generally divided into two parts, namely, the nuclei segmentation part and the tracing part. The nuclei segmentation part will analyze the 3D time-lapse images of C. elegans and detect the nuclei positions, size and intensity at every time point. Then the tracing part will build the lineage tree according to the nuclei information generated by the nuclei segmentation part. The main reason for StarryNite/MatLab StarryNite not being able to reach the 550 cell stage, with an acceptable error rate, is that too many cells could not be detected by the nuclei segmentation algorithm that works relatively well at or before the 350 cell stage. In fact, the detection accuracy of the nuclei segmentation part in StarryNite can only work for up to 350 cells at any one single time point (MatLab StarryNite, 500 cells). Therefore at the end of embryogenesis (cell number at 550), more than 200 cells cannot be detected correctly. Therefore, for every 3D image data set, even if the tracing part is error free, the number of errors would still be more than 7000, which means that one to two weeks of manual editing time are needed for generation of a 550-celled lineage. Therefore, development of a novel or optimization of the existing nuclei segmentation algorithm is the key to speed up the automatic lineaging throughout C. elegans embryogenesis.
In order to improve the performance of StarryNite, several techniques have been developed and implemented into StarryNite. Santella et al.  presents a blob-slice nuclei segmentation model to improve the correct nuclei detection rate and this model can be also used to analyze 3D images of other species, such as zebrafish and Drosopila . Aydin builds a new tracing system with a SVM classifier, which could greatly reduce the number of false positive errors that can be detrimental during the editing of the lineage tree . Richards et al.  increase the 3D time-lapse image resolution in the Z-axis with a resonance-scanning confocal microscope . In his work, the lineage tree generated by Matlab StarryNite with these new images could reach the 500 cell stage with manual editing. In all the above-mentioned methods, the 3D local maximum is a must for the nuclei segmentation, based on the principle that one particular 3D local maximum point in images corresponds to a certain nucleus. However, as the number of images increases, nuclei at the top of the embryo will be flattened and crowded. Consequently, one 3D local maximum point in images at this stage cannot always indicate one nucleus but two (or more) nuclei. Besides, in 3D local maximum, there are several other techniques available for 3D nuclei segmentation such as adaptive thresholds, [12, 13] mode finding, [14, 15] gradient flow tracking,  watershed [17, 18] and level set methods . These methods perform well when nuclei are widely spaced, but given conditions in which nuclei are crowded, as in the low resolution along the Z axis in the 3D image, all these methods show some disfigurement to various degrees.
Judging from these aspects, we speculate that under certain circumstances, all the above methods have difficulty in finding certain nuclei at some time points due to the low image quality. In view of this issue, we develop a technique called bi-directional prediction to reduce false negative errors and improve the precision and accuracy of nuclei segmentation. Due to the minimum position shift of a nucleus between two adjacent time points, if there is a nucleus that can hardly be detected based on the image of one single time point, we can predict it by studying the images before and after the current time point.
In this paper, a new nuclei segmentation method is attempted to substitute that for StarryNite, which enables a more accurate extraction of useful information from 3D time-lapse images. The results achieved by our new nuclei segmentation method demonstrate better performance in nuclei segmentation and reduced false negative errors. Currently, editing a lineage up to 350 or 450 cell stage requires roughly 0.5 or 8 hours respectively using our image series.
Error rate analysis
Error analysis of StarryNite
Total cell number
The main source of nuclei segmentation errors includes failure to identify nuclei (false negatives), misjudgments of non-existent nuclei (false positives), false diameter errors and dislocation errors. As illustrated in Table 1, false negative errors are the most common errors in the StarryNite/MatLab StarryNite output. What is more important is that nuclei segmentation errors will propagate to tracing errors. In fact, the majority of tracing errors in the StarryNite/MatLab StarryNite output which do not detect matching errors and false division identification errors, are caused by nuclei segmentation errors like dislocations or false diameters. According to the above analysis, it is the nuclei segmentation errors, especially the false negative errors, that become the biggest stumbling barrier in using the StarryNite/MatLab StarryNite system. For this reason we have built a new nuclei segmentation system which focuses on eliminating false negative errors.
There are three reasons why the multiplication of false negative errors occurs at the end of the embryonic development. First, the image resolution in the Z-axis is 0.71 microns, which means that one small nucleus (with nucleus diameter < 1.8 microns) can be clearly seen in only two images. Second, the increase in nuclei numbers makes the nuclei at the top of the embryos become flattened and crowded to such an extent that they are not always recovered by the iterative nuclei identification.4 This crowding becomes a significant source of errors after the 150-cell stage, when the nuclei first enter the very top of the embryo. Third, neighboring cells may have different levels of GFP expression. When two nuclei are close to each other, the local signal density of the weaker one may blend into the edge of the high signal from the brighter one and, thus, no longer obey the rule of local maximum. Accordingly, we designed our nuclei segmentation system for the purpose of avoiding false negative errors from the above three original sources of error formation.
Step1: Pre-processing with level-set method (LSM)
There are several limitations in the 3D time-lapse images of C. elegans embryo: First, the image resolution in the Z-axis is 0.71 microns, while the diameter of some nuclei at the 450 cell stage is around 1.8 microns, so each of these small nuclei can only be mapped in one or two images. These defects cause extraordinary difficulties in nuclei detection. Second, the GFP expression of some nuclei are quite low, so distinguishing the noises of the original images from a low GFP expression signal is a major issue. Third, the nuclei boundaries are not clear, because of the crowed nuclei at the 450 cell stage, thus a boundary enhancement method is needed for an accurate segmentation. In our work we overcome the above difficulties by pre-processing the images with the level-set method before segmentation.
Error analysis of Matlab StarryNite
Total cell number
Error analysis of before and after level set processing
130108NHR25p1 (before Level Set)
130108NHR25p1 (after Level Set)
Step2: 2D segmentation using the modified region growing method
In order to use the region growing method, the first and foremost step is to predetermine the initial ‘seed’ points, each of which will then grow into one 2D nuclei slice after being processed by a region growing algorithm. The method we use to search seed points is the 2D local maximization algorithm. Following our goal of obtaining precise segmentation of nucleus slices, we improve the design of the region growing algorithm in our nuclei segmentation system, so that it can perform well even under the condition in which the nuclei are crowded together. The traditional region growing algorithm is an iterative process, in which neighboring pixels of initial “seed points” are examined and one fixed threshold is used to determine whether the pixel neighbors should be added to a certain region. In our modified region growing method, the value of the threshold is determined by the distance between each pixel and seed point respectively, as a precaution against nuclei with lower expression values being overshadowed by brighter neighboring ones when they are overcrowded.
According to the nuclei diameter (1.8-7 microns) and the image resolution in the Z-axis (0.355 microns now, improved by Level Set), every nuclei will correspond to 4-18 2D nuclei slices, so even if some 2D nuclei slices of one nucleus cannot be recognized in some Z-planes, we can still identify them by the 2D local maximum points in the other Z-plane images. However, StarryNite/MatLab StarryNite uses only the method in which each 3D local maximum point is considered as one nucleus, therefore, one nucleus will be left out in the case where the corresponding 3D local maximum point is omitted. Thus it can be said in this sense that our method based on 2D nuclei slices is fault tolerant, which can reduce the false negative errors significantly.
Step3: Nuclei generation
where is the mean gradient of slice n on image Z, and is the mean gradient of slice n on image Z-1, the z-axial position of which is located just below image Z.
Due to the fact that the center slices contain data of the center position of each nucleus, we can work out which of the remaining slices are the constituent parts of one certain nucleus by calculating distances between the remaining slices and each center slice respectively. Each nucleus turns out to be composed of 4-18 slices in this case and information about exact position, radius, volume and brightness of every nucleus is achieved through this procedure.
Step4: Bi-directional prediction
After step 3, owing to the poor quality and low Z-axial resolution of images, some nuclei segmentation errors especially the false negative errors surface inevitably if we do nuclei segmentation based only on images of one single time point. During the process of embryogenesis, every new nucleus in the lineage tree comes from cell division. These new nuclei will divide into other new nuclei or die, in other words, the same nuclei may not be captured via single-time images. Further, the shift distance of every nucleus between two adjacent time points is very short, and the nuclei size will not change much within short-time intervals, so we can refer to the image before or after the current time point to forecast a nucleus in order to avoid the failure of capturing it at one single time point. This process effectively decreases the emergence of false negative errors. Considering the possibility that one nucleus may divide into two nuclei during embryogenesis, a bi-directional prediction instead of only a forward or backward prediction is necessary and effective in nuclei prediction.
where A and B are the symbolic representation of the nuclei; DS AE is the the dissimilarity degree between nuclei A and B; d is the shift distance between A and B; r A , r B are the radii of A and B respectively; and V A and V B are the volume of A and B respectively.
Both false positive errors and false negative errors can be decreased significantly through our bi-directional prediction method. Bi-directional prediction is capable of dealing with some problems that cannot be solved by simply improving image segmentation, and this is the most significant innovation of our method.
Step5: Nuclei filtering
Through the above processing, all the possible nuclei for each time point are recognized and saved in one nuclei list. However, at several consecutive time points, some false positive errors which cannot be eliminated by our bi-directional prediction method, still remain, so we need a filter to deal with this kind of error.
Then we select nucleus A or B to be deleted based on their volume and gray values.
Finally, an optimal nuclei list is obtained, representing the expected properties including nuclei ID, 3D center coordinate, nucleus volume, radius, 2D slice and gray value. With this list we can generate the information for lineage tracing. In fact for the existing tracing method, the necessary information is the 3D coordinate of the nuclei center and gray value. Other very important information can be used to improve the tracing algorithm in StarryNite.
We used our method on data sets obtained from our Leica SP5 microscope. [Additional files 2, 3, and 4] We also applied our nuclei segmentation results to the generation of a lineage tree with the tracing algorithm in StarryNite. The performance of the nuclei segmentation method can be evaluated by the following aspects: Nuclei number, accuracy rate and editing time needed for building a lineage tree.
All our data sets are sampled with 90 seconds time resolution, but a better result will be reached if the time resolution is increased to 60 seconds as in the image protocol in Aydin . This is because the nuclei will become shorter as the time resolution is increased from 90 seconds to 60 seconds, which will greatly reduce the difficulty in the following bi-directional prediction.
Editing time needed for one correct lineage tree is another important performance criterion. In Bao , StarryNite requires 4-10 hours to edit the tree at the 350 cell stage, and due to the aforementioned maximum-number limit of nuclei by the nuclei segmentation part in StarryNite, it is quite hard to edit the lineage tree beyond the 350 cell stage. One lineage tree, which was edited to the 450 cell stage using StarryNite, required more than one week of manual editing . Santella’s 2D blob-slice method performed much better than StarryNite, Richards  used Santella’s method and could reach the 550 cell stage with 8-16 hours manually editing. Their images have a double high resolution in the Z direction, that is there are 60 images in every time point, nearly twice the number of images used in our case. As the number of images double, the storage space required for image and computation time also doubles as a consequence. We performed our method on data sets 130108NHR25p1 and 130108PHA4p1, which requires only 30 minutes for a 350 cells lineage tree, and eight hours for a 450 cell lineage tree.
There have been several methods for automatic cell lineage tree generation. StarryNite  together with a tool named AceTree  provided an efficient system for cell recognition and annotation. Santella  and Richard  improved the nuclei segmentation with a different method. Santella  revised the nuclei segmentation with a 2D blob-slice method, and Richard  advanced the output of the tracing procedure by increasing the image resolution of the Z-axis with a new type of resonance-scanning microscopy. Due to poor image quality and the rapid multiplying of nuclei after the 350 cell stage, there are substantial nuclei segmentation errors, especially false negative ones, causing large numbers of unavoidable follow-up errors and long time manual editing of the cell lineage. Therefore a bi-directional prediction technique was developed by us as an efficient way to alleviate these issues. At the same time, other methods such as level-set, region growing and a modified gradient algorithm has also been adopted for the purpose of improving our nuclei segmentation system.
In our work, instead of using single time point images for nuclei segmentation, we take advantage of the bi-directional prediction technique, analyzing the images before and after the current time point to improve nuclei identification. As the interval between two adjacent sampling times is only 50-90 seconds, the coordinates, size, and expression value of the same nucleus in two adjacent time points remain nearly unchanged. This feature significantly reduces the segmentation errors in several continuous time points. One omitted nucleus could certainly be discovered from nuclei segmentation in another time point. We also apply level-set, 2D region growing and the gradient method to our whole nuclei segmentation system. First, the level-set method boosts the image resolution in the Z-axis, and ensures the integrity of the edge information by smoothing the image under the curvature. Moreover, noises existing in the images also get reduced by the level-set method. Second, our 2D region growing method can achieve an accurate 2D slice segmentation. It is mainly because the threshold for our region growing method varies inversely with the size of the nuclei slice, which guarantees that the slices will not be eclipsed by a brighter neighbor. Third, the modified gradient method is an effective way to search for the center slices, each of which represents one exclusive nucleus. As demonstrated in experiment results, when the error rate is low, the outcome of our nuclei segmentation method is much better in nuclei reorganization. But due to the existence of errors caused by the tracing part of StarryNite, it is still impossible to get a perfectly accurate lineage without manual editing. Future research may include work on improving the performance of the tracing algorithm. In fact, the weakest part of the tracing algorithm lies in the cell-division detection, as the algorithm is mainly based on the assumption that the parent cell is located at the center position of two daughter cells. However, at the 350 cell stage, a satisfactory result cannot be achieved via the present day cell-division detection algorithm, since the space distribution pattern of a cell becomes more complicated. At the same time, nuclei segmentation errors are common, which indicates that there is still room for much improvement of the StarryNite method through the optimization of the tracing algorithm.
In this paper, we present a novel nuclei segmentation system which can be used for automatic cell lineage tracing. Several new designs are introduced using our method, including level-set, region-growing, gradient tracing and bi-directional prediction. Among these techniques, the bi-directional prediction, which is capable of detecting nuclei based on backward and forward information, provides a major advance in the nuclei segmentation algorithm. The results using our image data also demonstrate that the proposed algorithm is accurate and practical.
There are still two aspects which need to be noted. First, our level-set preprocessing is time-consuming. This problem can be solved in the future by using more efficient algorithms and more powerful computers. Second, at a very late stage (at the 500 cells stage), some small and non-spherical nuclei may appear, which are hard to distinguish with the limitations in the current image resolution.
Subsequent research may focus on the combination of nuclei segmentation and tracing procedure. A real-time lineage tracing system may then be achieved with this kind of architecture.
The data sets we used are obtained from a Leica SP5 confocal microscope using Hybrid detector, and images were taken every 90 seconds with the same x and y resolution as described previously (ref). The data sets contain 41 image planes in every time point, the distance between two adjacent images planes is 0.71 microns (Additional file 1: Table S1).
The level set method (LSM) is a kind of numerical technique which can be used for tracking interfaces and shapes. The advantage of the level set method is that one can perform numerical computations involving curves and surfaces on a fixed Cartesian grid without having to parameterize these objects (this is called the Eulerian approach) . The level set method provides us with an easy way to smooth the 3D time-lapse image of C. elegans without losing important image details. The method has laid a foundation for 2D slice segmentation by region growing.
We have used a Matlab toolbox from http://barissumengen.com/level_set_methods/. The iteration number is 6, and at the 500 cell stage, the iteration number is 4.
Region growing is a kind of pixel-based image segmentation method . This approach to segmentation examines neighboring pixels of initial “seed points” and determines whether the pixel neighbors should be added to the region depending on a region membership criterion. In our method, this kind of growing process is iterated on, and every seed point will grow into a nucleus slice. These nuclei slices will be used in nuclei recognition.
G A and G B is the gray value of pixel A and B respectively, r is the max radius of all nuclei in previous time points, d BS is the distance between B and S.
This work is supported by the Hong Kong Research Grants Council (Project HKBU5/CRF/11G).
- Boyle TJ, Bao Z, Murray JI, Araya CL, Waterston RH: AceTree: a tool for visual analysis of Caenorhabditis elegans embryogenesis. BMC Bioinformatics. 2006, 7: 275-10.1186/1471-2105-7-275.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones TR, Kand IH, Wheeler DB, Lindquist RA, Papallo A, Sabatini DM, Golland P, Carpenter AE: CellProfiler Analyst: data exploration and analysis software for complex image-based screens. BMC Bioinformatics. 2008, 9 (482):Google Scholar
- Chen X, Zhou X, Wong STC: Automated segmentation, classification, and tracking of cancer cell nuclei in time-lapse microscopy. IEEE T Bio-med Eng. 2006, 53 (4): 762-766. 10.1109/TBME.2006.870201.View ArticleGoogle Scholar
- Bao Z, Murray JI, Boyle T, Ooi SL, Sandel MJ, Waterston RH: Automated cell lineage tracing in Caenorhabditis elegans. Proc Natl Acad Sci U S A. 2006, 103: 2707-2712. 10.1073/pnas.0511111103.PubMed CentralView ArticlePubMedGoogle Scholar
- Murray JI, Bao Z, Boyle TJ, Boeck ME, Mericle BL, Nicholas TJ, Zhao Z, Sandel MJ, Waterston RH: Automated analysis of embryonic gene expression with cellular resolution in C. elegans. Nat Methods. 2008, 5 (8): 703-709. 10.1038/nmeth.1228.PubMed CentralView ArticlePubMedGoogle Scholar
- Murray JI, Bao Z, Boyle TJ, Waterston RH: The lineaging of fluorescently-labeled Caenorhabditis elegans embryos with StarryNite and AceTree. Nat Protoc. 2006, 1: 1468-1476. 10.1038/nprot.2006.222.View ArticlePubMedGoogle Scholar
- Rembold M, Loosli F, Adams RJ, Wittbrodt J: Individual cell migration serves as the driving force for optic vesicle evagination. Science. 2006, 313: 1130-1134. 10.1126/science.1127144.View ArticlePubMedGoogle Scholar
- Bischoff M, Parfitt DE, Zernicka-Goetz M: Formation of the embryonic-abembryonic axis of the mouse blastocyst: relationships between orientation of early cleavage divisions and pattern of symmetric/asymmetric divisions. Development. 2008, 135: 953-962. 10.1242/dev.014316.PubMed CentralView ArticlePubMedGoogle Scholar
- Santella A, Du Z, Nowotschin S, Hadjantonakis A-K, Bao Z: A hybrid blob-slice model for accurate and efficient detection of fluorescence labeled nuclei in 3D. BMC Bioinformatics. 2010, 11: 580-10.1186/1471-2105-11-580.PubMed CentralView ArticlePubMedGoogle Scholar
- Aydin Z, Murray JI, Waterston1 RH, Noble WS: Using machine learning to speed up manual image annotation: application to a 3D imaging protocol for measuring single cell gene expression in the developing C. elegans embryo. BMC Bioinformatics. 2010, 11: 84-10.1186/1471-2105-11-84.PubMed CentralView ArticlePubMedGoogle Scholar
- Richards JL, Zacharias AL, Walton T, Burdick JT, Murray JI: A quantitative model of normal Caenorhabditis elegans embryogenesis and its disruption after stress. Dev Biol. 2013, 374 (1): 12-23. 10.1016/j.ydbio.2012.11.034.PubMed CentralView ArticlePubMedGoogle Scholar
- Harder N, Mora-Bermudez F, Godinez WJ, Wunsche A, Eils R, Ellenberg J, Rohr K: Automatic analysis of dividing cells in live cell movies to detect mitotic delays and correlate phenotypes in time. Genome Res. 2009, 19: 2113-2124. 10.1101/gr.092494.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Long F, Peng H, Liu X, Kim SK, Myers E: A 3D digital atlas of C. elegans and its application to single-cell analyses. Nat Methods. 2009, 6: 667-672. 10.1038/nmeth.1366.PubMed CentralView ArticlePubMedGoogle Scholar
- Li G, Liu T, Tarokh A, Nie J, Guo L, Mara A, Holley S, Wong ST: 3D cell nuclei segmentation based on gradient flow tracking. BMC Cell Biol. 2007, 8: 40-10.1186/1471-2121-8-40.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen Y, Ladi E, Herzmark P, Robey E, Roysam B: Automated 5-D analysis of cell migration and interaction in the thymic cortex from time-lapse sequences of 3-D multi-channel multi-photon images. J Immunol Methods. 2009, 340: 65-80. 10.1016/j.jim.2008.09.024.PubMed CentralView ArticlePubMedGoogle Scholar
- Li G, et al: 3D cell nuclei segmentation based on gradient flow tracking. BMC Cell Biol. 2007, 8: 40-10.1186/1471-2121-8-40.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin G, Adiga U, Olson K, Guzowski JF, Barnes CA, Roysam B: A hybrid 3D watershed algorithm incorporating gradient cues and object models for automatic segmentation of nuclei in confocal image stacks. Cytometry A. 2003, 56A: 23-36. 10.1002/cyto.a.10079.View ArticleGoogle Scholar
- Yan P, Zhou X, Shah M, Wong ST: Automatic segmentation of high-throughput RNAi fluorescent cellular images. IEEE Trans Inf Technol Biomed. 2008, 12: 109-117.PubMed CentralView ArticlePubMedGoogle Scholar
- Mosaliganti K, Cooper L, Sharp R, Machiraju R, Leone G, Huang K, Saltz J: Reconstruction of cellular biological structures from optical microscopy data. IEEE Trans Vis Comput Graph. 2008, 14: 863-876.View ArticlePubMedGoogle Scholar
- Murray JI, Boyle TJ, Preston E, Vafeados D, Mericle B, Weisdepp P, Zhao Z, Bao Z, Boeck ME, Waterston R: Multidimensional regulation of gene expression in the C. elegans embryo. Genome Res. 2012, 22 (7): 1282-1294. 10.1101/gr.131920.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Osher S, Sethian JA: Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations. J Comput Phys. 1988, 79: 12-49. 10.1016/0021-9991(88)90002-2.View ArticleGoogle Scholar
- Pratt WK: Digital Image Processing. 2007, Los Altos, California: John Wiley & Sons, Inc., 4View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.