Device for calculating positional data of standard points of photogrammetric target

[0002] 1. Field of the Invention

[0003] This invention relates to a calculation of camera position (photographing position) in a photogrammetric analytical measurement system, in which a survey map is produced on the basis of a set of photographed pictures obtained at two different photographing positions.

[0004] 2. Description of the Related Art

[0005] For example, photogrammetry is carried out at a traffic accident spot. The traffic accident spot is photographed by a camera in at least two different positions, and a survey map of the traffic accident spot is produced based on a set of photographed pictures obtained at the different positions.

[0006] In particular, a two-dimensional coordinate system is defined on each of the photographed pictures, and two-dimensional positions of the objects, which are recorded on each picture, are determined by the two-dimensional coordinate system. Then, a three-dimensional coordinate system is defined on the basis of the two sets of two-dimensional coordinate systems, and three-dimensional positions of the recorded objects are determined from the three-dimensional coordinate system. Accordingly, it is possible to produce a survey map of the traffic accident spot by drawing the objects on a sheet of paper in such a manner that the objects are projected on one of the three planes defined by the three-dimensional system.

[0007] Before accurately scaled distances and lengths can be reproduced on the survey map, a standard measurement scale must be recorded together with the objects in the photographed pictures. Also, a standard reference plane, on which the survey map should be drawn, must be defined in the photographed pictures.

[0008] Usually, in order to define the standard measurement scale and the reference plane, three respective cone-shaped markers, which are identical to each other, are positioned at suitable locations around a traffic accident spot. Namely, a distance between two apexes of the cone-shaped markers is measured, for example, with a measuring tape, and set as the standard measurement scale, and a plane, defined by the three apexes of the cone-shaped markers, is utilized as the reference plane. As the positioning of the cone-shaped markers and the measurement of the distance between the two apexes are performed by operators, a process of preparing for photographing is onerous and requires a lot of time.

[0009] In order to solve the above mentioned problems, photogrammetric analytic measurement, using a target which has a triangular frame instead of the three cone-shaped markers, is disclosed in Japanese Unexamined Patent Publication (Kokai) No. P10-221072. Standard point members are respectively mounted on each of three apexes of the target. A distance between the standard point members is set as the standard measurement scale, and a plane defined by the standard point members is utilized as the reference plane. A reflecting sheet, for example, is attached on the standard point members so that the standard point members can be viewed easily in photographed pictures.

[0010] The determination of the three-dimensional positions of the recorded objects from the three-dimensional coordinate system is carried out by iterating a series of approximate calculations, using a computer having a monitor on which the set of photographed pictures is displayed.

[0011] Initially, before a first series of approximate calculations is executed, the standard point members and a suitable point on one of the recorded objects in the set of photographed pictures displayed on the monitor, are selected with a cursor, by manipulating a mouse, whereby three sets of two-dimensional coordinates, corresponding to the standard point members, and a set of two-dimensional coordinates on the suitable point of the selected object are inputted to the computer. Namely, the execution of the first series of approximate calculations is based upon the inputted two-dimensional coordinates of the standard point members and the inputted two-dimensional coordinates of the suitable point on the selected object, thereby determining three-dimensional positions of the standard point members and a three-dimensional position of the suitable point of the selected object from the three-dimensional coordinate system.

[0012] Then, another suitable point on the selected object in the set of photographed pictures displayed on the monitor is indicated with the cursor, by manipulating the mouse, whereby a set of two-dimensional coordinates of the other suitable point on the selected object is inputted to the computer, and a second series of approximate calculations is executed on the basis of the inputted two-dimensional coordinates of the other suitable point on the selected object, thereby determining a three-dimensional position of the other suitable point on the selected object from the three-dimensional coordinate system. This procedure is continued until a sufficient number of points on the selected object are indicated, to thereby specify a three-dimensional profile of the selected object.

[0013] The aforementioned series of approximate calculations is executed based on positional data of the camera at which the pictures are photographed. The positional data of the camera includes a distance from the camera to the selected object, an angle of the camera against the object. The positional data of the camera is calculated based on the two-dimensional coordinates of the standard point members selected by the operator.

[0014] However, it takes a comparatively long time to select the standard point members in the set of photographed pictures displayed on the monitor, with the cursor, by manipulating the mouse. Therefore, when there are a lot of the photographed pictures, it becomes onerous for the operator to operate the above mentioned manipulation of the mouse. Further, the accuracy of the photogrammetric analytic measurement is dependent upon the skill of the operator.

[0015] Therefore, an object of the present invention is to provide a device for calculating positional data of standard points of a target of photogrammetric analytic measurement.

[0016] In accordance with an aspect of the present invention, there is provided a device, for calculating positional data of standard points of a target of photogrammetric analytic measurement, comprising an image data obtaining processor and a value of two-dimensional coordinates calculating processor.

[0017] The image data obtaining processor photographs an object of photogrammetric analytic measurement with a target that has a first standard point and a second standard point and a third standard point.

[0018] Each length of distance between the first standard point and the second standard point and the third standard point is predetermined. An angle defined by a first straight line connecting the first standard point and the second standard point, and a second straight line connecting the second standard point and the third standard point, is predetermined. A number of assistant points provided on the first straight line and a number of assistant points provided on the second straight line are different.

[0019] The value of two-dimensional coordinates calculating processor calculates each two-dimensional coordinate value of the first standard point, the second standard point and the third standard point in a photograph coordinate system that is a two-dimensional coordinate system of a photographed image, obtained by the image data obtaining processor, after determining each position of the first standard point, the second standard point and the third standard point in the photographed image, by comparing the number of the assistant points provided on the first straight line and the number of the assistant points provided on the second straight line.

[0020] The first standard point, the second standard point, the third standard point and the assistant points provided on the first straight line and provided on the second line are made of a reflecting material, such that the first standard point, the second standard point, the third standard point and the assistant points have a relatively higher brightness than other objects in the photographed image.

[0021] The value of two-dimensional coordinates calculating processor is provided with: a binarization processor that classifies each pixel composing of digital image data of the photographed image into a high brightness class and a low brightness class after comparing a brightness of each pixel with a predetermined threshold value; a group extracting processor that extracts groups in which a plurality of the pixels of the high brightness class sequentially exist in the digital image data, an area of which is within a predetermined range; a calculating processor that calculates a value of two-dimensional coordinates of a center of gravity of each of the groups, based on the brightness of each pixel composing the groups and a value of coordinates of the each pixel in the photograph coordinate system; a first points range extracting processor that extracts a first points range on which a plurality of the centers of gravity lie, a number of the plurality of the centers of gravity on the first points range being equal to a total point number including the first standard point, the second standard point and the number of the assistant points on the first straight line; a second points range extracting processor that extracts a second points range on which a plurality of the centers of gravity lie, a number of the plurality of the centers of gravity on the second points range being equal to a total point number including the second standard point, the third standard point and the number of the assistant points on the second straight line; a combination extracting processor that extracts a combination of the first points range and the second points range in which the first points range and the second points range share the center of gravity of an identical end point; and a standard point determining processor that determines a center of gravity of an end point of the first points range selected by the selecting processor, being opposite to the identical end point shared by the first points range and the second points range, as the first standard point, determines the center of gravity of the identical end point shared by the first points range and the second points range as the second standard point, and determines a center of gravity of an end point of the second points range selected by the selecting processor, being opposite to the identical end point shared by the first points range and the second points range, as the third standard point.

