Abstract Rorschach inkblots have had a striking impact on the worlds of art and science because of the remarkable variety of associations with recognizable and namable objects they induce. Originally adopted as a projective psychological tool to probe mental health, psychologists and artists have more recently interpreted the variety of induced images simply as a signature of the observers’ creativity. Here we analyze the relationship between the spatial scaling parameters of the inkblot patterns and the number of induced associations, and suggest that the perceived images are induced by the fractal characteristics of the blot edges. We discuss how this relationship explains the frequent observation of images in natural scenery.

Citation: Taylor RP, Martin TP, Montgomery RD, Smith JH, Micolich AP, Boydston C, et al. (2017) Seeing shapes in seemingly random spatial patterns: Fractal analysis of Rorschach inkblots. PLoS ONE 12(2): e0171289. https://doi.org/10.1371/journal.pone.0171289 Editor: Stefan Glasauer, Ludwig-Maximilians-Universitat Munchen, GERMANY Received: July 29, 2016; Accepted: January 19, 2017; Published: February 14, 2017 This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication. Data Availability: All relevant data are within the paper. Funding: This study was supported by the Australian Research Council, the Research Corporation for Science Advancement, and the W. M. Keck Foundation. Competing interests: The authors have declared that no competing interests exist.

Introduction In 1921, Herman Rorschach (1884–1922) published the images of ten inkblot patterns that were believed to unlock the hidden secrets of the human unconscious [1]. Promoted as a ‘psychological X-ray’, psychiatrists adopted the blots as a probe of mental health based on the phenomenon of pareidolia in which familiar patterns are perceived in stimuli when none actually exist. Rorschach patterns remain unrivaled in their application since quite possibly “no other psychological test has been administered to more millions of people throughout the world” [2]. Today, their role as a projective psychological tool has only historical value [3]. Instead, the percepts induced by the blots are interpreted as a signature of observers’ creativity [4]. In support of this view, artists ranging from the Surrealists to Andy Warhol have employed blot patterns as a vehicle to trigger their imaginations. The remarkable propensity of these simple shapes to induce such a rich variety of imagery, with up to 300 different percepts recorded for each blot [5], remains intriguing. If the pattern characteristics that stimulate this prolific visual activity can be identified and quantified, the resulting advances could potentially impact on diverse applications ranging from camouflage design to artificial vision. A previous empirical investigation of pareidolia [6] was motivated by the observation that humans readily perceive and identify meaningful images in many naturally occurring but largely unstructured configurations ranging from clouds, rocks, and cracks in the ground to the surface of the Moon. All of these structures are fractal, featuring patterns that repeat at increasingly fine magnifications [7]. The study explored the relationship between the scaling characteristics of fractal patterns, as quantified by their fractal dimension D [7], and the patterns’ ability to evoke the perception of namable objects. Using computer-generated stimuli, fractals with low D values were found to elicit a higher ability to evoke the perception of namable objects. This was the same regardless of whether the fractal patterns contained only contours or filled-in regions. A more recent study showed that incorporating left-right symmetry (a central feature of inkblots) into computer-generated fractal stimuli further increased pareidolia [8]. When observers experience pareidolia, investigations of neural responses reveal that regions of the brain associated with object recognition are activated [8, 9]. For a diverse variety of patterns, ranging from fractals and overhead satellite imagery to magazine covers, ratings of ‘conspicuity’ (the ease of noticing images) elicited by these patterns were remarkably consistent within and between observers and remained so over a period of at least one year [10]. These results strongly suggest that perceived conspicuity is driven by pattern-inherent sensory factors common to all observers. Consistent with this pattern-driven approach to perceived images, we present analysis showing that the number of induced namable images perceived in Rorschach patterns is related to fractal characteristics that occur at the edges of the blots. Our analysis, of course, does not address the idiosyncratic inter-individual variability in the type of images perceived in such patterns, but convincingly shows that the number of induced images can be related to the fractal scaling parameters. The results of our analysis also provide an explanation for the frequent appearance of recognizable imagery within natural objects such as clouds, rock faces and coastlines.

