JP2005078349A - Geographic information measuring system, method and program - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To provide a geographic information measuring system that allows easily obtaining necessary three-dimensional information without requiring a huge amount of work in the construction of the system. <P>SOLUTION: A two-dimensional map is stored in the storage part of a geographic data providing server as it is associated with latitudes and longitudes. Aerial images (front view image, image immediately below, and rear view image) are stored in the storage part of a computer terminal as they are associated with the latitudes, longitudes and altitudes of their respective central projection points. A user specifies a certain position on the two-dimensional map displayed in the window W1 of a display. The front view image, the image immediately below, and the rear view image of the position specified in the window W2 are displayed. The user designates one point (one pixel) on the image immediately below as a measuring point. The computer terminal detects the corresponding points to the measuring points on the front view image and the rear view image using area correlation method, thereby calculating the three-dimensional information of the measuring points through stereo matching. <P>COPYRIGHT: (C)2005,JPO&NCIPI

Description

  The present invention relates to a GIS (Geographic Information System).

The three-dimensional GIS disclosed in Non-Patent Document 1 is created based on pre-measured three-dimensional information.
Seiichi Hirata, “3D GIS”, Surveying 2003.7. Electronic Land, Electronic Land Editing Subcommittee-Surveying Editing Committee, p. 13-20 "Business Mapping-updated on 2003/8/4", Top page-MSR data product, MS-Map option data, [online], Mitsui Engineering & Shipbuilding System Engineering Co., Ltd. e-business division, [search August 13, 2003 ], <Http://www.msr.mes.co.jp/mapinfo/>

  The three-dimensional GIS based on pre-measured three-dimensional information has a wide geographical range targeted by the GIS, and therefore requires a huge amount of work to obtain highly accurate three-dimensional information in advance over the entire target geographical range. There was a problem.

  The present invention has been made in view of the above problems, and provides a geographic information measurement system and the like that makes it possible to easily obtain necessary three-dimensional information without requiring a huge amount of work for system construction. The purpose is to do.

  The geographic information measuring system of the present invention is a two-dimensional map display means for displaying a two-dimensional map on a display, and a two-dimensional image obtained by imaging each location on the two-dimensional map and stored in advance in a storage means. A two-dimensional image display means for picking up an image of a designated part on the two-dimensional map and displaying it on the display; and a point, line or region designated on the displayed two-dimensional image 3D to calculate 3D coordinates by stereo-matching the 2D image and another 2D image obtained by imaging the point, line, or region from an angle different from that of the 2D image And an information measuring means.

  In the geographic information measuring method of the present invention, a first step of causing a computer to display a two-dimensional map on a display, and a certain location on the two-dimensional map are specified, thereby allowing the computer to display the two-dimensional map on the two-dimensional map. A second step in which the designated part is picked up from the two-dimensional images stored in advance in the storage means obtained by imaging each part, and displayed on the display; By specifying a point, line, or region on the image, another two-dimensional image obtained by imaging the two-dimensional image and the point, line, or region from a different angle from the two-dimensional image. And a third step of calculating a three-dimensional coordinate of the point, line, or region by stereo-matching with an image.

  The geographic information measuring program of the present invention is obtained by imaging each location on the two-dimensional map with a computer, and a specific location on the two-dimensional map is selected from the two-dimensional images stored in the storage means in advance. The three-dimensional coordinates of a point, line, or region designated on the two-dimensional image selected as captured and displayed on the display are the two-dimensional image, and the point, line, or region is the two-dimensional image. The step of calculating by performing stereo matching with another two-dimensional image obtained by imaging from different angles is performed.

  In the present invention, since it is not necessary to measure highly accurate three-dimensional information in advance, an enormous amount of work is not required for system construction. The user of the present invention can easily display a stereo image (two-dimensional image) of a survey target site with reference to a two-dimensional map, and further specify a survey point, survey line, or survey area on the stereo image. Can perform highly accurate three-dimensional surveying.

  In the geographic information measuring program of the present invention, in the step, a point in a single two-dimensional image is designated, and a point on the other two-dimensional image corresponding to the designated point is detected by an area correlation method. Is preferable.

  The geographic information measurement program of the present invention is used to acquire a cross-sectional shape of a measurement object existing on the two-dimensional image. In the step, a line passing through the cross-section in a single two-dimensional image is used. The points on the other two-dimensional image corresponding to the plurality of feature points on the designated line are detected by the area correlation method, and the three-dimensional coordinates of the feature points are calculated. It is preferable to execute a fourth step of interpolating the three-dimensional coordinates of each of the feature points to generate an image showing the cross-sectional shape.

  In the geographic information measuring program of the present invention, in the step, a region in a single two-dimensional image is designated, and points on the other two-dimensional image corresponding to a plurality of feature points in the region are area correlated. In addition, the three-dimensional coordinates of each feature point are calculated, and the computer interpolates the three-dimensional coordinates of each feature point to generate an image showing the shape of the region. Is preferably executed.

  According to these, it becomes unnecessary for the user to specify corresponding points on the stereo image.

  Another aspect of the geographic information measurement program of the present invention is obtained by imaging each location on a two-dimensional map with a computer, and from the two-dimensional image stored in advance in the storage means, on the two-dimensional map. A first step of selecting an image of a point including a designated point, line or region, and detecting a corresponding point of the point, line or region in the selected two-dimensional image; Performing a second step of calculating coordinates by stereo-matching the two-dimensional image and another two-dimensional image obtained by imaging the corresponding portion from a different angle from the two-dimensional image It is characterized by making it.

  In this aspect of the present invention, the user can obtain the above-described operation of the present invention without referring to a stereo image.

  Another aspect of the geographic information measurement system of the present invention is that the three-dimensional coordinates of a measurement target specified on a two-dimensional image obtained by imaging a specific location on a two-dimensional map are represented by the two-dimensional image and the measurement target. 3D information measuring means for calculating by stereo-matching with another 2D image obtained by imaging an object from an angle different from that of the 2D image, and 3D calculated by the 3D information measuring means 3D information display means for displaying a coordinate or a three-dimensional shape of the measurement object obtained based on the three-dimensional coordinate on a display so as to be superimposed on the two-dimensional map.

  In this aspect of the present invention, since it is not necessary to measure highly accurate three-dimensional information in advance, an enormous amount of work is not required for system construction. The user of the present invention can perform highly accurate three-dimensional survey only by specifying a survey point, survey line, or survey area on a stereo image, and simultaneously refer to the result of the three-dimensional survey and a two-dimensional map. can do.

  According to the present invention, a precise survey can be performed by operating a computer without going to a survey target site. Moreover, the result of precise surveying can be reflected in GIS.

  Hereinafter, preferred embodiments of the present invention will be described in detail with reference to the accompanying drawings.

  FIG. 1 is a functional configuration diagram of a GIS to which the geographic information measurement program 3 of the present embodiment is applied. The GIS of the present embodiment includes a computer terminal 1 that executes a command of the geographic information measurement program 3, and a geographic data providing server 4 that is connected to the terminal 1 via a communication network.

  The storage unit of the geographic data providing server 4 stores a spatial index file f1 and a geographic data file f2. The geographic data file f2 stores a two-dimensional map of a certain area defined by latitude and longitude. The spatial index file f1 stores the latitude and longitude of the center of each section of the two-dimensional map, and the altitude of the ground surface at discrete points.

  The storage unit of the terminal 1 stores a TLS original image / GPS data file f3. The TLS original image / GPS data file f3 stores a captured image of the ground surface obtained by scanning the above-described area with a TLS (Three Line Scanner) mounted on an aircraft. The TLS is one in which three CCDs are mounted on the aircraft so as to be aligned on the flight direction line, and the three CCDs respectively connect the front-view image Pf, the direct-view image Pn, and the rear-view image Pb on the ground surface. Take an image.

  The TLS original image / GPS data file f3 stores the three-dimensional coordinates (latitude, longitude, altitude) of the projection center points Cf, Cn, Cb of the captured image at each imaging time point. Specifically, the projection center point Cf indicates the three-dimensional coordinates of the center point of the imaging surface of the forward-view image capturing CCD at the time when a certain forward-view image Pf is captured. Similarly, the projection center point Cn corresponds to the nadir view image Pn, and the projection center point Cb corresponds to the rear view image Pb. As a method for calculating the three-dimensional coordinates of the projection center points Cf, Cn, Cb, GPS (Global Positioning System) or the like can be used.

  When the geographic information measurement program 3 stored in the CD-ROM 2 is installed in the terminal 1, a two-dimensional map presentation module 31, a TLS image presentation module 33, a three-dimensional information calculation module 35, a cross-sectional shape / three-dimensional calculation module 37, The three-dimensional information presentation module 38 and the corresponding location detection module 39 respectively convert the terminal 1 into a two-dimensional map presentation unit, a TLS image presentation unit, a three-dimensional information calculation unit, a cross-sectional shape / three-dimensional calculation unit, a three-dimensional information presentation unit, and It functions as a corresponding location detection means.

  When a request for presenting a specific section of the two-dimensional map is input to the input device of the terminal 1, the two-dimensional map presenting means acquires image data of the section from the geographic data providing server 4, and the section is displayed on the display 12. Present the displayed window.

  The TLS image presenting means calls a TLS captured image at a location designated by the user in the map section displayed on the display 12 from the TLS original image / GPS data file f3, and presents a window on which the TLS captured image is displayed on the display 12 To do.

  The three-dimensional information calculation means calculates the three-dimensional coordinates of the points, lines, or regions on the TLS captured image by stereo matching the TLS captured image.

  The cross-sectional shape / three-dimensional calculation means calculates a cross-sectional shape and a three-dimensional shape by interpolating the three-dimensional coordinates of the points calculated by stereo matching.

  The three-dimensional information presenting means projects the calculated three-dimensional information onto a window of the two-dimensional map, or presents another window on which the three-dimensional information is displayed on the display 12.

  When the survey execution location is specified on the two-dimensional map, the corresponding location detection means detects the corresponding location of the survey execution location on the TLS captured image.

  Below, operation | movement of the terminal 1 according to the instruction | command of the geographic information measurement program 3 is demonstrated.

