Geospatial Modeling System Providing Building Roof Type Identification Features and Related Methods

Information

  • Patent Application
  • 20080021683
  • Publication Number
    20080021683
  • Date Filed
    July 20, 2006
    18 years ago
  • Date Published
    January 24, 2008
    16 years ago
Abstract
A geospatial modeling system may include a geospatial model database and a processor. The processor may cooperate with the geospatial database for identifying a building roof type defined by building roof data points as being from among a plurality of possible building roof types. This may be done based upon applying multi-directional gradient calculations to the building roof data points.
Description

BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is a schematic block diagram of a geospatial modeling system in accordance with the invention.



FIG. 2 is a flow diagram illustrating a geospatial modeling method in accordance with the invention.



FIG. 3 is a series of top views of flat building roof data points illustrating the flat building roof identification method steps of FIG. 2.



FIG. 4 is a series of top views of sloped building roof data points illustrating the sloped building roof identification method steps of FIG. 2.



FIGS. 5-8 are schematic block diagrams illustrating complex sloped roof rendering in accordance with the invention.



FIG. 9 is top view of domed roof building data points and corresponding histogram representations therefor illustrating the doomed roof identification method steps of FIG. 2.



FIG. 10 is a perspective view of buildings with flat and sloped roofs rendered in accordance with the present invention



FIGS. 11 and 12 are perspective views of buildings with complex sloped roofs rendered in accordance with the present invention.



FIG. 13 is a perspective view of a domed roof building rendered in accordance with the present invention.





DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

The present invention will now be described more fully hereinafter with reference to the accompanying drawings, in which preferred embodiments of the invention are shown. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and completer and will fully convey the scope of the invention to those skilled in the art. Like numbers refer to like elements throughout.


Referring initially to FIG. 1, a geospatial modeling system 20 illustratively includes a geospatial model database 21 and a processor 22 that may advantageously be used for identifying and rendering different types of building roofs. By way of example, the processor 22 may be a central processing unit (CPU) of a PC, Mac, or other computing workstation, for example. A display 23 may also be coupled to the processor 22 for displaying geospatial modeling data, as will be discussed further below. The processor 22 may be implemented using a combination of hardware and software components to perform the various operations that will be discussed further below, as will be appreciated by those skilled in the art.


By way of example, the geospatial data may be captured using various techniques such as stereo optical imagery, Light Detecting and Ranging (LIDAR), Interferometric Synthetic Aperture Radar (IFSAR), etc. Generally speaking, the data will be captured from overhead (e.g., nadir) views of the geographical area of interest by airplanes, satellites, etc., as will be appreciated by those skilled in the art. However, oblique images of a geographical area of interest may also be used in addition to (or instead of) the nadir images to add additional 3D detail to a geospatial model. The raw image data captured using LIDAR, etc., may be processed upstream from the geospatial model database 21 into a desired format, such as a digital elevation model (DEM), or this may be done by the processor 22.


Turning now to FIGS. 2 through 9, beginning at Block 30, the processor 22 may remove ground and vegetation data points from geospatial model data to select building roof data points for processing, at Block 31. By way of example, this may be done using the process set forth in the above-noted U.S. Pat. No. 6,654,690 to Rahmes et al., as well as with RealSite® or other suitable approaches. Once the building roof data points are isolated, the processor 22 may begin identifying the building roof type defined by the building roof data points as being from among a plurality of possible building roof types. By way of example, the building roof types may include flat roofs, sloped roofs, complex slopeds roofs, and domed roofs (or combinations thereof). Generally speaking, this identification is done based upon applying multi-directional gradient calculations to the building roof data points.


More particularly, in the example shown in FIG. 3, building roof data points 50 are from a building with a flat roof. The multi-direction gradient calculations may include Robert's cross calculations and, more specifically, calculations based upon an eight-way (i.e., 3×3) Robert's cross grid 51, for example. This is done by centering the Robert's cross grid 51 for each building data point 50, summing distnaces of points on opposite sides in the grid, and identifying building roof data points whose sum is within a threshold from zero as being localized peak data points, at Block 32. That is, any building roof data point 50 having a sum of the differences of the points on opposite sides of the grid that is substantially zero will have a substantially zero slope, as will be appreciated by those skilled in the art. This indicates that the given point is not on an inclined surface, i.e., it is either located at a localized peak or on a flat surface (both of which are referred to as localized peak points herein).


