1 Introduction
Visual analysis is a vital tool in life sciences. Scientists can automate experiments and capture large result sets of digital images Scientists10 . Automatic analysis of the enormous sets of biomedical images is a must. Algorithms and techniques of image analysis are highly parameterized and significant human input is intertwined to optimize parameter settings Ruddle11
. It is common to utilize machinelearning algorithms in minimizing user overhead. However, users are asked to determine objects of interest in some input images, then the algorithms create classifiers to recognize similar objects in other images. This approach highly relies on the experience of domain experts and it requires a great amount of object labeling to be done by hand.
Also, lacking a standard metavisualizationmeta11 approach is a challenge. Highlevel description of a visualization algorithm can not be gained by debugging, forcing developers of the algorithm to manually create highlevel illustrations to understand the workflow of the algorithm.
This paper proposes an alternative iterative approach for optimizing input parameters, saving time by minimizing the user involvement, and allowing for understanding the workflow of algorithms and discovering new ones. The main focus is on developing an interactive visualization technique that enables users to analyze the relationships between sampled input parameters and corresponding output. This technique is implemented in a prototype called Veni Vidi Vici , following the famous Latin sentence claimed to be said by Julius Caesar venividivici . It translates as ”I came, I saw, I conquered.” This strategy is inspired by the mathematical formulas of numbering computable functions Cutland80 and is developed atop ImageJ ImageJ , a scientific biomedical image analysis program.
Section 2 summarizes related work, before surveying ImageJ in Section 3. Proposed work is presented in Section 4. This includes explaining Veni Vidi Vici strategy (Section 4.1), presenting the underpinning mathematical foundations (Section 4.2), and implementing the framework as a custom plugin for ImageJ (Section 4.3). Evaluation of the proposed approach with case studies are shown in Section 5. Finally, presenting conclusion and outlining future work in Section 6.
2 Related Work
Many visualization techniques are proposed to facilitate investigating the relationships between parameters and outcomes. Linked charts are among the first attempts to use interactive visualization to parameter optimization Tweedie95 . The mapping , where are the parameters and are the outcomes, can be represented as n + m interactive histograms. This approach is then extended to prosection matrices Tweedie98 . An alternative approach is proposed for visualizing the parameter search process as a directed graph Ma99 .
Exploration graphs for history management are used in VisTrails visualization management system VisTrails05 . The exploration is still linear and time consuming. This can be partially improved by allowing users to write code snippets to perform iterations over parameter values VisTrails06 .
Many techniques are proposed for imagebased analysis and visualization as a structured gallery JankunKelly0 . An innovative technique is proposed based on preset based interaction with high dimensional parameter spaces Preset03 . A preset is a reference point in parameter space used to compute the outcome by calculating the linear combination of the inverse distances from the preset to the outcome.
Many successful analysis and visualization systems are available such as IRIS ExplorerIRIS95 . CellProfiler CellProfiler06 is an image analysis software for identifying and quantifying cell phenotypes. VolumeShop VolumeShop05 is an interactive system for direct volume illustration.
Vismon Maryam11 is a data analysis visual tool for fisheries to enable decision makers to quickly narrow down all possible management options to an agreeable small set. A risk assessment framework for chum salmon is developed to evaluate different management policies CPZ09 . A detailed tradeoff among the few chosen management options can be performed Maryam12 .
Most of previous systems are domainspecific. However, some are for a wide range of applicability. Tuner Tuner11 is an application for finding optimal parameter settings for complex algorithms in highdimensional parameter space. It supports sampling, resampling, and visually seeing the stability of certain parameter settings with respect to the objective measures.
For more wisdom on new trends in visualization pipeline, algorithmic stability, and the exploration of parameters, readers are encouraged to check IEEE VisWeek 2012 tutorial on ”Uncertainty and Parameter Space Analysis in Visualization” VisWeek12 .
3 ImageJ
ImageJ ImageJ is a free image processing program developed at the National Institutes of Health. Originally, it is designed to be used by biomedical researchers working with microscope images and videos.
Being based on Java which supports algorithmic digital image processingJava7 , ImageJ is designed with an open architecture that provides extensibility via Java plugins and recordable macrosImageJ2 . These plugins provide custom solution for many analysis and visualization problems such as:

threedimensional livecell imagingImageJ3 ,

to radiological image processingImageJ4 ,

multiple imaging system data comparisonsImageJ5 ,

automated hematology systemsImageJ6 , and

Visualizing multidimensional biological image dataImageJ7 .
The following four subsections show a comprehensive list of plugins.
3.1 Analysis

