US 20020035553 A1 Abstract The present invention relates generally to interpolation methods for automatic generation of an accurate digital elevation model, and more particularly relates to intelligent interpolation methods for accurate extraction of 3-dimensional digital elevation models from satellite images, aerial photographs, or land surveying. The present invention is composed of three basic processes: The first process (
102) of COG/ECI elimination, the second process (104) of hole-fill segmentation, and the third process (103) of noise-remove segmentation. The invention produces an accurate digital elevation model even for the areas with open boundaries such as the sea off the coast and over a river. Claims(6) 1. An interpolation method for automatically generating a digital elevation model, which removes incorrect elevation outside stereo-matched regions, comprising the steps of;
the first process ( 102) where the interpolated elevation from the result (101) of a conventional interpolation method is eliminated if the corresponding Center-of-Gravity (COG) or Empty-Center-Index (ECI) is greater than a threshold; the second process (hole-fill segmentation, 104) where the eliminated region after the first process (102) is segmented, and it is decided whether the eliminated elevation should be reconstructed based upon the size of the segment; and the third process (noise-remove segmentation, 103) where the interpolated region after the first process (102) is selected, and it is decided whether the interpolated elevation should be eliminated based upon the size of the merged segment. 2. In the method as defined in 102) wherein x_{o}, y_{o }are horizontal coordinates, max-dist is the maximum of the distance between (x_{o}, y_{o}) and one of the interpolants, r_{k }is the distance between (x_{o}, y_{o}) and an interpolant, unit grid size Δr_{unit }is a constant to prevent w_{eq }reaching singularity, w_{eqECI }is a weighting factor to equalize the distribution of interpolants with respect to r_{k}, and
(EQUATION 5)
3. In the method as defined in 102) wherein w_{eqCOG}(x_{k}−x_{o}) is a weighting factor to equalize the distribution of interpolants with respect to (x_{k}−x_{o}), δ_{anti-sigu }is a constant to prevent w_{eqCOG }reaching singularity, and
(EQUATION 6)
similarly for {overscore (Δy)}.
4. The method as defined in 104) comprises of the following three steps:
the (a) step (
418) which searches neighboring grid cells of no valid elevation when the current grid cell does not have valid elevation, assigns a segment identification number (ID) to the present cell if the neighboring empty cells have any segment ID, assign a new ID if no neighbor cells have an ID, repeats the above tasks from the top-left to the bottom-right, and completes one scan for the image segmentation; the (b) step (
419) which merges the fragmented segments by repeating the (a) step (418) starting from top-right, bottom-left, and bottom-right for each scan, and assigns a neighboring ID to the current cell when the current cell's ID is greater than the neighboring IDs; and the (c) step (
420) which recovers the removed elevation of the merged segment if the size of the segment is smaller than a threshold, and leaves it as a hole when the segment is larger than the threshold. 5. The method as defined in 103) comprises of the following three steps:
the (A) step (
517) which search neighboring grid cells when the current grid cell has valid elevation, assigns a segment identification number (ID) to the present cell if the neighboring cells have any segment ID, assigns a new ID if no neighbor cells have an ID, repeats the above tasks from the top-left to the bottom-right and completes one scan for the image segmentation; the (B) step (
511) which merges the fragmented segments by repeating the (A) step (517) starting from top-right, bottom-left, and bottom-right for each scan, and assign a neighboring ID to the current cell when the current cell's ID is greater than the neighboring IDs; and the (C) step (
512) which removes the interpolated elevation of the merged segment if the size of the segment is smaller than a threshold. 6. The method as defined in 517).Description [0001] 1. Field of the Invention [0002] The present invention relates generally to interpolation methods for automatic generation of an accurate digital elevation model, and more particularly relates to intelligent interpolation methods for accurate extraction of 3-dimentional digital elevation models from satellite images, aerial photographs, or ground surveying. [0003] 2. Description of the Related Art [0004] Digital elevation model (DEM) is the topographic data composed of the set of horizontal coordinates (x, y) and its elevation (z) for each grid cell in a given region. The process of DEM can be summarized in three steps; (1) image-matching step to find overlapped regions from stereo images, (2) camera-modeling step to obtain (x,y,z) from the results of the image-match with a camera model relating the photo plane with the absolute ground coordinates, and (3) interpolation step to change the uneven horizontal distribution resulted from camera-model into the uniform horizontal square distribution. In ground surveying, image-match and camera-model could be omitted. [0005] Interpolation is essential for an automatically generated DEM. For a manually generated DEM, interpolation is necessary, but the effect of interpolation is insignificant. Many studies on interpolation have progressed since 1970's, and focused on finding the optimal interpolation method from dozens of mathematical model such as Gaussian, Kriging, Nearest Neighbor, Moving Window Average, Multiquadric, Modified Shepard, Spline, and Minimum Curvature. We call these old methods as conventional interpolation methods (CIMs). [0006] Since the CIMs handle less than a few hundreds of input data as an experimental trial, following problems show up when the CIMs are applied to the results of the real data; artificial lands are created with the size of the interpolation radius near coast, breakwater, lake, and river. An intellectual interpolation method is introduced to deal with the real data (Krupnik, A., 1998, Automatic detection of erroneous areas in automatic surface reconstruction, Proceedings of ISPRS Commission III Symposium, Columbus, Ohio, p. 132-137). In this method, the edges of a lake or farmland are matched, and the interpolation is performed on the matching results. However, this method can not be applied to the areas with open boundaries such as the sea off the coast and a river. [0007] It is an object of the present invention to provide intelligent interpolation methods for automatic generation of an accurate digital elevation model, which can solve the above-mentioned problems near a coast or a river. [0008] The present invention is composed of three basic processes: The first process ( [0009] The invention produces an accurate digital elevation model even for the areas with open boundaries off the coast and over a river. [0010] Exemplary embodiments of the present invention will be described in conjunction with the drawings in which: [0011]FIG. 1 is an overall block diagram of the intelligent interpolation methods of the present invention; [0012]FIG. 2 represents definitions of Center-of-Gravity (COG) and Empty-Center-Index (ECI) of the present invention; [0013]FIG. 3 shows the number of interpolants as a function of the distance from the center of the interpolation and the figure is referenced in Center-of-Gravity calculation of the present invention; [0014]FIG. 4 shows the number of interpolants as a function of the distance from the center of the interpolation and the figure is referenced in Empty-Center-Index calculation of the present invention; [0015]FIG. 5 exemplifies fragmentation during segmentation of the present invention; [0016]FIG. 6 is a flow chart for the hole-fill segmentation scheme of the present invention; and [0017]FIG. 7 is a flow chart for the noise-remove segmentation scheme of the present invention. [0018] [0019] [0020] [0021] [0022] [0023] [0024] [0025] [0026] [0027] The present invention will be better understood with regard to the following description, appended claims, and accompanying figures. The present invention is composed of three basic processes: In the first process ( [0028] As mentioned earlier, when the CIMs are used, artificial lands are created with the size of the interpolation radius near the edge of the matched region. In order to remove these unwanted artifacts, the present invention utilizes COG/ECI elimination scheme. Since artificial lands show up due to the concentration of the interpolants near the edge of the interpolation disc, the artificial elevation is removed using COG criterion; the interpolated elevation is eliminated if the corresponding COG is greater than a threshold. However, when the interpolants are distributed evenly in all directions as shown in the right of FIG. 2, the corresponding COG is close to zero, and the COG criterion can not eliminate the false elevation even though the interpolants concentrate near the edge of the interpolation region. In this case, noting the emptiness of the center of the interpolation disc, we use ECI elimination. In ECI elimination, the interpolated elevation is removed if the corresponding ECI is greater than an effective threshold. [0029] The COG in FIG. 2 is defined as: [0030] (EQUATION 1)
[0031] similarly for {overscore (Δy)}. [0032] Here x [0033] This weighting factor in FIG. 3 is defined as: [0034] (EQUATION 2)
[0035] Here δ [0036]FIG. 3 shows the number of interpolants as a function of the distance from the center of the interpolation when the interpolation disc is filled with interpolants; left shows the number of interpolants before the equalization, and right shows after. [0037] The ECI in FIG. 2 and FIG. 4 is defined as: [0038] (EQUATION 3)
[0039] Here r [0040] This weighting factor in FIG. 4 is defined as: [0041] (EQUATION 4)
[0042] Here unit grid size Δr [0043]FIG. 4 shows the number of interpolants as a function of the distance from the center of the interpolation where the interpolants are located in the hatched area; center column shows the number of interpolants before the equalization, and rigth column shows after. [0044] After the first process ( [0045] In the second process ( [0046] In the (a) step ( [0047] In the (b) step ( [0048] In the (c) step ( [0049] In the (a) step ( [0050] Segmentation mentioned in the (b) step ( [0051] In the third process ( [0052] In the (A) step ( [0053] In the (B) step ( [0054] In the above mentioned processes ( [0055] As mentioned above, the present invention can handle the elevation near a breakwater, a coast, and a river with open edges, and the invention provides the intelligent interpolation methods for automatic generation of an accurate digital elevation model. Therefore, the invention makes the DEM more useful in the field of disaster prevention and coastal simulation. [0056] When the invention is introduced to a commercial software to produce a DEM off the coast, it can provide better results. The intellectual interpolation methods according to the invention require less than 5% of the calculating time compared to the CIMs and 100% additional computer memories. The cost of the memory is not so high, but the accuracy of interpolation is improved very much. These intelligent interpolation methods could be used with the CIMs, and applied to stereo-images with different horizontal resolutions, visible range images, and radar images. Referenced by
Classifications
Legal Events
Rotate |