[0022] The number of the assistant points provided on the first straight line is two, and the number of the assistant points provided on the second straight line is one, so that the first points range is a four points range in which four centers of gravity range on a same straight line in the photographed image, and the second points range is a three points range in which three centers of gravity range on a same straight line in the photographed image.

[0023] In accordance with another aspect of the present invention, there is provided a method for calculating positional data of standard points of a target of photogrammetric analytic measurement comprising: a first step that photographs an object of photogrammetric analytic measurement with a target that has a first standard point and a second standard point and a third standard point, by a camera; and a second step that calculates each two-dimensional coordinate value of the first standard point, the second standard point and the third standard point in a photograph coordinate system that is a two-dimensional coordinate system of a photographed image photographed by the camera.

[0024] Each length of distance between the first standard point and the second standard point and the third standard point is predetermined. An angle, defined by a first straight line connecting the first standard point and the second standard point, and a second straight line connecting the second standard point and the third standard point, is predetermined. A number of assistant points provided on the first straight line and a number of assistant points provided on the second straight line are different.

[0025] In the second step, each two-dimensional coordinate value of the first standard point, the second standard point and the third standard point in the photograph coordinate system, after determining each position of the first standard point, the second standard point and the third standard point in the photographed image based on a difference between the number of the assistant points provided on the first straight line and the number of the assistant points provided on the second straight line.

[0026] In accordance with another aspect of the present invention, there is a provided a recording medium in which a calculating program of standard points of a target of photogrammetric analytic measurement is stored.

[0027] The target has a first standard point and a second standard point and a third standard point. Each length of distance between the first standard point and the second standard point and the third standard point is predetermined. An angle, defined by a first straight line connecting the first standard point and the second standard point, and a second straight line connecting the second standard point and the third standard point, is predetermined. A number of assistant points provided on the first straight line and a number of assistant points provided on the second straight line are different.

[0028] The calculating program comprises: a values of two-dimensional coordinate calculating routine that calculates each two-dimensional coordinate value of the first standard point, the second standard point and the third standard point in a photograph coordinate system that is a two-dimensional coordinate system of a photographed image, photographed in such a manner that an object of photogrammetric analytic measurement is photographed with the target by a camera, after determining each position of the first standard point, the second standard point and the third standard point in the photographed image based on a difference between the number of the assistant points provided on the first straight line and the number of the assistant points provided on the second straight line.

[0029] The object and other objects of the present invention will be better understood from the following description, with reference to the accompanying drawings, in which:

[0030] FIG. 1 is a plan view of a target in an operational state, to which an embodiment, according to the present invention, is applied;

[0031] FIG. 2 is a side view of the target of FIG. 1;

[0032] FIG. 3 is a sectional view of the target of FIG. 1;

[0033] FIG. 4 is a plan view of a non-reflecting member viewed from a bottom side of a second bar;

[0034] FIG. 5 is an enlarged view of a control unit box of the target;

[0035] FIG. 6 is a sectional view of the control unit box of FIG. 5;

[0036] FIG. 7 is a conceptual view of a positional relationship between the target and a camera;

[0037] FIG. 8 is a photographed image photographed by the camera;

[0038] FIG. 9 is a conceptual view of a positional relationship between the photographed image and the target;

[0039] FIG. 10 is a conceptual view of a relationship between a scene coordinate system and a camera coordinate system;

[0040] FIG. 11 is a flowchart indicating a process of calculating a position of a camera;

[0041] FIG. 12 is a flowchart of procedure of calculating a center of gravity of a group;

[0042] FIG. 13 is a conceptual view of each pixel of a digital image data of the photographed image;

[0043] FIG. 14 is a first half of a flowchart of procedures of extraction combinations of the groups;

[0044] FIG. 15 is shows the latter half of the flowchart of the procedures of extraction combinations of the groups;

[0045] FIG. 16 is a flowchart of procedures of calculation of a camera position; and

[0046] FIG. 17 is a block diagram of a coordinate-calculator system, in which a photogrammetric analytic measurement is performed, according to the present invention.

[0047] FIG. 1 is a plan view of the target 10, with portions broken away for clarity, and FIG. 2 is a side view of the target 10. The target 10 has an L-shaped figure, comprising a first bar 12 and a second bar 14. The first and second bars 12, 14 are made of metal material. The first and second bars 12, 14 respectively have a shape of a quadratic prism, which is hollow. A non-reflecting sheet, which does not reflect light, is attached to the whole outer side surface of each of the first and second bars 12, 14. The width of each of the first and second bars 12, 14 is a length L_W, and the thickness of each of the first and second bars 12, 14 is a length L_H.

[0048] An adhesive is applied to the surface of the non-reflecting sheet, which contacts the surface of the bars 12, 14. The other surface of the sheet is colored black, and rough. Incident light on the black and rough surface is absorbed and diffused, so that a luminance amount of reflected light is extremely reduced. Note that, for example, a black flattering or matter agent can be applied on the outer surface of each of the first and second bars 12, 14, instead of utilizing the non-reflecting sheet.

[0049] A controlling unit box 20, which is shaped parallelepiped, is unitarily fixed at one end of the first bar 12. The controlling unit box 20 is made of metal material. The non-reflecting sheet is attached to the whole outer surface of the controlling unit box 20. The thickness of the controlling unit box 20 equals the thickness L_Hof the first bar 12. The width of the controlling unit box 20 is twice the width L_Wof the first bar 12. The controlling unit box 20 is positioned in such a manner that a side surface 20band a side surface 12bof the bar 12 lie on a same plane. A side surface 20cof the controlling unit box 20 is parallel to a side surface 12cof the bar 12, being offset in a direction opposite to the side surface 20b.

[0050] One end 14aof the second bar 14 is rotatably mounted on the side surface 20cof the controlling unit box 20, by a hinge 15. A side surface 14bof the second bar 14 and an end surface 20dof the controlling unit box 20, opposite to the end surface at which the first bar 12 is fixed, lie on a same plane, when the target 10 is utilized for the aforementioned photogrammetry.