(First operation example)
When the user inputs a request for presenting a specific part of the two-dimensional map, the two-dimensional map presenting means refers to the section center coordinate table stored in the storage unit of the terminal 1 and determines the latitude and longitude of the center of the part. get. As an example of the presentation request method, input of a place name of a section and scrolling of the screen of the presently presented section can be considered. In addition, the latitude / longitude may be specified directly. Subsequently, the two-dimensional map presentation means transmits a data transfer request for the presentation request section to the geographic data providing server 4 together with the latitude and longitude of the center of the section. The geographic data providing server 4 refers to the spatial index file f <b> 1, calls the corresponding map section data, and transfers it to the terminal 1. The two-dimensional map presentation means displays on the display 12 a window W1 on which the received map section data is displayed. In this way, by adopting a method of acquiring geographic data from the server, the amount of information stored in the storage unit of the terminal can be reduced, and updating of the geographic data is facilitated.

  FIG. 2 shows a display example of the window W1. As shown in FIG. 2, eight-way scroll buttons are displayed at the lower left portion of the window W1. By clicking the scroll button, the screen of the two-dimensional map displayed in the window W1 is scrolled. A scale adjustment bar is displayed in the lower left middle part of the window W1. By operating the scale adjustment bar, the scale of the two-dimensional map displayed in the window W1 can be increased or decreased.

  An operation tool section such as a layer selection section is displayed in the upper left middle part of the window W1. The two-dimensional map data is composed of a layer displaying a road, a layer displaying a building, a layer displaying a railway track, and a plurality of layers such as an ortho image, but the window W1 can be adjusted by adjusting the layer selection unit. You can select the layer displayed in In the example of FIG. 2, the uppermost layer is an ortho image occupying the center of the window. By displaying the ortho image superimposed, the user can obtain detailed information on the ground surface.

  A latitude / longitude / elevation display column is displayed at the upper left portion of the window W1. In the latitude / longitude / elevation display column, the latitude, longitude and altitude of the position where the mouse point is located are displayed. However, the elevation is displayed only when the mouse point is at a point where elevation data exists.

  The TLS image presentation means calls a TLS captured image obtained by capturing the center of the section displayed in the window W1 or a TLS captured image obtained by capturing a point designated by the mouse point from the TLS original image / GPS data file f3, and obtains the TLS image. A window W2 on which an image is displayed is presented on the display 12. As a specific example of the process of calling the TLS captured image, the latitude / longitude of the projection center point Cn is captured at the same time as the direct view image Pn closest to the center point of the display section or the designated point, and the direct view image Pn. There is a method of calling the forward view image Pf and the back view image Pb.

  FIG. 3 shows a display example of the window W2. As shown in FIG. 3, a front view image Pf, a direct view image Pn, and a rear view image Pb are displayed on the left side, the center, and the right side of the window W2, respectively. An ortho image of the corresponding portion is displayed at the upper end of the window W2.

  Next, an operation for calculating and displaying the three-dimensional coordinates of a point designated in the nadir view image Pn will be described. FIG. 5 is a flowchart showing the processing procedure of this operation.

