JP2015111248A - Fault diagram creation device, fault diagram creation method and fault diagram creation program - Google Patents

Abstract

PROBLEM TO BE SOLVED: To obtain a fault diagram creation device that easily reads a fault in which a width ranges over about several meters from an image, and enables a fault line diagram based on the reading result to be easily obtained.SOLUTION: A fault diagram creation device comprises: a red relief image creation unit 10; a fault search range setting unit 12; an imaging unit 14; a range definition unit 16; an existing/red relief read drawing unit 18; a γ-ray survey point/fault line drawing unit 24; and a line adjustment unit 26. A γ-ray local survey result linear line LGj (yellow) based on the local survey result is configured to indicate the γ-ray local survey result linear line LGj substantially in parallel to a line having a plurality of drawn tiny fault lines arranged in a red relief image RGi by reading from the red relief image RGi, and capable of being caused to easily correspond to a segment across tiny fault lines DLj.

Description

  The present invention relates to a tomographic map creation device capable of easily and accurately connecting tomographic lines based on local measurement results on an image of a region acquired by an aircraft.

  Inland earthquakes often occur, and it is important to understand the location and activity of faults. Active faults are the cause of inland earthquakes, but new faults that have not been known in the past are often discovered recently.

  Such determination of the position of the tomography is mainly based on a method in which a stereo photograph is taken with an aircraft and the stereo photograph is stereoscopically detected to detect the tomography.

  However, in areas where slopes covered with vegetation are widely distributed, it is difficult to recognize a fine fault displacement topography (hereinafter referred to as a fine fault) of 10 m or less.

  On the other hand, a map creation device that generates a DEM (Numerical Elevation Model) that emits a laser from an aircraft, obtains reflected light, and expresses the shape of the ground surface, and emphasizes the unevenness of the ground surface using the DEM Is disclosed in Patent Document 1.

  The topographic map creation device of Patent Document 1 emphasizes fine irregularities to represent the fine topography of the ground surface based on the elevation data of the digital elevation model. That is, the unevenness is emphasized by amplifying the height difference. The saturation of the unevenness is determined in proportion to the altitude. Further, the hue is determined by the altitude, and a color image obtained by combining these is obtained.

JP 2008-242298 A

  However, micro faults that are considered active faults are generally several meters wide. However, since the topographic map creating apparatus of Patent Document 1 is to emphasize and emphasize the difference in altitude of fine irregularities while maintaining a local average value of altitude, it can detect fine faults with a width of several meters. It is unsuitable.

  Moreover, even if the location of the fault is read from the photograph and drawn on the map, the fault is generally several kilometers or more. In addition, faults are complicated and complicated.

  For this reason, it is not easy to determine which fault is connected to which fault and how long it extends in which direction.

  Furthermore, the fault length is said to depend on the fault width and ground activity. For this reason, it is desirable to connect the fault lines in consideration of the fault width, activity, and the direction of the fine fault measured in the field.

  The present invention has been made in order to solve the above problems, and it is easy to accurately detect fault lines in the direction of a fault with a fault length reflecting a local measurement result on an image of a target area where a fault line is drawn. It is an object of the present invention to obtain a tomographic map creating apparatus capable of obtaining a tomographic map that can be well connected.

The tomographic map creating apparatus according to the present invention includes a first storage means for storing a digital elevation model obtained by emitting laser pulses from an aircraft to a region at intervals capable of detecting fine tomographic lines DLj;
A second storage means for storing the red three-dimensional image RGi;
The fault length L of the fine fault line DLj and the gamma ray measurement obtained based on the gamma ray intensity and the measurement interval of a plurality of measurement points at the gamma ray measurement point Pi in the area determined based on the red three-dimensional image RGi. A second storage means in which the direction of the fine fault line DLj at the point Pi is stored as the local γ-ray measurement information Fγi;

A tomographic detection range setting unit for displaying an input box for setting a search range in which the tomographic zone width FW of the fine tomographic line DLj can be detected on the screen of the display unit;
The red three-dimensional image RGi based on the numerical elevation model of the first storage unit and the search range set by the tomographic detection range setting unit is generated and stored in the second storage unit and displayed on the screen. A red three-dimensional image creation unit to be displayed;
A surface for drawing the fine tomographic line DLj is displayed on the red three-dimensional image RGi on the screen, and the movement trajectory of the cursor on this surface is defined as the fine tomographic line DLj, which is displayed in a predetermined color. Means,
On the surface for drawing the fine tomographic line DLj, a surface for designating the γ-ray measurement point Pi is displayed, and the point designated on this surface is used as the γ-ray measurement point Pi for the fine tomographic line. Means for displaying in a color different from DLj;
A surface for displaying the γ-ray field measurement result straight line LGj is displayed above or below the surface for drawing the fine fault line DLj, and is defined as a surface for designating the γ-ray measurement point Pi on this surface. The γ-ray measurement point Pi is defined, and a straight line of the direction and the fault length L of the local γ-ray measurement information Fγi of the second storage means is drawn from the γ-ray measurement point Pi, means for displaying the γ-ray field measurement result line LGj in a color different from the γ-ray measurement point Pi and the fine tomographic line DLj;
Means for assigning the direction to the γ-ray field measurement result straight line LGj;
When the γ-ray field measurement result line LGj is designated by a cursor, the direction assigned to the γ-ray field measurement result line LGj is read, and the γ-ray field measurement result is made to follow the movement of the cursor. And a means for moving the straight line LGj in the read direction.

As described above, according to the present invention, the fine tomographic line LGj is read from the red three-dimensional image RGi displayed on the screen, and γ-rays are measured at a plurality of points at the local point where the fine tomographic line LGj exists. The obtained local γ-ray measurement information Fγi is stored.
When the γ-ray measurement point Pi is specified on the red stereoscopic image RGii on the screen, the tomographic length L included in the local γ-ray measurement information Fγi is measured based on the γ-ray measurement point Pi. The direction obtained in the field included in the information Fγi is displayed, and the γ-ray field measurement result straight line LGj is moved in the opposite direction to the above-described direction and displayed.

  For this reason, the γ-ray field measurement result straight line LGj is displayed so as to be substantially parallel to each small fine tomographic line DLj displayed on the red three-dimensional image RGi. Can be easily determined whether the fault is connected by a single line.

1 is a schematic configuration diagram of a tomographic drawing creation apparatus according to Embodiment 1. FIG. It is explanatory drawing of DEM data. It is a flowchart explaining the operation | work procedure until it obtains a hybrid fault line. It is explanatory drawing of the red three-dimensional image RGi. It is explanatory drawing of a fault. It is explanatory drawing of the fine tomographic line DLj seen by the red three-dimensional image RGi. It is explanatory drawing of a fault. It is comparison explanatory drawing of a contour map and red three-dimensional image RGi. It is comparison explanatory drawing of a contour map and red three-dimensional image RGi. It is explanatory drawing of the 2nd process (2) of FIG. It is explanatory drawing of the 3rd process (3) of FIG. It is explanatory drawing of the local measurement condition. It is a measurement figure of the gamma ray measurement of the field measurement situation. It is explanatory drawing of the measuring method of the gamma ray measurement of a local measurement condition. It is explanatory drawing of the measuring method of the gamma ray measurement of a local measurement condition. It is explanatory drawing of the 4th process (4) of FIG. It is a flowchart explaining the process of the 5th process of FIG. It is a flowchart explaining the process of the 5th process of FIG. It is explanatory drawing of a display of the gamma ray measurement point Pi. 6 is an explanatory diagram of a dialog box BPi of a γ-ray measurement point / tomographic line drawing unit 24. FIG. It is explanatory drawing of the related information of the gamma ray measurement point Pi. It is explanatory drawing of the gamma ray field measurement result straight line LGj. It is explanatory drawing of a movement of the gamma ray field measurement result straight line LGj. It is explanatory drawing of preparation of the hybrid tomographic line HBLj. It is explanatory drawing of deletion of the gamma ray field measurement result straight line LGj. FIG. 6 is a schematic configuration diagram of a tomographic map creating apparatus according to a second embodiment. FIG. 12 is an explanatory diagram of a display screen of the tomographic map creating device according to the second embodiment. It is a principle explanatory view of the ground opening and the underground opening. It is explanatory drawing of the main patterns of the ground opening degree and underground opening degree. It is a three-dimensional explanatory diagram of the ground opening and the underground opening. It is explanatory drawing of the sample point and distance of a ground opening degree and an underground opening degree. FIG. 6 is a schematic configuration diagram of a tomographic drawing creation apparatus according to a third embodiment. It is explanatory drawing of the hierarchy of a layer order | rank setting part. 4 is an explanatory diagram of a data configuration of a memory 80. FIG. It is explanatory drawing of the process of obtaining the gamma ray field measurement result straight line LGj. It is explanatory drawing of a movement of the gamma ray field measurement result straight line LGj. It is explanatory drawing of the attribute information of the hybrid fault line HBLj. It is explanatory drawing of drawing of the hybrid tomographic line HBLj. It is explanatory drawing of the layer order | rank setting part at the time of making the 1st layer into the memory 76 for γ-ray field measurement result straight line working. It is explanatory drawing of the gamma ray field measurement result straight line LGj at the time of making the 1st layer into the memory 76 for gamma ray field measurement result straight line working.