Autocorrelation,

MRI t2 calculations,

Line Analyzer,

Image Correlator ,

Particle Remover,

Circularity ,

Modulation Transfer Function,

Specify ROI ,

Specify Line Selection,

16bit Histogram ,

Draw line or point grids ,

Moment Calculator ,

Batch Statistics ,

Cell Counter ,

Oval Profile Plot ,

Color Comparison ,

Radial Profile Plot ,

Microscope Scale ,

MRI Analysis Calculator ,

Sync Measure 3D ,

Hough Circles ,

Convex Hull, Circularity, Roundness ,

Fractal Dimension and Lacunarity ,

Measure And Label ,

Colocalization ,

Granulometry ,

Texture Analysis ,

Named Measurements ,

Cell Outliner ,

Grid Cycloid Arc ,

RGB Profiler ,

Colocalization Finder ,

Spectrum Extractor ,

Contact Angle ,

RG2B Colocalization ,

Color Profiler ,

Hull and Circle ,

MR Urography ,

Template Matching ,

Extract IMT from ultrasound images ,

ITCN (Imagebased Tool for Counting Nuclei) ,

Multi Cell Outliner ,

FRETcalc  FRET by acceptor photobleaching ,

JACoP (Just Another Colocalization Plugin),

FRET and Colocalization Analyzer ,

CASA (Computer Assisted Sperm Analyzer) ,

Radial Profile Plot Extended ,

Concentric Circles (nondestructive overlay),

Azimuthal Average ,

Slanted Edge Modulation Transfer Function,

Calculate 3D Noise ,

FWHM (analyze photon detector pinhole images),

SSIM index (calculate structural similarity),

Image Moments (image moments of nth rank) ,

MS SSIM index (multiscale structural similarity),

Colony Counter (count colonies in agar plates),

Levan (chromosome morphology) ,

EXTRAX (electron diffraction intensity extraction),

Fractal Surface Measurement ,

Foci Picker3D (finds local maxima in 2D and 3D images),

Diameter (measures the diameter of a blood vessels),

Graph Demo (creates particle adjacency lists) ,

Asymmetry Analysis (HRTEM image conditions) ,

2D NMR Analysis (integrates peaks in 2D NMR spectra) ,

MetaData and Intracellular Calcium Line Scan Analysis,

GHT (General Hough Transformation object recognition),

IntraCell (nanoparticle colocalization within cells),

Lemos Asymmetry Analysis (dental panoramic radiographs),

Merz Grid Macro (semicircular lines and points in overlay),

Stress Granule Counter (counts SGs in eucaryotic cells),

Vamp 2D and 3D (isolate puncta in 2D and 3D images) ,

Sampling Window (unbiased sampling window) ,

Map Bone Microstructure (histomorphometry parameters),

Results and Text ,

Comment Writer
3.2 Filters

Real Convolver

Fast Fourier Transform (FFT)

LoG Filtering

Background Subtraction and Normalization

Contrast Enhancer

Background Correction

Byte Swapper

Discrete Cosine Transform (DCT)

FFT Filter

FFTJ and DeconvolutionJ

Unpack 12bit Images

Deinterlace

2D Gaussian Filter

DualEnergy Algorithm

Anisotropic Diffusion (edgepreserving noise reduction)

Grayscale Morphology

2D Hybrid Median Filter

3D Hybrid Median Filter

Spectral Unmixing

Haar Wavelet Filter and Adaptive Median Filter

’A trous’ Wavelet Filter

Kuwahara Filter

Granulometric Filtering

WindowedSinc Filter (low pass time series filter)

Anisotropic Diffusion 2D (edgepreserving noise reduction)

Auto Gamma for gamma correction

Linearize Gel Data

Radon Transform (back projection, sinogram)

Correct X Shift of Confocal Images

Multi Otsu Threshold

Spectral Unmixing of Bioluminescence Signals

Lipschitz Filter

Float Morphology (erode, dilate, open, close)

X Shifter (correct pixel mismatch of confocal images)

Sigma Filter (edgepreserving noise reduction)

Rolling Ball Background Subtraction

Mean Shift Filter (edgepreserving smoothing)

Accurate Gaussian Blur

Add Poisson Noise

CLAHE (Contrast Limited Adaptive Histogram Equalization)

Floyd Steinberg Dithering

Polar Transformer (corrects radial and angular distortions)

Gaussian Blur 3D

Image Rotator (rotates image around ROI center of mass)

Mexican Hat (2D Laplacian of Gaussian)
3.3 Segmentation