The user displays the window W2 on the display 12 (S501). The user designates one point (one pixel) of the directly viewed image Pn displayed in the window W2 as the measurement point M _p (S503).

The three-dimensional information calculation means acquires a calculation window A centered on the measurement point M _p (S505). The calculation window is, for example, a 3 × 3 pixel matrix centered on the target point. The three-dimensional information calculation means calculates an epipolar line at the measurement point M _p in the front view image Pf and the back view image Pb (S507).

The three-dimensional information calculation means acquires calculation windows B _fk , B _bk (k = 1, 2, 3,...) Centering on one point on the epipolar line (S509), and the calculation windows B _fk , B A correlation value between each of _bk (k = 1, 2, 3...) and the calculation window A is calculated (S511). The processing of S509 and S511 is performed for each point extracted from the epipolar line at a constant interval, and calculation windows B _fk and B _bk showing the highest correlation value are detected. The centers of these calculation windows B _fk and B _bk are the corresponding points of the measurement point M _p (S513).

The three-dimensional information calculation means calculates a forward line-of-sight vector for the measurement point M _{p in} each of the forward-view image Pf, the direct-view image Pn, and the rear-view image Pb based on the position of the measurement point M _p in each captured image and the coordinates of the projection center point. A direct gaze vector and a back gaze vector are calculated (S515).

