While frequency-based methods for the characterization of fractals are popular and effective in many applications, they have limitations when applied to irregularly shaped images, such as nuclear images. The irregularity renders texture characterization by frequency domain methods, based upon Fourier transform, problematic. To address this situation, this paper presents an algorithm based upon the signal analysis in the spatial domain. An autocovariance function can be estimated regardless of the shape and size of regions where the image is defined. As in the continuous fractional Brownian motion (FBM) that results from inputting white noise into a specific fractional integral system, a discrete FBM can be related to white noise by a specific fractional summation system (FSS) that is linear, causal and shift-invariant. Although the method of direct sampling is not valid for converting a continuous fractional integral to a discrete fractional summation, discrete fractional summations similar to the sampled system functions can be obtained through an iterative process. While the continuous system function of a fractional integral is linear in the frequency domain when plotted in log-log scales, unfortunately, it is not true for the comparable discrete system function. The discrete system function is actually approximately linear in the log-log scales over a very limited range. The slope of the straight line that approximates the function curve in the mean-square-error (MSE) sense in a specific time range provides a description of the autocovariance function that reveals the statistical relations among the local textures. Applications to characterization of ovary nuclear images in groups of normal, atypical and cancer cases are studied and presented.