The following embodiments exemplify devices and methods for embodying the technical idea of the invention, and the technical idea of the invention is not limited to the following structure and arrangement. .
The technical idea of the invention can be variously modified within the technical scope described in the claims. It should be noted that the drawings are schematic and the configuration of the apparatus and system is different from the actual one.

  In general, a fine fault that seems to be an active fault has a fault length of several kilometers, and this fault length is said to depend on the fault width. The fault length is said to be affected by the ground activity. A gamma ray measuring device is used to detect this activity.

  That is, in the first embodiment, the red three-dimensional image RGi obtained from the aviation laser measurement data is displayed on the screen, the fine tomography is read, the read fine tomography is drawn, and the γ-ray measurement result measured on site Based on the length and direction of the fault using, the fine faults can be easily connected with high accuracy.

<Embodiment 1>
FIG. 1 is a schematic configuration diagram of a tomographic drawing creation apparatus according to the first embodiment. As shown in FIG. 1, the tomographic map creating apparatus of the first embodiment displays the hybrid tomogram HHGi on the screen of the display unit 200 with the following program configuration in the personal computer main unit 100.

  The personal computer main body 100 includes a red three-dimensional image creation unit 10, a tomographic search range setting unit 12, an imaging unit 14, a range definition unit 16, an existing / red interpretation / drawing unit 18, a γ-ray measurement point, A program configuration such as a tomographic line drawing unit 24 and a line adjustment unit 26 is provided.

  Further, as shown in FIG. 1, the PC main body 100 has a DEM database 30 storing a digital elevation model (DEM) and a memory 32 storing a red three-dimensional image RGi, as shown in FIG. The memory 39 of the imaging unit 14 and the fine tomographic line definition memory 34 for drawing a fine tomographic line having a tomographic width of 10 m or less (hereinafter referred to as a fine tomographic line DLj) read based on the red three-dimensional image RGi. And an existing fault line memory 36 in which an existing fault line (hereinafter referred to as an existing fault line KLj) is drawn, a memory 38 in which a hybrid fault line HBLj is drawn, and on-site measurement in which on-site γ-ray measurement information Fγi is stored. A result information memory 40 is provided. In the present embodiment, the fine fault line DLj has a width of the fault line of about 10 m or less.

  The digital elevation model stored in the aforementioned DEM database 30 is a desired grid interval d (for example, 0.2 m, 0.5 m, 1 m or 1,5 m, 2.5 m) in the entire measurement area. Etc.), and filtering to remove buildings and trees other than the ground surface from the elevation data measured mainly by the last returned pulse (laser pulse) of the laser reflected pulses. This is ground level altitude data obtained by interpolation.

  A specific example is shown in FIG. As shown in FIG. 2, the DEM data Di includes X coordinates (longitudes Xi) and Y coordinates (latitudes Yi) of the center points of the respective lattices assigned lattice numbers i (i = 1, 2,..., N). , Z coordinates (ground elevation value Zgi).

  As an example of the elevation interpolation method, create a contour map connecting the same elevation values, create an irregular triangle network (TIN) for this contour map, restore the ground, and the intersection of TIN and each grid There is a way to find the height of the.

The red three-dimensional image creation unit 10 sequentially sets sample grid points in the digital elevation model of the DEM database 30 and searches the tomographic range set by the tomographic search range setting unit 12 for each sample grid point. In this fault search range, three parameters of the slope, the ground opening, and the underground opening are obtained based on DEM (Digital Elavation Model) data Di, and the plane distribution is obtained as a gray scale image. save.
The tomographic search range setting unit 12 displays an input box (not shown) for setting a search range in which the tomographic zone width FW of the fine fault line DLj described above can be detected.

  Then, by putting gray in the difference image between the ground opening and the underground opening and creating a pseudo color image with the slope in the red channel, the ridges and mountain peaks are whitish and the valleys and depressions are darkened. Then, a red three-dimensional map RGi expressed in red as a portion with a steep slope is generated and stored in the memory 32. The red three-dimensional image creation unit 10 will be described later with reference to the drawings.

  The imaging unit 14 combines (superimposes) the various output images and displays them on the display unit 200. At this time, the color information added to the image is displayed on the screen.

  The range definition unit 16 reads the planar range (X, Y) of the red three-dimensional image RGi stored in the memory 32 via the red three-dimensional image creation unit 10 and reads these ranges into the fine tomographic line definition memory 34, The existing tomographic line memory 36 and the memory 38 are defined. The plane defined in the fine fault line definition memory 34 is referred to as a red interpretation result plane, the plane defined in the existing fault line memory 36 is referred to as an existing plane, and the plan view defined in the memory 38 is hybridized. It is called a plane for use.

  That is, the red interpretation result plane of the fine tomographic line definition memory 34 has the origin coordinates of the red three-dimensional image RGi and is defined as a mesh size having the resolution of the red three-dimensional image RGi.

  Further, the existing plane of the existing tomographic line memory 36 has the origin coordinates of the red three-dimensional image RGi and is defined to have a mesh size of the resolution of the red three-dimensional image RGi or a mesh size finer than this resolution. .

  Further, the hybrid plane of the memory 38 has the origin coordinates of the red three-dimensional image RGi and is defined as a mesh size having the resolution of the red three-dimensional image RGi.

  When there is an instruction to draw the fine tomographic line DLj on the red three-dimensional image RGi on the screen, the existing / red interpretation / drawing unit 18 sets the coordinates associated with the movement of the cursor on the fine tomographic line as the cursor moves on the screen. The line is defined by converting to the coordinates of the red interpretation result plane in the line definition memory 34, and the line defined on the plane is displayed as the fine tomographic line DLj by the imaging unit 14.

  Further, when there is an instruction to draw an existing tomographic line KLj on the red three-dimensional image RGi on the screen, the existing / red interpretation / drawing unit 18 stores the coordinates in the existing tomographic line memory as the cursor moves on the screen. The coordinates are defined by converting into the coordinates of 36 existing planes, and the lines defined in the planes are displayed by the imaging unit 14 as existing tomographic lines KLj.

  Further, the field measurement result information memory 40 has γ rays obtained based on γ ray intensities and measurement intervals of a plurality of measurement points at the γ ray measurement points Pi in the field determined based on the red three-dimensional image RGi on the screen. The fault length L of the local measurement result line LGj and the direction Lθi of the γ-ray local measurement result line LGj at the local γ-ray measurement point Pi are stored. These are collectively referred to as local γ-ray measurement information Fγi. In this local γ-ray measurement information Fγi, a line number LGBi of the γ-ray local measurement result straight line LGj is generated and associated with the γ-ray measurement point Pi.

  The γ-ray measurement point / tomographic line drawing unit 24 reads the point designated by the cursor in the red three-dimensional image RGi displayed on the screen of the display unit 200 as the γ-ray measurement point Pi, and this γ-ray measurement point Pi is read. This is defined by converting into the coordinates of the hybrid plane in the memory 38. Then, the line having the length of the fault length L in the direction Lθi included in the local γ-ray measurement information Fγi stored in the local measurement result information memory 40 is displayed as the γ-ray local measurement result straight line LGj. Then, the γ-ray field measurement result straight line LGj is designated by the cursor, and with the movement of the cursor, the coordinates of the cursor are converted to the coordinates of the red interpretation result plane in the fine tomographic line definition memory 34 and displayed. This movement is made in the direction opposite to the direction Lθi.

The line adjustment unit 26 compares the γ-ray field measurement result straight line LGj and the existing fault line KLj on the screen, and if there is an input of a creation instruction of the hybrid fault line HBLj, the coordinates of the cursor are displayed on the red interpretation result plane. Define the coordinates of.
As the cursor moves, the moving coordinates are converted to the coordinates of the red interpretation result plane in the fine tomographic line definition memory 34 and displayed as hybrid tomographic lines HBLj (pink).

  That is, when the fine tomographic lines DLj displayed on the screen are connected by lines while being traced, these lines and the connected fine tomographic lines DLj are displayed in different colors (pink) as the hybrid tomographic lines HBLj. Let