The three-dimensional information calculation means calculates the three-dimensional coordinates of the measurement point M _p from the line-of-sight vector (S517). As an example of a specific calculation method, an intersection approximation point is calculated for each of two sets of (front line of sight vector, direct line of sight vector) and (direct line of sight vector, back line of sight vector), and the difference between the two intersection approximate points is a threshold value. If the difference between the two intersection approximation points is equal to or greater than a threshold value, the middle point is determined as the position of the measurement point M _p in the three-dimensional coordinates. It is conceivable that a set of intersection approximate points including the line-of-sight vector corresponding to the window B for measurement is the position of the measurement point M _p . Here, the intersection approximate point of two line-of-sight vectors refers to the midpoint of both points when the distance between one point on one line-of-sight vector and one point on the other line-of-sight vector is the shortest. By using the three line-of-sight vectors in this way, it is possible to prevent a ranging error when the measurement point is hiding any obstacle.

The three-dimensional information presenting means displays the calculated altitude at the position of the measurement point M _{p on} the two-dimensional map displayed in the window W1 (S519).

In this operation example, the corresponding points in the front view image Pf and the back view image Pb of the measurement point M _p specified in the direct view image Pn are automatically calculated, but the user himself / herself himself / herself the front view image Pf and the back view image Corresponding points may be designated on Pb.

(Second operation example)
Next, an operation for calculating and displaying the vertical cross section (cross section passing through the line segment and perpendicular to the horizontal plane) designated by the nadir view image Pn will be described. FIG. 6 is a flowchart showing the processing procedure of this operation.

The user displays the window W2 on the display 12 (S601). As shown in FIG. 3, the user designates a line segment M ₁ (line segment AB in FIG. 3) in the nadir view image Pn displayed in the window W2 (S603).