[0051] As shown in FIG. 1, the side surface 14cof the second bar 14 and the side surface 20cof the controlling unit box 20 define an angle θ. Namely, the angle θ is a right angle made by an axis 12x(represented as a broken line) of the first bar 12 and an axis 14x(represented as a broken line) of the second bar 14. A stay 16, which is a fixing member, is connected to the first and second bars 12, 14, at the side of the angle θ. Rotational movement of the second bar 14 is prevented by the stay 16, so that the positional relationship between the first and the second bars 12, 14, shown in FIG. 1, is maintained. The width and the thickness of the stay 16 are respectively smaller than the length L_Wand L_H, of the first and second bars 12, 14. Further, a length along the longitudinal direction of stay 16 is shorter than the length of the longitudinal direction of the first and second bars 12, 14.

[0052] The stay 16 is rotatably connected to the first bar 12 by a stay hinge 92, being attachable to and removable from the second bar 14 by a lock hinge 94. When the target 10 is utilized in an operational position shown in FIG. 1, the stay 16 is positioned so as to make a predetermined angle with each of the first and second bars 12, 14, so that the first and second bars 12, 14 make a right angle.

[0053] On the top surface of the target 10, namely, on the top surfaces of the bars 12, 14 and the controlling unit box 20, three standard point members 31, 34 and 36, and three assistant point members 32, 33 and 35 are mounted. The standard point member 31 (first standard point member) is a circular plate, the diameter of which is smaller than the width L_Wof the first and second bars 12, 14. The standard point members 34 (second standard point member), 36 (third standard point member) and the assistant point members 32, 33 and 35 are identical to the standard point member 31.

[0054] The standard point member 31 and the assistant point members 32, 33 are mounted on the top surface 12eof the first bar 12. The standard point member 34 is mounted on the top surface 20eof the controlling unit box 20. The assistant point member 35 and the standard point member 36 are mounted on the top surface 14eof the second bar 14. The assistant point members 32 and 33 are positioned in such a manner that centers of the assistant point members 32 and 33 lie on a straight line (first straight line) parallel to the axis 12x, which connects a center of the standard point member 31 and a center of the standard point member 34. Also, the assistant point member 35 is positioned in such a manner that a center of the assistant point member 35 lies on a straight line (second straight line) parallel to the axis 14x, which connects a center of the standard point member 34 and a center of the standard point member 36.

[0055] As described above, there are two assistant point members (32, 33) on the straight line connecting the standard point members 31 and 34, and there is one assistant point member (35) on the straight line connecting the standard point members 34 and 36. Namely, the number of assistant point member on each of the straight lines is different from each other. Further, the distance between the center of the standard point member 31 and the center of the assistant point member 32, the distance between the center of the assistant point member 32 and the center of the assistant point member 33, and the distance between the center of the assistant point member 33 and the center of the standard point member 34, are identical to each other. The distance between the center of the standard point member 34 and the center of the assistant point member 35, and the distance between the center of the assistant point member 35 and the center of the standard point member 36, are identical to each other. Furthermore, the distance between the center of the standard point member 31 and the standard point member 34, and the distance between the center of the standard point member 34 and the center of the standard point member 36, are identical to each other.

[0056] The aforementioned reference plane for the photogrammetric analytical measurement is defined by the standard point members 31, 34 and 36, and the assistant point members 32, 33 and 35. The standard measurement length is defined by a side length of an isosceles triangle apexes of which correspond to the standard point members 31, 34 and 36. Namely, a length of the distance between the standard point members 34 and 36, and a length of the distance between the standard point members 36 and 31 can be used the standard measurement length, as they are predetermined.

[0057] Note that, the angle θ is not restricted to 90°, and it is unnecessary to make the distance between 31, 34 and the distance between 34, 36 equal to each other. Namely, it is required that each value of the angle θ, the length of the distance between the standard point members 31, 34 and the length of the distance between the standard point members 34, 36 are predetermined. Considering facility of calculation in the photogrammetric analytical measurement system, it is preferable that the value of the angle θ is predetermined to be 90° and the distance between the standard point members 31, 34 and the distance between the standard point members 34, 36 are equal to each other.

[0058] As is apparent from FIG. 1, on each of the equilaterals of the isosceles triangle, a number of the assistant point members is different. Accordingly, as a positional relationship between the target 10 and other objects in a picture can be easily judged, determination of the photographing position of each picture for the photogrammetric analytical measurement system is facilitated.

[0059] It is desirable that the standard plane defined by the standard point members 31, 34, 36 and the assistant point members 32, 33, 35 is parallel to the surface of a road on which the target 10 is positioned. As the controlling unit box 20 is mounted on the first bar 12, the first bar 12 is heavier than the second bar 14. Accordingly, if the second bar 14 is connected to the first bar 12 only by the hinge 15 at the end portion thereof, the second bar 14 rises from the surface of the road, so that the standard point member 31 and the assistant point members 32, 33 on the first bar 12 do not lie on the same plane on which the assistant point member 35 and the standard point member 36 on the second bar 14 lie. Also, the standard plane might not be parallel to the surface of the road.

[0060] Accordingly, the first and second bars 12, 14 are connected by the stay 16 at center portions thereof, in order to prevent the second bar 14 from rising from the surface of the road and to keep the standard plane parallel to the surface of the road for more accurate photogrammetric analytic measurement.

[0061] The mounting of the hinge 15 produces an interval between the side surface 20cof the controlling unit box 20 and the end surface 14aof the second bar 14. An elastic member 19 (see FIG. 5) is disposed in the interval. Accordingly, instability of the second bar 14 is prevented. The elastic member 19 is made of, for example, rubber and sponge, being attached on the end surface 14aor the side surface 20c. Note that, a spring can be utilized instead of the elastic member 19.

[0062] A reflecting sheet is attached on the standard point members 31, 34, 36 and the assistant point members 32, 33, 35. As the surface of the reflecting sheet is smooth, and colored white, the amount of reflection of incident light on the sheet is high. The standard point members 31, 34 and 36 are respectively surrounded by non-reflecting members 41, 44 and 46. Similarly, the assistant point members 32, 33, 35 are respectively surrounded by non-reflecting members 42, 43 and 45. On the non-reflecting members 41, 42, 43, 44, 45 and 46, the non-reflecting sheet is attached, respectively. Accordingly, discrimination of the standard point members 31, 34, 36 and the assistant point members 32, 33, 35 in the pictures is facilitated, so that the photogrammetric analytic measurement can be performed accurately.

[0063] The target 10 comprises a first tilt sensor 52 and a second tilt sensor 54. Tilt angles of the target 10 around the axes 12xand 14xare sensed by the first and second tilt sensors 52, 54. The first tilt sensor 52 is mounted in the first bar 12, being placed between the assistant point members 32 and the standard point members 34. The tilt angle to a horizontal plane around the axis 12xis sensed by the first tilt sensor 52. The second tilt sensor 54 is mounted in the second bar 14, being placed between the standard point member 34 and the assistant point member 35. The tilt angle of the horizontal plane around the axis 14xis sensed by the second tilt sensor 54.

[0064] The first and second tilt sensors 52, 54 are connected to the controlling box 20 by a cable 17. Data of the tilt angles is transmitted to the controlling box 20 by the cable 17.