(Detailed explanation)
FIG. 3 is a flowchart for explaining a work procedure until a hybrid fault line is obtained.

  As shown in FIG. 3, a laser is emitted to the target area by an aircraft, and a DEM having a mesh of 2.5 m or less (for example, 1 m) is obtained in the DEM database 30 (S100).

  Then, the red three-dimensional image RGi is generated by the red three-dimensional image generation unit 10 (S200). FIG. 4 shows a red three-dimensional image RGi generated by the red three-dimensional image creating unit 10. Generation by the red three-dimensional image creating unit 10 is also referred to as a first step (1).

  In the generation of the red three-dimensional image RGi, the search range is made small (about 50 to 100 m). This is a red three-dimensional image RGi that is a single image processed so that the terrain looks three-dimensional, and the image resolution is 1 m (for example), which is the same as the mesh size of aviation laser data. The key micro-topography can be fully read.

  Here, the reason why the fine tomographic line (fine tomographic line DLj (thick black dotted line)) can be detected from the red three-dimensional image RGi will be described. As shown in (a) to (d) of FIG. 5, a fault is a fault if it is shifted all at once along a line with a stepped terrain extending in a line or a terrain such as a ridge / valley. Suspected to be.

  However, even if the step terrain extends linearly, if the direction is parallel to the maximum slope direction or the river flow direction, there is a high possibility that the terrain is not a fault but an erosion terrain. In addition, stepped terrain may be formed even when the rocks forming the ground are different at a certain line and the erosion resistance of the two is different. If there is a terrain that is suspected to be a fault, imagine and compare several terrain formation processes including the surrounding area, and if it is easier to explain that the terrain factor is a fault, it is estimated as a fault To do.

  When the terrain suspected of being a fault is seen in the red three-dimensional image RGi, the red three-dimensional image RGi is calculated in a search range of about 50 to 100 m using, for example, a 1 m mesh DEM, as shown in FIG. Colors that are characteristically different from the surroundings are arranged in a linear fashion. This is because the amount of topography such as slope is different from the surroundings only in the fault topography. That is, as shown in FIG. 7, the cross section between A and B in FIG. 6 can detect a portion where the ground is displaced.

  Here, a description will be given in comparison with a contour map.

In the contour map, there is no information between the contour lines (see FIGS. 8B and 9B).
For this reason, it is difficult to follow the continuity of terrain such as steps, and as a result, the smaller the terrain scale (height and width of steps), the more difficult it becomes.

  In addition, when the vertical displacement due to the fault (less than 1 m when small) is smaller than the interval between contour lines (generally 10 m at 1/25000 and 1 m at 2500) (see FIGS. 8B and 9B) ) Cannot be expressed by contour lines.

  On the other hand, since the red three-dimensional image RGi is created in the search range of 50 to 100 m using a DEM of about 1 m mesh (which may be 0.2 m or 0.5 m), FIG. As shown in FIG. 9A, the fault line (fine fault line DLj) can be clearly read.

  Further, the existing material (map, geological map, etc.) is read (S300), and the fault is read in the red three-dimensional image RGi (S400).

  And what is known as a fault (existing fault line KLj) from the existing data is transferred onto the red three-dimensional image RGi (see FIG. 10A), while the fault (fine fault line) in the red three-dimensional image RGi. A portion that can be read as (DLj) is drawn with a line (see FIG. 10B) (S500). This is also referred to as a second step (2).

  That is, by using the red three-dimensional image RGi, a fine slice that cannot be interpreted by the conventional method (aerial photo interpretation) is interpreted, and the interpretation result is drawn.

  However, although it is possible to decipher a fault “possible terrain”, it may include non-fault lineaments (eg, terrain gaps developed along geological boundaries).

  For example, the fault continuity cannot be determined when it is covered with sediment or when the fault is later divided by erosion or artificial alteration of the mountain stream.

  For this reason, these are investigated on site by γ-ray measurement (S600). This is also referred to as a third step (3).

  Before investigating these by γ-ray measurement at the site, as shown in FIG. 11, for example, ridges, ridges, and rocks that are not covered with thick sediment are removed from the fault read in step (2). Select locations that are close by and have good approaches such as close to the roadway. However, in FIG. 11, contour lines are superimposed.

  In this way, for example, not only the existence, position, and direction of the fault can be estimated by reading the red three-dimensional map in step (2), but also suitable sites for γ-ray measurement (field survey) can be selected based on conditions such as topography and geology. be able to. It will also lead to more efficient field surveys. In addition, field surveys can be conducted at more points.

  In addition, by selecting an appropriate measurement location, such as a location where there is no thick sediment on the ground surface, the boundary position of the fault crush zone can be determined with an accuracy of the order of 10 cm.

  The measurement situation at the measurement location determined in this way will be described below. FIG. 12A is a photograph of a γ-ray measuring apparatus used at a local measurement location. FIG. 12B is a photograph of the measurement situation. There are a total of 6 γ-ray survey locations. The γ-ray measuring apparatus uses a scintillation survey meter (total counting method).

  The station interval is 3m, for example. Then, at a point where an abnormal value of γ-ray intensity appears, measurement is performed at an interval of 10 cm 2, and the γ-ray intensity abnormal value range (crush zone) is determined in detail.

  As shown in FIG. 12B, the γ-ray measuring device is installed on the ground parallel to the direction of the crushing band and perpendicular to the measuring direction.

  That is, a survey line is provided in a direction orthogonal to the direction of the estimated fault fracture zone. Next, strip off the surface soil. Next, the gamma ray intensity on the ground is measured while moving on the survey line, and recorded as shown in FIG.

  In FIG. 13, if the γ-ray intensity clearly changes suddenly in a certain section on the survey line compared to both sides, it can be estimated that the section is a fault fracture zone.

  And the parallel survey line is set, the same measurement is performed, and the position of a fault crush zone is specified. The width and direction of the crushing zone can be specified from the gamma ray measurement of the two survey lines.

  The measurement method will be specifically described with reference to FIGS. As shown in FIG. 14, in STEP “0”, a survey line A (rope) is set at the ground survey point determined in the third step (3).

  In STEP “1-1”, a gamma ray measuring device is arranged at a rough interval (for example, 5 m), and it is plotted whether it shows an abnormal value or a normal value.

  In STEP “1-2”, the interval is narrowed (for example, 0.5 to 2 m). In STEP “1-3”, the left and right boundaries of the abnormal value are narrowed to, for example, 0.1 m. That is, (a1) and (a2) in FIG. 14 are the boundaries of the abnormal section.

  Then, in “STEP 2-1” in FIG. 15, a survey line B (rope) is installed on the ground in parallel with the survey line A (rope).

  Then, in “STEP 2-2”, the γ-ray intensity at the measurement line B is measured (the same method as FIG. 14). Then, in “STEP 2-3”, the γ-ray abnormal boundary positions b1 and b2 in the measurement line B are determined. Next, in “STEP2-4”, the abnormal boundary positions a1 and a2 of the survey line A and the abnormal boundary positions b1 and b2 of the survey line B are connected. Then, in “STEP2-4”, the length of the line connecting the measurement line A and the measurement line B at a right angle is determined as the fault zone width FW.

  Then, the γ-ray measurement results are organized and analyzed using the fault zone width FW determined in this manner (S700). This is referred to as a fourth step (4) in the present embodiment.

  The fourth step (4) will be described with reference to FIG. In the fourth step (4), the fault zone width FW obtained in the field in the third step (3) is input to the personal computer to obtain the fault length L.

  For example, when the gamma ray intensity measurement results at the local survey lines A and B are drawn on a graph as shown in FIG. 16A, when the fault zone width FW is 8.8 m, FIG. As shown, the fault length L is obtained as L = 0.36FW + 0.62.

  In the present embodiment, since the fault zone width FW is 8.8 m, the fault length L is determined to be 3.8 Km. The coefficients “0.36” and “0.62” for obtaining the fault length L are examples.

  Then, these measurement data are input to a tomographic map creation device using the LrG method shown in FIG. 1, and a tomographic line (γ-ray field measurement result straight line LGj) is automatically created as shown in FIG. 3 (S800). This is referred to as a sixth step (6) in the present embodiment.