The three-dimensional information calculation means _detects feature points M _lpk (k = 1, 2, 3,...) On the line segment M _l (S605). Here, the feature point extraction processing will be described in detail with reference to FIG. FIG. 7 is a flowchart showing a processing procedure for extracting feature points on the line segment _Ml .

One end point of the line segment M _l is defined as a start point, and the other end point is defined as an end point. In this feature point extraction processing, the target point on the line segment M _l extracted at regular intervals toward the end from the starting point of the line segment M _l will, each interest point is determined whether or not the condition is satisfied feature points To do. Specifically, first, a point on the line segment M _l that is a fixed distance L away from the previous point of interest toward the end point is set as the current point of interest (S701). However, at the time of initial point of interest extraction, a point that is separated from the start point by a distance L is regarded as the point of interest, and the start point is regarded as the previous point of interest.

  When the difference between the pixel value (density) of the point of interest and the pixel value of the previous point of interest is equal to or greater than the threshold value and the distance from the feature point detected immediately before is equal to or greater than the threshold value, It is set as a feature point (S703, S705, S707). If any of these conditions is not satisfied, the current point of interest is not used as a feature point (S703, S705, S709). When the current point of interest becomes the end point (or when the end point elapses), the above process ends (S711).

The three-dimensional information calculation means detects the corresponding points in the forward view image Pf and the backward view image Pb of each feature point M _lpk by the area correlation method (S607), and the forward view vector and the direct view vector for each feature point M _lpk. Then, the backward line-of-sight vector is calculated (S609), and the three-dimensional coordinates of each feature point _Mlpk are calculated from the line-of-sight vector (S611). The processes of S607 to S611 are the processes of S505 to _S517 repeated for each feature point _Mlpk .

The cross-sectional shape / three-dimensional calculation means calculates the three-dimensional coordinates of the line segment M _l by interpolating the feature points M _lpk to obtain the vertical cross-sectional shape of the line segment M _l (S613).

The three-dimensional information presenting means displays the elevation of the feature point M _lpk on the two-dimensional map of the window W1, and displays the cross-sectional shape in another window W3 (S615). FIG. 4 shows an example of the calculated cross-sectional shape.

(Third operation example)
Next, an operation for calculating and displaying the three-dimensional shape of the region designated by the nadir view image Pn will be described. FIG. 8 is a flowchart showing the processing procedure of this operation.

The user displays the window W2 on the display 12 (S801). The user specifies an area M _a in nadir images Pn displayed in the window W2 (S803).

Three-dimensional information computing means detects the feature points of the region _{M a M apk (k = 1,2,3} ···) (S805). This feature point extraction processing is performed by repeating the processing of S701 to S711 for pixel rows at regular intervals.

The three-dimensional information calculation means detects corresponding points in the forward-view image Pf and the rear-view image Pb of each feature point M _apk by the area correlation method (S807), and the forward view vector and the direct view vector for each feature point M _apk. Then, the backward line-of-sight vector is calculated (S809), and the three-dimensional coordinates of each feature point _Mapk are calculated from the line-of-sight vector (S811). Processing S807~S811 are those repeated for each feature point M _apk processing S505～S517.

Sectional shape / three-dimensional calculation means obtains a three-dimensional shape by calculating the 3-dimensional coordinates of the area M _a by interpolating the feature point M _apk (S813).

The three-dimensional information presenting means displays the three-dimensional shape of the feature point _Mapk superimposed on the two-dimensional map of the window W1 (S815).

(Fourth operation example)
In the above operation example, a survey point, survey line, or survey area is specified on the nadir image displayed in the window W2, but instead, an operation that specifies a survey point, survey line, or survey area on the two-dimensional map. Examples are also possible. In this operation example, when the user designates a survey point, a survey line, or a survey area on the two-dimensional map of the window W1, the corresponding location detection means uses the horizontal plane two-dimensional coordinates (latitude / longitude) of the projection center point Cn as the survey. Call the direct view image Pn closest to or included in the horizontal plane two-dimensional coordinates (latitude / longitude) of the point, survey line or survey area from the storage unit of the terminal 1, and detect the corresponding point of the survey point in the direct view image Pn. . As a specific example, it is conceivable to perform texture analysis on the nadir view image Pn and detect a corresponding portion from a relative position with a divided region corresponding to a certain feature region on the two-dimensional map.

  The processing procedure of the operation for calculating and displaying the three-dimensional information of the area designated by the two-dimensional map will be described. FIG. 9 is a flowchart showing the processing procedure of this operation.

