$6.37 LIBRARY Qan State niversity Michi U This is to certify that the thesis entitled A COMPARISON OF USER PERFORMANCE ON SPECTRAL COLOR AND GRAYSCALE CONTINUOUS- TONE MAPS presented by MICHAEL D. HYSLOP has been accepted towards fulfillment of the requirements for the Master of degree in Geography Arts MN //a/"‘7/ Maidr'Professor’s Signature 7%)”? Date MSU is an affirmative-action, equal-opportunity employer -.-.--v-n-—-— —-—-—n-o-.-I-n-r—-o---o--.-----.--.-.--.-u—.-.-u-.-.—----o----o--.-- PLACE IN RETURN BOX to remove this checkout from your record. To AVOID FINES return on or before date due. MAY BE RECALLED with earlier due date if requested. DATE DUE DATE DUE DATE DUE 6/07 p:/ClRC/Date0ue.indd-p.1 A COMPARISON OF USER PERFORMANCE ON SPECTRAL COLOR AND GRAYSCALE CONTINUOUS-TONE MAPS By Michael D. Hyslop A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of MASTER OF ARTS Department of Geography 2007 ABSTRACT A COMPARISON OF USER PERFORMANCE ON SPECTRAL COLOR AND GRAYSCALE CONTINUOUS-TONE MAPS By Michael D. Hyslop Continuous surfaces such as elevation, precipitation, and temperature are often mapped using isarithms. These representations may be difficult for map readers to understand. Continuous-tone maps are an alternate way of depicting continuous surfaces. Three real surfaces were used to generate continuous-tone maps using spectral color and grayscale color schemes. To assess the effectiveness of these maps, fifty nine subjects were tested. Questions were designed to evaluate how the color scheme affected perception of the surface, and to compare performance on specific map reading tasks. Subjects were shown a grayscale and a spectral color representation of each of the three surfaces. Questions involved the location of surface extremes, elevation estimation, profile identification, landscape position, and interpretation of surface form. Subjects performed significantly better on the spectral color maps than on the grayscale maps. T o my family, for their support and encouragement iii ACKNOWLEDGEMENTS My thanks to the Department of Geography and the Center for Remote Sensing (now Remote Sensing & Geographic Information Science Research and Outreach Services) for financial assistance as a Teaching and Research Assistant while I was at MSU. My thanks also to Dr. Richard Groop and Dr. Judy Olson for their guidance throughout the course of my research. It would not have been completed without their help, suggestions, and answers to my many questions and problems. Dr. David Lusch served as third reader and provided invaluable comments and final edits to this manuscript. In addition, I must acknowledge my colleagues at Michigan Technological University. Their persistence in encouraging the completion of this thesis is appreciated. Finally, I would like to thank my parents for always supporting learning, especially my father for suggesting I “take a few Geography classes”; and special gratitude to Mary, Ian, Sean, Wes and Seth for putting up with my many late nights and weekends over the past year. iv TABLE OF CONTENTS LIST OF TABLES ...................................................................................................... vi LIST OF FIGURES ................................................................................................... vii Chapter I An Introduction To Quantitative Surface Maps ................................... 1 Chapter 11 Statement of Problem ......................................................................... 10 Chapter III The Map Test ...................................................................................... 13 Chapter IV Results ................................................................................................. 25 Chapter V Summary and Conclusions ................................................................. 52 References ........................................................................................... 55 Appendices Appendix 1 — Consent Form .................................................. 59 Appendix 2 — Test Booklet Part I ........................................... 60 Appendix 3 — Test Booklet Part 11 Version A ........................ 68 Appendix 4 - Test Booklet Part II Version B ........................ 75 LIST OF TABLES Table l. The Twelve Test Illustrations .......................................................................... 15 2. Order of Triads in Test Booklet Section I ....................................................... 18 3. Order of maps in test booklet Section II a ....................................................... 22 4. Order of maps in test booklet Section II b ....................................................... 22 5. Summary of responses, Test Section I ............................................................. 25 6. Kruskal-Wallis H test for differences between test groups ............................. 28 7. Linear mixed effects analysis to test for learning, aggregate ........................... 3O 8. Linear mixed effects analysis to test for learning, separate ............................. 31 9. Average scores for map types and the twelve maps, Section II ....................... 33 10. Chi-square analysis of test booklet Section I ................................................... 36 l 1. Chi-square tests of differences, color and grayscale maps .............................. 4O 12. Wilcoxon signed rank test comparing color and grayscale maps .................... 4] 13. Chi-square test of differences, color vs. gray maps, by question .................... 43 14. Chi-square test comparing profile magnitude selections ................................ 44 vi LIST OF FIGURES Figure 1. 2. 10. ll. 12. l3. 14. 15. 16. l7. l8. Stepped and smooth statistical surfaces ............................................................. 3 USGS isarithm map ........................................................................................... 4 Shaded relief map combined with elevation tint ............................................... 5 Perspective block diagram ................................................................................. 6 Isoline and continuous-tone grayscale maps of elevation ................................. 6 Groop and Smith hexagon snowfall map ........................................................... 7 Johnson dot-matrix map .................................................................................... 8 Lavin quasi-random dot map ............................................................................. 8 Visual sensitivity and spatial frequency curves ............................................... 10 Growing Degree Days for Lower Michigan, 1993 ........................................... l3 Elevations for a portion of Keweenaw County, Michigan ............................... 14 SEMCOG rainfall, 1993 ................................................................................... 14 Orientation and scheme of test illustrations ..................................................... 16 Sample grayscale and color map triads from test Section I ............................. 19 Sample maps from test Section II .................................................................... 20 Test questions from Section II ......................................................................... 21 Average map score differentiating color and grayscale maps ......................... 32 Maps B3 and B5 .............................................................................................. 47 vii Chapter I An Introduction to Quantitative Surface Maps Most maps are two-dimensional representations showing objects relative to a set of x and y coordinates. A typical example would be a standard highway map showing the locations of roads, railroads, cities, lakes, etc. Other examples include maps illustrating land use, forest cover type, soil series, or election results by state or county. Standard symbols may be used to show these features, and they tend to be conceptually straightforward. In the last 50 years, there has been growth in the production of quantitative maps. Quantitative maps are more complex than simple reference maps or maps showing what is where rather than how much of something is there. Quantitative maps Show the location of features measured on ordinal, interval, or ratio scales. Ordinal data provides the map user with information about rank or hierarchy, for example, populated places classified as a city, town or village. Interval and ratio data are continuous and provide more detailed and precise information. They employ a scale of measurement, often revealed in the map legend. Interval data have no natural zero and ratios are meaningless; temperature in degrees Fahrenheit is interval data (2 degrees is not twice as warm as 1 degree). Ratio data have a natural zero point and ratios make sense; rainfall in number of inches per year is ratio data (twenty inches per year is twice as much as ten inches per year). Ordinal, interval, and ratio data as a group are referred to as quantitative data, and may be contrasted with qualitative, or nominal data (Slocum, 1999). Qualitative data show categories that have no quantitative relationship to one another, such as the cities, roads, and lakes on the road map. Quantitative maps, and the surfaces they represent, require innovation and map reading skills beyond the mere visual recognition of standard symbols. The number of methods developed for mapping quantitative surfaces (discussed below) support this idea. A common feature shown on a quantitative map is elevation of the earth’s surface above or below a datum. Elevation is not the only phenomenon that may be thought of as a “volume” to be mapped, however. Any feature that changes quantitatively as it varies across space may be represented by a statistical surface (J enks, 1963). Depending on how data are enumerated, or on the type of feature being mapped, there may be discontinuities in the data. These gaps may be caused by abrupt changes at the boundaries of the mapping units, or by voids where the mapped feature does not exist. Surfaces such as these are known as stepped statistical surfaces (Figure 1, top maps). An example of a stepped surface is population mapped at the county level. Smooth statistical surfaces consist of data that vary continuously, not discretely, over geographic Space (Figure 1, lower maps). Some examples of smooth statistical surfaces include elevation above sea level, air pressure, ocean temperature, and precipitation. To map statistical surfaces, it is necessary to Show both position and change in magnitude. Figure 1. Stepped and smooth statistical surfaces, fiom Groop and Smith, 1982. Reprinted with permission from The American Cartographer 9, October 1982, 123-131. Cartographers have developed a number of ways to map smooth statistical surfaces. Perhaps the oldest and most often used method of portrayal is the isarithmic, or isoline, map (Figure 2). Though appropriate for many purposes, isarithmic maps can be difficult to read because they are two-dimensional representations of three-dimensional data (J enks, 1963). Studies have shown that map users may perform well on some tasks when using isarithmic maps, such as estimating the value at a location, but may have difficulty understanding the overall surface (Phillips et al, 1975; Griffin and Look, 1979; Phillips, 1979). The lines are not continuous symbols across the field of the map, so readers must interpolate values between them (Lavin, 1986). Because of these difficulties, cartographers have sought alternate methods for depicting smooth statistical surfaces. Figure 2. USGS topographic quadrangles include isarithms (“contours” in this case because they represent land elevation). Contours are brown on these maps. Hill shading is one alternative method for showing continuous surfaces. With hill shading, various tints Show the aspect of the statistical surface, with northwest-facing slopes illuminated and southeast-facing slopes in shadow. This combination of illumination and shadow gives a 3-d efi'ect to the surface being mapped. It is difficult to produce quality hill shadings manually, so automated (computer) methods have been . developed. Early automated methods could not match the quality of shadings produced manually (Y oeli, 1967; Brassel, 1979), but more recent improvements in computing power and software have made shaded relief maps not only faster and easier to produce, but also of extremely high quality (Thelin and Pike 1992, Lewis 1992). Hill shading gives a general impression of the form of the statistical surface, but as it shows aspect and not the values of the distribution, it is not useful for estimating values on the surface. Hill shading is sometimes combined with tints that represent value (Figure 3). Figure 3. Shaded relief (hill shade) map combined with elevation tints, Baraga County, MI. Perspective block diagrams, also known as “fishnet” diagrams, perspective traces, or transects, Show elevation by placing lines across the statistical surface (Figure 4). Tedious to produce by hand, perspective diagrams can be produced quickly using computers. They do have their limitations, however. Regions of high variability—and thus high relief—can Obscure lower areas behind higher zones, requiring rotation about the vertical and horizontal axes to minimize masking. On some block diagrams, obscured areas cannot be eliminated due to the complexity of the distribution. There is no “rule of thumb” for producing the ideal representation, which will vary from surface to surface. Figure 4. Perspective block diagram, Baraga County, MI. The Continuous-T one Method Another way of portraying smooth statistical surfaces is the continuous-tone map. On continuous-tone maps, gray tones or colors (hues) vary smoothly to illustrate changes in value from high to low. Unlike an isarithmic map, they Show value everywhere (Figure 5). Figure 5. Isoline and continuous-tone grayscale maps of elevation, Baraga County, MI. Groop and Smith (1982) developed a method that used graduated hexagons for creating virtually continuous-tone maps where hexagon size was plotted proportional to data value. Over a field such as a county, apparent darkness then represented the value (Figure 6). The procedures for this technique were more suited to maps plotted with a pen plotter than current output technologies, and the method used extensive and (then) expensive computer time. The plots had to be photographically reduced to page size. Figure 6. Groop and Smith hexagon snowfall map. Reprinted with permission fiom The American Cartographer 9, October 1982, 123-131. Johnson (1984) introduced a procedure that employed a dot matrix printer, driven by a BASIC program, to create continuous shadings (Figure 7). Testing proved the dot matrix maps as effective as isoline maps at conveying information, but test subjects did not like their coarse appearance. Figure 7. Johnson dot-matrix map. Lavin (1986) developed a technique that used dot-density shading to produce continuous tones. Dots were plotted quasi-randomly in a gridded data matrix to assure non-overlap (Figure 8). Output was produced on a TektronixTM plotter. Lavin concluded that dot density maps are one alternative method of smooth surface portrayal that can be used to supplement, but not replace, isoline maps. Figure 8. Lavin quasi-random dot map. Reprinted with permission fi'om The American Cartographer 13, April 1986, 140-150. Kumler (1988) created a hybrid method for producing color continuous-tone maps, which used both digital and photographic techniques. He developed a series of BASIC programs that displayed gridded data on a high-resolution color computer monitor in shades of red, green, and blue. He then photographed the monitor, one color at a time, and the color-separated negatives were process-printed to produce test maps. These procedures were chosen to suit the technologies of the time. Kumler’s subject testing indicated improved map reading performance on the color continuous-tone maps when compared to standard isoline maps for the questions he asked, which included locating and marking surface lows and highs, estimating the relative elevation of a pair Of points, estimating the elevation at a point, and interpreting the slope of a line between map points. Chapter II Statement of Problem On continuous-tone maps, color scheme choice and surface complexity can affect the characteristics of the distributions that map readers perceive. Questions arise about the efficacy of the maps and just what information people are able or likely to gain from them. For example, Rogowitz and Treinish (l 994) report that broad surface variations (low spatial frequency data) are more easily visualized when mapped using colors, but local surface variations (high spatial frequency, or fine details) are more evident when mapped using gray tones (Figure 9). Stated differently, the details (high spatial frequency components) A v Luminance —m:m—-< ~ 0 and "-" when Xa-Xb < 0. The signed ranks are then summed and used to calculate W. R version 2.41 was used for this calculation, and the resulting W statistic was 1066, with a p-value of 0.00014. This result, like the chi-square test, indicates there were significant performance differences between the color maps vs. the grayscale maps: subject performance was significantly better on the spectral color maps. The result of the Wilcoxon Signed Rank Test is listed in Table 12. Table 12: Wilcoxon Signed Rank Test comparing color and grayscale map performance Ho: Performance differences are evenly distributed (a = 0.05) W value Probablility Decision 1066 0.0001382 Reject Ho 41 Testing the Other Research Hypotheses, Test Section II Section II of the test was designed to assess how well subjects could interpret the three surfaces when shown in spectral colors vs. gray tones. Research hypotheses three through six addressed specific map reading tasks from Section II of the test: H4 H5 H6 Subjects will select the profile of correct shape with least vertical exaggeration as correct on grayscale maps. Subjects will select the profile of correct shape with higher vertical exaggeration as correct on spectral color maps. Subjects will estimate elevation more accurately on spectral color maps. Subjects will locate surface extremes more accurately on spectral color maps. A chi-square analysis was performed on the results from Section II of the test on a question-by-question basis to compare the performance differences between spectral color and grayscale maps. The results of this analysis are listed in Table 13. 42 Table l 3: Chi-Square Test Comparing Spectral Color and Grayscale Map Performance Ho: No difference in performance with color vs. grayscale maps (a = 0.05) Aggregate Scores Number Number correct correct ‘ Question Color Gray X 2 Probability Decision 1 171 170 0.0029 0.9568 Accept Ho 2 165 160 0.0769 0.7815 Accept Ho 3 117 132 0.9036 0.3418 Accept Ho 4 110 77 5.8235 0.0158 Reject Ho 5 99 80 2.0168 0.1556 Accept Ho 6 131 142 0.4432 0.5056 Accept Ho 7 130 95 5.4444 0.0196 Reject Ho 8 72 75 0.0612 0.8046 Accept Ho Significant chi-square values are in bold A brief assessment of the results of this chi-square analysis would suggest that subjects performed equally well on the spectral color and grayscale maps on six of the eight associated questions. Only questions four and seven produced chi-square values that were significant. These two questions (four and seven) dealt with the landscape position of a point and the form of a slope between two points. Subjects performed significantly better on the spectral color surfaces on questions four and seven. To assess research hypotheses three and four, which suggested that higher magnitude profiles would be chosen on spectral color maps and lower magnitude profiles on grayscale maps, the compiled answers to question seven were scrutinized. There were 130 correct answers selected on color maps: 59 of these were the profile of lesser magnitude, and 71 were higher magnitude. On the grayscale maps, 95 correct answers were recorded: 50 were the lesser-magnitude profile and 45 the higher magnitude. At first 43 glance, these results support the research hypotheses, as a majority of subjects preferred the higher-magnitude profile on color surfaces and the lesser-magnitude profiles on grayscale surfaces. To test these observations for significance, the correct answers—— which included both the low and high magnitude profiles—were subjected to chi-square analysis. Both the color and grayscale profile results did not have chi-square values large enough to be significant. Research hypotheses three and four were rejected and the null hypotheses of no difference was accepted. The results of the chi-square analysis comparing actual and high-magnitude profile selection are presented in Table 14. Table 14: Chi-Square Test Comparing Profile Magnitude Selections, ' Spectral Color and Grayscale Maps Ho: Performance differences are evenly distributed (or = 0.05) Number Number correct correct Scheme Actual Higher x 2 Probability Decision type Magnitude Color 59 71 . 1.1077 0.2926 Accept Ho Grayscale 50 45 0.2632 0.6080 Accept Ho Research hypothesis five predicted that subjects would be able to estimate elevation at a point more accurately on spectral color maps. Though performance was better on the color maps when compared with the grayscale maps, with 99 and 80 correct responses respectively, the difference in performance was not statistically significant. The resulting chi-square value was 2.01 , below the necessary significance threshold of 3.84. Research hypothesis five was, therefore, rejected and the null hypothesis of no difference was accepted. 44 Research hypothesis six proposed that test subjects would be better able to locate surface extremes—the high and low points—on spectral color maps than on grayscale maps. Test questions one and two dealt with marking the surface extremes on each map. With 171 (color) and 170 (gray) correct answers for question one, and 165 and 160 correct answers for question two, performance was again (marginally) better on the color maps when compared with the gray surfaces. The resulting chi-square values—0.0029 and 0.0769 for questions one and two—were far from significant enough to accept the hypothesis. Research hypothesis six was also rejected and the null hypothesis accepted. Summary of Results The test instrument discussed above was designed with two goals in mind. The first was to see if gray tones better represent high spatial fi'equency data, and spectral colors (hue) better represent low spatial frequency data. The second was to compare map reading performance on spectral color vs. grayscale representations of the same statistical surfaces. The first part of the test was designed to help assess how subjects saw spatial fi'equency. Research hypotheses one and two predicted that a majority of subjects would choose a residual map as most like an original surface when shown in grays, and a trend surface as most like the original surface when mapped in spectral colors. A chi-square analysis of the results fi'om Section I showed significant differences in the selections made by test subjects. As a result, both research hypotheses one and two were accepted. In Section II of the test, subjects were asked to answer a series of eight questions on each of six maps — three surfaces shown in two color schemes: once in grayscale and 45 once in spectral color. Differences in performance on these map reading tasks were tested for significance, to assess whether the color scheme type—spectral color or grayscale—affected performance. Subjects scored higher overall on the color surfaces when compared with the grayscale surfaces—about six percent better on average. A chi- square analysis of the subject performance showed that subjects performed significantly better on the spectral color maps than on the grayscale maps. When performance on selected indiVidual questions were subjected to the chi-square test, however, research hypotheses three through six, which dealt with profile selection, elevation estimation, and location of surface extrema, were rejected, as significant differences between the grayscale and spectral color surfaces were not evident in the subject response data. Discussion Some oversights in test design and production occurred while carrying out this research. In addition to the inverted color schemes used in triads six and seven, one other error was present in the test booklets. In Section II of test “B”, map five was rotated to help reduce the chance that it would be recognized as the same surface as map three, but the corresponding map points were not moved accordingly. Answers were adjusted so they were correct relative to the points as actually shown, but the comparability of the two maps was not as well controlled as in the other pairs. Maps B3 and B5 are shown in Figure 18. 46 Figure 18. Maps B3 (left) and B5 (right, map points unrotated). This omission reduced the effectiveness of directly comparing B5 with its paired map, B3. Like three of the other six map pairs, however, there was no significant difference in subject performance on map pair B3/B5, so the oversight may not have effected a significant change in answers. Conversations with some of the subjects after testing revealed an interest in the study and questions about how the maps might be used. One subject mentioned that he felt the maps were easy to use, but found the color maps easier than the grayscale, which were “too dark”. This may be only a preference due to the appearance of the bright spectral color scheme used, but it also suggests that the overall darkness of the grayscale maps may have been part of the reason for inferior performance with them. On the other hand, the values ran virtually the full gamut fiom white to black and any changes to lighten the maps would have meant less differentiation in the lighter end of the scale. The increased visual differentiation afforded by the spectral scheme is much more likely responsible for the better performance. It is interesting that all the research hypotheses involving specific questions were rejected (i.e., the null hypotheses were accepted). Since the one involving question 7 47 (selection of profile) did not deal with correct vs. incorrect answers, but whether subjects answering correctly selected the less exaggerated or more exaggerated profile, its rejection is not contrary to the overall hypothesis of better performance with color maps. The others, however, involving Questions 1, 2, 3, 5, and 6, each showed no difference in performance with the two maps. The significant overall difference, then, was attributable to other questions, and Table 13 indicates these are 4 and 7. Perhaps a different criterion for assessing questions 1 and 2 would have yielded different results, e.g., a smaller radius in which answers were marked correct. However, it is likely that marking the “blackest” or “whitest” regions on a grayscale map is no more difficult than marking the “reddest” or “bluest” locations on a spectral color map, and the non-sigrrificant results are, therefore, not surprising. It may be that the non-linear nature of the spectral color schemes used in the test affected the results for question 3, as performance was marginally, but not significantly, better on the grayscale maps. Question 5 dealt with estimating elevation at a point on the map, which had better, but not significantly so, results on the spectral color maps. Again, a different criterion, such as a smaller margin for correct answers, might have given a different outcome. Questions 6 and 8 were both about roughness of the surfaces, and non-significant differences were found again. Performance was slightly better on the grayscale maps, so “roughness” may be a characteristic better communicated by the luminance changes of a grayscale presentation. Questions 4 and 7 were the two with significant chi-square values, and these questions addressed landscape position and slope profile selection. It may be that the non- linear appearance of the spectral color scheme aided in assessing the location of point A 48 on question 4, e.g., the rapid change between red and green (through yellow) or betWeen blue and green (through cyan) gave better clues to the form of the landscape than the steady progression of grays in the grayscale maps. Perhaps the same clues in the color table assisted the subjects in selecting profiles for question 7. What does all this mean for mapmaking? It seems clear that spectral color allows better performance on at least some of the questions, but not on all. It is particularly interesting that it was form of the distribution on which subjects performed better with the spectral scheme because it has been the non-sequential nature of spectral colors that has long made cartographers reluctant to use them for representing quantitative data. One potential area for additional research would be to test monochromatic color schemes against grayscale, and even more important, a spectral scheme with a progression of value (say, light yellow through medium-value greens through dark blues and even darker reds), since progressive schemes resulted in good performance in Mersey's (1990) study and Brewer (1997) recommends them for diverging data. Would we, for example, see the shift of Visual emphasis from broad patterns to local features if we compared a light-to-dark spectral scheme to the li ght-to-dark grayscale? And for the questions posed in Section II of the test, would we continue to see superior performance (perhaps even more superior performance) with a light-to-dark spectral scheme? For the questions they asked, Brewer et al. (1997) had good results with a diverging spectral scheme (yellow for the middle category, with light oranges to highly saturated red in one direction and light greens to highly saturated blue in the other), and even better results for more conventional bi-color diverging schemes. The spectral scheme used in this thesis research does not have continuous-appearing values and is not arranged as a diverging 49 scheme. This means certain hues—such as yellows and cyans that occupy a narrow portion of the color spectrum compared to greens or blues—can give the impression of greater importance (because of greater contrast) than is warranted to the areas mapped in those colors. This phenomenon can be especially pronounced on some data distributions, where the combination of data values, surface form and color scheme can cause representations that are not true to the surface due to visual artifacts. Another area for further research is the perception of broad patterns on spectral maps and local variation on grayscale ones. The fourth-order surfaces showed strong differences in choices between trend and residual maps depending on whether the maps were in color or grayscale. But was it just because of the trial and error selection of the fourth-order that was based on our visual inspection? Is there something inherent in that level of trend surface? Rogowitz et a1 (1995) do not use trend surfaces to determine the appropriate color scheme to use for a given data set. Instead, their software uses what is, in effect, a low-pass filter to derive new data from the original. It then subtracts the original data values fi'om the filtered. The standard deviation of difference values is computed, then divided by the original data values. The resulting “normalized standard deviation” is used to assess the complexity of the surface. An image with a normalized standard deviation (N SD) of less than 0.1 is deemed “low frequency” and can use a spectral color map, whereas an image with a NSD of greater than 0.1 is flagged as “high frequency” and only monochromatic color maps are available. Is the Rogowitz method the best way to choose a color map? Not necessarily: they state that the method was chosen for its speed and simplicity, and it maylnot stand up to 50 rigorous scrutiny. It might not be the best for a data set that has both high and low frequency areas. This thesis research has added to favorable results with spectral color maps. It has also provided food for thought. 51 Chapter V Summary and Conclusions Mapping quantitative surfaces has been somewhat problematic for cartographers. A number of techniques have been developed to portray statistical surfaces, with varying degrees of success. Research has shown that many map readers have difficulty in perceiving quantitative surfaces when depicted with isolines, and numerous alternate techniques have drawbacks as well, both with reader perception and production by map makers. The purpose of this research was twofold: to develop a new technique for producing continuous-tone maps using both free and commercial software, and to test the effectiveness of these maps at representing quantitative surfaces by testing map readers. Moore’s Law, observed in 1965 by Gordon Moore of Intel, the CPU manufacturer, roughly states that computing power doubles every two years. This phenomenon has had several effects: computers that were once too large or too expensive for the average user have become ubiquitous over the last twenty years, and tasks that were once relegated to mainfi'ames or minicomputers can now be easily accomplished with widely available software on desktop computers. Continuous improvement has also occurred in output devices, and currently, hi gh-resolution color laser printers can be purchased for under $500. The map production technique outlined in this study was performed on desktop computers with software available commercially and for free, and the test illustrations were produced on a Hewlett-Packard color laser printer that is relatively stande equipment in modern computer labs and offices. The test instrument was designed to assess two things: first, whether high-spatial frequency data is best represented with grayscale color schemes and lower spatial 52 frequency data with spectral colors, and second, to compare subject map reading performance on spectral color and grayscale depictions of the same three quantitative surfaces. Three real surfaces fiom Michigan were used: elevation, rainfall, and growing degree days. To answer the color scheme vs. spatial frequency question, readers were shown a small version of each surface, along with two additional representations: a fourth-order polynomial trend surface derived from the original, and the residuals from the trend surface. These map triads were shown once in each of the four color schemes, for a total of twelve triads. A significant number of subjects chose the trend map as most like the original surface when shown in spectral colors, and the residuals map when shown in grays. These results indicate that a significant number of subjects perceive high spatial frequency features when surfaces are shown in grayscale, and they perceive broad trends when surfaces are shown in spectral colors. To compare map reading performance on spectral color and grayscale maps, readers answered questions about high and low spots on each map, the landscape position of points, the slope between map points, and surface roughness. Subjects saw each of the three surfaces twice, once in colors and once in grays. The maps were rotated to help reduce possible recognition of duplication, but map points were, with one inadvertent exception, in the same locations so direct comparisons from map to map could be made. Overall, the subjects performed significantly better on the color maps than on the grayscale representations. The highest average score was on maps shown in the blue-to- red spectral color scheme, and the worst performance was observed on the white-to-black 53 grayscale maps. Only two of the eight individual questions, however, showed significantly higher performance on the color maps. Mapping quantitative surfaces as outlined above is practicable. No unusual hardware or expensive software is necessary to produce the maps, and high-quality output may be obtained on an inexpensive color laser printer with ordinary paper. Continuous-tone maps are planimetric, and, unlike isopleth maps, present the map reader with a continuous-appearing symbology for interpretation. They represent the quantitative surface as faithfully as other maps produced from the same data. Previous studies have shown them to be more effective at conveying information about the mapped surface than isoline maps. The results from the fifty-nine subjects tested during the course of this research showed that spectral color, continuous-tone maps are easier for subjects to interpret than continuous-tone, grayscale maps. However, caution should be used when applying a spectral color scheme to a surface with potentially important components that are high in spatial frequency, as they could be masked. In these circumstances, regional components may look more important than is warranted, due to the non-linear nature of the spectral color scheme. 54 References Brassel, Kurt E. and Jack J. Utano. January 1979. “Design Strategies for Continuous- Tone Area Mapping” American Cartographer 6: 39-50. Brassel, Kurt. January 1974. “A Model for Automatic Hill Shading” The American Cartographer 1:1: 15-27. Brewer, Cynthia A., 1997. “Spectral Schemes: Controversial Color Use on Maps,” Cartography and Geographic Information Systems, 24(4): 203-220. Brewer, Cynthia A., Alan M. MacEachren, Linda W. Pickle, and Douglas Herrmann. 1997, “Mapping Mortality: Evaluating Color Schemes for Choropleth Maps,” Annals of the Association of American Geographers, 87(3): 411-438. Griffin, TLC. and BF. Lock. December 1979. “The Perceptual Problem in Contour Interpretation.” The Cartographic Journal 16: 61-71. Groop, Richard E. and Paul Smith. October 1982. “A Dot-Matrix Method of Portraying Continuous Statistical Surfaces.” The American Cartographer 9: 123-131. Jenks, George F. March 1963. “Generalization in Statistical Mapping” Annals of the Association of American Geographers 53: 15-26. Johnson, William F. 1984. “A Dot-Matrix Method for Representing Smooth Statistical Surfaces.” Unpublished Master's thesis, Michigan State University. Kumler, Mark. 1988. “Mapping Smooth Surfaces with Continuous Tone Maps.” Unpublished Master's Thesis, Michigan State University. Kumler, Mark P. and Richard E. Groop. October 1990. “Continuous-Tone Mapping of Smooth Surfaces.” Cartography and Geographic Information Systems 17: 279- 290. Kutner, Michael H. et al. 2005. “Applied Linear Statistical Models” 5th Ed. Duxbury Press. ISBN 0-07-238688-6 Lavin, Stephen. April 1986. “Mapping Continuous Geographical Distributions Using Dot Density Shading.” The American Cartographer 13: 140-150. Lavin, Stephen, Jay Hobgood and Paul Kramer. May 1986. “Dot-Density Shading: A Technique for Mapping Continuous Climatic Data.” Journal of Climate and Applied Meteorology 25: 679-690. Lewis, Peirce. 1992. "Introducing a Cartographic Masterpiece: A Review of the US. Geological Survey's Digital Terrain Map of the United States, by Gail Thelin and 55 all. Richard Pike," Annals of the Association of American Geographers, Vol. 82, No. 2 (June), 289-299. Mersey, Janet E. 1990. Colour and Thematic Map Design: The Role of Colour Scheme and Map Complexity in Choropleth Map Communication. Monograph 41, Cartographica 27 (3). Phillips, Richard J ., Alan DeLucia and Nicholas Skelton. 1975. “Some Objective Tests of the Legibility of Relief Maps.” The Cartographic Journal 12: 39-46. Phillips, Richard J. “An Experiment With Contour Lines.” December 1979. The Cartographic Journal 16: 72-76. Pike, Richard J. and Gail P. Thelin. 1992. "Visualizing the United State in Computer Chiarascuro," Annals of the Association of American Geographers, Vol. 82, No. 2 (June), 300-302. Rogowitz, Bernice E. and Lloyd A. Treinish. 1994. “Using Perceptual Rules in Interactive Visualization.” SPIE Proceedings Vol. 21 79 Human Vision, Visual Processing and Digital Display V: 287-295. Rogowitz, Bernice E. and Lloyd A. Treinish. 1995. “How NOT to Lie with Visualization.” Proceedings of the 1996 IBM Visualization Data Explorer Symposium. http://opendx.npaci.edu/cdsiproceedings96/pravda/truevis.html Slocum, Terry A. 1999. Thematic Cartography and Visualization. New Jersey: Prentice Hall. Smith, Paul. 1980. “Representing the Statistical Surface Using Continuous-Appearing Grey-Tone Variation.” Unpublished Master's Thesis draft, Michigan State University. Yoeli, Pinhas. 1967. “The Mechanisation of Analytical Hill Shading.” The Cartographic Journal 4: 82-88. 56 all. General References Bergman, Lawrence D., Bernice E. Rogowitz and Lloyd A. Treinish. 1995. “A Rule- based Tool for Colorrnap Selection.” Proceedings of the IEEE Conference on Visualization, 1 18-125. Brewer, Cynthia A. 1999. “Color Use Guidelines for Mapping and Visualization.” Chapter 7 from Visualization in Modern Cartography Castner, Henry W. and Roger Wheate. December 1979. “Re-assessing the Role Played by Shaded Relief in Topographic Maps.” The Cartographic Journal 16: 77-85. DiBiase, David and Thomas Paradis. 1994. “Weighted Isolines: An Alternative Method for Depicting Statistical Surfaces.” Professional Geographer 46(2): 218-228. Irnhof, Eduard. 1982. Cartographic Relief Presentation. New York: de Gruyter. Marchak, Frank M. et a1. 1993. “The Psychology of Visualization.” Proceedings of the IEEE Conference on Visualization: 351-354. Peterson, Michael P. 1999. “Cognitive issues in Cartographic Visualization.” Chapter 3 from Visualization in Modern Cartography. Rheigans, Penny. 1992. “Color, Change, and Control for Quantitative Data Display.” Proceedings of the IEEE Conference on Visualization: 252-259. Robertson, Philip K. 1990. “A Methodology for Scientific Data Visualization: Choosing Representations Based on a Natural Scene Paradigm.” Proceedings of the first IEEE Conference on Visualization: 114-123. Robinson, Arthur H. 1961. “The Cartographic Representation of the Statistical Surface.” International Yearbook of Cartography 1: 53-61. Rogowitz, Bernice E., Daniel T. Ling and Wendy A. Kellogg. 1992. “Task Dependence, Veridicality, and pre-Attentive Vision: Taking Advantage of Perceptually-Rich Computer Environments.” SPIE Proceedings Vol. 1666: Human Vision, Visual Processing and Digital Display 111: 504-513. Rogowitz, Bernice E. and Lloyd A. Treinish. 1993. “Data Structures and Perceptual Structures.” SPIE Proceedings Vol. 1913: Human Vision, Visual Processing and Digital Display IV: 600-612. Rogowitz, Bernice E. and Lloyd A. Treinish. 1993. “An Architecture for Rule-Based Visualization.” Proceedings of the IEEE Conference on Visualization: 236-243. 57 Tanaka, Kitiro. 1950. “The Relief Contour Method of Representing Topography on Maps.” Geographical Review 40: 444-456. White, Daniel J. 1974. “An Application of Continuous Tone Photography in Cartographic Design.” Unpublished Master's Thesis, Bowling Green State University. 58 PP!“ Appendix 1 Consent Form Map Study You are being asked to participate in a map perception study. The study is being conducted by Michael Hyslop as part of a master’s program in geography, under the direction of Drs. Richard Groop and Judy Olson of the Geography Department at Michigan State University. You will be asked to answer questions on maps provided in a test booklet. The length of time for participation is approximately 25 minutes. There are no foreseeable risks to you if you participate in this study. Your participation will not affect your grade in any course, and you are free to discontinue participation at any time without penalty. You are assured that you will remain anonymous in any reports of the research. Within these restrictions you may receive a copy of results upon request. You may receive further information about the study at the end of this session if you request. For further information about the research, you may contact: Michael Hyslop voice: 906-487-2308 MTU School of Forestry fax: 906-487-2915 101 UJ Noblet Forestry and Wood Products Bldg 1400 Townsend Drive Houghton, MI 4993 1-1295 01' Dr. Richard Groop or Dr. Judy Olson voice: 517/355-4649 Department of Geography fax: 517/432-1671 Natural Science Building Michigan State University East Lansing MI 48824 . For information about your rights as a subject in a research project, you may contact: David Wright voice: 355-2180 Univ. Comm. for Research Involving Human Subjects fax: 517/432-1 171 246 Hannah Administration Bldg. Michigan State University East Lansing, MI 48824 I understand the above statements and freely consent to participate in this study. Signed Date 59 Appendix 2 - Test Booklet Part I Test Booklet PLEASE DO NOT OPEN UNTIL INSTRUCTED. 60 Section I In this part of the test you will be shown maps in sets of three. Circle one of the two lower maps - the one that looks the most like the top map. Example: Which map, b or c. looks more like map a? (circle your choice) Go ahead and complete this section of the test 61 62 63 66 67 Appendix 3 — Test Booklet Part H Version A Section II In this part of the test you will be shown individual maps and asked a series of questions. Circle the correct answer, fill in the blank, or draw on the map as asked. Example: a t 25 50 '75 100 125 150 175 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map inc (1 line from point A to point B. This slope along most of this line is generally a) uphil b) downhill c) flat (circle one) 4) Point A is a) on a ridge or hilltop ”or depression c) on a slope benveen a ridge or hilltop r r . a ley or depression (circle one) 5) Estimate the elevation at point C 2 § 2 6) book at the circle to the right - 3 m . and imagine it centered around points T and U. The area around T is a) flatter , ;« lar than the area around U (circle one) V’N 8) Which quadrant is most irregular? (circle one) MI 68 1600 1800 2000 2200 2400 I) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from point A to point B. This slope along most of this line is generally a) uphill b) downhill c) flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and C. The area around T is a) flatter b) more irregular than the area around C (circle one) 7) Imagine a line from point C to point B. Circle the profile below that most closely matches its slope 8) / b) \ c) d) e) 0 8) Which quadrant is most irregular? (circle one) in \ / 9 MI 600 700 800 900 1000 1 100 1200 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from point A to point B. This slope along most of this line is generally a) uphill b) downhill c) flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and U. The area around T is a) flatter b) more irregular than the area around U (circle one) 7) Imagine a line from point C to point B. Circle the profile below that most closely matches its slope a) b) C) d) e) 0 J \r 8) Which quadrant is most irregular? (circle one) n n- 70 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from point A to point B. This slope along most of this line is generally a) uphill b) downhill c) flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and U. The area around T is a)fiatter b) more irregular than the area around U (circle one) 7) Imagine a line from point C to point B. Circle the profile below that most closely matches its slope 8) e) f) b) c) d M N“ .N m ~— -~ 8) Which quadrant is most irregular? (circle one) In M- 71 3200 3600 4000 4400 4800 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from point A to point B. This slope along most of this line is generally a) uphill b) downhill c) flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and C. The area around T is a) flatter b) more irregular than the area around C (circle one) 7) Imagine a line from point C to point B. Circle the profile below that most closely matches its slope 3) b) C) d) e) f) In 8) Which quadrant is most irregular? (circle one) u 72 I) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from pointA to point B. This slope along most of this line is generally a) uphill b) downhill c) flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and U. The area around T is a) flatter b) more irregular than the area around U (circle one) 7) Imagine a line from point C to point B. Circle the profile below that most closely matches its slope 8) b) c) d) e) f) 8) Which quadrant is most irregular? (circle one) II 1200 1300 1400 1500 1600 1700 1800 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from pointA to point B. This slope along most of this line is generally a) uphill b) downhill c) flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and U. The area around T is a) flatter b) more irregular than the area around U (circle one) 7) Imagine a line from point C to point B. Circle the profile below that most closely matches its slope a) b) e) t) «frxww 8) Which quadrant is most irregular? (circle one) u- Appendix 4 — Test Booklet Part H Version B Section II In this part of the test you will be shown individual maps and asked a series of questions. Circle the correct answer, fill in the blank, or draw on the map as asked. Example: 75 100 12 5 150175 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map ine a line from point A to point B. This slope along most of this line is generally a) uphil b) downhill c)flat (circle one) 4) Point A IS a) on a ridge or hiIIto b) m a valley or depr sion c) on a slope between a ridge or hilltop701 I i epression (circle one) 5) Estimate the elevation at point C 6) Look at u =- '- « right of the map and imagine it centered around points T and U. The area around T is a) flatter h more irregular an the area around U (circle one) . u - int C to point B. Circle the profile below that most closely matches its slope C) \/\\d )A/ e) \/\0 ~/\/ 8) Which quadrant' 15 most irregular? (circle one) I "m -m 7) Imagine a line a) 75 3200 3600 4000 4400 4800 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from point A to point B. This slope along most of this line is generally a) uphill b) downhill c) flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and C. The area around T is a)fiatter b) more irregular than the area around C (circle one) 7) lrnagine a line from point C to point B. Circle the profile below that most closely matches its slope a) \ b) c) d) e) o 8) Which quadrant is most irregular? (circle one) II 76 I) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from point A to point B. This slope along most of this line is generally a) uphill b) downhill c)flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and U. The area around T is a) flatter b) more irregular than the area around U (circle one) 7) Imagine a line from point C to point B. Circle the profile below that most closely matches its slope a) b) C) d) e) t) 8) Which quadrant is most irregular? (circle one) I u MI 77 600 700 800 900 1000 1100 1200 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from point A to point B. This slope along most of this line is generally a) uphill b) downhill c) flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and U. The area around T is a) flatter b) more irregular than the area around U (circle one) 7) Imagine a line from point C to point B Circle the profile below that most closely matches its slope a) ) c \W-r/d\/\ 8) Which quadrant is most irregular? (circle one) in 10 12 14 1B 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from point A to point B. This slope along most of this line is generally a) uphill b) downhill c) flat (circle one) 4) PointA is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and U. The area around T is a) flatter b) more irregular than the area around U (circle one) 7) Imagine a line from point C to point B. Circle the profile below that most closely matches its slope 8) b) C) d) e) f) ‘-\/\ l\/" W V\- -\/\ ~— 8) Which quadrant is most irregular? (circle one) an El 1200 1300 1400 1500 1600 1700 1800 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from point A to point B. This slope along most of this line is generally a) uphill b) downhill c) flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a valley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and U. The area around T is o)flatter b) more irregular than the area around U (circle one) 7) Imagine a line from point C to point B Circle the profile below that most closely matches its slope a) b) c) d) e) JM/“NLJ 8) Which quadrant is most irregular? (circle one) l“ ml 80 1600 1800 2000 2200 2400 1) Put an H at the highest point on the map 2) Put an L at the lowest point on the map 3) Imagine a line from point A to point B. This slope along most of this line is generally a) uphill b) downhill c)flat (circle one) 4) Point A is a) on a ridge or hilltop b) in a volley or depression c) on a slope between a ridge or hilltop and a valley or depression (circle one) 5) Estimate the elevation at point C 6) Look at the circle to the right of the map and imagine it centered around points T and C. The area around T is o) flatter b) more irregular than the area around C (circle one) 7) Imagine a line frbom point C to poicnt B. Circle the profile below that most closely matches its slope 3) 0 \W/\/\/ 8) Which quadrant is most irregular? (circle one) u- IIiii‘liifliiilflljljiilWill .........