The automatic creation processing of the tomographic line (γ-ray field measurement result straight line LGj) will be described later with reference to the drawings.

  In addition, when automatically creating a fault line (γ-ray field measurement result straight line LGj), a fault line (fine fault line DLj) is drawn at a position considered to be a fault in the red three-dimensional image RGi as shown in FIG.

  Then, the fine fault line DLj obtained in the sixth step is corrected based on the γ-ray field measurement result straight line LGj (S900), and finally the hybrid fault line HBLj is obtained (S1000).

  In the following, the tomographic line automatic creation process (γ-ray measurement point / tomographic line drawing unit 24) in the tomographic map creation apparatus using the LrG method performed in the fifth step will be described with reference to the flowchart of FIG.

  First, as shown in FIG. 17, initial setting is performed (S801).

  In this initial setting, the tomographic search range setting unit 12 sets a search range (for example, 50 m) in the red three-dimensional image creation unit 10, and the red three-dimensional image creation unit 10 converts the red three-dimensional image RGi shown in FIG. Displayed on the display unit 200.

  The reason why the search range is set to about 50 m is that, in a normal red three-dimensional map, the fine terrain and the large terrain are expressed so that they can be grasped at the same time.

  However, when the purpose is specialized in the interpretation of fault terrain, it is possible to make the image easy to recognize only the fine terrain and reduce the calculation time by setting the search range to about 50 m. is there.

  Further, it is assumed that a red interpretation result plane is defined by the range definition unit 16, an existing plane is defined in the existing tomographic line memory 36, and a hybrid plane is defined in the memory 38.

  Furthermore, it is assumed that the fault length L is obtained by the field survey and the direction Lθi of the fault length L is obtained.

  Further, in the second step, the existing tomographic line KLj (black) based on the existing material is drawn on the red three-dimensional image RGi as shown in FIG. 10A, and the red three-dimensional map RGi is shown in FIG. 10B. It is assumed that a fine fault line DLj (black dotted line) is drawn in FIG.

  In this state, as shown in FIG. 17, it is determined whether the γ-ray measurement point Pi is set on the red three-dimensional map RGi displayed on the display unit 200 by the γ-ray measurement point / tomographic line drawing unit 24. (S803).

  If it is determined in step S803 that the γ-ray measurement point Pi is set on the red three-dimensional map RGi, the γ-ray measurement point / tomographic line drawing unit 24 sets the γ-ray measurement point Pi as a hybrid plane in the memory 38. At the same time, as shown in FIG. 19, the γ-ray measurement point Pi is displayed in color (for example, yellow) (S805).

  Then, as shown in FIG. 20, the γ-ray measurement point / tomographic line drawing unit 24 displays a dialog box BPi for requesting a distance and an angle (S807). In the dialog box BPi, for example, the fault length L = 3800 m obtained from the calculation formula is input. Further, for example, “154 °” is input in which the direction of the fault (traveling N64 ° W) is converted into an angle that is counterclockwise with 0 ° as the east.

  Next, the γ-ray measurement point / tomographic line drawing unit 24 reads the distance (3800 m) and the angle (154 °) input to the dialog box BPi and passes them to the line tool (S809). That is, when the γ-ray measurement point Pi is defined on the hybrid plane of the memory 38, as shown in FIG. 21, the γ-ray measurement point / tomographic line drawing unit 24 sets the line tool and dialog box BPi to the γ-ray measurement point Pi. Is associated with the distance (3800 m) and the angle (154 °).

  The above-mentioned line tool reads the distance (3800 m) and the angle (154 °), and the γ-ray field measurement result straight line LGj at the distance (3800 m) at the angle (154 °) with the γ-ray measurement point Pi as the origin (reference). Is drawn on the hybrid plane of the memory 38 (S811), and the γ-ray field measurement result straight line LGj (black dotted line) is superimposed and displayed on the red three-dimensional image RGi by the imaging unit 14 as shown in FIG. 22 (S813). ).

  As shown in FIG. 22, since the γ-ray measurement point Pi is the origin, a γ-ray field measurement result line LGj (yellow) is drawn only in the west direction.

  Then, as shown in FIG. 18, the γ-ray measurement point / tomographic line drawing unit 24 determines whether or not the γ-ray field measurement result straight line LGj (black dotted line) is adjusted (S815).

  Specifically, when the γ-ray local measurement result line LGj is designated (clicked) with the cursor, it is determined that the γ-ray local measurement result line LGj is adjusted.

  If it is determined in step S815 that the γ-ray local measurement result straight line LGj is adjusted, the adjustment of the γ-ray local measurement result straight line LGj is accepted, and as shown in FIG. The local measurement result straight line LGj is moved (S817).

Specifically, the γ-ray field measurement result straight line LGj is moved through the γ-ray measurement point Pi until the click is made while maintaining the angle (154 °). That is, the γ-ray local measurement result line LGj is moved through the hybrid plane γ-ray measurement point Pi of the memory 38 until the click is made while maintaining the angle (154 °), and this is moved by the imaging unit 14. Display.
That is, when the fault plane is inclined, the line where the fault plane and the ground surface intersect = the fault line does not have a complete straight line shape. For this reason, “faults based on existing data” and “topography considered to be faults by interpretation” passing through the vicinity of the line segment and having good continuity are judged to be a series of faults, and the traced lines are drawn.

  Accordingly, as shown in FIG. 23, the γ-ray field measurement result straight line LGj in FIG. 22 is moved in parallel at an angle (154 °) through the γ-ray measurement point Pi. In FIG. 23, the fact that the entire 3800 m γ-ray field measurement result straight line LGj was not visible in the image RGi corresponds to all the fine tomographic lines DLj estimated to pass the γ-ray measurement point Pi as shown in FIG. Therefore, it is possible to easily grasp how the fine fault lines DLj are connected.

  Then, as shown in FIG. 24, the γ-ray field measurement result straight line LGj and the fine fault line DLj (shown by the dotted line: blue) are the fault lines that connect the respective fine fault lines DLj while visually confirming the relationship. If the operator determines that they are connected, the operator issues an instruction to create a hybrid fault line HBLj (pink) and connects the fine fault line DLj as shown in FIG. 24 (thick line: color is pink) (S819). ). That is, a thick line is defined on the hybrid plane of the memory 38.

  Then, as shown in FIG. 25, the γ-ray field measurement result straight line LGj is deleted (S821), and it is determined whether or not the process is finished (S823). If it is determined that the processing has ended, the imaging unit 14 stores the screen image in the memory 39 as a hybrid tomogram HHGi (S825).

  Therefore, since fault continuity has a great influence on activity evaluation, the fault continuity, which was unknown only by interpretation from the red three-dimensional map, should be supplemented by the γ-ray field measurement result straight line LGj from γ-ray measurement. Is done.

  In addition, it is possible to grasp the distribution of fault fracture zones that had not been seen so far and their characteristics (width, length, etc.) over a wide area. Based on the information obtained in this way, it is also possible to evaluate the activity of active faults (evaluation based on fault fault zone width / length and geological age).

  Therefore, it can be used for risk assessment of earthquakes and landslide disasters on infrastructure such as railways and roads, and prediction of occurrence and impact on landslide disasters such as deep collapse.

<Embodiment 2>
FIG. 26 is a schematic configuration diagram of the tomographic drawing creating apparatus according to the second embodiment. In FIG. 26, the description of the same reference numerals as those in FIG. 1 is omitted.

  As shown in FIG. 26, the tomographic map creating apparatus of the second embodiment stores a database 41 in which orthophoto images in the same range as the red three-dimensional image RGi are stored, and a contour map in the same range as the red three-dimensional image RGi. And a switching unit 44.

  When the orthophoto or contour map is designated by the operator, the switching unit 44 reads out the orthophoto image in the database 41 or the contour map in the database 42 and displays it via the imaging unit 14.

  For example, when the hybrid tomogram is displayed according to the first embodiment as shown in FIG. 27A, when the operator inputs an orthophoto, the switching unit 44 images the orthophoto image in the database 41. To the unit 14. When the orthophoto image is output from the switching unit 44, the imaging unit 14 stops reading the red three-dimensional image RGi from the memory 32 and replaces the red three-dimensional image RGi with the orthophoto image from the switching unit 44. Is displayed as shown in FIG.

  When the operator inputs a contour map when the hybrid tomogram is displayed according to the first embodiment as shown in FIG. 27A, the switching unit 44 converts the contour map of the database 42 into an imaging unit. 14 for output. When the contour map is output from the switching unit 44, the imaging unit 14 stops reading the red three-dimensional image RGi from the memory 32, and displays the contour map from the switching unit 44 instead of the red three-dimensional image RGi. It is displayed as shown in 27 (c).

  Therefore, the hybrid tomogram can be confirmed from the orthophoto image or the contour map.

<Red three-dimensional image>
A method for generating the red three-dimensional image RGi will be described below.

  In generating the red three-dimensional image RGi, the concept of opening is used. The degree of opening is a quantification of the extent that the point is protruding above the ground and the depth of penetration into the basement. That is, as shown in FIG. 28 (a), the ground opening represents the size of the sky that can be seen within the distance LR from the sample point of interest, and the underground opening is shown in FIG. 28 (b). As described above, when looking over the ground with a handstand, it represents the size of the underground in the range of the distance LR.

  The opening degree depends on the distance LR and the surrounding terrain. FIG. 29 shows the ground opening and underground opening for nine types of basic topography in an octagonal graph of the ground angle and underground angle for each direction.

  In general, the ground opening increases as the point protrudes higher from the surroundings, and takes a larger value at the summit and ridge and smaller at the depression and valley bottom. On the contrary, the underground opening becomes larger as the depth of penetration into the basement is greater, with a large value at the depression and valley bottom, and a small value at the summit and ridge. Actually, since various basic landforms are mixed even within the range of the distance LR, the octagonal graph of the ground angle and the underground angle is often deformed and the opening degree also takes various values.

  In addition, the opening degree map can extract information suitable for the size of the terrain by specifying the calculation distance, and can display without depending on directionality and local noise.

  In other words, it is excellent in extracting ridge lines and valley lines, and abundant terrain / geological information can be read. As shown in FIG. 30, on a certain range of DEM data (the ground surface: solid: ( In a)), an angle vector θi formed by a horizontal line and a straight line L1 connecting the point B that is the maximum vertex when one of the eight directions from the set point A is viewed is obtained.

  This angle vector is calculated over 8 directions, and the average of these is called the ground opening, and a solid (Fig. 3) with an air layer pressed onto a certain range of DEM data (ground surface: solid). The point C (corresponding to the deepest point) when viewing either one of the eight directions from the point A in the inverted DEM data (FIG. 30C) with the (30) (b) turned upside down. The angle formed by the straight line L2 connecting the horizontal line) and the horizontal line. Obtaining and averaging this angle over eight directions is referred to as underground opening ψi.

  That is, the ground opening is generated by generating a topographic cross section for each of the eight directions on the DEM data included in the range from the point of interest to a certain distance, and connecting each point and the point of interest (see (a) in FIG. 30). The maximum value of the slope of L1) (when viewed from the vertical direction) is obtained. Such processing is performed for eight directions.

  The angle of inclination is the angle from the zenith (90 degrees for flats, 90 degrees or more for ridges and peaks, and 90 degrees or less for valleys and depressions).

  Further, in the range from the point of interest of the inverted DEM data to a certain distance, a topographic section is generated for each of the eight directions, and the maximum value of the slope of the line connecting each point and the point of interest (the ground surface in FIG. 30A). When the L2 is viewed from the vertical direction, the minimum value) is obtained. Such processing is performed for eight directions.

  In the three-dimensional view of the ground surface in FIG. 30 (a), the angle when viewing L2 from the vertical direction (underground opening ψi) is 90 degrees if flat, 90 degrees or less at the ridge or mountain peak, and 90 degrees or more at the valley or depression. It is.