Mixture Modeling Thresholding ,

Otsu Thresholding ,

Watershed Segmentation Maximum Entropy Thresholding ,

MultiThresholder (IsoData, MaxEntropy, Otsu, etc),

Multi Otsu Threshold ,

SIOX (Simple Interactive Object Extraction),

RATS (Robust Automatic Threshold Selection) ,

Densitometry ,

Blob Labeler (labels connected blobs of pixels)
3.4 Collections

UCSD Confocal Microscopy Plugins ,

MBF ImageJ for Microscopy Collection
4 Proposed Framework
4.1 Veni Vidi Vici
Let us discuss the general idea behind the proposed new iterative approach for optimizing input parameters and allowing for understanding the workflow of algorithms and discovering new ones. The main focus is on developing an interactive visualization technique that enables users to analyze the relationships between sampled input parameters and corresponding output. This technique is implemented in a prototype called Veni Vidi Vici , following the famous Latin sentence claimed to be said by Julius Caesar venividivici . It translates as ”I came, I saw, I conquered.”
It falls into three main parts:

Veni: corresponds to analysis and user interaction,

Vidi: corresponds to visualization, and

Vici: corresponds to segmentation or similar tasks.
Figure 1 shows the schematic diagram of the proposed strategy.
Note that rectangles represent files and rounded rectangles represent processes. Veni takes the original image and the outcome of previous runs as input. The dashed line is not presented in the first run. User interaction is minimized to Veni. Vidi is sandwiched between the other two processes. Vici corresponds to a desired task to perform on the image such as segmentation or counting objects. The final outcome is then fedback to Veni. These three processes are optional ,i.e. any one can be skipped, with at least one applied.
4.2 Mathematical Foundations
Recall that the over all process involves a set of n parameters: . Assume that these n parameters are divided among the three processes of Veni, Vidi, and Vici. Assume that and are the parameter shares of Veni, Vidi, and Vici respectively. So, . Veni, Vidi, and Vici are descried mathematically in table 1.
Operation  Summary 