Images perceived in Rorschach inkblots A key question in our analysis concerns the relationship between the fractal properties of the blot boundaries and the number of images perceived when observing the blots. To address this, we re-visited two original empirical investigations concerning the number of responses to individual Rorschach blots. The first was the compilation of responses for each of the Rorschach blots conducted by psychologist Marguerite Hertz [20]. As the founding member of the Rorschach Institute, Hertz pioneered standardized scoring of the inkblot tests in the 1930s and produced the Frequency Tables for Scoring Rorschach Responses. These tables were created by meticulous cataloguing of different percept types (e.g. an image of a bat, a person, etc.) for each of the ten Rorschach blots generated by a cumulative sample of N = 1050 subjects aged 11–19 [5]. For the purpose of our analysis, we simply counted the number of different percept types reported for the five black blots and plotted these empirically determined values n against their respective D values obtained from our analysis (Fig 6, red symbols). Note that, just as D quantifies the boundary of the whole blot, so too n quantifies the number of percepts types induced by the whole boundary, and does not include percepts induced by component regions of the blots. The data reveal a decreasing trend indicating the importance of D in observing the number of induced percepts. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 6. Number of percepts, n, induced by the five black blots plotted as a function of the D value of their boundaries. The inset tabulates the D values of the five blots along with the corresponding data from Hertz [5] and Wirt and McReynolds [21] respectively. For the purpose of direct comparison, the raw n scores have been standardized as z-scores and are shown in parentheses. The horizontal dotted line represents the mean score for each of the two sets of data. The red and blue lines are linear fits to the respective data and the dashed lines represent the boundaries of 95% confidence intervals. Note that three of the blue data points (for Blots 1, 5 and 6) are obscured by the equivalent red points. https://doi.org/10.1371/journal.pone.0171289.g006 A strikingly similar trend can also be observed with the second independent empirical sample of percepts for each Rorschach blot, reported by Wirt and McReynolds in 1953 [21]. The total number of responses is one of the standard variables in an administration of the Rorschach test (the so called Rorschach variable “R”) and Wirt and McReynolds were concerned with the extent to which this measure is reliable, reproducible and consistent from one testing occasion to another. For that purpose, they determined the mean and standard deviation of the total number of responses for each blot in the samples of normal (N = 76), neurotic (N = 32) and schizophrenic (N = 50) subjects. The blots were then divided into two groups of five blots to determine whether the responses to the two groups were similar. Correlations between the scores for the two groups of blots equaled 0.88, 0.91 and 0.77 for the normal, neurotic and schizophrenic sample respectively, indicating high reliability (a score of 1 would represent a perfect correlation). While the relatively high reliability of the total number of responses is encouraging in itself, we took the opportunity to compare the reported number of total responses for each of the five black blots to those estimated based on the responses compiled by Hertz [5]. We note that the values recorded by Wirt and McReynolds represent the average number of percepts one person would report per blot in one administration of the test. More precisely, each value refers to the average number of responses given to the question “What might this be?” for each individual blot. On average, one person would see from one to three shapes in each blot, totaling 20 to 30 for the entire test. Thus, the n values for the Wirt and McReynolds test are much smaller than those from Hertz data, which are based on the cumulative number of percept types reported by more than 1000 participants. The Hertz numbers indicate the variability of responses on a group level while the Wirt and McReynolds numbers are to some extent an indicator of variability on the level of one individual. Despite these differences in recording n, the standardized z-scores for the average number of total responses reported by Wirt and McReynolds, plotted as blue symbols in Fig 6, closely follow the same inverse relationship between the n and D values found in the Hertz data.

Images perceived in computer generated fractals Naturally, the boundary’s fractal structure is not the only cause of the large number of percept induced by the black blots–other likely factors include the blot’s (left-right) symmetry [8] and also shading caused by variations in the ink’s opacity in the regions enclosed by the boundaries. One way to isolate the contribution of variations in fractal dimension is to investigate the perceived images in computer generated fractal patterns varying solely in D. A second motivation for using computer generated stimuli is that the blots study was inevitably limited to just ten stimuli. By extending the studies to computer stimuli, we investigated responses to 24 additional images. Participants For this study, 23 first year Psychology students from the University of New South Wales were used (N male = 10, N female = 13) with an average age of 18.9 years. The students received extra course credit in exchange for their participation and had normal or corrected to normal vision. Prior to the start of the experiment, all participants signed an informed consent form. All procedures and protocols were approved by the University of New South Wales Human Research Ethics Approval Panel (HREAP C). Materials Three different ‘seed’ patterns were used to generate fractal images using a Fourier transform technique reported elsewhere [22]. Each seed set featured 8 images varying in D from 1.05 to 1.95. All images were 512 by 512 pixels wide and had the same mean luminance and contrast. One set of images is illustrated in the top section of Fig 7. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 7. The standardized number of reported percepts, n, induced by the computer-generated fractal patterns plotted as a function of their D value. One set of computer-generated images is illustrated above the graph. The dotted lines represent the boundaries of 95% confidence intervals. The inset tabulates the D values of the images along with the corresponding average number of reported perceived shapes. https://doi.org/10.1371/journal.pone.0171289.g007 Procedure Prior to data collection, participants received printed versions of images showing the variations in D and performed 4 practice trials. The study was conducted in a room with controlled lighting and participants were placed 60cm away from the monitor. All stimuli subtended a visual angle of 12 by 12 degrees and were centered against a uniform grey background. The presentation of each of the images was randomized and each image was presented once and for a duration of 10 seconds. Before each trial, a fixation cross appeared for 500 milliseconds. After each image was shown, participants were asked how many shapes they could discern in the image. The response options included “none”, “1 or 2 shapes”, “3 or 4 shapes”, “5 or 6 shapes” or “7 or more.” Results Before analysis, the answers were converted to numerical values by taking the average of the numbers in each choice except for “none” which was given “0” and “7 or more” which was given “7.5”. The reported number of perceived shapes, averaged across the three different sets of images, is shown in Fig 7. We then applied the standard statistical procedure, Analysis of Variance (ANOVA), to the data. One-way repeated measures ANOVA revealed the significant effect of D on the number of reported shapes (F 2.925,64.34 = 3.17, p = 0.03). The relationship between D and n observed for the Rorschach blots and our computer-generated fractal patterns is also consistent with a previous study of computer-generated fractal patterns [6] and confirms that low D fractal boundaries provide the best stimulus for induced associations with namable objects. Our study confirms the effect for an extended range of D and for much finer interval steps between different D values. More crucially, the previous study simply displayed groups of four images with different D values and recorded the percentage of times each image was chosen over the others based on which image induced more precepts. In contrast, our experiment recorded the number of percepts induced by each image, allowing a direct comparison with the original blot tests.

