We present a vector field method for obtaining the spatial organization of three-dimensional patterns of gene expression based on gradients and lines of force obtained by numerical integration. The convergence of these lines of force in local maxima are centers of gene expression, providing a natural and powerful framework to characterize the organization and dynamics of biological structures. We apply this methodology to analyze the expression pattern of the enhanced green fluorescent protein (EGFP) driven by the promoter of light chain myosin II during zebrafish heart formation.