[0065] The angle of the standard plane to the horizontal plane is obtained by sensing the tilt angles around the axes 12xand 14xwhich are perpendicular to each other.

[0066] The target 10 comprises three legs 18. The legs 18 are mounted on the bottom surface opposite to the top surface on which the standard point members 31, 34, 36 and the assistant point members 32, 33, 35 are mounted. Each of the legs 18 respectively correspond to the standard point members 31, 34 and 36. In a spot for photogrammetric analytic measurement, the target 10 is supported above a road with the interval corresponding to the height of the leg 18 therebetween. Accordingly, the target 10 is positioned so as to be parallel to general inclination of the road, without being effected by roughness of the surface of the road.

[0067] With reference to FIGS. 3 and 4, the constructions of the assistant point member 35 and the non-reflecting member 45 are described below. FIG. 3 is a sectional view of the target 10, taken in the direction of the arrows substantially along the line IX—IX of FIG. 1. FIG. 4 is a plan view of the non-reflecting member 45, shown from the side of the second bar 14.

[0068] A magnet holder 62 is mounted on the top surface 14eof the second bar 14. A magnet 64, which is ring-shaped, is held in the magnet holder 62. The outer diameter of the magnet holder 62 approximately equals the width L_Wof the second bar 14. The magnet 64 and the magnet holder 62 are integratedly fixed to the second bar 14 by a screw 66. A reflecting sheet 68 is attached to a head 67 of the screw 66. The assistant point member 35 comprises the magnet holder 62, the magnet 64, the screw 66 and the reflecting sheet 68.

[0069] The non-reflecting member 45 comprises a circular plate 72, made of a material through which an electric (radio) wave transmits, for example, resin or rubber material. If the circular plate 72 is made of rubber material, breakage of the circular plate 72 is avoided when the circular plate 72 is accidentally dropped. A non-reflecting sheet 74 is attached to one surface of the circular plate 72. The diameter of the non-reflecting member 45 is approximately seven times the diameter of the head 67 of the screw 66. The thickness of the non-reflecting member 45 is slightly smaller than the thickness of the head 67.

[0070] An engagement hole 76 is formed at the center of the non-reflecting member 45. The diameter of the engagement hole 76 approximately equals the diameter of the head 67. An iron ring 78 abuts around the engagement hole 76, at the surface opposite to the surface on which the non-reflecting sheet 74 is attached. The inner diameter of the iron ring 78 approximately equals the diameter of the engagement hole 76, and the outer diameter of the iron ring 78 approximately equals the outer diameter of the magnet holder 62.

[0071] The non-reflecting member 45 is attachable to and removable from the assistant point member 35. When the target 10 is utilized for the photogrammetric analytic measurement, the non-reflecting member 45 is mounted on the assistant point member 35 in such a manner that the head 67 of the screw 66 is engaged with the engagement hole 76. In accordance with the engagement of the head 67 and the engagement hole 76, the iron ring 78 is fixedly attached to the magnet holder 62 by magnetic force of the magnet 64, so that the iron ring 78 is securely engaged with the head 67. As is apparent from FIG. 3, when the non-reflecting member 45 is mounted on the assistant point member 35, the reflecting sheet 68 and the non-reflecting sheet 74 lie on a same plane. On the other hand, when the target 10 is not utilized, the non-reflecting member 45 is removed from the assistant point member 35 by a user, by resisting the magnetic attraction between the magnet 64 and the iron ring 78.

[0072] Note that, the standard point members 31, 34, 36 and the assistant point members 32, 33 are identical to the assistant point member 35. Further, the non-reflecting members 41, 42, 43, 44 and 46 are identical to the non-reflecting member 45. Namely, each of the non-reflecting members is attachable to or removal from the corresponding standard point member. As each of the non-reflecting members is removable from the corresponding standard point member, the target 10 is handy to carry.

[0073] Further, when the target 10 is utilized, the non-reflecting members are respectively mounted on the corresponding standard point members, so that each of the reflecting sheets (68) is surrounded by the non-reflecting sheet (74). Photogrammetric analytic measurement may be performed under a situation in which an amount of luminance is not ideal for photographing, for example, in rain or at nighttime. Also, the measurement may be performed at a spot, in which a surface of a road reflects incident light. However, even if photographing is performed under such adverse conditions, areas of the standard point members 31, 34, 36 and the assistant point members 32, 33, 35 can be viewed clearly in photographed pictures, as each of the reflecting sheets is surrounded by the non-reflecting sheet.

[0074] Note that, the ratio between the diameter of the assistant point member 35 and the diameter of the non-reflecting member 45, namely an amount of area of each of the reflecting sheet 68 and the non-reflecting sheet 74, is not restricted to this embodiment. It is just required that the reflecting sheet 68 is large enough for being able to be viewed clearly in photographed pictures. Further, the assistant point member 35 and the non-reflecting member 45 may have any shape other than a circle.

[0075] FIG. 5 shows an enlarged view of the controlling unit box 20 and any other members close to the box 20, with portions broken away for clarity. FIG. 6 is a sectional view of the controlling unit box 20, taken in the direction of the arrows substantially along the line XII—XII of FIG. 5. In FIG. 6, the construction of the controlling unit box 20 is simplified.

[0076] The controlling unit box 20 comprises a battery room 83. The battery room 83 is positioned at the side of the side surface 20d. A battery 87, which supplies an electric power to the target 10, is installed in the battery room 83. The battery room 83 has an opening in the side of the side surface 20d. The opening is closed by a cover 83a. A switch 85 is unitarily mounted on the side surface 20d. The power supply to the target 10 is stopped and started by manipulation of the switch 85.

[0077] An opening 81 is formed in the top surface 20eof the controlling unit box 20. The opening 81 is closed by a cover 82. The cover 82 is made of a material through which an electric (radio) wave can transmit, for example, resin material. An antenna 88 is mounted on an inner surface of the cover 82, in such a manner that the antenna 88 is wound along the periphery of the cover 82. In the controlling unit box 20, a control board 84 and an azimuth sensor 86 are mounted. The azimuth sensor 86 and the tilt sensors 52, 54 are connected to the control board 84, whereby the operation of the sensors is controlled by the control board 84.

[0078] As is apparent from FIG. 5, the azimuth sensor 86 is placed between the assistant point member 33 and the standard point member 34. Namely, when the distance between the assistant point member 33 and the standard point member 34 is L_A, the distance between the center of the azimuth sensor 86 and the assistant point member 33 is L_A/2.