Veni(g, )  Apply Veni on images g, setting its a parameters to x’s values 
Vidi(g, )  Apply Vidi on images g, setting its b parameters to y’s values 
Vici(g, )  Apply Vici on images g, setting its c parameters to z’s values 
The proposed strategy for exploring parameter space is inspired by the mathematical formulas of numbering computable functions Cutland80 . Let us define some functions that allow us to map any set of Veni, Vidi, and Vici to
a unique code, or Gödel number godel .
Mapping ordered pairs to N
A function such that maps . That is, maps the ordered pair (x,y) to a single number .
To see how this mapping is done, note that is a number such that
has x factors of 2 and a remaining odd number 2y + 1. As an example, let
= 103935. You can factories + 1 into 9 factors of 2 and you are left over with . Hence, yielding x = 9 and y = 101. Since every number has a different factorization into powers of two, is a bijective function mapping. The inverse of is defined such that (z) returns the exponent of = 2 in the prime factorization of z, and (z) returns the exponent of = 3.Mapping Ordered Triples to N
A function maps an ordered set of three natural numbers to a single number. To define this mapping, we can use the function , and write which maps the first parameters m and n to one number using . The resulting value and the third number q is then mapped to one number again using . Hence, the final result is one natural number. The purpose in subtracting 1 from the parameters m,n and q is to include zero in the mapping.
Mapping a Finite Sequence to N
The concept of converting bases can be used to devise a bijection from a finite sequence of numbers to one natural number.
The binary representation is simple yet powerful. A decimal number can be converted into a binary representation by repeatedly subtracting the next largest possible power of two from the previous remainder, until no other powers exist.
Consider the conversion of 44 decimal to binary: . By putting 1 for every existing power of 2, the binary representation of 44 is 101100. Let us define a function mapping from a finite sequence of numbers to one natural number. Let the function be defined as As an example, let us perform the calculations for
4.3 Implementation
Let us propose the following algorithm.
The algorithm starts with loading default settings from setting file containing the values on parameters. The parameters are divided among Veni Vidi Vici operations. The algorithm also asks for the exploration range. The smaller it is, the fewer possibilities the algorithm examines. Then the algorithm load an image or a set of images, and perform successive Veni Vidi Vici operations. Each operation emits a set of files, and has a code calculated using the underpinning mathematical foundations presented in Section 4.2. The total code of Veni Vidi Vici operations is then calculated in the same way.
The total code allows for enumerating outcomes, one by one along the exploration range. Each outcome has a corresponding decoded settings. In 17^{th} line, user interaction is required to determine which settings suits his or her needs.
The framework can be implemented in many ways such as a recordable macro or a plugin as shown in figure 3.
An editor pop ups to modify the run method. After compiling the plugin, the created plugin is shown in figure 4. Note that the name of the plugin must contains underscores to be automatically loaded when ImageJ starts.
5 Evaluation
To load an image, just dragdrop it on ImageJ. By clicking the menu item for the plugin, it starts working.
Many trails were done to test the proposed framework. This was done on images provided on ImageJ site, shown in figure 5.
Many Veni operations can be done on the loaded image such as analysis plugins listed in section 3. Figures 6,7,8 shows the application of edge detection, plot profile, and surface plot.
Many Vici operations can be done such as segmentation plugins listed in section 3. Figure 9 shows the application of Otsu thresholding technique otsu79
. This technique divides the histogram in two classes and then the interclass variance is minimized. This plugin outputs a thresholded image with the selected threshold.
6 Conclusion and Future Work
This paper proposes an alternative iterative approach to optimize input parameters and save time by minimizing the user involvement. This strategy is developed as custom plugin in ImageJ.
The main focus was on developing an interactive visualization technique that enables users to analyze the relationships between sampled input parameters and corresponding output. This technique is implemented in a prototype called Veni Vidi Vici. It provides users with a visual overview of parameters and their sampled values. To find optimal parameter settings, users can select best possible configuration.
Not only the proposed framework facilitates better algorithm understanding, but also it can be used to explore the parameter space.
Future work have many dimension. One dimension is to consider the application of the proposed framework on temporal data such as feature point tracking and trajectory analysis for video imaging mosaic .
Another dimension is extending the theoretical foundation of the proposed framework to be able to answer some questions regarding the parameters of Veni, Vidi, and Vici. Questions such as : How does a slight parameter change modify the result? How stable is a parameter? In which range is a parameter stable?
Another dimension is to enhance the accuracy of the framework by incorporating transfer learning
Sinno , to be able to generalize what learnt from one case to another. In many machine learning, there is an assumption that the training and future data must be in the same feature space and have the same distribution. However, in many realworld applications, this assumption may not hold. Knowledge transfer, if done successfully, would greatly improve the performance of learning by avoiding much expensive datalabeling efforts.References
 (1) T. Walter, D. W. Shattuck, R. Baldock, M. E. Bastin, A. E. Carpenter, S. Duce, J. Ellenberg, A. Fraser, N. Hamilton, S. Pieper, M. A. Ragan, J. E. Schneider, P. Tomancak, and J.K. H erich e. Visualization of image data from cells to organisms. NatureMethods Supplement, 7(3):S26 S41, 2010
 (2) A.J. Pretorius, M.A.P. Bray, A.E. Carpenter and R.A. Ruddle, Visualization of parameter space for image analysis, IEEE Transactions on Visualization and Computer Graphics, vol. 17, no. 12, pp. 24022411, 2011.
 (3) Peter Sikachev and Artem Amirkhanov and Robert S. Laramee and Gabriel Mistelbauer. Interactive Algorithm Exploration using Meta Visualization, Technical report, Institute of Computer Graphics and Algorithms, Vienna University of Technology, 2011.
 (4) http://www.unrv.com/fallrepublic/venividivici.php.
 (5) N. J. Cutland. Computability: An Introduction to Recursive Function Theory, Cambridge University Press, 1980, pp. 7284.
 (6) Collins TJ. ImageJ for microscopy. BioTechniques, 2007, volume 43, Suppl 1, pp.25 30.
 (7) L. Tweedie, B. Spence, H. Dawkes, and H. Su. The influence explorer. In Proceedings of the Conference on Human Factors in Computing Systems, pages 129 130, 1995.
 (8) L. Tweedie and R. Spence. The prosection matrix: a tool to support the interactive exploration of statistical models and data. Computational Statistics, 13(1):65 76, 1998.
 (9) K.L. Ma. Image graphs a novel approach to visual data exploration. In Proceedings of the IEEE Conference on Visualization, pages 81 88, 1999.
 (10) L. Bavoil, S. P. Callahan, P. J. Crossno, J. Freire, and H. T. Vo. VisTrails: enabling interactive multipleview visualizations. In Proceedings of the IEEE Conference on Visualization, pages 135 142, 2005.
 (11) S. P. Callahan, J. Freire, E. Santos, C. E. Scheidegger, C. T. Silva, and H. T. Vo. Managing the evolution of dataflows with VisTrails. In Proceedings of the International Conference on Data Engineering Workshops, pages 71 75, 2006.
 (12) T. J. JankunKelly and K.L.Ma. A spreadsheet interface for visualization exploration. In Proceedings of the IEEE Conference on Visualization, pages 69 76, 2000.