To determine whether the identified points are part of a flat roof, the processor 22 determines a standard deviation for the localized peak data points over a selected building roof area 52, at Block 33. The building roof area 52 may be selected based upon an aggregation of localized peak data points Selection is automatic based on a low (near zero slope) value provided by Robert's Cross Operator. If the selected building roof area 52 is greater than a threshold portion of a total building roof area, and if the standard deviation is below a threshold standard deviation, then the processor 22 identifies the building roof type as a flat roof type, at Blocks 34-35. By way of example, the threshold portion of the total building roof area may be greater than about 65%), and the threshold standard deviation may be less than about one meter. Of course, different thresholds may be used in different embodiments. The processor 22 may accordingly substitute one or more geometric objects/shapes in place of the building data points that approximates the flat roof building, such as the building 100 in FIG. 10, and optionally display the shape on the display 23 with or without accompanying terrain data (i.e., ground, foliage, etch).


If either of the above conditions is not met, however, then the processor 22 will proceed to determine whether the building has one or more sloped roof portions. Referring more particularly to FIG. 4, the processor 22 removes groups of contiguous localized peak data points being less than a threshold in number, at Block 36. For example, a group of less than three contiguous points may be excluded from consideration in determining whether the building roof data points correspond to a sloped roof. In the illustrated example, two sets of building roof data points 60a, 60b for buildings each having a standard sloped roof are shown. After performing calculations based upon the 3×3 Robert's cross grid 61 as described above, three groups 62a, 63a, 64a and 62b, 63b, 64b of localized peak data points remain for the buildings 60a, 60b, respectively.


If a given length of a contiguous localized peak data point group remains, at Block 37, and this length has an area that is significantly less than the entire area of the entire roof but is significantly greater than the lengths of any other remaining groups, then the processor 22 identifies the roof as being of a single sloped roof type (Block 39) if no more than one such group is present, at Block 38. In the illustrated example, the groups 63a and 63b meet this criteria and therefore represent the apex of the sloped roof structure. However, if there are more than four corners, this identifies the roof as a complex sloped roof (Block 40). That is, a complex sloped roof has multiple sloped roof features with multiple apexes.


To construct a sloped or complex sloped roof, the processor 22 generates sloped roof/building shapes 65 (FIG. 5). For a single sloped roof, the processor 22 fits the sloped roof shape 65 to the width, length, and height of the building corresponding to the building data points. Exemplary single sloped roof buildings 101, 102 rendered using the sloped roof shape 65 are shown in FIG. 10. For rendering a complex sloped roof which includes more than one sloped roof shape 65, the complex sloped roof can be segmented into a plurality of intersecting sloped roof shapes, as shown in FIG. 5, which are merged to provide a single complex sloped roof/building structure 66.


More particularly, the process of rendering different sloped roof shapes is further illustrated in FIGS. 6 through 8. The processor 22 first verifies that the polygon that approximates the building shape is segmentable. That is, the processor 22 checks to make sure that all of the angles of the building are 90° or 270° (i.e., negative 90 degrees) angles (FIG. 6(b)). Then, the shape is graphically skeletonized (FIG. 6(c)) and reduced to nodes 67 and edges by removing points with exactly two neighbors, and then merging any adjacent points (FIG. 6(d)).


The edge points 67 are then classified based upon the number of their respective neighbors. That is, an edge point having one neighboring edge point is classified as an end point 67a. An edge point with three neighboring edge points, two of which are end points 67a, is classified as a sloped end point 67b. Moreover, an edge point with three neighboring edge points, one of which is an end point 67a, is classified as a corner point 67c. Further, an edge point with three neighboring edge points and no adjacent end points 67a is classified as a T-junction point 67d. Additionally, an edge point having four neighboring edge points is classified as a cross-junction point (i.e., it is located at the center of a “+” shaped complex sloped building roof structure), which is not shown in the example illustrated in FIG. 6.


