JP3535597B2 - Graphic processing unit - Google Patents

Description

DETAILED DESCRIPTION OF THE INVENTION [0001] BACKGROUND OF THE INVENTION The present invention relates to a road (block) map
This shows the status of the installation of various facilities such as electricity, water, gas, etc.
A computer that superimposes and manages equipment maps
Data mapping system, especially the alignment between multiple drawings.
The present invention relates to a graphic processing device suitable for obtaining a combination. [0002] 2. Description of the Related Art In recent years, telegraph telephones, electric power, gas, water,
Management of facilities that occupy roads, such as water and subways
Of facilities such as occupied properties and underground
Map database for installation status, etc.
Including the terrain, roads, blocks, structures, names, etc.
New computer by integrating with
・ A mapping system is being constructed. [0003] Each business operator is responsible for underground burial such as gas and water pipes.
When constructing facilities and utility poles, etc.,
In order to avoid damage beforehand, place
It is necessary to understand the relationship. Each company has its own
Between the generated map or map database, surveying,
Same due to errors caused by assumptions such as plotting and numerical data conversion
The location coordinate data of the point is shifted, and the equipment of another
Data of the installation status of
Even if output is superimposed on the
It is not possible to accurately grasp the positional relationship with [0004] The positional relationship between roads and equipment or different types of equipment, etc.
To determine the location of each facility,
Included in figure database (first figure database)
Convert equipment data to a common basic map database
Method shown on map data of (second graphic database)
Has been requested. Basic map database (2nd figure
Database) is an institution with a center system
Are provided in the form of a database. However, the equipment map database (first figure)
Shape database) and basic map database (second figure)
Means to match data in the form of
Is not effectively functioning, so each operator
Extract the necessary data from the shape database and
Print on non-shrinkable transparent paper and overlay these transparent papers
The location of map data and equipment data
I was letting it. In other words, the basic map database (the second figure
Roads, blocks, etc. listed on the map
To match the equipment map database (first figure
Move the equipment listed on the database) map
Trace, and use this as input manuscript paper,
Data had to be entered. This involves non-
Because of the constant labor and expense required, each operator
Participation in building a computer mapping system
A problem of hesitation. Therefore, a plurality of computers are automatically
To match the data contained in the graphic database
Various methods described below have been attempted. (A) JP-A-62-278679,
A method described in JP-A-64-15873, etc.
Is an affine transformation based on the four corner point data
And other parameters. (B) As described in JP-A-1-240983.
What is listed, that is, polygons for houses etc.
And calculate the center of gravity and inclination from the vertex coordinates of
To convert and input using point coordinates. (C) The operator operates the equipment on the display screen.
Correction position designation is now automated without performing for each equipment
As an example, a figure correction device (Japanese Unexamined Patent Application Publication No.
5) ". This "Graphic correction device" has two types of maps
If there is a misalignment between
The area enclosed by the lines indicating the world is automatically extracted and both maps are used.
Between each area and equipment data within the area
And translate the other map as a reference (background)
Is converted to data and corrected.
Excellent in that. (D) "Graphic processing device (Japanese Patent Laid-Open No. 5-614)
07), “graphic processing device (JP-A-5-6292)
6) "," graphic processing device (Japanese Patent Laid-Open No. 5-62927) ",
"Map database verification method (Japanese Unexamined Patent Publication No.
8) "," Figure accuracy test method "
9) "," Figure processing method "
0) "," Graphic processing method (Japanese Patent Laid-Open No. 6-251121) "
Has been filed by the present applicant. In these ways
First compares and overlays the data of the two maps, (d1)
The deviation distance between the corresponding line segments included in each map and
Direction, (d2) of corresponding points included in each map,
(D3) Line segment included in one map
Between the point on the corresponding line segment included in the other map and
Find the distance and direction. Thus, (d1), (d2),
(D4) Any one of the displacement distances and directions is set to a number of sets.
For the affine transformation
I was looking for a meter. Between the corresponding points in the two drawings,
Affine transformation using X- and Y-direction shift values
Including the method of obtaining parameters for
Adjustment method (Japanese Patent Laid-Open No. 5-174127) "
Has been filed. In general, a displacement between two types of maps is used.
The reasons for this are that the topographic survey is flat
Differences in surveying methods such as quantity, reference point survey error,
Measurement errors, trace errors, and input when creating a database
This is due to a force error or the like. Ministry of Construction public survey work regulations
According to the 1/500 scale survey, the total standard error was set to 0.
It is within 5mm, but work due to the same measurement error
At different times, 2-3mm error seems to have occurred
is there. Considering these together, the map database
Is less than 10mm for a 1m map, ie
It may be considered that it is 1/100 or less. Therefore,
Small movement, small angle principle, linear handling is acceptable
Will be done. This is a statistically determined parameter
Implies that affine transformation can be performed. [0012] SUMMARY OF THE INVENTION However, the above-described conventional
Each of the techniques has the following problems. In the method (a), if the surveying is different,
The relative position between the map and the map in the picture frame is shifted, so use
It turned out not to work. The method (b) is extended to the whole drawing and considered.
As a result, post and reference points are used as point coordinates of map data.
There are many things that are commonly described on both maps
It is recommended that the equipment map database owned by each
(The first graphic database) includes roads, town blocks, sidewalks,
The map is centered on line data such as the median strip.
There are very few points common to both maps.
The available point is the intersection of the above lines
There are only intersections. However, intersections have complex shapes
Therefore, the surveying accuracy around the area is poor, and the intersection
Conversion using point coordinates (intersection coordinates) because there are many changes
The problem with the method is that the conversion accuracy must be poor.
was there. In the method (c), the areas are associated with each other by a road.
Lines, block lines, etc.
For example, road expansion work and corner cutting at intersections
Even if road lines and block lines change over time due to
Area extraction and association,
Executes parallel movement of equipment data.
The displacement of the data may increase
There was a problem. Further, the extraction and the correspondence of the areas are automatically performed.
In doing so, create a program with fewer exceptions.
Difficult to achieve, and the development period is long,
There was a problem that it had to be expensive. In the method (d), a geometric
The relationship between the relationship and the affine transformation becomes complicated and
There was a problem that it was difficult to be. In addition, correspondence of two maps
To determine the angle of deviation between the segments
The sign of the required angle changes depending on the direction.
Takes great care when programming with care
There was a problem. In order to specify a required point on the map,
Extreme precision is required for the operator to specify data
Problem. The present invention has been made in view of such a problem.
The purpose is to use the facilities on one map
Affine transformation of data position with high accuracy and efficiency
Efficiently to match the map data etc.
Ability to determine accurate affine transformation parameters
To provide a graphic processing device. [0020] [MEANS FOR SOLVING THE PROBLEMS] To achieve the above-mentioned object.
In order to achieve this, the present invention uses a map
Position coordinates and location of equipment such as gas pipes and water pipes
A first graphic database storing mark data and a map data
A second graphic database storing the position coordinate data;
(A) the first figure
At least one point in a predetermined area of the shape database;
The same point as the corresponding point in the second graphic database
A coordinate value providing means for providing a standard value; and (b) the first figure
Included in the predetermined area on the databaseBlock line
Line segment dataFrom a plurality of sets of a first line segment and a second
First line segment designating means for designating a line segment;
, A straight line including the first line segment and the second line
To calculate a first intersection with a straight line including the line segment
Point calculation means; and (d) the first figure database
On the second graphic database corresponding to the first line segment
A third line segment and a second line segment on the first graphic database;
A fourth on the second graphic database corresponding to the line segment
Second line segment designating means for designating a line segment;
A straight line including the third line segment, and a straight line including the fourth line segment.
A second intersection calculating means for obtaining a second intersection of (f)
An X direction between the first intersection and the second intersection;
A shift amount calculating means for calculating a shift amount with respect to the Y direction;
(I) Least squares method using a plurality of sets of the shift amounts
Parameter determinator for determining affine transformation parameters
Steps and (k) the location of the equipment in the first graphic database
The coordinate data is stored in the position of the map in the second graphic database.
The parameters must be modified to correspond to the coordinate data.
Parameters of the affine transformation determined by the data determination means
Affine transformation means for performing affine transformation using
A graphic processing apparatus characterized in that: [0021] The first figure for a predetermined common target area
The first block line data extracted from the database and the second block line data
The second block line data extracted from the figure database of
Are distinguished by color, etc.
In the state where the display is folded and displayed, first,
The step allows a small number of objects within the target area in the first graphic database.
At least one point and the corresponding point in the second graphic database
And the same coordinate value. Next, the first line segment designating means
Specify first and second line segments from block line data
Then, the first intersection calculating means calculates the first line segment and the second line segment.
Is calculated. Then, the second line segment designator
Steps correspond to the first line segment from the second block line data
The fourth line segment corresponding to the third line segment and the second line segment
And the third line segment and the fourth line are calculated by the second intersection calculating means.
Calculate the second intersection with the minute. Next, the deviation amount calculating means
Thus, the X direction between the first intersection and the second intersection,
The shift amount from the Y direction is calculated. From this line designation,
By repeating the process up to the calculation of the deviation amount,
Data. Next, a plurality of sets of
Parameter for affine transformation by the least squares method
Data. Affine transformation means is a parameter determinant
Using the parameters determined by the step, the first graphic
Affine transformation of position coordinate data of
The position coordinate data of the map in the second graphic database.
Modify to correspond to the data. [0024] BRIEF DESCRIPTION OF THE DRAWINGS FIG.
Will be described. FIG. 1 illustrates the entire implementation environment of this embodiment.
FIG. 2 is a diagram showing an example of a system configuration according to the present embodiment.
It is. First, the outline of the components will be described. Figure
1. As shown in FIG. 2, a graphic processing apparatus 1 according to the present embodiment
Is the host computer 3 and the host computer
External storage device 5 connected to line 3 by terminal, terminal device group
7, an electrostatic plotter 9 is provided. The host computer 3 has an internal storage
A work area in the internal storage device 11
13 are provided. The external storage device 5 has already
The first graphic database 21 which is a database
2 graphic database 23, application program
And the like are stored. The first diagram in the internal storage device 11
The first data extracted from the shape database 21
The graphic work file 27 and the second graphic database
The second graphic workflow that extracts necessary data from
File 29 is stored in the first graphic work file 27.
A third word for converting, modifying and saving equipment data.
Prepare a backup file 31. Later, a third workf
File 31 is integrated as necessary, and the converted
The third graphic database 33 including preparatory data is
It is built in the storage device 5. Each hardware and host in the graphic processing apparatus 1
Is connected to the remote computer 3 by a communication line 35 or the like.
The terminal device group 7 includes a tablet 41, a digitizer
43, display (graphic display) 4
5, PC 47, hard copy machine 49, display
Pointer for moving the position of the cursor mark 51 on 45
53, light pen, digital camera connected to tablet 41
The cursor 57 and the keyboard 59 connected to the
Have. Also, as shown in FIG.
Of the center system 71
And the host computer 73 via a communication line 75.
Continued. Center Host computer 73 center
-The external storage device 77 has a center facility map database.
And the basic map database 83
Is done. Center host of center system 71
The computer 73 includes various items other than the graphic processing device 1.
Various user systems 91 provided by traders
a, 91b... are also connected to the respective host computers.
You. The graphic processing device 1 includes a tablet
41 is used when searching the graphic database, and
Various instructions are given to the strike computer 3. On this occasion,
A menu with the contents of various operations (processes) for search printed
Paste the display sheet on the tablet 41, and
After specifying the four corners of the menu sheet frame with
A desired menu is selected and designated with the light pen 55. The digitizer 43 is provided with, for example, the position data of a gas pipe.
Data (coordinate data), gas pipe material, thickness, etc.
It is used to input graphic data.
Use the menu sheet to specify input contents. However
The selection of a menu or the like is performed by the cursor 57. Display (graphic display)
B) 45 is each figure data stored in the external storage device 5
It displays information on the base in a graphic format.
Color display of various data based on the indicated search content
Etc. and selectively displayed. As needed
Display method for partially enlarging or displaying the progress of processing, etc.
Change and display. The indicator 53 is located on the display 45
Move the cursor mark 51 on the display 45
By specifying the shape of the data, the data displayed as the shape
Used to specify The personal computer 47 is provided with a name such as a town or a house.
Used when inputting character data other than standard data.
You. The hard copy machine 49 has a display
Copy the graphic data displayed in b45 as a drawing on paper.
The electrostatic plotter 9 is used for external storage.
The information on the graphic database stored in the device 5 is used as a drawing.
Used for drawing on paper. The center basic map database 83 is
The second storage stored in the external storage device 5 of the shape processing device 1
This is a higher-level database including the graphic database 23.
That is, from the center basic map database 83,
Extract data of specified contents for the area of
The result is the second graphic database 23. C
The center equipment map database 81 is stored in the graphic processing device 1.
Various things including the built first graphic database 21
It is composed of original equipment map data created by the contractor
I have. In FIG. 1, the center system 71
And the graphic processing device 1 or the user systems 91a, 9
1b are connected by a communication line 75,
Center equipment map database 81 or center basic map
Extract and record only predetermined data from database 83
First graphic database with recorded magnetic tape, etc.
21 or the second graphic database 23
It is also possible to give and receive. FIG. 3 shows the processing performed by the graphic processing apparatus 1.
FIG. 6 is a diagram for explaining a method of preparing data for use in the data processing. This implementation
In the example, rent a magnetic tape and read it,
Of the basic map database 83
Stored in the external storage device 5 as the shape database 23
It is explained as if it is. In addition, about corresponding area
First figure data consisting of map data and equipment data
The base 21 is also already constructed and stored in the external storage device 5
It is assumed that Next, the block lines are obtained from the two graphic databases.
Retrieve position coordinate data. The data retrieved here is
Topographic data such as sidewalks and center lines may be used,
Block lines are used because there is little change. The following is an example of a block line
I will explain. The first graphic database 21 stores the first
All block line data 101 and facility data 103 are stored.
In the second graphic database 23, the second all block lines
Data 105 is stored. Here, the first whole block
Coordinates of end points of each block line data included in the line data 101
And a corresponding mark included in the second all block line data 105
Generally, there is no space between the coordinates of the end point of each block line data.
This exists. The host computer 3 has an external storage device 5
1st all block line data of the first graphic database 21
From the data 101, the first block line data included in the predetermined area
The data 107 is searched and extracted, and the first
In the graphic work file 27 of FIG. Also, the same prescribed
From the second graphic database 23,
Search and extract the block line data 109 of No. 2 and store it in the internal storage
In the second graphic work file 29 in the storage 11. The host computer 3 has
The terminal device group 7 processes each block line data in the storage device 11.
After converting to format data in a format that
It transmits to the specified terminal equipment group 7. Terminal equipment group 7 transmits
Converts the input data into internal commands and displays
Display of color, thickness, line type, blinking, etc. specified for 45
Display according to the method. Thus, the first block line day
Data 107 and the second block line data 109 are displayed.
A) Superimposed and displayed on the screen 45, distinguished by color etc.
In the state, the operator uses the graphic processing apparatus 1 of the present embodiment to
Simple position specification work and automatic affine transformation
Is performed continuously. FIG. 4 is a functional block diagram of the graphic processing apparatus 1.
is there. As shown in FIG. 4, the graphic processing device 1
The first graphic work file is output by the
27 and at least one second graphic work file
The same coordinate value is given to the corresponding point in 29. Next, the first line segment designating means 303
From the first block line data 101, the first line segment and the second line segment
Is designated, and the first intersection calculating means 305 specifies the first line segment.
The first intersection of the line segment with the second line segment is calculated. Continued
The second line segment designation means 307 uses the
From the data 103, a third line segment corresponding to the first line segment,
Specify the fourth line segment corresponding to the second line segment, and calculate the second intersection
The second intersection of the third line segment and the fourth line segment
calculate. Next, the first amount is calculated by the shift amount calculating means 309.
The X direction and the Y direction between the intersection and the second intersection,
Calculate the shift amount. From designation of this line segment to calculation of deviation
Is repeated to obtain a plurality of pieces of deviation amount data. Next, the parameter determination means 413
For affine transformation by least squares method using several sets of shift amounts
Determine the parameters. The affine transformation means 415
Using the parameters determined by the meter determination means,
First graphic work file 27 or first graphic data
Affiliate the position coordinate data of the equipment stored in the base 21
Of the map in the second graphic database 23
Correct so that the positional relationship matches the positional coordinate data. FIG. 5 shows a plurality of maps by the graphic processing apparatus 1.
9 is a flowchart illustrating a matching operation between data. First
The operator of the graphic processing apparatus 1 pastes the
Select and enter numbers from the attached menu sheet
This specifies the mesh number. The specified message
Of the host computer 3 via the communication line 35.
Transmitted to Then, the host computer 3 sends the second
The second all block line data 105 in the graphic database 23
The second block line data corresponding to the specified mesh number.
Data 109, convert the format, and
Converted to the spray 45 internal command, the second block line 1
Displayed on the display 45 in the designated color as 21
(Step 501). At this time, the display 3
The display content is, for example, a display example of the second block line 121 in FIG.
Is shown by a dotted line in FIG. Similarly, the host computer 3 sends the first
Whether it is the first all block line data 101 in the graphic database 21
Then, the corresponding first block line data 107 is extracted. The reference of the first graphic database 21
The coordinates are different from the reference coordinates of the second graphic database 23.
For example, if the
Shape Processing Engineering (Nikkan Kogyo Shuppan, June 10, 1981)
Line) ”on page 81.
The position coordinates of the first graphic database 21 are
Position coordinates based on the reference coordinates of the graphic database 23
(Step 502). The first block line data 107 or the first town
Regarding the data obtained by converting the position coordinates of the section line data 107,
To further convert the format, and
To the first block line 123 and the second block line
121 and displayed on the display 45 (step
503). FIG. 7 shows the first block 123 and the second town.
FIG. 4 is a diagram in which a ward line 121 is distinguished by a solid line and a dotted line.
However, on the display 45, different colors such as red and yellow
If the colors are displayed in different colors, the distinction becomes easier. As shown in FIG. 7, the second block line 121
There is a gap with the first block line 123,
From step 504, in step 512, this deviation is minimized.
An affine transformation parameter such as First, e
The strike computer 3 determines how many sets of deviation amounts have been calculated.
Initialize the indicated parameter i to “0 (zero)”
(Step 504). Next, the first deviation amount is obtained.
Therefore, the value of the parameter i is increased by “1” (the
Step 505). Next, the first line segment designating means 403
The operator selects an arbitrary block line among the first block lines 123
Nearby point Ga_iMove the cursor mark 51 to
When the tool 53 or the like is clicked, the coordinate value of the point Gai is entered.
(Step 506). Here, the operator
You do not need to specify the exact location of the
You just need to specify your body position. Next, the position coordinates of the cursor mark 51 are
Sent to the host computer 3 through the communication line 35,
The host computer 3 stores the first graphic database 21
Searches the block location coordinate data and sends the cursor marker
Segment data of the block line that is closest to the position coordinates of block 51
Search for. From the position coordinates of the cursor mark 51, for example,
The method of calculating the distance to all line segments was September 1986.
Sep. 10th issue of Bit magazine “Computational Geometry and Geographic Information Processing”
On page 40, etc. Also searched
The shortest line segment and the position of the cursor mark on that line segment
The position coordinates of the intersection P1 with the perpendicular drawn from the mark are
It can be easily obtained from any school. From cursor mark 51, etc.
If there are multiple block lines in the distance, the display
Each block line is distinguished and displayed on the display 45, and the indicating tool 53 and the like are displayed.
Each time you click, the block line blinks one by one in order
For example, the operator can easily extract the target block line.
To do. Thus, one block line is converted to the first line segment ga
_iAs blue, distinguished from other block line colors, for example, blue
Is displayed on the display 45 (step 50).
7). Subsequently, the first line segment designating means 403 operates
Data is another point Gb on the display 45_iSpecify
And the first intersection calculating means 405 determines whether the first block line 123
To another block line, the second line segment gb_iExtracted as the first
Line segment ga_iAnd a second line segment gb_iIncluding
The intersection with the line is calculated, and this is defined as the first intersection G_iShow as
(Step 508). Here, the first line segment ga_iWhen
The second line segment gb_iAnd if they actually intersect
If no intersection is found and they do not intersect,
Or the intersection G on the extension of both line segments_iIs required
You. Next, the second line segment designation means 407
The operator moves the first line segment ga from the second block line 121.
_iThe third line segment fa corresponding to_iAnd the second line segment gb_iTo
The corresponding fourth line segment fb_iTo specify the third
Line segment fa_iFa in the vicinity of_iAnd the fourth line segment fb_iNear
Nearby point Fb_iIs designated, the second intersection calculating means 40
9, the third line segment fa_iAnd a fourth line segment f
b_iIntersection F with a straight line containing_iIs calculated and displayed
(Step 509). Next, the shift calculating means 411 calculates the two intersections
G_iAnd F_iΔX in the X direction of the coordinates_i, Y direction difference
Minute δY_iIs calculated (step 510). This corresponding
There are two methods to calculate the corresponding intersection by specifying the block line
To find and specify the corresponding point on the existing database
The corresponding block line itself has a cursor mark 51
Each of the corresponding block lines without having to specify exactly
Since it is only necessary to specify the neighborhood, the burden on the operator is reduced.
You. Therefore, much easier for corresponding sets of points
Parameter for more accurate affine transformation
You can ask for a meter. Steps 505 to 510 are performed.
And the two intersections G_i, F2_iSet of
Specify the positions differently, and select them according to the predetermined criteria.
If you have finished specifying the line segment
When designation of up to a predetermined number of sets, for example, N sets, is completed,
It is assumed that the operation of specifying the block line by the operator has been completed.
You. Block line finger by operator for any reason
It is determined whether or not to end the regular work (step 511),
Affine transformation, where N is the value of parameter i at that time
It is used when calculating parameters for use. Specify in step 511
Is completed, the parameter determination unit 413 determines that N
Two intersections G in a set_i, F_iX direction calculated for the coordinates of
From the data of each difference between the Y direction and the Y direction,
The parameters are calculated (step 512). Next,
The conversion unit 415 is a unit in the work area 13 before conversion.
Data 111 is used for the obtained affine transformation parameters.
Affine transformation (step 513). Affine strange
The converted equipment data 13 obtained by the conversion is the third graphic data.
Is stored in the database 33 (step 514). The affine transformation is applied to all the meshes.
After the execution, the first
Can be deleted. sand
That is, the database in the external storage device 5 of the graphic processing device 1
As shown in FIG. 3, before the conversion process,
The user map database 19 created by the vendor
First graphic database 21 including equipment data as background
And the center basics provided by the center system 71
The second figure data with the map in the map database 83 as the background
Database 23 and all meshes
After the conversion process is completed, as shown at the bottom of FIG.
The second graphic database 23 and the second graphic database
FIG. 3 is a third diagram showing the positions of the facilities with the map of FIG.
Only the shape database 33 is stored. Next
In step 512 of FIG.
Regarding the method of obtaining the parameters, the first intersection G in FIG.
_iAnd the second intersection F_iDecision processing and affine transformation
This will be described in detail with reference to FIG. FIG. 8 (a)
Intersection G_iAnd the intersection F_i8 (b).
Is the difference δX in the X direction obtained in step 510._iand
Difference δY in Y direction_i, Deviation δL_iIs shown. Toes
F on the display 45_iCoordinates of (X_i+ Δ
X_i, Y_i+ ΔY_i). The shift δ shown in FIGS.
L_iIs the difference δX in the X direction._iAnd difference in Y direction
Minute δY_iFrom the following data. [0060] (Equation 1) N sets of differences δX in the X direction_iAnd the difference δY in the Y direction_iof
Least squares separately for X and Y directions using data
The affine transformation parameters by the
512). That is, in step 315, generally, the first
Each coordinate point (X, Y) of the graphic work file 27 of
2 (x
^*, Y^*) x^*= − (AX + bY + c) + X y^*=-(DX + eY + f) + Y For the entire first graphic work file 27.
Affine transformation parameters a, b,
It is necessary to determine c, d, e, and f. For each data set, the first figure
I-th point G of work file 27_iCoordinates (X_i, Y
_i) To the coordinates of the second graphic work file 29
Coordinates (x^* _i, Y^* _iAffine transformation to correct to)
To x^* _i=-(AX_i+ BY_i+ C) + X_i y^* _i=-(DX_i+ EY_i+ F) + Y_i Then, obtained from the second graphic work file 29,
Point G of the first graphic work file 27_iPoint F corresponding to
_iCoordinates (x_i, Y_i) And the point G after the affine transformation_iof
Coordinates (x^* _i, Y^* _i) And the error ε in the X direction_xi, Y direction
Error ε_yiAre represented by the following equations. [0063] ε_xi= X_i-X^* _i = Δx_i-(AX_i+ BY_i+ C) ε_yi= Y_i-Y^* _i = Δy_i− (DX_i+ EY_i+ F) At this time, the error ε in the X direction_xi, The error ε in the Y direction_yi
From the error magnitude δl after affine transformation_iIs [0064] (Equation 2) Can be obtained by Figure 9 shows before and after affine transformation
FIG. 6 is a diagram for explaining a positional relationship and the like of each point of FIG._iCoordinates
(X_i, Y_i), Point G_iPoint F corresponding to_iCoordinates
(X_i, Y_i), Point G_iAnd point F_iDistance δL_i, X direction
Direction difference δX_iAnd the difference δY in the Y direction_iAffine strange
By exchange, point G_iTo the coordinates of the second graphic work file 29
Coordinates (x^* _i, Y^* _i), Affine strange
Error magnitude δl after conversion_i, X-direction error ε_xi, Y direction
Error ε_yiShows the relationship. An error in the X and Y directions after the affine transformation
Difference ε_xiAnd ε_yiTo minimize the sum of squares of
The conditions to be satisfied are [0066] (Equation 3)It is. Therefore, the error ε_xiAnd ε_yiThe sum of squares of
Affine transformation parameters a, b,
To find c, d, e, and f, [0067] (Equation 4) Solve the simultaneous equations of To solve this system of equations
Is, for example, “Scientific calculation
Routine Library "(Keigaku Shuppan P84-86)
If you prepare for the program, the host computer 1
Can be calculated. The parameter for affine transformation
A is the expansion and contraction rate in the X direction, b + d is the flattening rate, and c is X
E is the amount of expansion and contraction in the Y direction, f is the movement in the Y direction
It is about quantity. Next, the equipment data 111 before the conversion and the conversion
The later equipment data 113 will be described. For example, figure
As shown in FIG.
The installation line that is the equipment data 111 before conversion
When the first block line 123 is indicated by a thin solid line with a wavy line
If there is no
When the block line 121 of FIG.
(The facility data 111 before conversion) and the second block line 121
It can be seen that the relative positional relationship between the two greatly deviates.
However, due to affine transformation, it is arranged as shown in FIG.
When the line is a thick broken line, the second diagram shown by a thin broken line
With the second block line 121 included in the shape work file 29,
It can be understood that the relative positional relationship is appropriately displayed. This coordinate conversion is performed in the first graphic database 2
1 for all equipment data 103 (arrangement lines)
This is performed for the constituent line segments. As described above, the surveying accuracy is high,
By using the line data of the map with little change over the years,
Determining highly accurate affine transformation parameters
it can. Also, if you use line data, the line data
Internal processing even when a point near the specified point is specified without specifying
Can be used to specify line data.
Data input operation is also facilitated. Moreover, affine transformation para
The meter is calculated using a statistical method called least squares.
Therefore, the operator specifies a slightly inappropriate block line etc.
Does not significantly affect the conversion accuracy and the operator's attention
Since a stable value that is not affected by force can be obtained,
Highly accurate conversion can be performed. Next, a second embodiment will be described. Second
Is an example of the type of block line used to determine the intersection,
For example, depending on national roads, city roads, ward roads, city roads, private roads, sidewalks, etc.
And affine transformation parameters are determined.
For example, national roads, city roads, ward roads, city roads, private roads, sidewalks
If one block line is specified for each,
These are 6, 5, 4, 3, 2, 1 finger each
Determine affine transformation parameters as specified
You. In this way, it is possible to set
Affine transformation is performed with higher accuracy as the position coordinate data
It becomes. As described above, the surveying accuracy is high,
Two figures using the line data of the map with little year change
From the database, a number of two line segments corresponding to each other
Find the set and calculate the amount of shift between the intersections of one set of two line segments.
And process the data by statistical methods.
Determines the parameters for affine transformation with high accuracy
can do. Also, if line data is used, line data
If you do not specify the data exactly and specify a point near it
Can also specify line data by internal processing
Thus, the input operation of the operator becomes easy. And Affi
Parameters are calculated using a statistical method called the least squares method.
Block line that is somewhat inappropriate for the operator because it is calculated more
Even if the user designates, etc., the error that occurs is negligible. [0073] As described above in detail, according to the present invention,
Can be used to efficiently and accurately determine the affine transformation
By converting equipment data based on parameters
Equipment data created based on a predetermined map.
Accurately and efficiently match other maps
A graphic processing device can be provided.

BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a view for explaining an entire implementation environment of the present embodiment. FIG. 2 is a view showing an example of a system configuration of a graphic processing apparatus 1. FIG. FIG. 4 is a diagram illustrating an example of a data preparation method. FIG. 4 is a functional block diagram of the graphic processing apparatus 1. FIG. 5 is a flowchart showing a matching operation between a plurality of map data. FIG. 7 is a diagram showing a display example. FIG. 7 is a diagram showing an example of simultaneous display of the second block line 121 and the first block line 123. FIG. 8 is a diagram showing a first intersection G _i and a second intersection F _i . FIG. 9 is a diagram illustrating a determination process and an affine transformation. FIG. 9 is a diagram illustrating a positional relationship and the like of each point before and after the affine transformation. FIG. FIG. 11 is a diagram showing an example of simultaneous display with facility data 111 of FIG. 11. FIG. Diagram showing an example of simultaneous display with replaced equipment data 111 [Explanation of reference numerals] 1... Graphic processing device 3... Host computer 5... External storage device 7. ... internal storage device 13 ... work area 21 ... first figure database 23 ... second figure database 27 ... first figure work file 29 ... second figure work file 33 ... ... Third graphic database (including converted equipment data) 45 Display (graphic display) 53 Indicator 71 Center system 81 Center equipment map database 83 Center Basic map databases 91a, 91b... User system 101... First whole block line data 103... Facility data 105... Data 107 First block line data 109 Second block line data 111 Facility data before conversion 113 Converted facility data 401 Coordinate value assigning means 403 ... First line segment specifying means 405... First intersection point calculating means 407... Second line segment specifying means 409... Second intersection point calculating means 411. ... Parameter determination means 415 ... Affine transformation means