[0079] The sensing of the azimuth sensor 86 is affected by a magnetic material placed close to the sensor 86, for example, the controlling unit box 20 and the battery 87. Further, as the non-reflecting members 41, 42, 43, 44, 45 and 46 are respectively attached to the members 31, 32, 33, 34, 35 and 36 by magnetic force, by which the sensing of the azimuth sensor 86 may be affected. If the sensing of the azimuth sensor 86 is affected by the magnetic force, it is necessary to compensate azimuth data sensed by the azimuth sensor 86. However, as described above, the azimuth sensor 86 is intermediately placed between the assistant point member 33 and the standard point member 34. Namely, at the portion where the sensor 86 is placed, an influence of the magnetic force generated around the assistant point member 33 and the standard point member 34 is at a minimum. Accordingly, the sensing of the azimuth sensor 86 is not significantly affected.

[0080] When the switch 85 is turned ON to start the power supply, the sensors 52 and 54 measure the tilt angles and the azimuth sensor measures the azimuth at regular intervals in accordance with a control pulse output from the control board 84. Data measured by the sensors 52, 54 and 86 is output to the control board 84. After the data is subjected to predetermined operations, for example, compensation, in the control board 84, the data is transmitted to a receiver from the antenna 86 by wireless. The receiver is mounted on an external device, for example, a digital camera (not shown).

[0081] As described above, the antenna 88 is placed on the inner surface of the cover 82. As the cover 82, the non-reflecting member 43, 44 are made of the material through which an electric (radio) wave can transmit, the electric wave output from the antenna 88 is transmitted to the receiver, without being interrupted by the cover 82 and the non-reflecting member 43, 44.

[0082] As described above, if the receiver is mounted in a digital camera, the tilt angle data, the azimuth and image data photographed by the digital camera are saved to a recording medium. By inputting such data to a computer from the recording medium, processing of the image data can be performed in the computer more rapidly, so that a survey map can be made more accurately.

[0083] A spot of the photogrammetric analytic measurement, with the target 10 placed at a predetermined position on the spot, is photographed by a camera,.

[0084] A procedure for calculating positional data of a camera using the target 10 in this embodiment will now be explained. FIG. 7 shows a positional relationship between the target 10 and a camera 100. The target 10 is photographed by the camera 100 at a camera position M. The camera position M is defined as a back principal point of a photographing lens system of the camera 100. An optical axis of the photographing lens system is represented by a broken line O. As described above, the target 10 is provided with the three standard point members and the three assistant point members thereon. However, only the three standard point members are used for clarity of explanation. In FIG. 7, a plane, defined by the standard point members 31, 34 and 36, is the standard plane. The length of the distance between the standard point members 31, 34, and the length of the distance between the standard point members 34, 36 are the standard length L.

[0085] Note that, as is apparent from FIG. 7, an IC card 24, a readable recording medium, is mounted in the camera 100. The IC card 24 is removable from the camera 100, and attachable to other devices. Image data of the photographed pictures is recorded in the IC card 24. When the photogrammetric analytic measurement is performed, the image data recorded in the IC card 24 is loaded in a memory of the device of the photogrammetric analytic measurement as described below, after mounting the IC card 24 in the device of the photogrammetric analytic measurement.

[0086] FIG. 8 shows a photographed image at the camera position M. A two-dimensional coordinate system (Xp, Yp) is set in the photographed image. An origin of the two-dimensional coordinate system is a photographing center point CP. As apparent from FIG. 2, in the photographed image, image points of the standard point members 31, 34 and 36 respectively correspond to two-dimensional coordinates P1(Xp₁, Yp₁), P2(Xp₂, Yp₂) and P3(Xp₃, Yp₃). Note that, the two-dimensional coordinate system of the photographed image is referred to as a “photograph coordinate system”.

[0087] FIG. 9 shows a positional relationship between the photographed image photographed by the camera 100 and the target 10. A hatched area is the standard plane defined by the standard point members 31, 34 and 36. In FIG. 9, a three-dimensional coordinate system of the camera position M is shown as (Xc, Yc, Zc). An origin of the three-dimensional coordinate system corresponds to the back principal point of the photographing lens system of the camera 100. Further, the Zc-axis corresponds to the optical axis O of the photographing lens system at the camera position M. The Xc-axis is parallel to the Xp axis of the two-dimensional coordinate system of the photographed image shown in FIG. 8, and the Yc-axis is parallel to the Yp axis of the two-dimensional coordinate system. In this specification, the three-dimensional coordinate system at the camera position M is referred to as a “camera coordinate system”.

[0088] A three-dimensional coordinate system (Xs, Ys, Zs) shown in FIG. 10, which is a right-handed coordinate system, is set in order to determine the camera position M based on the photographed image. An origin of the three-dimensional coordinate system (Xs, Ys, Zs) is accorded to the standard point member 34 of the target 10. A Zs-axis is a direction from the standard point member 34 to the standard point member 36. An Xs-axis is perpendicular to the Zs-axis, lying on the standard plane. A Ys-axis is an axis on which the standard point member 34 lies, being perpendicular to the standard plane, namely being perpendicular to the sheet of drawing. As described above, the length of the distance between the standard point members 31, 34 and the length of the distance between the standard point members 34, 36 are L. Accordingly, image points of the standard point members 31, 34 and 36 respectively correspond to three-dimensional coordinates Ps1(−L, 0, 0), Ps2(0, 0, 0) and Ps3(0, 0, L). In this specification, the three-dimensional coordinate system of FIG. 10 is referred to as a “scene coordinate system”.

[0089] The camera position M is defined by a relationship between the scene coordinate system and the camera coordinate system. Namely, the camera position M is defined by a relationship of a movement distance of the origin of the camera coordinate system from the origin of the scene coordinate system to, a rotational angle of the Xs-axis around the Xc-axis, a rotational angle of the Ys-axis around the Yc-axis and a rotational angle of the Zs-axis around the Zc-axis.

[0090] A relationship between Psi(Psxi, Psyi, Pszi) (i=1, 2, 3) in the scene coordinate system and Pci(Pcxi, Pcyi, Pczi) (i=1, 2, 3) in the camera coordinate system is stated in an equation (1) as below.

[0091] $\begin{array}{cc}Pci=R\left(Psi-\Delta \right):& \left(1\right)\\ Pci=\left(\begin{array}{c}Pcxi\\ Pcyi\\ Pczi\end{array}\right);Psi=\left(\begin{array}{c}Psxi\\ Psyi\\ Pszi\end{array}\right);& \\ R=\left(\begin{array}{ccc}\mathrm{Cos}\beta \mathrm{Cos}\gamma & \mathrm{Cos}\alpha \mathrm{Sin}\gamma +\mathrm{Sin}\alpha \mathrm{Sin}\beta \mathrm{Cos}\gamma & \mathrm{Sin}\alpha \mathrm{Sin}\gamma -\mathrm{Cos}\alpha \mathrm{Sin}\beta \mathrm{Cos}\gamma \\ -\mathrm{Cos}\beta \mathrm{Cos}\gamma & -\mathrm{Cos}\alpha \mathrm{Cos}\gamma -\mathrm{Sin}\alpha \mathrm{Sin}\beta \mathrm{Sin}\gamma & \mathrm{Sin}\alpha \mathrm{Cos}\gamma +\mathrm{Cos}\alpha \mathrm{Sin}\beta \mathrm{Sin}\gamma \\ \mathrm{Sin}\beta & -\mathrm{Sin}\alpha \mathrm{Cos}\beta & \mathrm{Cos}\alpha \mathrm{Cos}\beta \end{array}\right);& \\ \Delta =\left(\begin{array}{c}\Delta X\\ \Delta Y\\ \Delta Z\end{array}\right)\left(i=1,2,3\right)& \end{array}$