The user designates an area M _a1 on the two-dimensional map (S901). As a designation method, a method of tracing a region on a two-dimensional map with a mouse point, a method of designating a feature region (building, station, parking lot, etc.) registered in association with the two-dimensional map can be considered.

Corresponding portion detecting means, call the nearest nadir images Pn on the latitude and longitude of the center point of the latitude and longitude of the projection center point Cn is region M _a1 from TLS original image / GPS data file f3 (S903). Corresponding location detection means texture-analyze the nadir view image Pn and divide it into a plurality of regions (S905). As an example of texture analysis. There are area expansion method and level rising method.

Corresponding part detection means detects the divided area M _a2 of the directly viewed image Pn corresponding to the feature area overlapping the area M _a1 (S907). As an example of a specific method, a line where a corresponding region can exist in the nadir image Pn is calculated from the two-dimensional coordinates of the feature region and the three-dimensional coordinates of the projection center point Cn, and the divided region and the feature region on the line are calculated. It is conceivable to compare the shapes.

The corresponding part detecting means detects the feature point M _a2pk (k = 1, 2, 3,...) _Of the divided region M _a2 (S909), and each feature point M in the forward view image and the backward view image by the area correlation method. _a2pk (corresponding point) is detected (S911), the line-of-sight vector in the forward view image, the nadir view image and the back view image is calculated (S913), and the three-dimensional coordinates of each feature point are calculated from the line-of-sight vector (S915), A three-dimensional shape of the divided region M _a2 is calculated by interpolating the three-dimensional coordinates of each feature point M _a2pk (S917). The processing of S909 to S917 is the same as the processing of S803 to S811. The three-dimensional shape of the divided area _Ma2 calculated by fitting the structure template, clustering the structure segment, or the like may be shaped.

Three-dimensional information presenting means, and displays superimposed three-dimensional shape of the divided region M _a2 in the two-dimensional map of the window W1, and presents a three-dimensional shape of the divided region M _a2 in a separate window (S919).

  FIG. 10 is a diagram showing a modification of the above embodiment. 1 and 10, the same reference numerals are assigned to the same components. FIG. 1 shows an example in which a geographic information measurement program is composed of six modules, but the present invention is not limited to such a geographic information measurement program having a module configuration. For example, as shown in FIG. 10, a combination of a two-dimensional geographic information program 3b and a three-dimensional geographic information measurement program 3a is used, and position information and image information are exchanged between the two programs as a whole. It is also possible to achieve the function of the geographic information measurement program 3. In the user's computer, both programs are processes that operate independently, but the exchange between the programs is performed using an inter-process communication technique prepared in the operating system of the user's computer. In particular, in the modified example of FIG. 10, since the two-dimensional geographic information program 3b and the three-dimensional geographic information measurement program 3a are combined at the time of execution, a program on one hand can be used, or two-dimensional geographic information can be obtained from different vendors. It is preferable that the program 3b and the three-dimensional geographic information measurement program 3a are provided or that only one of them is upgraded, so that operational flexibility can be provided. The method of dividing the program is not limited to that shown in FIG. Further, a distributed system in which a plurality of divided programs are executed by several computers may be used. In short, various embodiments are possible in accordance with the spirit of the present invention.

FIG. 1 is a functional configuration diagram of a GIS to which the geographic information measurement program 3 of the present embodiment is applied. FIG. 2 shows a display example of the window W1. FIG. 3 shows a display example of the window W2. FIG. 4 shows an example of the calculated cross-sectional shape. FIG. 5 is a flowchart showing a processing procedure of an operation for calculating and displaying the three-dimensional coordinates of a point designated in the nadir view image Pn. Figure 6 is a flowchart illustrating a processing procedure of an operation of calculating and displaying a vertical cross-section of a line segment M _l specified in nadir image Pn. FIG. 7 is a flowchart showing a processing procedure for extracting feature points on the line segment _Ml . FIG. 8 is a flowchart showing a processing procedure of an operation for calculating and displaying a three-dimensional shape of a region designated by the nadir view image Pn. FIG. 9 is a flowchart showing a processing procedure of an operation for calculating and displaying the three-dimensional information of the region designated by the two-dimensional map. FIG. 10 is a diagram illustrating a modification of the embodiment.