(13)
J. J. van Wijk and C. W. A. M. van Overveld.
Preset based interaction with high dimensional parameter spaces
. Data Visualization: The State of the Art, pages 391 406, 2003.
 (14) D. Foulser. IRIS Explorer: a framework for investigation. SIGGRAPH Computer Graphics, 29:13 16, 1995.
 (15) A. E. Carpenter, T. R. Jones, M. R. Lampbrecht, C. Clarke, I. H. Kang, O. Friman, D. A. Guertin, J. H. Chang, R. A. Lindquist, J. Moffat, P. Golland, and D.M. Sabatini. CellProfiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biology, 7(R100), 2006.
 (16) Bruckner S., Gröller M. E. VolumeShop: An interactive system for direct volume illustration. In Proceedings of IEEE Visualization 2005 (October 2005), pp. 671 678.
 (17) Maryam Booshehrian, Torsten Möller, Randall M. Peterman, and Tamara Munzner. Vismon: Facilitating Risk Assessment and Decision Making In Fisheries Management, Technical report TR 201104, School of Computing Science, Simon Fraser University, Burnaby, BC, Canada, September 2011.
 (18) Maryam Booshehrian, Torsten Möller, Randall M. Peterman, and Tamara Munzner. Vismon: Facilitating Analysis of TradeOffs, Uncertainty, and Sensitivity In Fisheries Management Decision Making, Computer Graphics Forum (Proceedings of Eurographics Conference on Visualization 2012 (EuroVis 2012)), vol. 31, no. 3, pp. 12351244.
 (19) Collie J. S., Peterman R. M., Zuelke B., A Risk Assessment Framework for AYK Salmon. Tech. rep., Project completion report for the Alaska Sustainable Salmon Fund, Nov. 2009.
 (20) Thomas TorsneyWeir, Ahmed Saad, Torsten Moller, HansChristian Hege, Britta Weber, and JeanMarc Verbavatz. Tuner: Principled Parameter Finding for Image Segmentation Algorithms Using Visual Response Surface Exploration. IEEE Transactions on Visualization and Computer Graphics, volume 17, number 12, pp. 18921901, 2011.
 (21) http://www.cg.tuwien.ac.at/research/publications/2012/VisWeekTutorial2012Uncertainty/
 (22) Burger W, Burge M. Digital Image Processing: An Algorithmic Approach Using Java. Springer. ISBN 1846283795, 2007
 (23) Girish V, Vijayalakshmi A. Affordable image analysis using NIH Image/ImageJ. Indian J Cancer vol. 41, no. 1, pp. 47, 2004.
 (24) Eliceiri K, Rueden C. Tools for visualizing multidimensional images from living specimens. Photochem Photobiol vol. 81 , no. 5, pp. 1116 1122, 2005.
 (25) Barboriak D, Padua A, York G, Macfall J. Creation of DICOM Aware Applications Using ImageJ. J Digit Imaging , vol. 18, no. 2, pp. 9199, 2005.
 (26) Rajwa B, McNally H, Varadharajan P, Sturgis J, Robinson J. AFM/CLSM data visualization and comparison using an opensource toolkit. Microsc Res Tech vol. 64, no. 2, pp. 176 184, 2004.
 (27) Gering E, Atkinson C. A rapid method for counting nucleated erythrocytes on stained blood smears by digital image analysis. J Parasitol vol. 90, no. 4, pp. 879 881, 2004.
 (28) Rueden CT, Eliceiri KW. Visualization approaches for multidimensional biological image data. BioTechniques vol. 43, Suppl. 1, pp. 33 36, 2007.
 (29) Davis, Martin D. Computability and Unsolvability. Dover Publishing, 1983
 (30) N. Otsu, A threshold selection method from gray level histograms, IEEE Trans. Systems, Man and Cybernetics, vol. 9, pp. 6266 1979
 (31) I. F. Sbalzarini and P. Koumoutsakos. Feature Point Tracking and Trajectory Analysis for Video Imaging in Cell Biology, Journal of Structural Biology vol. 151, no. 2, pp. 182195, 2005
 (32) Sinno Jialin Pan and Qiang Yang, A Survey on Transfer Learning, Knowledge and Data Engineering, IEEE Transactions on, 2010, volume 22, number 10, pages 13451359.