That is, the ground opening and the underground opening are considered as two sample points A (iA, jA, HA) and B (iB, jB, HB) as shown in FIG. Since the sample interval is 1 m, the distance between A and B is
P = {(iA − iB) 2 + (jA − jB) 2} 1/2… (1)
It becomes.

(A) of FIG. 31 shows the relationship between A and B of the sample points with an altitude of 0 m as a reference. The elevation angle θ of the sample point A with respect to the sample point B is θ = tan −1 {(H B −HA) / P
Given in. The sign of θ is positive when (1) HA <HB, and negative when (2) HA> HB.

A set of sample points within the distance LR in the direction D from the sample point of interest is described as DSL, and this is referred to as a “DL set of sample points of interest”. here,
DβL: Maximum value of the elevation angle for each element of DSL at the sample point of interest
DδL: The following definition is made as the minimum value of the elevation angles for each element of DSL at the sample point of interest (see FIG. 31 (b)).

Definition 1: The ground angle and underground angle of the D-L set of sample points of interest are each
DφL = 90− DβL
as well as
DψL = 90 + DδL
Means.

    DφL means the maximum value of the zenith angle at which the sky in the direction D can be seen within the distance LR from the sample point of interest. Generally speaking, the horizon angle corresponds to the ground angle when L is infinite. Further, DψL means the maximum value of the nadir angle at which the ground in the direction D can be seen within the distance LR from the sample point of interest. When L is increased, the number of sample points belonging to DSL increases, so that DβL has a non-decreasing characteristic, while DδL has a non-increasing characteristic. Therefore, both DφL and Dψ1 have non-increasing characteristics with respect to L.

  High angle in surveying is a concept defined with reference to a horizontal plane passing through the sample point of interest, and does not exactly match θ. In addition, if the ground angle and the underground angle are to be strictly discussed, the curvature of the earth must be taken into consideration, and definition 1 is not necessarily an accurate description. Definition 1 is a concept defined on the assumption that terrain analysis is performed using DEM.

  The ground angle and the underground angle are the concepts of the specified direction D, but the following definition is introduced as an extension of this.

Definition II: The ground opening and underground opening of the distance LR of the sample point of interest are ΦL = (0φL + 45φL + 90φL + 135φL + 180φL + 225φL + 270φL + 315φL) / 8
And ψL = (0ψL + 45ψL + 90ψL + 135ψL + 180ψL + 225ψL + 270ψL + 315ψL) / 8
Means.

  The ground opening represents the size of the sky that can be seen within the range of the distance LR from the sample point of interest, and the underground opening represents the size of the underground in the range of the distance LR when looking up underground with a handstand. (See FIG. 30).

  Then, the DEM data is meshed into a square, and the average slope of the square surface adjacent to the target point on the mesh is obtained. There are four adjacent squares, and any one of them is a square of interest. Then, the altitude and average inclination of the four corners of this square of interest are obtained. The average slope is the slope of the surface approximated from 4 points using the least square method.

  And by providing the first gray scale for expressing the ridge and valley bottom with brightness by emphasizing the convex part, the ground opening (average angle when viewing the range of L from the point of interest in 8 directions: high place Each time an index for determining whether or not the brightness (lightness) corresponding to the value of the ground opening θi is calculated.

  For example, when the value of the ground opening is within a range of about 40 degrees to 120 degrees, 50 degrees to 110 degrees is associated with the first gray scale and is assigned to 255 gradations.

  That is, the ridge portion (convex portion) has a larger value of the ground opening, so the color becomes white.

  Then, in a mesh area having a point of interest (coordinates) (when contour lines connecting the same Z values of DEM data are meshed with a square (for example, 1 m), and any of the four corners of this mesh is the point of interest) The color data based on the first gray scale is allocated and stored. Next, the ground opening layer obtained by inverting the color gradation of the stored ground opening image data by gradation correction is stored. That is, the ground opening layer adjusted so that the ridge is white is obtained.

  Next, a second gray scale for expressing the valley bottom and ridge with brightness is provided, and the brightness corresponding to the value of the ground opening is obtained every time the underground opening ψi (average in eight directions from the point of interest) is obtained. Is calculated.

  For example, when the value of the underground opening is within a range of about 40 degrees to 120 degrees, 50 degrees to 110 degrees is associated with the second gray scale and is allocated in 255 gradations. That is, since the value of the underground opening is larger in the valley bottom portion (concave portion), the color becomes black.

  Then, read the underground opening image data, mesh the mesh area with the point of interest (coordinates) (contour lines connecting the same Z value of the DEM data with a square (for example, 1 m), and select one of the four corners of this mesh. ) Is assigned, and the color data based on the second gray scale is allocated and stored. Next, the color tone of the underground opening image data is corrected by the tone correction process.

  If the color becomes too black, the tone curve is corrected. This is called the underground opening layer and saved.

  Next, a third gray scale for expressing the degree of inclination according to brightness is provided, and every time the inclination (average in four directions from the point of interest) is obtained, the value of this inclination is set. The brightness (lightness) of the corresponding third gray scale is calculated.

  For example, when the value of the inclination αi falls within the range of about 0 to 70 degrees, 0 degree to 50 degrees corresponds to the third gray scale and is assigned to 255 gradation. That is, 0 degrees is white and 50 degrees or more is black. The color becomes darker as the slope α (also referred to as slope αi) is larger.

  Then, a difference image between the underground opening image data and the ground opening image data is stored as an inclination image.

  At this time, a mesh region having a point of interest (coordinates) (when the contour line connecting the same Z values of DEM data is meshed with a square (for example, 1 m), and any of the four corners of this mesh is used as the point of interest) Are assigned color data based on the third gray scale. Next, red processing enhances R with the RGB color mode function. That is, a tilt-emphasized image in which red is emphasized as the tilt is large is obtained.

  Then, a composite image obtained by multiplying the ground opening layer and the underground opening layer is obtained. At this time, both balances are adjusted so that the valley portion is not crushed.

The above-mentioned “multiplication” is a term of a layer mode on a photoshop, and is an OR operation in numerical processing.

  In this balance adjustment, the distribution of the values of the ground opening and the underground opening is cut out from a ground surface having a certain radius (L / 2) around a certain point.

  When the entire sky has a uniform brightness, the size of the sky looking up from the ground surface gives the brightness of the ground.

  That is, the ground opening is bright. However, the value of the underground opening should also be taken into account when the light goes around.

  Depending on the ratio between the two, the ridge of the terrain can be emphasized or changed arbitrarily. If you want to emphasize the topography in the valley, increase the value of b.

Brightness index = a × ground opening−b × underground opening where a + b = 1
That is, a composite image of gray gradation expression obtained by multiplying the ground opening layer (ridge is emphasized white) and the underground opening layer (bottom is emphasized black) is obtained.

  On the other hand, a red three-dimensional map RGi in which the ridge combined with the tilt-emphasized image and the composite image is emphasized in red is obtained and displayed on the display unit 200.

  That is, it is possible to generate a three-dimensional image in which the ridges and mountain peaks are whitish, the valleys and depressions are dark, and the steeper slopes are red. For this reason, it is possible to grasp the degree of unevenness and the degree of inclination at a glance.

  In other words, for each direction from the DEM data of the point of interest on the ground surface to which the DEM data is connected, the angle vector formed by the DEM data that is the maximum vertex within a certain range and the horizontal line is obtained and averaged. Degree is obtained, and a ground opening degree image is obtained in which a bright color is assigned as much as the value of the ground opening degree.

  Then, from the DEM data of the point of interest of the inverted DEM data obtained by reversing the solid that presses the air layer on the DEM data in a certain range, the DEM data that becomes the maximum vertex up to the certain range and the horizontal line are formed for each of a plurality of directions. Each angle vector is obtained and averaged to obtain the underground opening, and an underground opening image in which a darker color is assigned to the magnitude of the value of the underground opening is obtained.

  Then, the ground opening image and the underground opening image are weighted and synthesized, and a synthesized image expressed in gradation according to this value is obtained as a red three-dimensional image RGi.

<Embodiment 3>
FIG. 32 is a schematic configuration diagram of the tomographic drawing creating apparatus according to the third embodiment. In FIG. 32, the same reference numerals as those in FIG.
As shown in FIG. 32, the tomographic map creating apparatus according to the third embodiment includes a mode determination unit 50, a memory definition unit 52, an existing / red interpretation tomographic line drawing unit 54, a layer order setting unit 56, and a γ-ray. A program configuration including a measurement point defining unit 58, an estimated tomographic line drawing unit 60, a hybrid tomographic line drawing unit 62, an image composition output unit 64, and the like is provided.

  The imaging unit 14 includes an image memory (not shown), and synthesizes (superimposes) various images stored in the image memory and displays the images on the display unit 200. At this time, color information, line information, etc. added to the image are displayed on the screen.

The mode determination unit 50 displays a dialog box (not shown) for inputting a screen to be displayed on the screen, and activates each corresponding unit based on the information input in the dialog box.
This dialog box includes, for example, a button for displaying the red three-dimensional map RGi in the memory 32, a button for displaying the existing fault line KLj in the existing fault line memory 36, and a fine fault line DLj in the fine fault line definition memory 34. A button for displaying and a button for displaying the hybrid tomographic line HBLj are provided.

In addition, a button for starting a mode for defining a γ-ray measurement point Pi, a button for starting a mode for creating a γ-ray field measurement result straight line LGj, and a mode for moving a γ-ray field measurement result straight line LGj. A button, a button for starting a mode for creating a hybrid tomographic line HBLj, and the like are provided.
Furthermore, a button for deleting an existing fault line KLj in an existing fault line working memory 71 described later, a button for deleting a fine fault line DLj in a fine fault line working memory 73 described later, and the like are provided.

The memory definition unit 52 reads the origin of the red three-dimensional image RGi upon activation, the fine tomographic line working memory 73 having this origin, the existing tomographic line working memory 71, and the γ-ray local measurement point defining working memory 75. The γ-ray field measurement result straight line working memory 76 and the hybrid tomographic line working memory 78 are generated in the working memory unit 70, and each defines the XY plane having the above-mentioned origin.
That is, the memory definition unit 52 defines the origin coordinates of the red three-dimensional image RGi in the range definition unit 16 on each plane.

  The plane defined in the fine fault line working memory 73 is referred to as a red interpretation result fine fault line working plane Ha, the plane defined in the existing fault line working memory 71 is referred to as an existing fault line working plane Hb, and γ The plane defined in the line local measurement point definition working memory 75 is referred to as a γ-ray field measurement point definition working plane Hc, and the plane defined in the γ-ray field measurement result straight line working memory 76 is defined as the γ-ray field measurement result straight line. The drawing working plane Hd is referred to as a plane defined in the hybrid tomographic line working memory 78 as a hybrid tomographic drawing working plane He.

The existing red-reading tomographic line drawing unit 54 displays the red three-dimensional image RGi stored in the fine tomographic line definition memory 34 on the screen as the γ-ray fine tomographic line working memory 73 is generated, and displays the tomographic image. The fine tomographic lines DLj (DL1, DL2,...) Determined to be copied are copied (copied) to the fine tomographic line working memory 73.
Further, the existing / red-reading tomographic line drawing unit 54 generates the existing fault line KLj (drawn on the existing plane stored in the existing fault line memory 36 as the existing fault line working memory 71 is generated. KL1, KL2) are copied to the existing tomographic line working memory 71.

  The layer order setting unit 56 sets the order of image layers (layers) to be displayed on the screen based on the input order in a memory (not shown), and causes the imaging unit 14 to display the images in the order of the layers.

  As shown in FIG. 33, in this hierarchy, for example, the first layer (topmost layer) is set to a working memory 75 for defining a gamma ray local measurement point (a plane Hc for defining a gamma ray local measurement point) (specifically, the memory) The second layer is made to be a γ-ray field measurement result straight line working memory 76 (γ-field field measurement result straightness drawing working plane Hd) (specifically, a memory address is preferably associated). The third layer is used as a hybrid tomography working memory 78 (hybrid tomography drawing working plane He) (specifically, it is preferable to associate a memory address), and the fourth layer is used for γ-ray fine tomography working. Memory 73 (red interpretation result fine tomographic line working plane Ha) (specifically, it is preferable to associate the memory address), the fifth layer The existing fault line working memory 71 (existing fault line working plane Hb) is used (specifically, a memory address is preferable).

The sixth layer is a red three-dimensional image RGi (32) (specifically, a memory address is preferable).
Upon activation, the γ-ray measurement point definition unit 58 displays a screen (not shown) for inputting the color and symbol type (for example, ◯) of the γ-ray measurement point Pi, and the identification number ( (Hereinafter referred to as γ-ray measurement point number PBi), and these are associated and stored in the memory 80 (see FIG. 34).

  Then, the γ-ray measurement point definition unit 58 reads the memory address of the γ-ray local measurement point definition working memory 75 into the layer order setting unit 56 and allocates the memory at the memory address from the working memory unit 70. The γ-ray local measurement point definition working plane Hc defined in the γ-ray local measurement point definition working memory 75 is output to the imaging unit 14 to display a screen (transparent when the γ-ray measurement point Pi is not defined). ).

Then, when the cursor is positioned on the screen and is clicked, for example, it is determined that the γ-ray measurement point Pi is determined, and the position of the cursor on the screen is determined as the coordinates of the γ-ray local measurement point definition working plane Hc. And is defined (calculated) as on-site measurement point coordinates PZi.
That is, the γ-ray measurement point definition unit 58 reads the point designated by the cursor in the red three-dimensional image RGi in the γ-ray measurement point / tomographic line drawing unit 24 as the γ-ray measurement point Pi, and this γ-ray measurement point Pi is read. The function of defining the gamma ray local measurement point definition working plane Hc is performed.
At this time, the color, symbol type (for example, ◯), γ-ray measurement point number PBi, etc., input to the local measurement point coordinates PZi are assigned (see FIG. 35).

Then, the γ-ray measurement point definition unit 58 informs the imaging unit 14 of the local measurement point coordinates PZi, the color (for example, yellow) of the γ-ray measurement point Pi defined on the γ-ray local measurement point definition working plane Hc, and the size. , Symbol types (for example, ◯) etc. are output as a set.
Accordingly, as shown in FIG. 19, the γ-ray measurement point Pi is displayed on the red three-dimensional image RGi on the screen. FIG. 19 shows an example in which the existing fault line KLj and the fine fault line DLj are already displayed. Further, in this example, the working memory 75 for defining γ-ray local measurement points is set at the highest level.

At this time, the field measurement point coordinates PZi are stored in the memory 80 in association with the γ-ray measurement point number PBi (see FIG. 34).
The estimated tomographic line drawing unit 60 displays a screen (not shown) for inputting a color, a line type (thickness), and the like for drawing the γ-ray field measurement result straight line LGj upon activation, and displays the γ-ray field A line number LGBi (including a unique code: year, month, day, time, and coordinates) of the measurement result straight line LGj is generated and stored in the memory 80 in association with the γ-ray measurement point number PBi (see FIG. 34).

The estimated tomographic line drawing unit 60 then stores the local measurement point coordinates PZi associated with the γ-ray measurement point number PBi in the γ-ray local measurement result straight line working memory 76 (γ-ray local measurement result straight line drawing working plane Hd). Define (see FIG. 35).
Then, the local γ-ray measurement information Fγi having the γ-ray measurement point number PBi or the line number LGBi of the γ-ray measurement point Pi stored in the memory 80 is read from the local measurement result information memory 40, and this local γ-ray measurement information Fγi is read. The direction Lθi, the fault length L, the width, and the like included in are read.

Then, the estimated tomographic line drawing unit 60 reads the memory address of the γ-ray local measurement result straight line working memory 76 (γ-ray local measurement result straight line drawing working plane Hd) set in the layer order setting unit 56. The memory at the memory address is allocated from the working memory unit 70.
Then, the local measurement point coordinate PZi of the γ-ray measurement point Pi on the γ-ray local measurement result straight line drawing working plane Hd of the γ-ray local measurement result straight line working memory 76 is used as a starting point, the direction Lθi, the fault length L (for example, 3800 m ) Draw (calculate) a γ-ray field measurement result line LGj having a width W (see FIG. 35). At this time, the line number LGBi, the direction Lθi, the fault length L (for example, 3800 m), the width W, the line type, the color, the local measurement point coordinates PZi, and the like are associated with the γ-ray local measurement result straight line LGj.

  The estimated tomographic line drawing unit 60 then transmits the direction Lθi, the fault length L (for example, 3800 m) and the width associated with the γ-ray local measurement result straight line LGj drawn on the γ-ray local measurement result straight line drawing working plane Hd. W, line type, color, local measurement point coordinates PZi, and the like are output to the imaging unit 14 as line attribute information ZLi. Therefore, it is displayed as shown in FIG.

  Further, in the mode in which the γ-ray local measurement result straight line LGj is moved, the estimated tomographic line drawing unit 60 reads the memory address where the γ-ray local measurement result straight line LGj is defined from the layer order setting unit 56. That is, the memory address of the γ-ray local measurement result straight line working memory 76 (γ-ray local measurement result straight line drawing working plane Hd) is read. Then, the line attribute information ZLi that is related to the γ-ray field measurement result straight line LGj on the γ-ray field measurement result line drawing working plane Hd is read.

Then, the estimated tomographic line drawing unit 60 moves the γ-ray field measurement result line LGj in the direction opposite to the direction Lθi of the line attribute information ZLi read in accordance with the movement of the cursor (see FIG. 36). At this time, the γ-ray local measurement result straight line LGj of the γ-ray local measurement result straight line drawing working plane Hd is output to the imaging unit 14 at regular intervals.
That is, the estimated tomographic line drawing unit 60 displays the γ-ray local measurement result line LGj in the γ-ray measurement point / tomographic line drawing unit 24, and moves the γ-ray local measurement result line LGj in the direction Lθi as the cursor moves. The function to move in the opposite direction is performed.

  Therefore, as shown in FIG. 23, the γ-ray field measurement result straight line LGj is displayed.

  The hybrid tomographic line drawing unit 62 displays a screen (not shown) for inputting a color (for example, pink), a line thickness, and the like for drawing the hybrid tomographic line HBLj with the activation. Then, an identification number of the hybrid fault line HBLj (hereinafter referred to as hybrid fault line number HBbi) is generated, and these are associated and stored in the memory 90 (see FIG. 37). Then, the memory address where the hybrid tomographic line HBLj is defined in the layer order setting unit 56 is read. That is, the hybrid tomographic line working memory 78 (hybrid tomographic line drawing working plane He) is allocated.

Then, this hybrid tomographic drawing working plane He is displayed by the imaging unit 14 (transparent).
Then, along with the movement of the cursor on the screen (on the fine tomographic line DLj), the coordinates of this cursor are converted into the coordinates of the hybrid tomographic line drawing working plane He and defined (calculated) as the hybrid tomographic line drawing working plane He. Then, the coordinates, color (pink), thickness, and the like are output to the imaging unit 14 (see FIG. 38).
That is, the hybrid tomographic line drawing unit 62 performs a function of displaying as the hybrid tomographic line HBLj (pink) in the line adjusting unit 26.

  Therefore, as shown in FIG. 24, the fine tomographic lines DLj are connected by lines while being traced.

  The image composition output unit 64 adds the red interpretation result fine fault line working plane Ha of the γ-ray fine tomography working memory 73 and the existing fault line working of the existing fault line working memory 71 to the red three-dimensional map RGi of the memory 32. Plane Hb, γ-ray field measurement point definition working plane Hc of the γ-ray field measurement point definition working memory 75, γ-ray field measurement result line-working memory 76 γ-ray field measurement result line drawing work plane Hd, hybrid fault Any or all of the images of the hybrid tomographic line drawing working plane He in the line working memory 78 are superimposed and output to a printer (not shown) for printing.

(Hybrid tomogram HHGi)
Hybrid tomogram HHGi is a selection of a button for deleting an existing fault line KLj, a button for deleting a fine fault line DLj, a button for deleting a γ-ray measurement point Pi, and a button for deleting a γ-ray field measurement result straight line LGj Get by.

  That is, in the γ-ray measurement point Pi deletion mode, the γ-ray measurement point definition unit 58 allocates the γ-ray local measurement point definition working memory 75 and designates the γ-ray measurement point Pi on the screen. The γ-ray measurement point Pi defined on the local measurement point definition working plane Hc is deleted.

  The existing / red-reading tomographic line drawing unit 54 allocates the existing tomographic line working memory 71 in the deletion mode of the existing tomographic line KLj. Then, as the existing fault line KLj on the screen is traced with the cursor, the existing fault line KLj on the existing fault line working plane Hb of the existing fault line working memory 71 corresponding to this coordinate is deleted.

  Further, the existing / red interpretation tomographic line drawing unit 54 allocates the γ-ray fine tomographic line working memory 73 in the deletion mode of the fine tomographic line DLj. Then, as the fine tomographic line DLj on the screen is traced with the cursor, the fine tomographic line DLj on the fine tomographic line working plane Ha in the γ-ray fine tomographic line working memory 73 corresponding to this coordinate is deleted. To do.

Further, the estimated tomographic line drawing unit 60 allocates the γ-ray local measurement result straight line working memory 76 in the deletion mode of the γ-ray local measurement result straight line LGj. Then, as the γ-ray local measurement result line LGj on the screen is traced with the cursor, the γ-ray local measurement result line LGj in the γ-ray local measurement result line working memory 76 corresponding to this coordinate is deleted.
The hybrid tomogram HHGi shown in FIG. 25 is obtained by such processing.

  Further, in the layer order setting unit 56, as shown in FIG. 39, when the first layer is the γ-ray local measurement result straight line working memory 76 (γ-ray local measurement result straight line drawing working plane Hd), FIG. As shown, the γ-ray field measurement result line LGj passes on the γ-ray measurement point Pi.

  In the above-described embodiment, the hybrid tomographic line is displayed in an overlapping manner on the orthophoto image and the contour map.

Further, in each of the above-described embodiments, a fine fault has been described, but a fault having a width of 15 m, 20 m,... May be used.
Moreover, although the color of each line is arbitrary, it is preferable that each is a different color, thickness, or line type.

DESCRIPTION OF SYMBOLS 10 Red map creation part 12 Fault search range setting part 14 Imaging part 16 Range definition part 18 Existing and red interpretation drawing part 24 γ-ray measurement point / fault line drawing part 26 Line adjustment part

The tomographic map creating apparatus according to the present invention includes a first storage means in which a digital elevation model obtained by emitting laser pulses from an aircraft to a region at intervals capable of detecting fine tomographic lines DLj,
A second storage means for storing the red three-dimensional image RGi;
Fault length L and the γ-ray of the fine fault line DLj determined on the basis of the measurement interval and γ-ray intensities of a plurality of measurement points in the γ-ray measurement point Pi in the area determined based on the red three-dimensional image RGi Third storage means in which the direction of the fine tomographic line DLj at the measurement point Pi is stored as the local γ-ray measurement information Fγi;

A tomographic detection range setting unit for displaying an input box for setting a search range in which the tomographic zone width FW of the fine tomographic line DLj can be detected on the screen of the display unit;
The red three-dimensional image RGi based on the numerical elevation model of the first storage unit and the search range set by the tomographic detection range setting unit is generated and stored in the second storage unit and displayed on the screen. A red three-dimensional image creation unit to be displayed;
A surface for drawing the fine tomographic line DLj is displayed on the red three-dimensional image RGi on the screen, and the movement trajectory of the cursor on this surface is defined as the fine tomographic line DLj, which is displayed in a predetermined color. Means,
On the surface for drawing the fine tomographic line DLj, a surface for designating the γ-ray measurement point Pi is displayed, and the point designated on this surface is used as the γ-ray measurement point Pi for the fine tomographic line. Means for displaying in a color different from DLj;
A surface for displaying the γ-ray field measurement result straight line LGj is displayed above or below the surface for drawing the fine fault line DLj, and is defined as a surface for designating the γ-ray measurement point Pi on this surface. The γ-ray measurement point Pi is defined, and a straight line of the direction and the fault length L of the local γ-ray measurement information Fγi of the third storage means is drawn from the γ-ray measurement point Pi, means for displaying the γ-ray field measurement result line LGj in a color different from the γ-ray measurement point Pi and the fine tomographic line DLj;
Means for assigning the direction to the γ-ray field measurement result straight line LGj;
When the γ-ray field measurement result line LGj is designated by a cursor, the direction assigned to the γ-ray field measurement result line LGj is read, and the γ-ray field measurement result is made to follow the movement of the cursor. And a means for moving the straight line LGj in the read direction.

Claims (21)

First storage means for storing a digital elevation model obtained by emitting laser pulses from an aircraft to an area at intervals capable of detecting fine tomographic lines DLj;
A second storage means for storing the red three-dimensional image RGi;
The fault length L of the fine fault line DLj and the gamma ray measurement obtained based on the gamma ray intensity and the measurement interval of a plurality of measurement points at the gamma ray measurement point Pi in the area determined based on the red three-dimensional image RGi. Second storage means in which the direction of the fine fault line DLj at the point Pi is stored as local γ-ray measurement information Fγi, and an input for setting a search range in which the fault zone width FW of the fine fault line DLj can be detected A tomographic detection range setting unit for displaying a box on the screen of the display unit;
The red three-dimensional image RGi based on the numerical elevation model of the first storage unit and the search range set by the tomographic detection range setting unit is generated and stored in the second storage unit and displayed on the screen. A red three-dimensional image creation unit to be displayed;
A surface for drawing the fine tomographic line DLj is displayed on the red three-dimensional image RGi on the screen, and the movement trajectory of the cursor on this surface is defined as the fine tomographic line DLj, which is displayed in a predetermined color. Means,
On the surface for drawing the fine tomographic line DLj, a surface for designating the γ-ray measurement point Pi is displayed, and the point designated on this surface is used as the γ-ray measurement point Pi for the fine tomographic line. Means for displaying in a color different from DLj;
A surface for displaying the γ-ray field measurement result straight line LGj is displayed above or below the surface for drawing the fine fault line DLj, and is defined as a surface for designating the γ-ray measurement point Pi on this surface. The γ-ray measurement point Pi is defined, and a straight line of the direction and the fault length L of the local γ-ray measurement information Fγi of the second storage means is drawn from the γ-ray measurement point Pi, means for displaying the γ-ray field measurement result line LGj in a color different from the γ-ray measurement point Pi and the fine tomographic line DLj;
Means for assigning the direction to the γ-ray field measurement result straight line LGj;
When the γ-ray field measurement result line LGj is designated by a cursor, the direction assigned to the γ-ray field measurement result line LGj is read, and the γ-ray field measurement result is made to follow the movement of the cursor. Means for moving the straight line LGj in the read direction;
A tomographic map creating apparatus characterized by comprising: A surface for connecting the fine tomographic lines DLj with lines is displayed on the surface for drawing the fine tomographic lines DLj, and the fine tomographic lines DLj and γ are followed by movement of the cursor on the surface. A line measuring point Pi, a line having a color different from the color of the γ-ray field measurement result line LGj as a hybrid tomographic line HBLj, and a means for drawing on a surface for connecting the fine tomographic lines DLj with lines;
The tomographic map creating apparatus according to claim 1, wherein: The fault length L is
L = 0.36FW + 0.62
The tomographic map creating apparatus according to claim 1, wherein:   The tomographic map creating apparatus according to claim 1, wherein a width of the fine tomographic line DLj is 10 m or less. The digital elevation model is
The tomographic map creating apparatus according to claim 1, wherein the mesh size is 2.5 m or less. The red three-dimensional image creation unit
Set a sample point on the grid of the digital elevation model, find the ground opening, underground opening and slope to the target grid existing within the search range set from this sample point for each setting, The red three-dimensional image RGi is assigned by assigning a brighter color as the ground opening value is larger, a darker color as the underground opening value is larger, and a red-enhanced color as the slope value is larger. The tomographic map creating device according to claim 1, wherein the tomographic map creating device is generated. The search range is
The tomographic map creating apparatus according to claim 1, wherein the tomographic map creating apparatus has a length of about 5 to 10 times the tomographic zone width FW.   The tomographic map creating apparatus according to claim 7, wherein the search range is about 50 to 100 m. First storage means for storing a digital elevation model obtained by emitting laser pulses from an aircraft to an area at intervals capable of detecting fine tomographic lines DLj;
A second storage means for storing the red three-dimensional image RGi;
The fault length L of the fine fault line DLj and the gamma ray measurement obtained based on the gamma ray intensity and the measurement interval of a plurality of measurement points at the gamma ray measurement point Pi in the area determined based on the red three-dimensional image RGi. A second storage means in which the direction of the fine tomographic line DLj at the point Pi is stored as the local γ-ray measurement information Fγi,
Computer
A tomographic detection range setting unit for displaying an input box for setting a search range in which the tomographic zone width FW of the fine tomographic line DLj can be detected on the screen of the display unit;
The red three-dimensional image RGi based on the numerical elevation model of the first storage unit and the search range set by the tomographic detection range setting unit is generated and stored in the second storage unit and displayed on the screen. Creating a red three-dimensional image to be displayed;
A surface for drawing the fine tomographic line DLj is displayed on the red three-dimensional image RGi on the screen, and the movement trajectory of the cursor on this surface is defined as the fine tomographic line DLj, which is displayed in a predetermined color. Steps,
On the surface for drawing the fine tomographic line DLj, a surface for designating the γ-ray measurement point Pi is displayed, and the point designated on this surface is used as the γ-ray measurement point Pi for the fine tomographic line. Displaying in a color different from DLj;
A surface for displaying the γ-ray field measurement result straight line LGj is displayed above or below the surface for drawing the fine fault line DLj, and is defined as a surface for designating the γ-ray measurement point Pi on this surface. The γ-ray measurement point Pi is defined, and a straight line of the direction and the fault length L of the local γ-ray measurement information Fγi of the second storage means is drawn from the γ-ray measurement point Pi, displaying the γ-ray field measurement result straight line LGj in a color different from the γ-ray measurement point Pi and the fine tomographic line DLj;
Assigning the direction to the γ-ray field measurement result straight line LGj;
When the γ-ray field measurement result line LGj is designated by a cursor, the direction assigned to the γ-ray field measurement result line LGj is read, and the γ-ray field measurement result is made to follow the movement of the cursor. Moving a straight line LGj in the read direction;
A method for creating a tomographic map, characterized in that: The computer is
A surface for connecting the fine tomographic lines DLj with lines is displayed on the surface for drawing the fine tomographic lines DLj, and the fine tomographic lines DLj and γ are followed by movement of the cursor on the surface. Drawing a line measuring point Pi and a line having a color different from the color of the γ-ray field measurement result line LGj as a hybrid tomographic line HBLj on a surface for connecting the fine tomographic lines DLj with lines. The tomographic map creating method according to claim 9. The fault length L is
L = 0.36FW + 0.62
The tomographic map creation method according to claim 9, wherein:   The tomographic map creation method according to claim 9, wherein the width of the fine fault line DLj is 10 m or less. The digital elevation model is
The tomographic map creation method according to claim 9, wherein the mesh size is 2.5 m or less. The red three-dimensional image creation step includes
Set a sample point on the grid of the digital elevation model, find the ground opening, underground opening and slope to the target grid existing within the search range set from this sample point for each setting, The red three-dimensional image RGi is assigned by assigning a brighter color as the ground opening value is larger, a darker color as the underground opening value is larger, and a red-enhanced color as the slope value is larger. The tomographic map creating method according to claim 9, wherein the tomographic map creating method is performed. The search range is
The tomographic map creating method according to claim 9, wherein the tomographic map width is about 5 to 10 times the tomographic zone width FW.   The tomographic map creation method according to claim 15, wherein the search range is about 50 to 100 m. First storage means for storing a digital elevation model obtained by emitting laser pulses from an aircraft to an area at intervals capable of detecting fine tomographic lines DLj;
A second storage means for storing the red three-dimensional image RGi;
The fault length L of the fine fault line DLj and the gamma ray measurement obtained based on the gamma ray intensity and the measurement interval of a plurality of measurement points at the gamma ray measurement point Pi in the area determined based on the red three-dimensional image RGi. Using the second storage means (and the direction of the fine tomographic line DLj at the point Pi is stored as the local γ-ray measurement information Fγi,
Computer
A tomographic detection range setting means for displaying an input box for setting a search range in which the fault zone width FW of the fine tomographic line DLj can be detected on the screen of the display unit;
The red three-dimensional image RGi based on the numerical elevation model of the first storage means and the search range set by the tomographic detection range setting means is generated and stored in the second storage means and displayed on the screen. Red three-dimensional image means to display,
A surface for drawing the fine tomographic line DLj is displayed on the red three-dimensional image RGi on the screen, and the movement trajectory of the cursor on this surface is defined as the fine tomographic line DLj, which is displayed in a predetermined color. means,
On the surface for drawing the fine tomographic line DLj, a surface for designating the γ-ray measurement point Pi is displayed, and the point designated on this surface is used as the γ-ray measurement point Pi for the fine tomographic line. Means for displaying in a color different from DLj;
A surface for displaying the γ-ray field measurement result straight line LGj is displayed above or below the surface for drawing the fine fault line DLj, and is defined as a surface for designating the γ-ray measurement point Pi on this surface. The γ-ray measurement point Pi is defined, and a straight line of the direction and the fault length L of the local γ-ray measurement information Fγi of the second storage means is drawn from the γ-ray measurement point Pi, means for displaying the γ-ray field measurement result line LGj in a color different from the γ-ray measurement point Pi and the fine tomographic line DLj;
Means for assigning the direction to the γ-ray field measurement result straight line LGj;
When the γ-ray field measurement result line LGj is designated by a cursor, the direction assigned to the γ-ray field measurement result line LGj is read, and the γ-ray field measurement result is made to follow the movement of the cursor. Means for moving the straight line LGj in the read direction;
Tomographic map creation program to execute the function as. The computer,
A surface for connecting the fine tomographic lines DLj with lines is displayed on the surface for drawing the fine tomographic lines DLj, and the fine tomographic lines DLj and γ are followed by movement of the cursor on the surface. Means for drawing a line having a color different from the color of the line measurement point Pi and the γ-ray field measurement result straight line LGj as a hybrid tomographic line HBLj on a surface for connecting the fine tomographic line DLj with a line;
The tomographic map creation program according to claim 17 for executing the function as. The fault length L is
L = 0.36FW + 0.62
The tomographic map creation program according to claim 17, characterized in that:   The tomographic map creation program according to claim 17, wherein a width of the fine fault line DLj is 10 m or less. The digital elevation model is
The tomographic map creation program according to claim 17, wherein the mesh size is 2.5 m or less.