| Title: | Furrow Segmentation |
|---|---|
| Description: | Image feature data and analysis codes for the Guglielmi, Barry et al. paper describing the application of an optogenetics tools to disrupt Drosophila embryo furrowing. |
| Authors: | Joseph Barry |
| Maintainer: | Joseph Barry <[email protected]> |
| License: | Artistic-2.0 |
| Version: | 1.40.0 |
| Built: | 2026-10-10 11:25:13 UTC |
| Source: | https://git.bioconductor.org/packages/furrowSeg |
Calculates dimensions of box for at a given DV position. Ensures that box does not exceed dimensions of image.
constructBox(dvPos, Lx=100, Ly=50, w=512, mid=NA)constructBox(dvPos, Lx=100, Ly=50, w=512, mid=NA)
dvPos |
Pixel location along DV of box center. |
Lx |
Half of box width in pixels. |
Ly |
Half of box height in pixels. |
w |
Image width in pixels. |
mid |
Location of midpoint along AP in pixels. If not specified defaults to half of the image width. |
A vector with locations of box corners. Nonclemature is 'xleft', 'ybottom', 'xright' and 'ytop'.
Joseph Barry, 2014
if (interactive()) vignette(topic="genPaperFigures", package="furrowSeg")if (interactive()) vignette(topic="genPaperFigures", package="furrowSeg")
An example movie on which furrowSeg segmention can be performed.
exampleFurrowMovieexampleFurrowMovie
A 4D array.
data(exampleFurrowMovie, package="furrowSeg") dim(exampleFurrowMovie)data(exampleFurrowMovie, package="furrowSeg") dim(exampleFurrowMovie)
Identifies furrowing line by identifying DV position of minimum area.
identifyFurrowPosition(x, nbinsExclude=3, h=100, plot=FALSE, myCex=1.4, w=512, px=0.293)identifyFurrowPosition(x, nbinsExclude=3, h=100, plot=FALSE, myCex=1.4, w=512, px=0.293)
x |
Feature table. |
nbinsExclude |
Number of pixel columns to exclude at the DV edges of the image. |
h |
Smoothing bandwidth, passed to locfit. |
plot |
Logical specifying whether or not to plot data and fit. |
myCex |
Size of axis labels. |
w |
Width of image in number of pixels. |
px |
Pixel dimensions in microns (assumed isotropic). |
The pixel index along DV indicating the furrowing position.
Joseph Barry, 2014
if (interactive()) vignette(topic="genPaperFigures", package="furrowSeg")if (interactive()) vignette(topic="genPaperFigures", package="furrowSeg")
Identifies time point where the cell areas attain a minimum.
identifyTimeMinArea(x, h=2, px=0.293, plot=FALSE, myCex=1.4)identifyTimeMinArea(x, h=2, px=0.293, plot=FALSE, myCex=1.4)
x |
Feature table. |
h |
Smoothing bandwidth, passed to locfit. |
px |
Pixel dimensions in microns (assumed isotropic). |
plot |
Logical specifying whether or not to plot data and fit. |
myCex |
Size of axis labels. |
Returns the time at which the tissue invaginates ('tstar') and the index of the corresponding time point ('tindex').
Joseph Barry, 2014
if (interactive()) vignette(topic="genPaperFigures", package="furrowSeg")if (interactive()) vignette(topic="genPaperFigures", package="furrowSeg")
Checks if a number is odd or adds one to make it odd. Useful for constructing filters.
isOdd(x) makeOdd(x)isOdd(x) makeOdd(x)
x |
An integer. |
A logical indicating if number is odd or an odd integer.
Joseph Barry, 2014
isOdd(seq(1:10))isOdd(seq(1:10))
Subsets feature table to include only cells whose center are in the interior of the specified box dimensions.
isolateBoxCells(x, box)isolateBoxCells(x, box)
x |
Feature table containing centroid positions as 'x.0.m.cx' and 'x.0.m.cy'. |
box |
Coordinates of box corners, specified as 'xleft', 'ybottom', 'xright' and 'ytop'. |
A subsetted 'x' containing box cells.
Joseph Barry, 2014
if (interactive()) vignette(topic="genPaperFigures", package="furrowSeg")if (interactive()) vignette(topic="genPaperFigures", package="furrowSeg")
Table containing all cell feature data for optogenetically perturbed samples and controls. Contains the following columns:
sampleUnique sample identifier referring to the .rda object from which the image analysis was loaded.
tInteger index of time point.
zInteger index of z-stack.
x.0.m.cxx position (along anterior-posterior axis) of cell center in number of pixel lengths.
x.0.m.cyy position (along dorsal-ventral axis) of cell center in number of pixel lengths.
x.0.m.majoraxisLength of major axis of the cell.
x.0.m.thetaAngle between the major axis of the cell and the anterior-posterior axis of the embryo.
x.0.s.areaArea of the cell in number of pixels.
x.0.s.perimeterPerimeter length of cell in number of pixel lengths.
x.0.s.radius.meanMean radius of cell in number of pixel lengths.
x.0.s.radius.maxMaximum radius of cell in number of pixel lengths.
e.xFirst component of anisotropy vector. Referred to as AP anisotropy in the paper.
e.ySecond component of anisotropy vector. Referred to as DV anisotropy in the paper.
dtTime between frames in seconds
pxSide length of a (square) pixel in microns. Note that the z-stack spacing is longer.
conditionFactor identifying which experimental condition cell is associated with.
optoopto
A data table.
data(opto, package="furrowSeg") head(opto)data(opto, package="furrowSeg") head(opto)
Plots mean and standard deviation of area and elongation features over time.
plotFeatureEvolution(x, dt=32.6/60, tMax, myTitle="", cex=1.4, cex.axis=1, px=0.293, mar=c(5.1, 5.1, 4.1, 4.1), legend=TRUE, line=2.5)plotFeatureEvolution(x, dt=32.6/60, tMax, myTitle="", cex=1.4, cex.axis=1, px=0.293, mar=c(5.1, 5.1, 4.1, 4.1), legend=TRUE, line=2.5)
x |
A feature table, as supplied by |
dt |
Timestep in minutes (numeric). |
tMax |
Latest time point to plot in minutes (numeric). |
myTitle |
Plot title (string). |
cex |
Label size. |
cex.axis |
See help for |
px |
Pixel width in microns. |
mar |
See help for |
legend |
A logical. Should figure legend be displayed or not? |
line |
Determines placement of right-hand axis label. See help for
|
Nothing is returned from this function.
Joseph Barry, 2014
if (interactive()) vignette(topic="genPaperFigures", package="furrowSeg")if (interactive()) vignette(topic="genPaperFigures", package="furrowSeg")
Converts area in pixels to microns squared and vice versa.
px2area(x, px) area2px(x, px)px2area(x, px) area2px(x, px)
x |
A vector of numbers. |
px |
Side-length of a pixel in microns. |
A vector of areas in new units.
Joseph Barry, 2014
# pixels side-length half a micron, square of 10x10 pixels px2area(x=10*10, px=0.5)# pixels side-length half a micron, square of 10x10 pixels px2area(x=10*10, px=0.5)
Converts length in pixels to microns and vice versa.
px2microns(x, px) microns2px(x, px)px2microns(x, px) microns2px(x, px)
x |
A vector of numbers. |
px |
Side-length of a pixel in microns. |
A vector of lengths in new units.
Joseph Barry, 2014
# map a contiguous block of 8 pixels to position in microns (here pixel side-length is half a micron) px2microns(x=seq(1:8), px=0.5)# map a contiguous block of 8 pixels to position in microns (here pixel side-length is half a micron) px2microns(x=seq(1:8), px=0.5)
Contains names of the images used in study, and assigns them to their respective experimental groupings. The time interval between frames is listed in seconds and the (isotropic) pixel dimensions in microns.
sampleTablesampleTable
A data table.
data(sampleTable, package="furrowSeg") head(sampleTable)data(sampleTable, package="furrowSeg") head(sampleTable)
Performs segmentation on furrow images using smoothing, adaptive thresholding and watershed algorithms.
segmentFurrowAllStacks(x, L=17, filterSize=3, threshOffset=0.001, closingSize=3, minObjectSize=2^5, maxObjectSize=2^10)segmentFurrowAllStacks(x, L=17, filterSize=3, threshOffset=0.001, closingSize=3, minObjectSize=2^5, maxObjectSize=2^10)
x |
A 4-dimensional image with dimensions x, y, z, t |
L |
The characterisitc diameter of a cell in pixels. |
filterSize |
The size of the filter for gaussian smoothing. |
threshOffset |
The offset value for the adaptive thresholding algorithm that is used to segment cytoplasmic fluorescence signal. |
closingSize |
The size of the brush that is used to perform a closing operation that smooths the cytoplasmic mask after the adaptive thresholding. |
minObjectSize |
Determines the threshold below which objects in the cytoplasmic mask are removed. |
maxObjectSize |
Determines the threshold above which objects in the cytoplasmic mask are removed. |
A list with items.
x |
A smoothed version of the original image array |
mask |
Cell masks |
hs |
An image showing highlighted segmentation of the cell masks |
Joseph Barry, 2014
if (interactive()) vignette(topic="exampleFurrowSegmentation", package="furrowSeg")if (interactive()) vignette(topic="exampleFurrowSegmentation", package="furrowSeg")