[0092] In the equation (1), α is the rotational angle of the Xs-axis around the Xc-axis, β is the rotational angle of the Ys-axis around the Yc-axis, γ is the rotational angle of the Zs-axis around the Zc-axis. Further, R is a rotational matrix of α, β and γ, and Δ is a vector from the scene coordinate system to the camera coordinate system.

[0093] A relationship between Pci (Pcxi, Pcyi, Pczi) and Di (Dxi, Dyi) (i=1, 2, 3) is stated in an equation (2) as below. Note that, Pci is the three-dimensional coordinates of the standard point Pi in the camera coordinate system, and Di is the two-dimensional coordinates of the standard point in the photograph coordinate system. Futher, in the equation (2), “f” is a focal distance of the photographing lens system of the camera 10.

[0094] Dxi=f·Pcxi/Pczi Dyi=f·Pcyi/Pczi}  (2)

[0095] (i=1,2,3)

[0096] As described above, the two-dimensional coordinates of each of the standard points in the photograph coordinate are calculated from the three-dimensional coordinates in the scene coordinate system based on the equations (1) and (2). On the other hand, the two-dimensional coordinates of each of the standard points are automatically extracted from the photographed image. Accordingly, with respect to the standard points, by comparing the values of the two-dimensional coordinates calculated from the equations (1) and (2) with the values of the two-dimensional coordinates extracted from the photographed image, the values defining the camera position, namely ΔX, ΔY, ΔZ, α, β and γ can be calculated.

[0097] FIG. 17 shows a block diagram of a coordinate-calculator system, in which the photogrammetric analytic measurement, as mentioned above, is performed on the basis of the digital image data stored in the IC card 24.

[0098] As shown in FIG. 17, the coordinate-calculator system is constructed as a computer system comprising a computer 184, which includes: a central processing unit (CPU) 184A; a read-only-memory (ROM) 184B having an operating system program, constants, etc. stored therein; a random-access-memory (RAM) 184C for storing temporary data, temporary constants, etc.; a hard disk 184D for storing data resulting from calculations executed by the CPU 184A; and an input/output interface (I/O) 184E.

[0099] The computer system also comprises an IC memory card reader 186 connected to the computer 184, via the I/O 184E. The IC memory-card reader 186 is provided with a slot for receiving the IC memory card 24, and includes an IC card driver 188 for reading a given frame of digital image data and other information data.

[0100] The computer system also comprises a CD-ROM reader 197 connected to the computer 184, via the I/O 184E. The CD-ROM reader 197 is provided with a tray for receiving a CD-ROM 198, and includes a CD-ROM driver unit 199. A calculating program of the standard points of the target and positional data of the camera, as described below, is stored in the CD-ROM 198. The calculating program is read from the CD-ROM 198, and installed to the hard disk 184D.

[0101] The computer system further comprises a monitor 190 for reproducing a photographed picture based on the frame of digital image data read from the IC memory card 24 and a survey map produced by the computer 184, a keyboard 192 for inputting various command signals and various data to the computer 184, a mouse 194 for manipulating a cursor displayed on the monitor 190, and a printer 196 for printing the survey map on a sheet of paper, if necessary.

[0102] FIG. 11 shows a flowchart indicating a process of the calculating program of the camera position, according to this embodiment, executed by the CPU 184A.

[0103] After the target 10 is situated on a predetermined position in a spot for photogrammetric analytic measurement, the spot is photographed by the camera 100. In step S1000, digital image data of the spot is loaded to a system controller of a photogrammetric analytic measurement device to which this embodiment is applied.

[0104] In step S1002, a gray scale operation is performed. By the gray scale operation, brightness information is calculated based on color information of one frame color image data. After the brightness information, which is represented by 256 intensity levels, with respect to each pixel composing the digital image data, the brightness information is stored to a variable Img of array type. Namely, the brightness information of each pixel is represented by Img (j, k). Note that, “j ” and “k” correspond to the values of the two-dimensional coordinates (the photograph coordinate system) of each pixel in the digital image data.

[0105] In step S1004, binarization of the brightness information of each pixel is performed. In this embodiment, the binarization is performed using thresholding. The brightness information stored in the variable Img is compared with a predetermined threshold. If the brightness information is higher than the predetermined threshold, “1” is stored to a variable Bin of array type. If the brightness information is not higher than the predetermined threshold, “0” is stored to the variable Bin. Namely, binary information of the brightness information of each pixel is represented by Bin (j, k).

[0106] As described above, the reflecting sheet is attached on the standard point members 31, 34, 36 and the assistant point members 32, 33, 35. The standard point members 31, 34 and 36 are respectively surrounded by the non-reflecting members 41, 44 and 46 which are black. Likewise, the assistant point members 32, 33, 35 are respectively surrounded by the non-reflecting members 42, 43 and 45 which are black. Accordingly, the binary information of pixels corresponding to the standard point members 31, 34, 36 and the assistant point members 32, 33, 35 is represented by “1”, and the binary information of pixels corresponding to the non-reflecting members 41, 42, 43, 44, 45, 46 is represented by “0”.

[0107] Then, in step S1006, labeling is performed. By the labeling, a set of pixels, binary information of which are “1”, is extracted as a group. Namely, regions, brightness information of which are high, are extracted from the digital image data.

[0108] In step S1008, an area of each group is calculated. Namely, with respect to each group, a number of pixels composing the group is calculated. The number of pixels of each group and two-dimensional coordinates corresponding to the pixels are stored to predetermined variables.

[0109] Then, in step S1010, it is judged based on the number of pixels of each group if each group corresponds to the standard point members or the assistant point members. It is checked that the number of the pixels is between a greatest lower bound Th1 and a least upper bound Th2. The values of Th1 and Th2 are determined, such that a range of the pixel numbers defined by Th1 and Th2 is able to contain the pixel number of the group corresponding to the standard point members and the assistant point members, with the range being as minimal as possible.

[0110] Accordingly, a group corresponding to noise of the digital image, a pixel number of which is one or two, and a group corresponding to a wide region of high brightness, for example, a painted region of bright color on a road of the spot, are deleted.

[0111] In step S1012, a center of gravity (spectral center) of each group is calculated based on the brightness information of each group. FIG. 12 shows a flowchart of a routine for calculating the center of gravity. In step S1020, one group is selected. In step S1022, variables SumX, SumY and SumImg are initialized. A continued product of the brightness information of each pixel composing the selected group and an X component of the two-dimensional coordinate of the each pixel is stored to the SumX. A continued product of the brightness information of each pixel composing the group and a Y component of the two-dimensional coordinate of the each pixel is stored to the SumY. A sum total of the brightness information of each pixel composing the group is stored to the SumImg.