For each edge without an end point, a rectangle is generated as shown in FIG. 7 based upon the given geometric building shape (e.g., “I” shaped, “L” shaped, “T” shaped, “U” shaped, “S” shaped, “J” shaped, or “H” shaped). For each rectangle, the slopeds are then projected. For each sloped point not between two endpoints, this point is replaced with the nearest intersection between its sloped segment and any other sloped segment, as illustrated in FIG. 8, to render the final complex sloped roof shape Exemplary complex sloped roof buildings 121 through 125 rendered using the above-described approach are shown in FIGS. 11 and 12.


If it is determined that the roof type is neither flat nor a single/complex sloped roof, the processor 22 then determines whether the building roof data points define a circular shape, at Block 41, as will be appreciated by those skilled in the art. Referring to FIG. 9(a), a set of circular building roof data points 90 is shown As used herein, “circular” is meant to include not only perfect circles but also variations thereof, including ovals, ellipses, semicircular shapes, etc. If the building roof data points 90 do define a circular shape, then the processor 22 generates an elevation histogram for the heights of the points, at Block 42 (FIG. 9(b)). Any statistical outlying points are then removed (FIG. 9(c)), and if the standard deviation of the remaining histogram points is below a threshold (e.g., it is low relative to its mean), then the roof is identified as a domed roof, and the processor 22 substitutes an appropriate dome shape therefor. An exemplary domed roof building 130 rendered using this approach is shown in FIG, 13.


Once the appropriate building roof type is identified and rendered (or if the roof type does not fall into one of the above-described categories), the illustrated method is concluded (Block 45). Of course, in an actual implementation the process may then be repeated to identify other building roof types for different sets of building rood data points.


It should be noted that while the above-described exemplary approach discusses an order of determining roof shapes beginning with flat roof shapes, then sloped roof shapes, and lastly doomed roof shapes, this order may be changed in different embodiments. Generally speaking, however, if a large urban area is being classified then it is desirable to start with the most common building roof type, as this may conserve processing resources. For example, in a downtown area, many of the buildings may have flat roofs, so it may be advantageous to first check if a given set of building roof data points corresponds to a flat roof. In a suburban area, many of the buildings may be houses, which are more likely to have sloped roof structures, and therefore in such applications it may be desirable to start with sloped roof identification.


Moreover, the processor 22 need not check for every possible roof type in every embodiment. For example, it may be desirable to identify only buildings that have a desired type of roof in some examples, for example. Also, the actual shapes of the buildings need not necessarily be rendered for display That is, the above-described approach may be used to automatically identify building roof types with or without subsequent automatic roof/building rendering.


Many modifications and other embodiments of the invention will come to the mind of one skilled in the art having the benefit of the teachings presented in the foregoing descriptions and the associated drawings Therefore, it is understood that the invention is not to be limited to the specific embodiments disclosed, and that modifications and embodiments are intended to be included within the scope of the appended claims.