Discussion Fractal dimension plays a crucial role in determining a pattern’s visual appearance. Returning to the scaling behavior of Fig 3, D corresponds to the gradient of the scaling plot. A high D value is therefore a signature of a large N(L) value at small L and reflects the fact that many small boxes are being filled by fine structure. This can be seen, for example, for the two computer-generated fractal boundaries shown in Fig 4. The fine features play a more dominant role for the D = 1.77 pattern than for the D = 1.26 pattern. The fine feature content for the D = 1.26 pattern is much closer to that of Blot Five (D = 1.23). Because of this relationship between D and fine structure content, D is a well-established tool for quantifying fractal complexity [7, 14]. Traditional measures of visual patterns quantify complexity in terms of the ratio of fine structure to coarse structure. D goes further by quantifying the relative contributions of the fractal structure at all the intermediate magnifications between the coarse and fine scales. Previous psychophysical experiments performed on fractal patterns confirm that raising the D value increases its perceived complexity [23–25]. Thus, the increase in number of induced percepts between the different blots might be linked to a reduction in their fractal complexity. In addition to the dependence of n on D, the importance of the fractal boundary for inducing percepts is further highlighted by the fact that the non-fractal, multi-colored blots induce fewer percepts (n = 140 to 170) than their fractal, black counterparts (n = 170 to 300) [5]. We note, however, that there may be additional causes for this drop in n. In particular, for the multi-colored blots, the component regions (each with a distinct color) might visually dominate over the whole blot. The boundary shapes of these component regions would then predominantly determine the percepts rather than the whole blot’s boundary. This would result in a lower n because n quantifies only percepts induced by the whole blot. This effect is consistent with findings suggesting the powerful role of color on image segmentation [26]. Although the n values are smaller for the colored blots than their black counterparts, they nevertheless still have significant magnitudes. In our current study, we focused on the box-counting dimension D because of the range of previous psychophysical experiments that linked the visual properties of fractals to D [6, 17, 22–25, 27–37]. However, as noted earlier, future studies should explore if other scaling parameters revealed by a multi-fractal analysis can explain the percepts induced by the colored blots. It is interesting that the fractals play such an important role given that the structure is limited to size scales smaller than approximately 2.5cm (patterns larger than 2.5cm weren’t generated by the fractal ink diffusion process but by Rorschach’s ink smearing actions prior to pressing). In Fig 8, we therefore demonstrate the visual importance of the fractal structure by eliminating this structure from the boundary of Rorschach Blot Seven by Fourier transforming the blot image and removing spatial frequencies corresponding to the fractal scaling regime (2-25mm). This process clearly impacts on the visual perception of the blot. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 8. (A) Rorschach Blot Seven, (B) Rorschach Blot Seven with the fractal features removed. https://doi.org/10.1371/journal.pone.0171289.g008 The human visual system’s ability to detect fractal characteristics based on only a limited magnification range has been highlighted in the previous psychophysical experiments [6, 17, 22–25, 27–37]. This remarkable sensitivity to fractals is perhaps not surprising given that many of the physical fractals that define our daily visual environment exhibit fractal characteristics over limited magnification ranges [38]. This relationship with nature highlights another crucial factor for inducing perceived images—the statistical quality of the fractal pattern. Fig 9 shows a fractal pattern based on the Koch curve, which has a similar D value (1.26) to the Rorschach Blot Five of Fig 2 and the low D computer-generated fractal of Fig 4. In this figure, we morph the ‘exact’ fractal (where the patterns repeat exactly at different magnifications) into the more natural-looking ‘statistical’ fractal (in which the statistical qualities of the fractal repeat at different magnifications) by introducing random variations into the fractal pattern as follows. The probability p for a ‘spike’ on the curve to be pointing up or down changes for the three curves in Fig 9. For the traditional Koch curve p = 0, corresponding to zero probability of having a spike pointing down. For the second image p = 0.25, so most of the spikes are pointing up. For the third curve, p = 0.5, corresponding to a 50% chance of pointing up or down. In each case, the spatial distribution of the up and down spikes is random. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 9. Koch curves quantified by D = 1.26. The exact fractal (top trace) is morphed into the statistical fractal (bottom trace). See the text for details. https://doi.org/10.1371/journal.pone.0171289.g009 These random variations preserve the fractal scaling properties such as D while removing the artificial appearance of the exact Koch pattern, allowing the emergence of a more ‘organic’ visual character that we hypothesize enhances the pattern’s ability to induce associations with namable objects and/or recognizable imagery. Integration of random variations with fractal scaling of D = 1.1 appears to be optimal. We note that purely random white noise also has the capacity to induce percepts [8]. However, because white noise is characterized by D = 2 [39], the data of Figs 6 and 7 suggest that the lower D values (corresponding to pink noise [39]) will induce more percepts than purely random patterns. Natural objects whose fractal characteristics might be expected to induce recognizable images include rock faces, coastline patterns, clouds and craters on the Moon. In particular, many clouds have fractal boundaries quantified by D = 1.3 (similar to the Rorschach blots) and are well-known for inducing percepts in our daily lives. Biological examples of fractals that induce associations with namable objects include animal markings [40], a phenomenon that inspired research of fractal camouflage [41]. Fractal-induced imagery is also evident in artworks. In particular, the Surrealists developed several techniques to produce patterns that induce imagery. A prime example is Oscar Dominguez’s technique, decalcomania [42], for which he painted a surface, pressed a sheet of paper down on that surface and then pulled it off. The infusion of air into the paint as the two layers were separated generated a fractal process called viscous fingering [43]. Dominguez described his resulting patterns as “unequalled in [their] power of suggestion”, emphasizing the propensity of these simple fractal patterns for triggering a striking variety of imagery. The Abstract Expressionist Jackson Pollock’s poured paintings are composed of fractals with D values that increased from 1.1 to 1.7 over the decade 1943–1952 [32, 44]. Intriguingly, Pollock seems to have been aware that his drive towards higher complexity paintings would reduce the number of induced precepts: “I try to stay away from any recognizable image; if it creeps in I try to do away with it… I don’t let the image carry the painting… It’s extra cargo—and unnecessary” [45].