[0112] In step S1026, with respect to one pixel composing the group, a product of the brightness information (Img (j, k)) and the X component (j) is calculated to be stored to the SumX after being added to a current value stored in the SumX. Then, in step S1028, with respect to the same pixel composing the group, a product of the brightness information (Img (j, k)) and the Y component (k) is calculated to be stored to the SumY after being added to a current value stored in the SumY. In step S1030, the brightness information of the same pixel is added to a current value stored in the SumImg to be stored to the SumImg. The above mentioned processes from step S1026 through step S1030 are performed repeatedly, with respect to each pixel composing the group.

[0113] Namely, a sum total of weighted brightness information based on X component of each pixel composing the group is stored in the SumX, and a sum total of weighted brightness information based on Y component of each pixel composing the group is stored in the SumY. As described above, the sum total of the brightness information of the pixels of the group is stored in the SumImg.

[0114] In step S1032, it is checked if all pixels composing the group have been subjected to the above mentioned procedures. If the above mentioned procedures have been performed with respect to all pixels composing the group, the process goes to step S1034. In step S1034, an X component of the center of gravity is calculated by dividing the sum total of weighted brightness information based on X component of each pixel by the sum total of the brightness information of each pixel. In step S1036, a Y component of the center of gravity is calculated by dividing the sum total of weighted brightness information based on Y component of each pixel by the sum total of the brightness information of each pixel.

[0115] Then, in step S1038, it is checked if all groups in the digital image data have been subjected to the calculation of the center of gravity. If there are any groups which have not been subjected to the calculation, the procedures from step S1020 through step S1036 are repeatedly performed. If all groups have been subjected to the calculation of the center of gravity, the process goes to step S1014 of FIG. 11.

[0116] FIG. 13 is a conceptual view of each pixel of the digital image data after the calculation of the center of gravity. In FIG. 13, pixels of low brightness information are hatched by oblique lines. As is apparent from FIG. 13, groups from G1 through G13 have been extracted. The centers of gravity of the groups are respectively represented by GC1 through GC13.

[0117] In step S1014, combinations of the groups, which satisfy a positional relationship between the standard position members 31, 34, 36 and the assistant position members 32, 33, 35 of the target 10, are extracted. FIG. 14 shows the first half of a flowchart of an extraction routine of the combinations of the groups. FIG. 15 shows the latter half of the flowchart. The procedures of the flowchart will be explained with reference to FIG. 13.

[0118] In step S1040, two groups are optionally selected. In step S1042, a straight line, which lies on the centers of gravity of the selected two groups, is calculated. Then, in step S1044, groups which have a center of gravity between the centers of gravity of the selected two groups, lying on the straight line, are detected.

[0119] In step S1046, a number of the groups detected in step S1044 is checked. If there are four groups, which have the center of gravity on the straight line, including the selected two groups in step S1040, the process goes to step S1048. If there are three groups, which have the center of gravity on the straight line, including the selected two groups in step S1040, the process goes to step S1050. Otherwise, the process goes to step S1052.

[0120] For example, if the group G1 and the group G2 are selected in step S1040, there is no other group which has a center of gravity between the centers of gravity GC1 and GC2 on a straight line on which the center of gravity GC1 and the center of gravity GC2 lie. Accordingly, the process goes to step S1052. If the groups G3 and G6 are selected, there are two centers of gravity, GC4 of a group G4 and GC5 of a group G5, on a straight line on which a center of gravity GC3 of the group G3 and a center of gravity GC6 of the group G6 lie. Accordingly, the process goes to step S1048.

[0121] Also, if the groups G10 and G13 are selected, there are two centers of gravity, GC11 of a group G11 and GC12 of a group G12, on a straight line on which a center of gravity GC10 of the group G10 and a center of gravity GC13 of the group G13 lie. Accordingly, the process goes to step S1048.

[0122] If the groups G1 and G3 are selected, there is one center of gravity, GC2 of the group G2, on a straight line on which the center of gravity GC1 of the group G1 and a center of gravity GC3 of the group G3 lie. Accordingly, the process goes to step S1050.

[0123] If the groups G7 and G9 are selected, there is one center of gravity, GC8 of a group G8, on a straight line on which a center of gravity GC7 of the group G7 line. Accordingly, the process goes to step S1050.

[0124] In step S1048, the four groups are registered as a four points range which lie on a same straight line. In step S1050, the three groups are registered as a three points range which lie on a same straight line. For example, a set of the groups G1, G2, G3 and a set of the groups G6, G7, G8 are respectively registered as three points ranges. Also, a set of the groups G4, G5, G6, G7 and a set of the groups G10, G11, G12, G13 are respectively registered as four points ranges.

[0125] In step S1052, it is checked if all combinations of two groups in the digital image data have been subjected to the procedures from step S1042 through step S1050. If there are some combinations of two groups which have not been subjected to the procedures, the process returns to step S1040 and the procedures from step S1040 through step S1050 are repeatedly performed. If there is no combination of two groups which have not been subjected to the procedures, the process goes to step S1054 of a flowchart of FIG. 15. This process may be considered to be somewhat flexible. That is, the process for detection of gravity on a straight line between the centers of gravity of the two selected groups should be considered to detect centers of gravity substantially on the straight line.

[0126] In step S1054, one combination of the registered three points range and the registered four points range is selected. Then, in step S1056, it is checked if a center of gravity of a group corresponding to end point of the registered three points range and a center of gravity of a group corresponding to end point of the registered four points range are in accord, with respect to the combination selected instep S1054. Namely, it is checked if there is an identical end point shared by the registered three points range and the registered four points range, of the combination. If the end points accord, the process goes to step S1058. If the end points do not accord, namely, the registered three points range and the registered four points do not share an identical end point, the process returns to step S1054 to select another combination of the registered three points range and the registered four points range. In FIG. 13, the group G3 is an endpoint of the registered three points range composed by the groups G1, G2, G3, and also an end point of the registered four points range composed by the groups G3, G4, GS, G6. Accordingly, if the registered three points range of the groups G1, G2, G3 and the registered four points range of the groups G3, G4, G5, G6 are selected in step S1054, the process goes to step S1058.

[0127] In step S1058, two-dimensional coordinate of a center of gravity of the group corresponding to another end point, of the registered four points range, opposite to the end point which accords with the end point of the registered three points range, is stored to a variable A1 (Ax₁, Ay₁). Also, two-dimensional coordinate of the center of gravity of the group corresponding to the end point shared by the registered three points range and the registered four points range is stored to a variable A2 (Ax₂, Ay₂). Further, two-dimensional coordinate of a center of gravity of the group corresponding to another end point, of the registered three points range, opposite to the end point which accords with the end point of the registered four points range, is stored to a variable A3 (Ax₃, Ay₃). Namely, two-dimensional coordinates of the center of gravity GC6 are stored to the variable A1, two-dimensional coordinates of the center of gravity GC3 are stored to the variable A2 and, two-dimensional coordinates of the center of gravity GC1 are stored to the variable A3.