──────────────────────────────────────────────────続 き Continued on the front page (58) Field surveyed (Int.Cl. ⁷ , DB name) G06T 11/60-11/80 G06F 17/50

Claims (1)

(57) [Claims 1] A first figure database storing map position coordinate data and position coordinate data of equipment such as gas pipes and water pipes, and map position coordinate data are stored. Second
(A) giving the same coordinate value to at least one point in a predetermined area of the first graphic database and a corresponding point in the second graphic database (B) a plurality of sets of a first line segment and a second line segment from the line segment data of block lines included in the predetermined area on the first graphic database; (C) calculating a first intersection of a straight line including the first line segment and a straight line including the second line segment for the plurality of sets; (D) a third line segment on the second graphic database corresponding to a first line segment on the first graphic database;
Second line segment designating means for designating a fourth line segment on the second graphic database corresponding to a second line segment on the first graphic database; (e) the third line A second intersection calculating means for obtaining a second intersection between the straight line including the second segment and the straight line including the fourth segment; and (f) a second intersection between the first intersection and the second intersection. A shift amount calculating unit that calculates a shift amount between the X direction and the Y direction; and (i) a parameter determining unit that determines an affine transformation parameter by a least square method using a plurality of sets of the shift amounts. (K) affine transformation parameters determined by the parameter determination means so as to correct the position coordinate data of the equipment in the first graphic database so as to correspond to the position coordinate data of the map in the second graphic database. Affine transformation using affine transformation Graphics processing apparatus characterized by comprising: a stage, a.