Conclusions Analysis of Rorschach inkblots provides an appealing framework for understanding the rich variety of visual associations induced by fractal patterns spanning psychology, art and nature. In each case, the repetition of structure at increasingly fine magnifications generates the visual complexity necessary to induce the perceived images. Within this fractal model, low D fractals provide the optimal distribution of spatial frequencies to induce the images. In contrast, the dominance of fine structure in high D fractals appears to reduce the perception of recognizable images. The ability to tune the number of percepts by adjusting the D value of computer-generated fractals has huge potential for both the visual arts and visual sciences. Given the recent emphasis on the connection between the images induced by the blots and the observer’s creativity [4], it is intriguing to note that self-reported creative people prefer to look at fractals with higher D values [46]. It would therefore be interesting to investigate how the n vs D dependence identified in Figs 6 and 7 is affected by observers’ creativity. Another potential focus for future investigations concerns the relationship between percepts and aesthetics. We note that a previous study of the fractal boundaries of blots focused on their aesthetics [47]. However, the blots were analyzed after undergoing filtering for red, blue and green colors. This difference between their images and our unfiltered ones prevents a comparison between the two studies.

Author Contributions Conceptualization: RPT. Formal analysis: TPM APM CB BCS MSF JHS RDM. Methodology: RPT BS. Project administration: RPT. Software: TPM APM CB BCS MSF RDM JHS. Validation: RPT. Writing – original draft: RPT BS. Writing – review & editing: RPT BS TPM APM CB BCS MSF JHS RDM.