Input Data Options
	special purpose data requiring pre-processing (0)
	tomographic data (1)
	segmented data (2)
	burn data (3)
	medial axis data (4)
	throat data (5)
	fluid data (6)
Enter choice: 2

Data Processing Options
	plot image (1)
	resize data (2)
	disconnected volume distribution (3)
	clean up segmented files (4)
	burn and compute LKC medial axis/surface (5)
	moment of inertia of disconnected components (6)
	convert data to/from ascii format (7)
	2-point covariance function (8)
	tomographic/segmented data covariance comparison (9)
	compress data in z direction (10)
	count voxels in spherical/cylindrical shells (11)
	set fiducial polygon exterior to grain phase (12)
Enter choice: 5
LKC computation of
	0) medial axis
	1) medial surface
Specify request (0,1): 0

Methods for handling the medial surface/axis computation at the fiducial polygon
or volume edges. We recognize two possibilities for the region exterior to the
edge, either i) it is of material type opposite to that for which we are
constructing the medial axis or ii) the composition of the exterior is
unknowable.  For case ii) edges we implement `absorbing' boundary conditions.
Current options for treating edges are
	0) all edges are case i)  [original LKC algorithm]
	1) all edges are case ii)
	2) top/bottom slice are case ii), all other edges are case i)
Specify choice (0,1,2): 1

Methods for examining/deleting voxels during erosion
	1) exmine/delete during each of the six directional sweeps
						[original LKC algorithm].
	2) delete only after all six directional examination sweeps done.
Specify choice (1,2): 1

Specify connectivity (6,26(dflt)): 26


    Files assumed to be labelled  basename.ext,
    basename  is limited to 255 characters maximum,
    ext       is a numerical designator lying between 000 and 999
    It is assumed that the files have consecutively numbered extensions
        ie.  000 -> 056,  or   021 -> 049

Enter basename for segmented files: ../c_seg/oil
Are files compressed? [y,n]: y


Enter first and last slice of data to use: 1 100 

The material/void boundary can be lightly smoothed.
Available options are
	0) no conversion
or convert those boundary voxels having
	1) exactly one neighbor of the same type
	2) less than a majority of neighbors of the same type
Enter choice (0(dflt),1,2): 0


Isolated clusters of grain and/or pore voxels up to a specifed size
can be assumed to be misidentified and converted to the opposite material type.
Convert isolated grain clusters? (y,n): n
Convert isolated pore  clusters? (y,n): n


Input data can be inverted to compute burn/medial axis of
grain space rather than void space.
Invert data (y,n(dfnt)): n

Exclude the first or last slices to compensate
for lack of information beyond the end slices?
[(n)one, (f)irst, (l)ast, (b)oth]: n
Enter basename for burned files: ../burn/oil
Are files to be compressed? [y,n]: y


Fiducial polygon generation methods
	NONE      (N)
	MANUAL    (M)
	AUTOMATIC (A)
Enter method: N

Enter basename for LKC ma files: ../ma/oil
Are files to be compressed? [y,n]: y