implemented cell automaton spreading only in direction of colony boundary, as determined by bwdist using Euclidean distance. Issues were, for calculation of 2d indices from bwdist index output, subtracting minus one from indices before modding or dividing, then adding it back, and apparently the bwdist index output is strictly integer, and must be cast to double, multiplying by 1.0 does not work. Resulting colonies are diamond-shaped due to use of Euclidean distance, must find a way around that. Also discovered that matrixToPlotScale must be initialized as uint16 matrix, had done this previously by calling uint16 on matrixToPlot, but had to discard that when using channel-specific scaling. If matrixToPlotScale is normal matrix, simply putting uint16 values in each channel will give solid colors whenever a nonzero value. Finally, to get realistic simulation of 20mm by 20mm Cong image, made dx=.1, sideLength=200, tripled DL and DR (need to check whether this might just be physically realistic), and multiplied biomass by 100 from OD scaled value (but really need to check whether the original OD scaling is valid). Also modified XThresh to assume 1 micrometer square bacterium, from Dave's advice, also need to check whether this fits given OD scaling