Claims
  • 1. A geospatial modeling system comprising: a geospatial model database; anda processor cooperating with said geospatial model database for identifying a building roof type defined by building roof data points as being from among a plurality of possible building roof types based upon applying multi-directional gradient calculations to the building roof data points.
  • 2. The geospatial modeling system of claim 1 wherein said processor identifies the building roof type as being from among a flat roof type, a sloped roof type, a complex sloped roof type, and a domed roof type.
  • 3. The geospatial modeling system of claim 1 wherein the multi-direction gradient calculations comprise Robert's cross calculations.
  • 4. The geospatial modeling system of claim 3 wherein the Robert's cross calculations comprise eight-way Robert's cross calculations.
  • 5. The geospatial modeling system of claim 1 wherein said processor further removes ground and vegetation data points from geospatial model data to select the building roof data points.
  • 6. The geospatial modeling system of claim 1 wherein the multi-direction gradient calculations comprise using a grid centered for each building data point, summing distances of opposite sides in the grid, and identifying building roof data points whose sum is within a threshold from zero as being localized peak data points
  • 7. The geospatial modeling system of claim 6 wherein the multi-direction gradient calculations further comprise determining a standard deviation for localized peak data points over a selected building roof area, and if the selected building roof area is greater than a threshold portion of a total building roof area, and if the standard deviation is below a threshold standard deviation, then identifying the building roof type as a flat roof type.
  • 8. The geospatial modeling system of claim 6 wherein the multi-direction gradient calculations further comprise removing groups of contiguous localized peak data points being less than a threshold in number.
  • 9. The geospatial modeling system of claim 8 wherein the multi-direction gradient calculations further comprise determining if a length of contiguous localized peak data points remains, and, if so, then identifying the building roof type as a sloped roof type.
  • 10. The geospatial modeling system of claim 9 wherein the multi-direction gradient calculations further comprise determining a plurality of intersecting lengths of contiguous localized peak data points, and, if so, identifying the building roof type as a complex sloped roof type.
  • 11. The geospatial modeling system of claim 1 wherein said processor determines whether the building roof data points define a circular shape; and wherein the multi-direction gradient calculations comprise calculating at least one elevation histogram for the roof building data points.
  • 12. The geospatial modeling system of claim 11 wherein the multi-direction gradient calculations further comprise calculating elevation histogram data for the roof building data points, and, if a standard deviation of the elevation histogram data is below a threshold, then indicating the building roof type as a domed roof type
  • 13. The geospatial modeling system of claim 1 wherein said geospatial model database stores at least one respective building roof shape for each of the plurality of possible building roof types; and wherein said processor cooperates with said geospatial model database to substitute in place of the building roof data points the at least one respective building roof shape for the identified building roof type.
  • 14. A geospatial modeling system comprising: a geospatial model database for storing at least one respective building roof shape for each of a plurality of possible building roof types; anda processor cooperating with said geospatial model database for identifying a building roof type defined by building roof data points as being from among the plurality of possible building roof types based upon applying multi-directional gradient Robert's cross calculations to the building roof data points, andsubstituting in place of the building roof data points the at least one respective building roof shape for the identified building roof type.
  • 15. The geospatial modeling system of claim 14 wherein said processor identifies the building roof type as being from among a flat roof type, a sloped roof type, a complex sloped roof type, and a domed roof type.
  • 16. The geospatial modeling system of claim 14 wherein said processor further removes ground and vegetation data points from geospatial model data to select the building roof data points.
  • 17. The geospatial modeling system of claim 14 wherein the multi-direction gradient calculations comprise using a grid centered for each building data point, summing differences of opposite sides in the grid, and identifying building roof data points whose sum is within a threshold from zero as being localized peak data points.
  • 18. A geospatial modeling method comprising: storing building roof data points in a geospatial model database; andidentifying a building roof type defined by building roof data points as being from among a plurality of possible building roof types using a processor based upon applying multi-directional gradient calculations to the building roof data points.
  • 19. The method of claim 18 wherein identifying comprises identifying the building roof type as being from among a flat roof type, a sloped roof type, a complex sloped roof type, and a domed roof type.
  • 20. The method of claim 18 wherein the multi-direction gradient calculations comprise Robert's cross calculations.
  • 21. The method of claim 18 wherein the multi-direction gradient calculations comprise using a grid centered for each building data point, summing distances of opposite sides in the grid, and identifying building roof data points whose sum is within a threshold from zero as being localized peak data points.
  • 22. The method of claim 21 wherein the multi-direction gradient calculations further comprise determining a standard deviation for localized peak data points over a selected building roof area, and if the selected building roof area is greater than a threshold portion of a total building roof area, and if the standard deviation is below a threshold standard deviation, then indicating the building roof type as a flat roof type.
  • 23. The method of claim 21 wherein the multi-direction gradient calculations further comprise removing groups of contiguous localized peak data points being less than a threshold in number.
  • 24. The method of claim 23 wherein the multi-direction gradient calculations further comprise determining if a length of contiguous localized peak data points remains, and, if so, then identifying the building roof type as a sloped roof type.
  • 25. The method of claim 24 wherein the multi-direction gradient calculations further comprise determining a plurality of intersecting lengths of contiguous localized peak data points, and, if so, identifying the building roof type as a complex sloped roof type.
  • 26. The method of claim 18 further comprising determining whether the building roof data points define a circular shape using the processor; and wherein the multi-direction gradient calculations comprise calculating elevation histogram data for the roof building data points, and, if a standard deviation of the elevation histogram data is below a threshold, then identifying the building roof type as a domed roof type.
  • 27. The method of claim 18 wherein the geospatial model database stores at least one respective building roof shape for each of the plurality of possible building roof types; and further comprising using the processor to substitute in place of the building roof data points the at least one respective building roof shape for the identified building roof type.