Explanation of symbols

DESCRIPTION OF SYMBOLS 1 ... Computer terminal, 12 ... Display, W1, W2, W3 ... Window, Pw ... Ortho image, Pf ... Front view image, Pn ... Direct view image, Pb ... Rear view image, 2 CD-ROM.

Claims (8)

2D map display means for displaying a 2D map on a display;
An image obtained by imaging each location on the two-dimensional map and picked up from the two-dimensional image stored in advance in the storage means and imaging the designated location on the two-dimensional map is selected and displayed on the display. Two-dimensional image display means for displaying;
It is obtained by capturing the three-dimensional coordinates of a specified point, line, or region on the displayed two-dimensional image from an angle different from that of the two-dimensional image and the point, line, or region. A geographic information measuring system comprising: three-dimensional information measuring means for calculating by stereo matching with another two-dimensional image. A first step of causing a computer to display a two-dimensional map on a display;
By designating a certain location on the two-dimensional map, the designated location is obtained from the two-dimensional images stored in advance in the storage means obtained by imaging each location on the two-dimensional map on the computer. A second step of selecting a captured image and displaying it on the display;
By specifying a point, line, or region on the displayed two-dimensional image, the computer can obtain the two-dimensional image and the point, line, or region from a different angle from the two-dimensional image. And a third step of stereo-matching with another two-dimensional image to calculate the three-dimensional coordinates of the point, line or region. On the computer,
Obtained by imaging each location on the 2D map, selected from the 2D images stored in advance in the storage means as an image of a specific location on the 2D map, and displayed on the display Another point obtained by imaging the three-dimensional coordinates of a point, line, or region designated on the two-dimensional image from a different angle from the two-dimensional image and the point, line, or region. A geographic information measurement program characterized by causing a step of calculating by stereo matching with a two-dimensional image. On the computer,
An image obtained by imaging each location on the two-dimensional map and imaging a location including a point, line, or region designated on the two-dimensional map from the two-dimensional images stored in the storage means in advance. A first step of selecting and detecting corresponding points of the points, lines or regions in the selected two-dimensional image;
The three-dimensional coordinates of the corresponding part are calculated by stereo-matching the two-dimensional image with another two-dimensional image obtained by imaging the corresponding part from a different angle from the two-dimensional image. 2. A geographic information measuring program characterized in that two steps are executed. 4. In the step, a point in the single two-dimensional image is designated, and a point on the other two-dimensional image corresponding to the designated point is detected by an area correlation method. The described geographic information measurement program. Used to obtain the cross-sectional shape of the measurement object present on the two-dimensional image,
In the step, a line passing through the cross section in a single two-dimensional image is designated, and points on the other two-dimensional image corresponding to a plurality of feature points on the designated line are detected by an area correlation method. In addition, the three-dimensional coordinates of each feature point are calculated,
The geographic information measurement program according to claim 3, wherein the computer is caused to execute a second step of interpolating the three-dimensional coordinates of the feature points to generate an image showing the cross-sectional shape. In the step, a region in a single two-dimensional image is designated, and points on the other two-dimensional image corresponding to a plurality of feature points in the region are detected by an area correlation method. 3D coordinates of feature points are calculated,
The geographic information measurement program according to claim 3, wherein the computer is caused to execute a third step of interpolating the three-dimensional coordinates of the feature points to generate an image indicating the shape of the region. The three-dimensional coordinates of the measurement object specified on the two-dimensional image obtained by imaging a specific location on the two-dimensional map are imaged from the angle different from the two-dimensional image and the measurement object. 3D information measuring means for calculating by stereo matching with another 2D image obtained by
3D information display means for displaying the 3D coordinates calculated by the 3D information measuring means or the three-dimensional shape of the measurement object obtained based on the 3D coordinates on the display so as to overlap the 2D map; Geographic information measurement system characterized by having.

Images