Detecting subsurface objects and voids using thermal inertia

11243329 · 2022-02-08

    Inventors

    Cpc classification

    International classification

    Abstract

    A method for gathering reference data for use in planning and interpreting infrared surveys for the purpose of detecting and locating underground features, such as tunnels, voids, or manmade devices. Measurements, images, or observations at a site having known underground features are recorded. Recorded details include a combination temperatures at or near a soil surface at multiple points across the site in addition to above surface factors such as shading, weather conditions, and objects or foliage. Analysis of the details recorded from the site having known underground features yields quantitative estimates of the effects of various above and below surface factors on temperatures at or near the soil surface.

    Claims

    1. A method including: recording measurements of temperature from thermal probes distributed near a soil surface across a first site throughout a monitoring period, the first site having known subsurface features and the probes being arranged in discrete groups distributed across the first site, each discrete group including at least one probe below the soil surface and at least one probe above the soil surface; monitoring weather conditions at the first site and video footage of the site throughout the monitoring period; identifying interactions between above and below soil surface temperatures at particular locations within the first site, the monitored weather conditions, visible above surface features and conditions within the first site, and known subsurface features of the first site to produce a database for discovering unknown subsurface features at a second site by observing weather, above surface features and conditions, and temperature contrasts across multiple soil surface locations at the second site at multiple times of day.

    2. The method of claim 1, further including: recording measurements of lighting and weather conditions occurring at the first site during the monitoring period; and quantitatively estimating effects of the lighting and weather conditions on temperatures of the soil surface at the particular soil surface locations.

    3. The method of claim 1, further including: estimating effects of above surface features of the first site on temperatures of the soil surface at the particular soil surface locations by comparison of temperatures measured at particular soil surface locations having varying proximity to the above surface features.

    4. The method of claim 3, wherein at least one of the above surface features is a shadow.

    5. The method of claim 1, wherein the identifying step includes identification of thermal inertia signatures created by the known subsurface features at the soil surface.

    6. The method of claim 5, further comprising identifying a threshold wind gust speed, wherein subsurface features having certain characteristics that could otherwise be discoverable by monitoring soil surface level temperature at the second site on a given day will not be discoverable by monitoring soil surface level temperature following the occurrence of a wind gust exceeding the threshold wind gust speed on the given day.

    7. The method of claim 6, wherein the certain characteristics include depth.

    8. The method of claim 6, further comprising determining a time of day during which a wind gust must exceed the threshold wind gust speed to render subsurface features having the certain characteristics undiscoverable by monitoring soil surface level temperature on the given day.

    9. A method of detecting features below a soil surface of a survey site, the method comprising: compiling a database according to the method of claim 1; interpreting survey data of the survey site with reference to the database, the survey data being derived from an infrared survey including a first pass and a second pass, the second pass conducted at a later time than the first pass, and the database including measurements of temperature from thermal probes distributed near a soil surface across a surrogate site, the surrogate site having known subsurface features and the measurements of temperature being recorded at regular intervals of less than ten minutes; wherein the reference to the database includes comparison of the survey data to correlations between temperature changes over time at the particular soil locations within the surrogate site and the known subsurface features of the surrogate site.

    10. The method of claim 9, wherein the first pass and the second pass are scheduled by: selecting a target depth; and consulting the database to determine times of year and times of day wherein a known subsurface feature of the surrogate site at the target depth creates detectable signatures in temperature change over time at the surrogate site.

    11. The method of claim 9, wherein the database includes estimates of effects of at least one type of above surface feature on temperatures measured at the soil surface of the surrogate site, and the interpreting step includes accounting for effects of an above surface features of the same type at the survey site by reference to the estimates.

    12. The method of claim 11, wherein the type of above surface feature at the surrogate and survey site is shade.

    13. The method of claim 11, wherein the estimates are derived by comparison of temperatures at particular soil locations at the survey site having varying proximity to the above surface feature at the surrogate site.

    14. The method of claim 9, wherein the infrared survey includes use of a first infrared camera and a second infrared camera, the first infrared camera configured for shorter wavelengths than the second infrared camera.

    15. The method of claim 9, wherein the infrared survey includes use of a cold filtered infrared camera.

    16. The method of claim 1, wherein the second site is within the same geographic region as the first site.

    17. The method of claim 1, wherein the database includes soil surface level thermal inertia signatures of columns of earth containing manmade anomalies as contrasted to columns of earth free of manmade anomalies.

    18. The method of claim 17, wherein the database further includes times of day and times of year wherein thermal inertia contrasts between columns of earth containing manmade anomalies and columns of earth free of manmade anomalies will be most apparent from soil surface level temperature measurements.

    19. The method of claim 18, comprising: identifying, from the database, minimum soil surface level temperature change over time contrasts between columns of earth containing manmade anomalies and columns of earth free of manmade anomalies necessary for the thermal inertia signatures to be detectable from successive surveys of soil surface level temperature at the second site; and identifying, from the database, wind conditions that would reduce the soil surface level temperature change over time contrasts below the minimum on dates on which the soil surface level temperature change over time contrasts would otherwise meet or exceed the minimum.

    Description

    BRIEF DESCRIPTION OF THE DRAWINGS

    (1) FIG. 1 is a schematic plan view of an arrangement of measuring devices at a surrogate site according to an arrangement of the present disclosure.

    (2) FIG. 2 is a schematic elevation view of a group of thermal probes according to the arrangement.

    (3) FIG. 3 is a schematic elevation view of a survey process according to the present disclosure.

    (4) FIG. 4 is a flowchart of an imaging method according to the present disclosure.

    (5) FIG. 5 is a flowchart of an example of an outcome of an imaging method.

    DETAILED DESCRIPTION

    (6) As shown in FIG. 1, a surrogate site 10 is populated with an array of thermal probe groups 14 arranged at, near, and/or slightly above or below a soil surface of the surrogate site, a weather station 18, and cameras 22. Though not illustrated, foliage, rocks, debris, height variations, and other surface variations are also distributed across the surrogate site 10. Generally, more sensitive and accurate thermal probes provide more useful information. In some examples, the thermal probes used in the probe groups 14 are accurate to two tenths of a degree Centigrade or better. In further examples, the thermal probes used in the probe groups 14 are accurate to one tenth of a degree Centigrade or better. Environmental factors and sources of clutter can reduce the effective accuracy of thermal probes used in the field, so thermal probes may be chosen to be more precise or accurate than necessary for the desired data to be gathered. For example, where data accurate to within two tenths of a degree Centigrade is desired, thermal probes accurate to within one tenth of a degree Centigrade may be used to compensate for environmental factors and clutter. The weather station 18 has apparatus for measuring various weather conditions, such as wind speeds, ambient air temperature, cloud cover, and/or any other measurable weather conditions. The cameras 22 are oriented toward the probe groups 14 so as to enable observation and recording of visible conditions affecting the probe groups 14. In the illustrated example, the cameras 22 are on opposite sides of the surrogate side and oriented to face each other. The cameras are spaced apart by a suitable distance to put all probe groups 14 within the field of view of both cameras, such as between 10 and 20 meters, or 15 meters. Further, a single camera 22 or a greater number of cameras 22 may be used. Though not shown in FIG. 1, an infrared camera may also be placed at the surrogate site 10 for use in cooperation with the probe groups 14, weather station 18, and cameras 22.

    (7) Measurements taken by the probe groups 14, weather station 18, and cameras 22 are recorded continuously or at regular intervals over a given period of time. In one example, temperature measurements by the probe groups 14 are taken every five minutes, or at any regular interval below 10 minutes within an observation period. In various further arrangements, the regular interval is 8 minutes, 6 minutes, 4 minutes, less than 8 minutes, less than 6 minutes, less than 5 minutes, less than 4 minutes, or approximately 8, 6, 5, or 4 minutes. The observation period is a given time of day, such as from noon to 5 PM, or an entire 24 hour period, and is limited to a single day, or extends across multiple days that may or may not be consecutive. Typically, temperature contrasts and differences in thermal inertia between various above and below surface features are greatest in the spring and fall of a given year. Optimum dates for detecting deep below surface features are therefore in the spring or fall in most regions. Measurements from the weather station 18 and images or video from the cameras 22 are recorded on the same interval as the measurements from the probe group 14 or on a different interval. In some examples, average or maximum wind speed, or wind gusts, meaning sudden changes in wind velocity, measured by the weather station 18 are recorded for every three hour period within the observation period. Use of thermal probes at the surrogate site 10 enables recordation of accurate temperature measurements. Particularly, thermal probes will not measure infrared clutter that would be indistinguishable from temperature variations from an infrared survey alone, such as clutter resulting from reflected earth or sky spectral radiance, bodies of water, restrahlung rays, vegetation, pockets of hot or cool air, or reflective above ground objects.

    (8) An exemplary visible condition detectable by the cameras 22 is a shade line 26, which separates a shaded area 30 from an unshaded area 34. Video or sequential images recorded by the cameras 22 therefore shows which probe groups 14 are in the shaded area 30 or unshaded area 34 throughout an observation period and a particular times of day. The cameras 22 also enable observation of above surface features near each probe group 14. Above surface features such as rocks, foliage, and debris are observed such that their influence on temperature measurements taken by the probe groups 14 can be accounted for.

    (9) The probe groups 14 are distributed near the surrogate site's 10 soil surface to measure near-surface conditions so that near-surface conditions typical of known subsurface features at the surrogate site 10 can be observed. In the illustrated example, the known subsurface features of the surrogate site 10 include solid ground and an anomaly. Solid ground here means a vertical column of earth having a typical composition for the region and free of significant natural or manmade irregularities at least to a measurable depth. The measurable depth is a function of the temperature measurement equipment available, and in some examples is 15 meters, 20 meters, or greater than 20 meters. Exemplary subsurface anomalies include a portion of loosened or disturbed soil, or a void such as a cavern or manmade tunnel. An anomalous area 38 of soil surface at the surrogate site is above a subsurface anomaly, and a solid ground area 42 of soil surface at the surrogate site is above solid ground. The anomalous area 38 is bounded by anomaly borders 46 corresponding to an edge of the subsurface anomaly below ground.

    (10) It should be understood that the arrangement of measurement devices, visible phenomena, and subsurface features illustrated in FIG. 1 and described above is merely exemplary. In various other examples, a differing number and arrangement of measurement devices such as probe groups 14, cameras 22, or weather stations 18 of various types are used. In further examples, the measurement devices are used at a surrogate site 10 having different shading and different above surface and below surface features than illustrated in FIG. 1. In yet further examples, a surrogate site 10 has multiple subsurface anomalies, and the subsurface anomalies are of various shapes, sizes, and types.

    (11) FIG. 2 shows an arrangement of thermal probes 50 within a probe group 14. Within the group 14, a single probe 50 is positioned above a soil surface 54, possibly with appropriate support structures (not illustrated), and three probes 50 are positioned in soil 58 below the soil surface 54. All probes 50 within the group 14 are positioned close to the soil surface 54. For example, all probes 50 are positioned within one meter of the soil surface 54. The probes 50 beneath the soil surface 54 are positioned at equal or various depths. All probes 50 within the probe group 14 are arranged relative to one another as appropriate for the surrogate site. In various examples, each probe 50 in the group 14 is arranged to effectively measure temperature at a respective point 10 centimeters or less from the soil surface. Similarly, in various arrangements, each probe 50 located above the soil surface 54 is located to measure temperature of air at a point 10 centimeters or less above the soil surface 54, and each probe 50 located below the soil surface 54 is located to measure temperature of soil at a point 10 centimeters or less below the soil surface 54. The illustrated arrangement of probes 50 within the probe group is merely exemplary, and in other examples probe groups 14 include differing total numbers of probes 50, differing numbers of probes 50 above and below the soil surface 54, and arrangements of probes 50 within a group that differ from the arrangement illustrated in FIG. 2. In some examples, probe groups 14 include an equal number of probes above and below the soil surface 54. In various further examples, probe groups 14 in a single surrogate site each have a same number and arrangement of probes 50, or have differing numbers and arrangements of probes 50.

    (12) The combined use of probe groups 14, one or more cameras 22, and one or more weather stations 18 at a surrogate site 10 having known subsurface features as generally described above enables the influence of various factors on temperature conditions near the soil surface 54 at the surrogate site 10 to be given quantitative estimates. The temperature conditions tracked in various examples include, mean daily temperatures, temperatures at specific times of day, and rates or magnitudes of temperature change throughout a particular time of day or over a course of multiple days. Such temperature conditions measured by each probe group 14 are compared in view of the factors affecting each probe group to yield the quantitative estimates. For example, differences in magnitudes of temperature change measured by probe groups 14 over a certain period of time are compared between probe groups 14 in a solid ground area 42, in an anomalous area 38, and on or near an anomaly boundary 46. Such a comparison reveals how various subsurface features affect thermal inertia characteristics of a column of soil. The thermal inertia characteristics so derived are useable to identify subsurface features at other sites by observing temperature at the other sites over time and matching the observed temperature changes over time to the thermal inertia characteristics derived for known subsurface features. Such thermal inertia characteristics acquired from comparison of multiple successive surveys enable identification of below surface features of the survey site 70 and removal of noise generated by clutter at the survey site 70.

    (13) Further analysis of the results of the comparison in view of the cameras 22 and weather station 18 data reveals how certain above surface factors affect the measurable differences between portions of soil above various subsurface features. Such further analysis enables quantification of, for example, how shade or wind diminishes differences in thermal inertia characteristics between a solid ground area 42 and an anomaly area 38. The quantified effects of the above surface features are useable when interpreting survey data of sites with unknown subsurface features. Further analysis in some examples includes averaging the temperature measurements recorded within a 24 hour period from each probe group 14 to obtain a daily mean temperature at each probe group's 14 location. Differences between daily mean temperatures corresponding to each probe group 14 within a day and between different days is calculated, and those differences are attributed variously to subsurface features below each probe group 14, the amount of time each probe group 14 was shaded, weather conditions such as wind speeds at various times of day, above surface objects near each probe group 14, or interactions between any of the foregoing factors.

    (14) Such further analysis also reveals what transient conditions render a survey for unknown subsurface features unproductive. Specifically, the further analysis can identify weather conditions so unfavorable as to warrant rescheduling a live survey. For example, it may be revealed that certain subsurface features tend to be undetectable on days on which wind gusts above a certain speed are detected during a certain time. In a more specific example, it may be revealed that wind gusts above 15 miles per hour occurring between 9 AM and noon on a given day are likely to render subsurface voids 10 meters below the soil surface 54 or deeper undetectable by measurement of temperatures or thermal inertia at the soil surface 54. Wind gusts above 15 miles per hour occurring between 9 AM and noon would therefore be an unfavorable weather condition that warrants rescheduling a planned live survey if the planned live survey is intended to discover subsurface features than 10 meters below the soil surface 54 or deeper. Similar conclusions may be drawn based on other weather events such as daily high or low air temperatures, cloud cover, or precipitation. Generally, calmer and warmer weather conditions are required to discover deeper or smaller subsurface features. Thus, more cloud cover, higher windspeeds, cooler weather, and more precipitation are tolerable when only large or shallow subsurface features are sought. In some examples, a decision of whether to reschedule a live survey takes into account considerations of the depth or size of the subsurface features sought and of the probability of better weather conditions occurring later in the same year.

    (15) The foregoing measurements and analyses are used to compile and apply a database for use in connection with live surveys of a survey site 60, as illustrated in FIG. 3. The survey site 60 includes unknown subsurface features or anomalies 64, and a live survey of the survey site 60 includes use of a survey device 68 to collect data from the soil surface 54 of the survey site 60. The survey device 68 is any device equipped with sensors for collecting imageable temperature data from the soil surface 54, such as one or more infrared cameras. In some examples, the infrared camera is a cold filtered infrared camera. In some examples, the survey device 68 is further equipped with an optical camera or any other sensor for gathering data from the survey site 60. Example survey devices 68 include manned or unmanned aerial or terrestrial vehicles including automobiles, drones, and airplanes, handheld devices, or any moveable apparatus positionable to survey an area of land with any of the above described equipment. In further arrangements, the survey device 68 may be one or more thermal probes 50 deposited directly within the survey site 60.

    (16) In some arrangements, the live survey includes two or more surveys of the survey site 60 at different times. Measurements taken in each survey may be analyzed consecutively, and the difference between measurements taken in successive surveys may be used to determine the thermal inertia characteristics of specific points throughout the survey site 60. Those thermal inertia characteristics, in view of other observations of the survey site 60 such as the topography, shading, weather, foliage, or any other above ground features, may be compared to the database compiled from measurements previously taken at the surrogate site 10.

    (17) FIG. 4 represents a method 72 of imaging unknown subsurface features by application of the above described concepts. A monitoring step 76 includes monitoring a surrogate site 10 as generally set out in FIGS. 1 and 2 and the associated description. A database compiling step 80 includes compiling and analyzing the data from the monitoring step 76. The compiling step 80 involves analysis of how various above and below surface factors interact to affect temperatures measured at a soil surface 54 in the surrogate site.

    (18) Particularly, the compiling step 80 includes identifying “signatures” in the soil surface heat change over time data created by the subsurface features of the surrogate site 10. Such signatures are identified by correlating locations and characteristics, such as depth, volume, and type, of the known subsurface features of the surrogate site to differences in temperature change over time measured at particular locations at the soil surface of the surrogate site of varying proximity to the known subsurface feature. For example, if higher daytime peak temperatures and lower night time trough temperatures are measured at a particular soil location that is above a known subsurface anomaly void than at a particular soil location above a column of soil containing no subsurface anomalies, the differences in temperature over time characteristics are a signature created by the known subsurface void. In some arrangements, the characteristics of the known subsurface features are used to determine what signatures would be created by hypothetical subsurface features having slightly differing characteristics. For example, the signature that would be created by a 12 meter deep void may be estimated by extrapolating from known signatures created by a 10 meter deep void and a 15 meter deep void. In yet further arrangements, identification of the signatures includes accounting for the effects of above surface features of the surrogate site, weather conditions, or both, on soil surface temperatures measured at the surrogate site 10. For example, in some arrangements, temperatures measured at particular soil surface locations having varying proximity to an above surface feature are compared to estimate effects of the above surface feature on soil surface temperature at various distances from the above surface feature. In a more specific example, in some arrangements, comparison between temperatures measured at shaded and unshaded locations in the surrogate site 10 is used to estimate how much shade affects soil surface temperatures over time and at various specific times throughout a day. Such estimates of the effect of above surface features are considered in the identification of signatures of subsurface features. With the effects of above surface features accounted for in the identification of the signatures, unknown subsurface features can be detected and described more accurately from survey data acquired at survey sites 60 having different combinations and arrangements of above surface features than the surrogate site 10.

    (19) In some arrangements, the estimating or quantification of effects of above surface factors including above surface features such as shade, rocks, foliage, or reflective objects, and weather such as wind or precipitation is achieved by comparing mean measurements from two subsets of the data that show differing degrees of effect by any given above surface factor. For example, in some arrangements, temperature measurements from a first and second probe group 14 are averaged over a period of time to obtain a first mean and a second mean, respectively, wherein the first probe group is shaded, and the second probe group is unshaded, but the two probe groups are otherwise similarly situated. The period of time may be any suitable amount of time, and in specific examples, may be an hour, three hours, or 24 hours, such as from midnight to midnight. The difference between the first and second mean is obtained to provide a quantitative estimate of the effect of shade on soil surface temperature. In other examples, the first and second probe group, or a first and second sensor or sensor group of a different type, differ in proximity to other above surface features such as foliage, reflective surfaces, or local high or low elevation points of the soil surface. The difference between the first and second mean is accordingly a quantitative estimate of the effect of the proximity of the above surface feature whose proximity differentiates the two subsets of data. In another example, infrared radiation measurements of two locations at differing proximity to a reflective object may be averaged over time to provide the first and second mean, and the difference between the first and second mean provides a quantitative estimate of the effect of the proximity of the reflective object on infrared radiation. Similarly, in various additional arrangements, measurements from a single sensor or group of sensors, such as thermal probes or infrared cameras, are averaged on days having different weather conditions to obtain a quantitative estimate of the effect of the weather condition. For example, in some arrangements, measurements from a group of thermal probes are averaged over a period of time on a first day to obtain a first mean, and measurements from the same group of thermal probes may be averaged over the same period of time on a second day. In such arrangements, a difference between the first mean and the second mean provides a quantitative estimate of the effect of any weather condition, such as windspeed, timing of wind gusts, precipitation, that differs between the first day and the second day.

    (20) In yet further arrangements, the compiling step includes identification of weather conditions that can obscure signatures of subsurface features during times of day or times of year when those signatures would usually be detectable within soil surface temperature measurements over time. For example, it may be determined that wind gusts above a certain speed at certain times of day will generally obscure the signatures of subsurface anomalies having certain characteristics. In a more specific example, it may be determined that wind gusts above 15 miles per hour before noon on a given day during the monitoring period rendered the signatures of subsurface voids located 10 meters or further below the soil surface at the surrogate site 10 undetectable within the temperature data acquired by the thermal probes 50. The above surface features and weather conditions described in the foregoing examples are for illustrative purposes, and a skilled person would appreciate that the compiling step may be include similar consideration of any observable features and conditions at the selected surrogate site 10, and that various factors such as soil composition and local climate will make the interaction of above surface features, weather conditions, and subsurface feature characteristics on surface temperature data unique to any selected geographic region. Further features or conditions such as, but not limited to, foliage, foliage density, protruding rocks, topography, cloud cover, and precipitation may be considered in the compiling step.

    (21) A planning step 84 includes use of the compiled database to determine what times of year and what times of day are likely to have the best conditions for using soil surface 54 temperature data to image unknown subsurface features 64 at a survey site 60 in the same geographic region as the surrogate site 10. In some examples, the planning step 84 includes selection of multiple viable days for conducting a live survey, and in further examples the planning step 84 includes rescheduling a live survey from a preferred viable day to a backup viable day if weather conditions on the preferred viable day are unfavorable. A live survey of the survey site 60 is conducted in a surveying step 88 at a time or times selected in the planning step 84. Finally, results from the surveying step 88 are used to image unknown subsurface features 64 at the survey site 60 in a processing step 92.

    (22) In some examples, the processing step 92 includes known imaging methods, such as the DIRECT method described in U.S. Pat. No. 8,787,619, which is hereby incorporated by reference in its entirety. In some examples, an infrared survey in a dual time arrangement includes at least a first pass over the survey site with one or more infrared cameras and a second pass over the survey site with one or more infrared cameras, with the second pass conducted at a later time than the first pass. In various examples, the second pass is conducted with the same infrared wavelength or wavelengths. The processing step 92 further involves use of quantitative estimates of the effects of above surface factors derived in the compiling step 80 to declutter, or otherwise improve, imaging derived from infrared data acquired in the surveying step 88.

    (23) In some examples, the effects of shadows or other above ground factors visible in optical images of the survey site 60 are accounted for in generating images from the infrared data acquired in the surveying step 88. In further examples, the accounting includes adjusting for an estimated temperature difference or for reflected infrared radiation attributable to the above ground features within a single survey. In another example, the accounting involves using data acquired from multiple surveys to filter out temperature or infrared anomalies that have thermal inertia profiles inconsistent with subsurface features of interest. In a specific example, in a single set of infrared survey data acquired from an airborne survey device, an above ground feature such as a bush presents an infrared anomaly similar to that of a subsurface void. However, the bush or other above surface feature has differing thermal inertia characteristics than a subsurface void, and therefore exhibits a differing rate of temperature change over time. Successive surveys of the same location therefore reveal which temperature or infrared anomalies have thermal inertia characteristics more consistent with above surface features or clutter than with subsurface voids. Anomalies thereby identified as above surface features could be disregarded or factored into interpretation of other survey data. In various other examples, other above surface sources of clutter including, e.g., shadows, clouds, tracks, stains, disturbed terrain, holes, vegetation or foliage, foreign objects, foreign materials, foreign soils, water, cool air pools and roughness variations, are differentiated from subsurface anomalies by their thermal inertia characteristics as observed from at least two successive surveys.

    (24) This additional process of accounting for above ground factors based on interactions identified in the compiling step 80 makes the imaged contours of unknown subsurface features 64 more accurate and enables discovery and imaging of subsurface features 64 disposed further below the soil surface 54 than could be located using infrared imaging alone.

    (25) FIG. 5 shows a thermal map 110 generally consistent in appearance with a product of the above described infrared survey methods. Underneath thermal map 110 in FIG. 5 is an aerial photo 120 of the same area imaged in thermal map 110. Linear pipe voids and gaps that are not are visible in a circled area 130 of thermal map 110 that are not apparent within circled area 130 in aerial photo 120. Arrows in thermal map 110 and aerial photo 120 point to respective locations of bushes, deep crossed pipe voids, construction gaps, foot traffic, and metal tape markers on cinder blocks. However, aerial photo 120 includes thermal probe wires and a yardstick not visible in thermal map 110. Comparison of thermal map 110 to aerial photo 120 therefore suggests which features may be included or excluded from an image by selecting an image type.

    (26) The above methods are described generally with reference to location of below ground anomalies, but other applications are available. For example, similar results may be obtained by application to manmade objects such as bridges, decks, building foundations, or anything else constructed from concrete.

    (27) Although the invention herein has been described with reference to particular embodiments, it is to be understood that these embodiments are merely illustrative of the principles and applications of the present invention. It is therefore to be understood that numerous modifications may be made to the illustrative embodiments and that other arrangements may be devised without departing from the spirit and scope of the present invention as defined by the appended claims.