[0128] In step S1060, an angle θ, by which a vector from the variable A2 to the variable A3 is rotated in the counterclockwise direction from a vector from the variable A2 to the variable A1, is calculated. It is checked that the angle θ is larger than a greatest lower bound Th3, being smaller than a least upper bound Th4. Note that, in this embodiment, 0 degrees is set to the greatest lower bound Th3, and 180 degrees is set to the least upper bound Th4. If the angle θ is between the greatest lower bound Th3 and the least upper bound Th4, the selected combination of the registered three points range and the registered four points range meets a positional relationship between the standard point members 31, 34, 36 and the assistant point members 32, 33, 35. Accordingly, the selection of the combination is ended and the process goes to step S1016 of the flowchart in FIG. 11.

[0129] If the angle θ is smaller than the greatest lower bound Th3 or larger than the least upper bound Th4, the process goes to step S1062. In step S1062, it is checked if all combinations of the registered three points range and the registered four points range have been subjected to the procedures from step S1056 through step S1060. While there are some combinations which have not been subjected to the procedures from step S1056 through step S1060, the procedures are repeatedly performed.

[0130] Namely, until a combination of the registered three points range and the registered four points range, which meets the positional relationship between the standard point members 31, 34, 36 and the assistant point members 32, 33, 35, is found, the procedures from step S1054 through step S1062 are repeatedly performed. When one combination, which meets the positional relationship between the standard point members and the assistant point members, is found, the procedures from step S1054 through step S1062 cease to be performed and the selection of combination of the registered three points range and the registered four points range is ended.

[0131] As described above, with reference to FIG. 13, a set of the groups G1, G2, G3 is selected as the three points range, and a set of the groups G3, G4, G5, G6 is selected as the four points range, and the two-dimensional coordinates of GC6 is stored to the variable A1, the two-dimensional coordinates of GC3 is stored to the variable A2 and, the two-dimensional coordinates of GC1 is stored to the variable A3. As is apparent from FIG. 13, an angle by which the vector from the variable A2 to the variable A3 is rotated in the counterclockwise direction from the vector from the variable A2 to the variable A1, is between the greatest lower bound Th3 and the least upper bound Th4. Accordingly, the combination of the groups G1, G2, G3 and the groups G3, G4, G5, G6 is extracted as a combination which meets the positional relationship between the standard point members 31, 34, 36 and the assistant point members 32, 33, 35.

[0132] In step S1016, a calculation of a camera position is performed. FIG. 16 shows a flowchart of a calculation routine of the camera position. In step S1070, variables ΔX, ΔY, ΔZ, α, β and γ of the aforementioned equation (1) are initialized.

[0133] In step S1072, (−L, 0, 0), (0, 0, 0), (0, 0, L) are respectively substituted for Ps1, Ps2, Ps3 of the aforementioned equation (1), and suitable values are respectively substituted for the variables ΔX, ΔY, ΔZ, α, β and γ. Then, three-dimensional coordinates of the standard point members 31, 34 and 36 in the camera coordinate system, namely Pc1(Pcx₁, Pcy₁, Pcz₁), Pc2(Pcx₂, Pcy₂, Pcz₂) and Pc3(Pcx₃, PCy₃, Pcz₃) are calculated.

[0134] In step S1074, three-dimensional coordinates of Pc1, Pc2, Pc3 calculated in step S1072, are substituted for Pc1, Pc2 and Pc3 of the aforementioned equation (2). Then, two-dimensional coordinates of D1(Dx₁, Dy₁), D2(Dx₁, Dy₁), D3(Dx₁, Dy₁) in the photograph coordinate system are calculated.

[0135] In step S1076, two-dimensional coordinates of A1, A2, A3 in the photograph coordinate system, which are automatically extracted from the digital image data, are compared with D1, D2, D3 calculated from the equations (1) and (2). Namely, an absolute value of difference between D1 and A1, an absolute value of difference between D2 and A2, an absolute value of difference between D3 and A3 are calculated, and the sum S of the absolute values is calculated.

[0136] In step S1078, it is checked if the sum S is at a minimum. If the sums S is not at a minimum, the procedures from step S1072 through step S1076 are repeated after each value of ΔX, ΔY, ΔZ, α, β, γ is changed. Namely, values of ΔX, ΔY, ΔZ, α, β, γ, by which the sum S is at a minimum, are calculated by least squares method.

[0137] If it is judged that the sum S is at a minimum in step S1078, the process goes to step S1080. Instep S1080, the values of ΔX, ΔY, ΔZ, α, β, γ are registered as parameters which define the positional data of the camera, and the calculation of the positional data of the camera is ended.

[0138] The above mentioned procedures of the flowcharts shown in FIGS. 11, 12, 14, 15 and 16 is the calculation of the two-dimensional coordinates of the standard points of the target and the positional data of the camera, in the photogrammetric analytic measurement. Namely, the above mentioned procedures correspond to, for example, procedures of flowcharts from step 1201 through step 1211 shown in FIG. 12 of U.S. application Ser. No. 09/017,237. Due to the above mentioned procedures, calculation of three-dimensional coordinates of reference points (correspond to the standard points in this embodiment) of step 1211 can be accurately performed without being affected by the skill of an operator, as designation of the reference points of step 1208 is automatically performed.

[0139] As described above, according to this embodiment, the number of the assistant point members between the standard point members 31 and 34, and the number of the assistant point members between the standard point members 34, 36, are different. Further, the reflecting sheet is attached on the standard point members and the assistant point members. Accordingly, it is easy to determine the positions of the standard point members and the assistant point members, in the photographed picture, so that two-dimensional coordinates of the standard point members in the digital image are automatically calculated.

[0140] Further, according to this embodiment, as the length between the standard point members 31, 34, the length between the standard point members 34, 36, and the angle defined by the first bar 12 and the second bar 14 are determined, the three-dimensional coordinates of the standard point members in the scene coordinate system are predetermined. Further, the three-dimensional coordinates of the standard point members in the photograph coordinate system are calculated based on the three-dimensional coordinates in the scene coordinate system. Accordingly, the positional data of the camera is automatically calculated by comparing the two-dimensional coordinates of the standard point members in the photograph coordinate system with the two-dimensional coordinates of the standard point members in the digital image.

[0141] As described above, according to the present invention, the positional data of a camera is automatically calculated in the photogrammetric analytic measurement.

[0142] The present disclosure relates to subject matter contained in the following Japanese Patent Application No. 10-334320 (filed on Nov. 25, 1998), which is expressly incorporated herein, by reference, in its entirety.