METHOD AND APPARATUS FOR PHOTOGRAMMETRIC ASSESSMENT OF BIOLOGICAL TISSUE

Abstract: A remote wound assessment method and apparatus includes forming an oblique image (32) of a wound (34) and a target object such as a plate (36) containing a rectangle (37) and placed near the wound. Using a novel method of determining vanishing points (V1, V2) where photographic images of parallel lines on the target object intersect, coordinate transformations are calculated which map the oblique image (36I) of the rectangle (37I) into a normal image (36N) thereof. Using the same coordinate transformations, an oblique image (34I) of a wound (34) adjacent to the target plate (36) is mapped into a normal view (34N) thereof, allowing precise measurement of wound features. By forming two separate images (32L, 32R) of a wound (34) and target plate (36) at two different oblique inclinations, three-dimensional features of a wound may be measured.

Description: MFiTPP" . Λ N APPARATUS FOR PHOT Qg ^KElRlg Background of the XnΛanLlan A. Field of the Invention The present invention relates to methods and apparatus for visually assessing the vitality of biological tisnue. More particularly, the invention relates to mβthodα and apparatus for remote photograminetric assessment of he extent of wound damage to external tissues of a human or animal . B . B&flrcr1pti.pn of fiac- rfrunri Art Treatment of external wounds to the bodieo of humano and animals usually includes visual evalυatione of the heπling process by a medical doctor, veterinarian, nurøe or clinical specialist, Such evaluations are necessary to determine whether the wound treatment regimen io successful in edu ing the volume or damage to necrotic tissue, indicating a satisfactory rate of healing of the wound, The expertise of experienced doctors or clinicians in visually assessing changes in color, texture or area size of damaged biological tissue is an important asset in determining the ef ectiveness of medical procedures used to treat wounds. However, a problem exists in measuring the si2e of wounds?, to determine whether the area of damaged biological tieeue ie increasing or decreasing. Such quantitative assessments are particularly important in monitoring the growth rate of oun " infected by gangrene or "flesh-eating" bacteria, aince. the rapid unchecked growth of these infections can neoesαitate amputation, or result in death. A number of devices have been disclosed that relate generally to the field of assisting in visual evaluation ol human body parts, some of which utilize optical means, These ir.elulf: t e fo3lowing U.S. PatentB; K9yahaE&....3.617.635. June IB. 1974. Device For Measuring Thι» Δstua! Dimension Of An ob j ect At The H∑ ΩS_££&£: / / / an object in the field of view of an endoscope, the forward end portion of which is adapted to be inserted into a hollow portion of a living body or the like for the inspection thereof. The endoscope is provided with a fixed focal length objective lens system, a differential mechanism, and a sizing means along with adjusting means for the objective lens system, and a distance measuring means. One sizing means disclosed comprises a caliper having measured spaced pointers. A second sizing means disclosed comprises a pair of blades pivotally mounted in symmetrical relation to each other, the blades provided with cut out portions adapted for size measurement . A third sizing means disclosed comprises an iris diaphragm. The actual size of an object appearing in the field of view is measured directly by differentially coupling the differential mechanism with the sizing means and either the adjusting means for the objective lens system, or the distance measuring means for measuring the size of the object appearing in the field of view. An indicating means coupled with the differential mechanism indicates directly the actual dimension of the object regardless of the variation in the distance between the object and the forward end portion of the endoscope . Mizumoto. 4,278.077. July 14. 1981, Medical Camera System: Discloses a capsule-shaped miniature camera comprising at least one permanent magnet, an induction coil, a lamp serially connected to the induction coil and a shutter device. The induction coil induces an electromotive force when a magnetic field generated by electromagnets outside the camera acts on it. The electromotive force turns on the lamp and drives the shutter device. Landwehr. 4.370.039. January 25. 1983. Method and Apparatus For Photographically Measuring An Object: Discloses a method in which a composite photographic picture is generated through double-exposure, whereby the object, as for instance a person, is photographed on one half of the picture whereas a grid pattern is photographed on the other half, for the second exposure, person and grid pattern change position In each instance of exposure, a line pattern is pro ected onto the object under utilization of one movable or two fixed, overhung projectors. Equipment is enclosed, permitting easy change of scenery. Zoltan. 4.535.782. August 20 , 1985. Method For Determining Wound Volume : Discloses a non-contacting volume determination method for use on wounds, m which a known pattern of lines optically projected onto the volume to be determined. The image of the projection viewed from an angle other than the projection, axis, along with the image of a reference volume located near the volume to be determined, are used to accurately determine the unknown volume. Hallous. 4.564.295. January 14. 1986. Apparatus And Method For Projection Moire Topography: Discloses an apparatus and method for obtaining an improved moire fringe pattern image of an object. A beam of incoherent light is projected at the object. A grating is disposed m the path of the beam projected at the object, this grating being referred to as a projection grating. Means are provided for focusing the beam reflected from the object to obtain an image at an image plane This movable grating is referred to as a reference grating. Means are provided for recording the image at the image plane, the recording means being, for example, a photographic or video camera In accordance with an important feature of the invention, means are provided for moving the projection grating and the reference grating in synchronism. In a preferred embodiment, the projection and reference gratings are mounted in spaced relation in a movable member, and the synchronized motion of the gratings is implemented by moving the member. In a form of the disclosure, a double-pro ection set-up is used wherein two projection beams are directed at the object through first and second movable projection gratings, respectively. In this embodiment, the first and second projection gratings are moved in synchronism with the reference grating. This is preferably achieved by mounting the three gratings in appropriate spaced relationship in a common member, and moving the member. Gormley. 4.569.358. February 11. 1986. Optical Micrometry Of Skin Surfaces : Discloses an apparatus and method for measuring changes in the topology of a skin surface in which a series of casts are taken of the skin surface, using an elastomeric material such as silicone rubber in a liquid state which is later polymerized, the casts being taken of a skin surface over a period of time and then measured by means of an automatically focusing microscope to produce a quantitative model of the changes in the topology of the surface. Landwehr. 4.786.925. November 22. 1988. Photographic Contour Measurement : Discloses a method for measuring contours of a human body in which horizontal lines are projected onto a person from a plurality of overhead projectors, each projecting at a 45 degrees angle, all of the projectors have parallel optical axes, the person being photographed with the projected raster superimposed. Steinhauer et al .. 4.996.994. March 5. 1991. Apparatus For Photogrammetrically Measuring The Human Head: Discloses a process device for the photogrammetrical measurement of the human head, especially of the middle region of the face with the eye sockets, the nose, the cheek bones and the brows, which is comprised of a projector for projecting a pattern image onto the face and two cameras which can simultaneously take two pictures from two different directions of the face and head with the pattern image projected on it. The two cameras and the projector are supported by an apparatus carrier, which can be moved in a vertical direction relative to a base permanently connected in the device and in at least one horizontal direction. This minimum of one direction coincides with the optical axes of the cameras. The device will facilitate bringing both the pro ector and the cameras on the one hand and the human head on the other into the required relative position to one another as necessary to take the pictures. Of the prior art references cited above, those most applicable to the problem of optically assessing external body parts include Landwehr, 4,370,039, which discloses a method and apparatus for measuring body contours that uses a projected grid pattern on the body, and a double photographic exposure, with the body in different positions for the two exposures, Zoltan, 4,535,782, which discloses a non-contactmg volume determination method for use on wounds, in which a grid pattern is projected onto the wound, and the wound and pattern photographed, Hallous, 4,564,295, which discloses a method and apparatus for remotely measuring surface contours of an object that uses one or more grid patterns produced by projecting light through a grating onto the object surface, and a reference grating that produces a Moire pattern with light reflected off of the object, the patterns to be photographed, and Landwehr, 4,786,925, which discloses a method for remotely measuring body contours by projecting a grid pattern onto the body, and photographing the reflected pattern. Steinhauer et al . 4,996,994, discloses an apparatus for making photogrammetric measurements of the human head using two cameras photographing a grid pattern projected on the head. None of the prior art references known to the present inventors is particularly well-suited to either qualitative or quantitative assessment of wound characteristics, because the projected grid patterns or Moire patterns used in prior art methods obscure certain wound features. Moreover, none of the prior art provides means for quantitatively establishing the color boundary of a wound. The present invention was conceived of to provide a method and apparatus for photogrammetric assessment of biological tissue that overcomes limitations of the prior art. / / / / / / Ob ects of the Invention An object of the present invention is to provide a method and apparatus for photographically assessing areas of damaged biological tissue. Another object of the invention is to provide a photogrammetric wound assessment method and apparatus that does not require physical contact with a wound. Another object of the invention is to provide a photogrammetric wound assessment method and apparatus which requires minimal or no contact with a patient's body. Another object of the invention is to provide a photogrammetric wound assessment method and apparatus that provides repeatably quantifiable photographic views of a wound, in spite of varying locations of an optical image forming device such as a camera employed to form the views. Another object of the invention is to provide a photogrammetric wound assessment method and apparatus that provides means for remotely and accurately measuring the area of a wound. Another object of the invention is to provide a photogrammetric method and apparatus for wound assessment that may optionally provide a three dimensional image and measurements of wound depth in addition to wound area. Another object of the invention is to provide a photogrammetric method and apparatus for wound assessment that provides means for transforming an oblique image of a wound formed by a camera having its line of sight inclined from a normal to a wound surface, to a normal image. Various other objects and advantages of the present invention, and its most novel features, will become apparent to those skilled in the art by perusing the accompanying specification, drawings and claims. It is to be understood that although the invention disclosed herein is fully capable of achieving the objects and providing the advantages described, the characteristics of the invention described in this specification are merely illustrative of the preferred embodiment. Accordingly, we do not intend that the scope of our exclusive rights and privileges in the invention be limited to details of the embodiments described. We do intend that equivalents, adaptations and modifications of the invention reasonably inferable from the description contained herein be included within the scope of the invention as defined by the appended claims. Summary of the Invention Briefly stated, the present invention comprehends a method and apparatus for assessing the condition of biological tissue, particularly tissue which has been damaged by a wound, by remote optical means. According to the present invention, an optical imaging apparatus such as a still or video camera is used to photograph, at an oblique, arbitrary angle, a wound and target plate located physically near to the wound. The target plate contains at least one pair of lines that intersect at a known angle, and preferably contains two pairs of parallel lines that are mutually perpendicular, forming a rectangle. When photographed at an arbitrary oblique angle, the image of the target rectangle is in general a quadrilateral . The method of the present invention includes locating on the photograph intersection points of lines on the target plate, either manually or, preferably, using a combination of image search algorithms and contrast edge discerning techniques . Using a novel method of determining a pair of vanishing points where photographic images of lines that are parallel in the real world intersect, coordinate transformation are calculated which map the photographic image into its real-world shape. These coordinate transformations are then used to map the intersecting lines of a quadrilateral image into an image normal to the plane of the target plate, thus mapping the quadrilateral image into the rectangular "real world" shape of the target plate. Using the same coordinate transformations, the oblique image of an adjacent wound is mapped into a normal view thereof . Measurement of the mapped image of the target plate allows precise two- dimensional measurements of the wound features, such as its perimeter, thereby permitting precise repeatable, remote measurement of the wound size and other characteristics. In a variation of the method of the present invention that permits three-dimensional measurements of features such as wound depth, two separate photographs of a wound are made at two different oblique inclinations of the camera's focal axis with respect to the wound. Preferably, the two viewing angles are separated by about 45°. Mathematically combining two-dimensional coordinates of the two photographs thus obtained permits locating the Z-axis or depth coordinates of wound features, as well as their X-Y coordinates in the plane of the wound. Brief Description of the Drawings Figure 1 is a partially diagrammatic view of a two- dimensional wound image acquisition apparatus and method according to the present invention. Figure 2 is a view similar to Figure 1, but showing the imaging acquisition apparatus therein located at two separate oblique viewing positions relative to a wound, thereby providing three-dimensional photogrammetric wound data. Figure 3 is a perspective diagrammatic view showing the geometrical relationship between the plane of a wound image formed according to Figure 1, and an imaging control plane constructed from the image of a target object located near the wound of Figure 1. Figure 4 is a diagrammatic view of a coordinate system based upon a perspective planar view of a picture image obtained with the apparatus shown in Figure 1, with the image oriented parallel to the plane of the figure. Figure 5 is an orthogonal view of the coordinate system of Figure 4, taken in the direction 5-5. Figure 6 is a perspective view of a coordinate system having an origin selected by an operator of the apparatus according to the present invention. Figure 7 is a block diagram of an apparatus used to make photogrammetric assessments of damaged biological tissue according to the present invention. Figure 8 is a simplified flow chart showing operation of the present invention. Figure 9 is a flow chart showing a method of locating features of a target plate, according to the present invention. Figure 10 is a perspective view of a target plate according to the present invention, showing the target plate placed adjacent to an ulcer on the right leg of a patient. Figure 11 is a planar image of the damaged tissue and target plate shown in Figure 10, as imaged at an arbitrary, possibly oblique angle, by the image forming apparatus shown in Figure 1. Figure 12 is a re-oriented planar image of the damaged tissue and target plate similar to that of Figure 11, but showing the damaged tissue and target plate viewed normally thereto, the re-orientation being performed by the method and apparatus of the present invention. Figure 13 is a perspective view similar to that of Figure 10, but taken at a later date. Figure 14 is a planar view similar to that of Figure 9, but taken at the same date as Figure 13. Figure 15 is a re-oriented planar view similar to that of Figure 12, but taken the same date as Figure 13. Figure 16 is a perspective view similar to that of Figure 13, but taken at a still later date. Figure 17 is a planar view similar to that of Figure 14, but taken at the same date as Figure 16. Figure 18 is a re-oriented planar view similar to that of Figure 15, but taken at the same date as Figure 16. Figure 19 is a perspective view similar to that of Figure 16, but taken at an even later date. Figure 20 is a planar view similar to that of Figure 17, but taken at the same date as Figure 19. Figure 21 is a re-oriented view similar to that of Figure 18, but taken at the same date as Figure 19. Figure 22 is a perspective view of an alternate embodiment of a target object according to the present invention, in which the target object has an open area which may be positioned to overlie a wound. Figure 23 is a partly schematic view showing the geometric relationship between camera views, control rectangle and a box representing a limb being viewed by the method of the present invention. Figure 24A is a planar camera image showing the relationships between dimensions of the control rectangle and box of Figure 23. Figure 24B is a view similar to that of Figure 24A, but showing a different camera image. Figures 25 and 26 are schematic views illustrating object size calculations for the images of Figures 23 and 24. Description of the Preferred Embodiments Figures 1-26 illustrate photogrammetric apparatus, articles and methods according to the invention for assessing conditions of biological tissues, such as wounds in humans or animals . Referring first to Figure 1, a method of acquiring a visual image of biological tissue using an apparatus and article according to the present invention is shown. As will be made clear later in the description of the invention, visual images acquired as shown in Figure 1 are processed by methods according to the present invention to quantify various parameters of tissues imaged. As shown in Figure 1, a biological tissue image acquisition apparatus 30 according to the present invention includes a visual imaging device 31, which may be a still photographic film camera such as 35 mm camera, a film motion picture camera, video camera, or a digital still camera. As shown in Figure 1, camera 31 is used to form an image 32 at -li¬ the focal plane 33 of the camera In the example of Figure 1, the biological tissue imaged on focal plane 33 of camera 31 includes a wound 34 on the limb 35 of a patient. According to the present invention, the image acquisition method and apparatus depicted in Figure 1 includes a target plate 36 placed on or near the surface of limb 35, adjacent wound 34. As shown in Figure 1, target plate 36 has a generally rectangular, preferably square shape, and has a rectangular central area 37 concentric with the perimeter 38 of the target plate. Central area 37 of target plate 36 is preferably of a different color or darkness than the remainder of the target plate Thus, as shown in Figure 1, central area 37 of target plate 36 may be white, while the remainder of the target plate may be black. As shown in Figure 1, target plate 36 is preferably placed on or near the surface of limb 35 at a location longitudinally displaced from wound 34, but in approximate longitudinal alignment with the wound. Preferably, the longitudinal axis of target plate 36 is aligned with a longitudinally disposed medial line 39 of wound 37 Also, target plate 36 is so positioned as to make the plane of the target plate approximately tangent to the upper surface of wound 34, at a longitudinal medial line 39 on the wound As will be made clear in the ensuing description of the invention, precise alignment of target plate 36 relative to wound 34 and limb 35 is not required to obtain adequate photogrammetric measurements of the wound, according to the present invention. Preferably, however, the plane of the target plate is inclined no more than about ± 10 degrees with respect to a tangent plane to the wound. This degree of alignment is readily accomplished by visual means alone. After imaging device 31 has recorded a still or video frame containing an image frame 32 including an image 32W of a wound 34 and 36T of a target plate 36, the images are digitally processed according to the present invention, as shown in Figures 7 and 8. When imaging device 31 is a photographic still camera, employing photographic film, processing of image frame 32 requires that the exposed film be developed in a conventional fashion The developed negative image may then be digitally processed according to the present invention. Alternatively, a positive print may be made of the negative, and digitally processed in a manner to be described below. Whether a film positive or negative image frame 32 is utilized, the method of the present invention preferably utilizes an optical scanner 40 and an analog-to-digital converter (A/D) 41 to convert the analog two-dimensional image into a digital image . An X-Y scanner and A/D converter that was found satisfactory for the present invention is the Hewlett Packard brand Scan Jet TIOX scanner, which has a resolution of 600 DPI and a gray scale resolution of 24 bits As was previously stated, imaging device 31 may be a video camera rather than a still camera In this case, optical scanner 40 and A/D converter 41 are replaced by a digital frame grabber 42 which converts optical images 32 into digital electronic images that are digitally processed according to the method and apparatus of the present invention. We have found that a suitable frame grabber for use m the present invention is a MRT brand video port model, having a resolution of 780 X 480 pixels, and a gray scale resolution of 24 bits. This device is manufactured by the MRT company, located in Norway. Whether digitized image 32 is obtained by means of a still camera 31S or video camera 31V, the digitized image 43 is processed by novel methods according to the present invention, in a manner which will now be described. As shown in Figure 7, the photogrammetric assessment method according to the present invention utilizes a computer 44 operated by a user 45 to make measurements 46 of certain parameters of wound image 34. As will now be explained m detail, the method includes re-orienting images of wound 34 and target plate 36 formed at arbitrary, oblique inclinations of the optical axis of camera 31 to limb 35 into normal (perpendicular) views of the target plate and wound, thereby allowing repeatable, accurate measurements of wound features, such as wound area. The method according to the present invention of re¬ orienting or mapping an oblique camera view of wound 34 into a normal view, perpendicular to the wound, may be best understood with reference to Figures 3-7. Mapping or re¬ orienting an image of a wound 34 photographed at an "arbitrary" oblique angle includes the step of mapping an oblique view of a target rectangle 36, which oblique view is in general a quadrilateral, into a rectangle. Now target rectangle 36 is assumed to have been located in a plane tangent to the wound when the composite image 32 of the wound and target plate was originally formed, as shown in Figure 1. Therefore, mapping or re-orienting the quadrilateral image into a rectangular image results in re-oriented, normal (perpendicular) "real world" views of the wound as well as the target rectangle. Also, since the size of the target rectangle is known, the dimensions of the normal image of the wound are of a known scale. These two factors permit accurate, repeatable remote measurements of wound dimensions such as the perimeter and area of the wound, independent of camera optical axis inclination relative to the wound, and independent of camera lens focal length and distance between the camera and the wound. The mapping or re-orientation method according to the present invention and referred to above is achieved by utilizing coordinate transformations which mathematically project the camera image into a "control plane" in which the target plate quadrilateral image is transformed into a rectangle having the same "real-world" scaled dimensions and angles as the actual target plate. This operation is shown in Figure 3, which depicts a perspective view of the mapping geometry. The mapping step of the wound assessment method by photographic measurements (photogrammetric) according to the present invention may be best understood by reference to Figures 3-6 Coordinate Systems The photogrammetric wound assessment method according to the present invention preferably utilizes four coordinate systems, such as those shown m Figures 3-6. In the coordinate systems of Figures 3-6, the focal point is defined as the point of convergence of all rays of light passing both through an oblique picture plane image and a re-oriented, normal image in the control plane. The x.y.z coordinate system has its origin at the bottom left corner of the picture The positive z axis is normal to the picture plane and directed towards the focal point, the positive y axis is parallel to the vertical edge of the picture and directed toward the top left corner of the picture. The x axis is normal to the yz plane The focal axis is the axis which is normal to the picture plane and passes through the focal point The point where the focal axis intersects the picture plane is the focal center (x n , y , 0) The x' , y' . z' coordinate system is defined by finding the objects m the picture that represent the real-world pairs of parallel, mutually perpendicular lines of the target plate For each of the two pairs of parallel lines, their points of intersection are determined and called vanishing points with coordinates (x 1( y l t 0) and (x , γ-, , 0) . The horizon line is defined as the line that passes through the two vanishing points . As shown in Figure 3 , the x' . y' . z' coordinate system has its origin at the focal center The positive z' axis is normal to the picture plane and directed toward the focal point, the positive y' axis is perpendicular to and directed at the horizon line and the x' axis is normal to the y'z 1 plane . / / / Referring still to Figure 3, the X' . Y' . Z' coordinate system may be seen to have its origin at the point where the focal axis intersects the control plane. The positive Z' axis is normal to the control plane and directed toward the camera focal point, while the positive Y' axis is perpendicular to the line defined by the intersection of the picture plane and the control plane and is directed away from the camera. The X' axis is normal to the Y'Z' plane. As shown in Figure 3, the X. Y. Z coordinate system has its origin at a user selected point in the control plane. The positive Z axis is normal to the control plane and directed toward the camera focal point, while the positive Y axis is directed toward a second user selected point. The X axis is normal to the YZ plane . Procedure For Calculating Camera Orientation From A Picture Calculate the vanishing point's coordinates Referring now to Figure 4 , the vanishing points are determined by the intersection of the two sets of parallel lines using the equations: X xvp = (B : »C 1 -B 1 »C 2 )/(B 1 »A 2 -B,*A 1 ) Eqn. (1) Y xvp = • ( C 1 +A 1 «X xvp ) /B 1 ( 2 ) X yvp = ( B 4 « C, - B 3 « C 4 ) /B 3 »A, - B, »A 3 ) ( 3 ) Y yvp = • ( C 3 +A 3 »X yvp ) /B 3 ( 4 ) where A 2 , B λ , (ft and A 2 , B , C-, are the coefficients of one set of parallel lines and A 3 , B 3 , C 3 and A , B 4 , C, are the coefficients of the other set of parallel lines. Calculate the rotation angle, θ The rotation angle equals the angle between the x-axis and the vector running from the y-axis vanishing point (yvp) to the x-axis vanishing point (xvp) . θ = Tan-1( (Y xvp -Y yvp ) / (X xvp -X yvp ) (5) Convert x. y coordinates to x' . y' coordinates: Referring still to Figure 4 and defining the angle between the positive x axis and the horizon line as the rotation angle θ, then to transform a point in the x, y coordinates system to the x', y' coordinates system. x' = (x-x 0 ) «cosθ+ (y-y 0 ) «sinθ (6) y' = (y-y 0 ) *cosθ- (x-x 0 ) »sinθ (7) Calculate the distance from the focal center to the horizon line, h: h = y 1 of v, (8) or h = y' of v, (9) Calculate the focal length, f, in the picture coordinate system: If xft, yft and z' 2 , y' 2 are the coordinates of the two end points of one of the known lengths (segment A of target rectangle) then, U 12 , which is defined as: (x' (h-y' , ) -x ' ? / (h-y' ft ) (10) V 12 , which is defined as: (y' ,/ (h-y' 3 ) -y' 2 / (h-y' 2 ) ) r (11) If x' 3 , yft and x' 4 , yft are the coordinates of the two end points of the other known length (segment B of target rectangle) then, υ 3A , which is defined as: (x' ,/ (h-y' ft -x' (h-yft) ) (12) V j4 , which is defined as: (y' 3 /(h-y'ft - y\/ (h-y ■ ) fl3) If D 12 is the length of segment A and D t4 is the length of segment B in the camera coordinate system, then And the focal length in the picture coordinate system is determined by: f = h- ( (U 12 -R«U 34 )/(R»V 34 -V 1 )-1) 4 ' (15) Calculate the distance from the camera to the control plane, Z in the real world coordinate system: Z. = ( (f»h»D ]2 ) 2 /H 2 «h »U 12 +H 4 «V 12 ) ) (16) Where H = (h 2 +f 2 ) Calculate the pitch angle, α in the real world coordinate svstem: The pitch angle a is defined as the angle between the focal axis and the control plane and is determined by: α = Tan(f/h) (17) Converting x' , y' coordinates to X' . Y' coordinates: X' = M«x'/(h-y') where M = Z.«H/f (18) Y- = N»y'/( -y') where N = Z.«H /(f-h) (19) / / / Calculate the yaw angle Φ: The yaw angle φ is defined as the angle of rotation between the X' , Y' , Z' coordinate system and the X, Y, Z coordinate system. Calculate the X' , Y' coordinates of two points in the picture and determine the equation of the line that passes through both points from the equation: A«X'+B»Y'+C = 0 (20) where A = Y' -Yft, B = Xft-Xft and C = X ' .»Y' x -X' ! »Yft If the line's equation in X, Y, Z space is defined as: D«X' +E»Y' +F = 0, (21) then the angle between the two lines and thus the yaw angle is : φ = Tan <A«E-D»B)/(A«D+B»E) ) (22) Converting X' . Y' coordinates to X. Y coordinates: X = (X' -Xft) »cosφ+ (Y' 0 -Y' 0 ) »sinφ (23) Y = (Y' -Y' 0 ) «cosφ- (X' -X' f ) »sinφ (24) Note that the Z component is equal to the Z' component and is defined as 0. Calculate the camera position XnYr. in the real world coordinate system: X ' = 0 (25) Y ' = -f»Z c /h (26) X. = X 1 0 «cosφ+(Y T -Y'ftsin (27) Y. = (Y '-Yft) 'cosφ-Xftsinφ (28) 3D Application: A variation of the photogrammetric method according to the present invention allows measurement of depth as well as surface dimensions of wound features. Thus, as shown in Figure 2, two pictures may be taken from two separate oblique positions (X.,, Y. : , Z.ft and (X._, Y„, Z. 2 ) viewing a limb 35, preferably separated by between about 15° and 165°. The three dimensional "3-D" coordinates of a point (Point A in Figure 2, for example) visible in both pictures, can then be determined by triangulation, according to the following procedure. Procedure: First, the X, Y Coordinates of the Point A as projected onto the control plane from the point of view of both pictures are located and identified as (Xj, Y : , 0) and (X t , Y-, 0) . Second, a line is constructed from focal point 1 to the projected point 1. Similarly, a line is constructed between focal point 2 and projected point 2. The coordinates of the Point A are the coordinates of the intersection of the two lines. Because of errors in calculations, the two lines may not intersect . The point where they come closest to each other would then be considered the coordinates of Point A. Proofs Proof for equation 19: Define Y' as a function of its projected point in the picture (x' , y') . See Figure 5. Y' = Z 0 # (tan (α+φ) -tanα) Y' = Z t ,« ( (tanα+tanφ) / (l-tanc.»tanφ) -tanα) ) Since tanα = f/h, tanφ = y'./f: Y 1 = z n »( (f/h+y'/f)/(l-y'/h) -f/h) Y' = Z 0 « ( (f+y »h/f)/(h-y' ) - (f-f»y' /h) /h-y' ) ) Y' = z 0 »y • ( f+f) / ( (h-y 1 ) »f»h) Define H = Define N - Z.»H 2 /(f»h) : yielding equation 19 Y' (x' , y ) = N»y'/ (h-y' ) Proof for equation 18: Define X' as a function of its projected point in the picture (x' , y' ) : X'/Z o = x'/z (38) z = f»cosα-y' *sinα (39) Since cosα = h/H, sinα = f/H: (40) X' = Z 0 «H«x'/ ( (h-y' ) »f) (41) Define M = Z 0 »H/f: (42) yielding equation 18 X' (x' , y' ) = M*x'/ (h-y' ) (43) Camera Orientation From A Picture Proof for equation 16: Define the height of the camera off the control plane, Z 0 : (X' 1 -X' 2 ) 2 +(Y' 1 -Y' 2 ) r = D 12 2 (44) (M«x' 1 / (h-y' ; ) -M«x ' 2 / (h-y' : ) ) 2 + ( N• y ' / (h-y 1 ) i - N»Yft/ (h-y'ft ) ; = D 12 2 (45) M « (x' (h-y' ft -x' 2 /(h-y' . ,) ) ^ +N » (y ' (h-y ' ft -y ' 7 (h-y 1 2 ) ) - = D ]: 2 (46) Define U 12 = (x ' ,/ (h-y' γ ) -x ' r / (h-y' 2 ) ) : : (47) Define V^ = y' 7 (h-y ' , ) -y ' ,/ (h-y ' 2 ) ) : (48) M »U 12 +N 2 «V 12 = D 12 (49) Substituting for M and N: (Z 0 »H/f) 2 *U 1 +(Z 0 «H 2 / (f«h) ) 2 *V 1 = D ιr (50) Reducing: (H «h 2 »U 12 +H"»V l;) ) / (f «h ; ) = D 17 /Z n :' (51) Z () 2 = (f»h«D 1 )7 (H-»h J »U 12 +H 4 »V 1? ) (52) yielding equation 16 Z c = ( (f oh »D 2 ) 2 /H 2 -h 2 »U 12 +H 4 »V 12 )) 1/2 (53) Proof for equation 15 : Define focal-length f: Equating the Z 0 2 equation for D 12 and D J4 : (f•h*D 1? ) 7 (H »h «U 12 +H"«V, 2 ) = (f•h»D 34 ) 7 (H 2 *h 2 »U 3 ,+H '! »V( 3 5&) Reducing: (D 12 2 /D j4 2 ) • (H ;' »h »U 3 ,+H' ; «V 34 ) = H 2 »h 2 «U 12 + H 4 »V 12 (55) Define R = (D ]2 /D 34 2 ) : Substituting for H and reducing: R » (h 2 »U 4 + (h 2 + fft »V 3 ft = h «U 12 + (h + f 2 ) »V 12 (56) R » (h 2 »U 34 +h :' «V + f «V 3 J = h 2 »U ι; +(h-«V 12 + f 2 ) •V 12 (57) f « (R»V 3< -V 12 ) = h 2 » (U 12 +V 1; -R«U 34 -R»V J4 ) (58) f = h» ( (U 12 +V 1 -R»U 4 -R»V 3 ,)/(R»V 34 -V 12 ) ) h (59) yielding equation 15 f = h* ( (U 12 -R»U 3 / (R«V 34 -V 12 ) -1)" (60) Operation of the Invention Referring now to Figures 7 through 9, in addition to Figure 1, a still or video picture of a wound 34 is photographed along with an adjacent target plate or control rectangle 36. Preferably, the camera focal length and distance from the wound are adjusted to maximize the size of the wound image in the camera focal plane, while still maintaining the control rectangle within the field of view of the camera . Also, the photograph of the wound preferably is time and date stamped, and preceded or succeeded with still or video photographic frames showing the patient's face, identification number, etc., to provide positive documentary evidence of the wound measurement procedure. After a photograph of a wound has been taken as shown m Figure 1 and described above, the photograph is scanned, digitized and input to computer 44. Then, as shown m Figures 7 and 8, an operator 45 selects a wound image 341 as shown in Figure 10, for display on monitor 47. As shown m Figures 3-6 and described above, a pair of intersecting lines forming a first corner of image 371 of target plate rectangle 37 in digitized image 43 is located visually on the screen of monitor 47, whereupon a cursor is placed on the corner by an operator, and the X-Y coordinates of the first corner in the monitor frame of reference entered into computer memory by the operator. These coordinates locate the origin of the X-Y coordinate system of the control plane. In an exactly similar manner, the coordinates of the remaining three corners of control rectangle image 371 are entered into computer 44. Then, the coordinates of the control rectangle in the picture plane are transformed into the "real world" coordinates of the rectangle in the control plane, by the coordinate transformation indicated m equations 6, 7, 18, 19, 23, 24, thus forming a normal, perpendicular view of the control rectangle . With the entry into computer 45 of the known lengths of the sides of control rectangle 37, the image of the rectangle in the control plane will be of known scale. Transforming a quadrilateral image of target plate rectangle 37 such as image 371 in Figure 11, into a "real world," i.e., normally viewed rectangle 37N, as shown in Figure 12, of known dimensions in the control plane, entails derivation of rotation angle θ, camera focal length f, distance from the camera to the control plane, pitch angle α, and yaw angle φ, according to equations (5) , (15) , (16) , (17) and (22) . These four parameters 34 are then stored in memory of computer 44, and used to calculate the actual shape 34N and dimensions of a wound 34 from the coordinates of its oblique image 341 in the picture plane. To perform this operation, operator 45 moves a cursor onto the desired wound features (Figure 12) , and enters the monitor or picture plane coordinates of the features into the computer, whereupon equations (23) and (24) are used to calculate the "real world" coordinates of the features. Thus, by tracing the perimeter of the wound image with the computer monitor cursor, a normal view of the wound perimeter may be calculated and displayed according to the method of the present invention, as shown in Figure 12. Also, since the scale of the normal view of a wound, such as the wound shown in Figure 12, is known, the outline or perimeter of the normal wound image may be used to calculate the exact area of the wound, using conventional area calculation algorithms such as the polygon algorithm described below. Area of polygon P 1 ,P : ,»»«O n Area = A(x,•y.+χ 2 »y 3 +««•x r , -1 »y n +x Il »y 1 -y,»x . -y.»x,-»»» -y n -lXn-Yn β Xl) In the preferred embodiment, the novel method of the present invention utilizes a contrast edge tracking algorithm to automatically locate the corners of target rectangle 36 in the picture plane. This algorithm, depicted in Figure 9 functions as follows: 1. Locate the steepest edge in entire picture and call this location start edge. 2. Set the current edge to start edge. 3. From the current edge's neighboring locations find the steepest edge and set the current edge to that location. 4. Repeat No. 3 until current edge equals start edge. Figures 13, 16 and 19 show photographic images 32 of different wounds, analogous to Figure 10. / / / Figures 14, 17, and 20 show digitized oblique planar images of the wounds of Figures 13, 16 and 19, respectively, analogous to Figure 11. Figures 15, 18, and 21 show digitized planar images similar to those of Figures 14, 17, and 20, but with the images re-oriented to normal views, similar to that of Figure 12. Example of 2D Measurement From A Picture Given A digital picture of resolution 580X480 (See Figure 24A) . The picture coordinates (x, y) of the four corners of a 1 inch square as it appears in the picture. 1. 528.657, 173.343 2. 464.753, 145.661 3. 421.317, 192.106 4. 483.920, 218.080 Calculate the line equation coefficients Using the general equation of a line A*x+B*y+C=0 (61) where A = Yι-Y 2 B = χ 2 -x,, C = Xι # y;-yi # x 2 (62) determine the four lines that form the quadrilateral . 1-2. A = 173.343-145.681=27.662, B = 464.753-528.657=- 63.904, C = 528.657*146.681-173.343*464.753=-3546.399 2-3. A = -46.425, B = -43.436, C = 27903.958 3-4. A = -25.974, B = 62.603, C = 1083.124 4-1. A = 44.737, B = 44.737, C = 31405.374 Calculate the vanishing points' coordinates (See Figure 26A) The vanishing points are determined by the intersection of line 1-2 with line 3-4 and the intersection of line 2-3 with line 4-1 using equations 1 and 2: x = (Bo-Ci-Bj β Cy) / (B,«A,-B 2 »A 1 ) (63) y = (d+Aι»x) /B, (64) Intersect line 1-2 and line 3-4 x = (62.603*-3546.399- [-63.904] •-1083.124) / (-63.904*- 25.974-82.603*27.882) = 4051.535 y = - (-3546.399+27.662*4051.535) /-63.904 = 1698.284 Interest line 2-3 and line 1-4 X = (44.737*27903.958- [-43.436] *-31405.374) / (- 43.436*44.737-44.737*-46.425) = -885.880 y = - (27903.958+ [-46.425] *-665.880) /-43.436 = 1587.880 Reorder the points (See Figure 26A) The control rectangle corner points must be ordered such that the first point is the point that forms the largest triangle with the vanishing points and the other points are ordered counter-clockwise around the rectangle. By visual inspection the correct order is 1. 464.753, 145.681 2. 528.657, 173.343 3. 483.920, 218.080 4. 421.317, 192.106 Calculate the rotation angle (See Figure 26A) Define the x-axis vanishing point as the intersection of line 1-2 and line 3-4 and define the y-axis vanishing point as the intersection of line 2-3 and line 4-1. Intersection of line 1-2 and line 3-4 (x-axis vanishing point) x xv . = 4051.535 y xvp = 1698.284 Intersection of line 2-3 and line 4-1 (y-axis vanishing point) x yvp = -885.880 y yvp = 1567.880 The rotation angle equals the angle between the x-axis and the vector running from the y-axis vanishing point (y^,) to the x-axis vanishing point (x vp ) . Using equation 5: θ = Tan "1 ((y X vp-y V vp)/(Xχvp-Xyvp) ) θ = Tan _1 ( (1698.254-1567.880. /4051.535- [-885.880] ) ) θ = 0.027 radians = 1.547 degrees Calculate the distance from the focal center to the horizon line, h (See Figure 26B) / / / Using equation 8: h = yftv Using equation 7: yft vp = (y xvp -y 0 ) β cosθ- (x xvp -x 0 ) *sinθ X n , y 0 = center of the picture = 560/2,480/2 = 280.240 h = (1698.284-240*cos (0.027) - (4051.535-280) *sin(0.027) = 1355.933 Calculate the focal length, f Using equations 6 and 7: x = (x-x 0 ) *cosθ+ (y-y n ) *sinθ y = (y-yo) # cosθ- (x-x 0 ) *sinθ x = (x-280) *cos (0.027)+(y-240)*sin(0.027) y = (y-240) *cos (0.027) - (x-280)*sin (0.027) x = (464.753-280) *-0.39862+ (145.651-240) *-0.91711 y = (145.681-240) *-0.39862- (145.681-280)* (0.91711 X = 182.139 y = -99.272 X = (528.657-280) *-0.39862+ (173.343-240*-0.91711 y = (173.343-240) *-0.39862- (528.657-280) *-091711 X = 246.767 y = -73.346 X = (483.920-280) *-039862+ (218.080-240) *-0.91711 y = (218.080-240) *-039862- (483.920-280) *-0.91711 X = 203.254 y = -27.417 X = (421.317-280) *-0.39862+ (192.106-240) *-0.91711 y = (192.106-240) *-0.39862- (421.317-280) *-091711 X = 139.973 y = -51.692 Using equation 10: U 13 = (x'7(-y ^-.^/(h-yft) ) 2 U 13 = (182.139/(1355.933- [-199.272] ) -203.254/ (1355.933- [-27.417] ) ) 2 = 0.00047 Using equation 11: v 13 = (y7(h-y 1 )-y 3 /(h-y 3 )) 2 V 13 = (-99.272/(1355.933- [-99.272] )- [-27.417]/ (1355.933- [-27.417] ) ) 2 = 0.00234 Using equation 12 : U,, = (xft/(h-y ft -xft/(h-yft,) ) 2 U v . = (246.757/ (1355.933- [-73.346] ) -139.973/ (1355.933- [-51.692] ) )- = 0.00536 Using equation 13 : V 24 = (y ft -yft/ (h-y 1 ,) )' V 4 = (-73.346/ (1355.933- [-73.346] - [-51.692] /(1355.933- [-51.692] ) )" = 0.00021 D 3j , D 24 is the real-world length of the diagonals of the control rectangle = 1.414. therefore, using equation 14, R = = 1 Using equation 15 the focal length is then f = h* ( (U 13 -R*U. 4 ) / (R*V, 4 -V 13 ) -1) f = 1355.933* ( (0.00047-1*0.00536) / (1*0.00021-0.00234) I) 4 = 1543.487 Calculate the distance from the camera to the control plane, Z o Using equation 16 : Z 0 = ( (f*h*D 1 )7(H 2 *h ^ *U 13 +H 4 *V 13 ) 7 where H - (h 2 +f 7 H = (1355.933 2 +1543.487 2 7 = 2054.484 Z 0 = ( (1543.487*1355.933*1.414)7 (2054.484 *1355.933 * 0.00047+2054.484''*0.00234) 7 Z 0 = 13.898 Calculate the pitch angle, α The pitch angle is the angle between the picture plane and the control plane and is determined by using equation 17: α = Tan ~] (f/h) α = Tan "3 (1543.487/1355.933) α = 0.850 radians = 48.701 degrees Converting x' , y' coordinates to X' . Y' coordinates Using equation 18 : X' = M*x' / (h-y' ) where M = Z 0 *H/f M = 13.698*2054.484/1543.487 = 18.499 X' = 18.499*x'/ (1355.933-y' ) Using equation 19: Y' = N * y' / (h-y' ) where N = Z 0 *H 2 /(f * h) N = 13.898*2054.484 : / (1543.487*1355.933) = 28.030 Y' = 28.030*y'/ (1355.933-y' ) Calculate the yaw angle Φ (See Figure 26D) Using equations 18 and 19, with appropriate substitutions : X' = 18.499*x' / (1355.933-y' ) Y' = 28.030*y'/ (1355.933-y' ) Xft = 18.499*182.139/ (1355.933- [-99.272] ) = 2.315 Yft = 28.030*-99.272/ (1355.933- [-99.272] ) = -1.912 Xft = 18.499*246.767/ (1355.933- [-73.346] ) = 3.194 Y'. = 28.030*-73.346/ (1355.933- [-73.346] ) = -1.438 Using equation 22 : φ = Tan -1 ( (Yft-Yft) /Xft-X'ft ) φ = Tan "1 ( (-1.438- [-1.912] ) /3.194-2.315) ) = 0.495 radians = 28.361 degrees Converting X 1 . Y' coordinates to X, Y coordinates X = (X- -X' ft *cosφ+ (Y' -Y' ft *sinφ Y = (Y' -Y' ,) *cosφ- (X' -X 1 ft *sinφ X = (X' -2.315) *cos (0.495+ (Y' - [-1.912] ) *sin(0.495) Y = (Y' - [-1.912] ) *cos (0.495) - (X 1 -2.315) *sin(0.495) Calculate the camera's position X f , Y : X' (l = 0 Using equation 26 : Y'o = -f*Z 0 /h = 1543.487*13.898/1355.933 = -15.820 X 0 = (0-2.315*cos(0.495) + (-15.820- [-1.912] ) *sin(0.495) -8.644 Y 0 = (-15.820- [-1.912] ) *cos (0.495) - (0-2.315) *sin(0.495) -13.338 Converting x. y coordinates to X. Y coordinates (see Figure 26) Now that we have the camera's position and orientation we can calculate the locations of any point in the picture as projected onto the control plane relative to the control rectangle. For an example we will calculate the locations of the control rectangle ' s corners . Using equations 6 and 7 with appropriate substitutions: x' = (x-280) * cos (0.027) + (y-240) * sin(0.027) y' = (y-240) *cos (0.027) - (x-280) *sin(0.027) then, using equations 18 and 19 with appropriate substitutions : X' = 18.499*x' / (1355.933-y' ) Y' = 28.030*y'/ (1355.933) -y 1 ) then, using equations 23 and 24 with appropriate substitutions : X = (X 1 -2.315) *cos(0.495) + (Y 1 - - 1.912 ) *sin ( 0.495 ) Y = (Y 1 - -1.912) *cos (0.495) - (X' -2.315) *sin(0.495) x' j = (464.753-280) *cos (0.027) + (145.681-240) *sin(0.027) = 182.139 yft = (145.681-240) *cos (0.027) - (484.753-280) *sin(0.027) = -99.272 X 1 ! = 18.499*182.139/ (1355.933- [-99.272] ) = 2.315 Yft = 28.030*-99.272/ (1355.933- [-99.272] ) = -1.912 X ! = ( 2.315 - 2.315 ) *cos ( 0.495 ) + ( - 1.912 - [ - 1.912] ) *sin(0.495) = 0 Y, = ( - 1.912 - [ - 1.912] ) *cos ( 0.495 ) - ( 2.315 - 2.315) *sin(0.495) = 0 xft = (528.657-280) *cos (0.027) + (173.343-240) *sin(0.027) = 246.767 yft = (173.343-240) *cos (0.027) - (528.657-280) *sin(0.027) = -73.346 X' 2 = 18.499*246.767/(1355.933- [-73.346] ) = 3.194 Yft = 28.030*-73.346/ (1355.933- [-73.346] ) = -1.435 X 2 = (3.194 -2.315) *cos ( 0.495 ) + ( - l .438- [ - 1.912] ) •sin(0.495) = 0.999 Y 2 = (-1.438- -1.912) «cos (0.495) - (-1.912- 2.315) *sin(0.495) = 0.000 xft = (483.920-280) *cos (0.027) + (216.080-240) *sin(0.027) - 203.254 yft = (218.080-240) *cos(0.027) - (483.920-280) *sin(0.027) - 27.417 Xft = 18.499*203.254/ (1355.933- [-27.417] ) = 2.718 Yft = 28.030*-27. 17/ (1355.933- [-27.417] ) = -0.556 X 3 = (2.718-2.315) *cos (0.495) + ( - 0.556 - [ - 1.912] ) *sin(0.495) = 0.999 Y 3 = ( -0.558- [-1.912] ) *cos (0.495) - (2.718- 2.315) *sin(0.495) = 1.002 xft = (421.317-280) *cos(0.027)+ (192.106- 240) •sin (0.027) = 139.973 yft = (192.106-240) *cos (0.027) - (421.317-280) *sin(0.027) = -51.692 X' 4 = 18.499*139.973/ (1355.933- [-51.692] ) = 1.840 Yft = 28.030*-51.692/ (1355.933- [-51.692] ) = -1.029 X 4 = (1.840-2.315) *cos ( 0.495 ) + ( - 1.029 - [ - 1.912] ) *sin(0.495) = 0.001 Y„ = ( -1.029- [-1.912] ) *cos ( 0.495) - ( 1.840- 2.315) *sin(0.495) = 1.003 X ? , Y 2 = 0.999, 0.000 X 3 , Y, = 0.999, 1.002 X 4 , Y 4 = 0.001, 1.003 The expected values for the above four coordinates are X 2 , Y 2 = 1, 0 X 3 , Y 3 = 1, 1 x 4 . Y. = o, l The differences between the expected result and the actual result is due to the inaccuracy in determining the picture coordinates of the control rectangle corners, assumptions made about the location of the focal center, and image distortions. Example of 3D Measurement From Two Pictures Given: Two digital pictures of the same subject (a control rectangle and a box) (See Figures 24A and 24B) Calculate the camera's position and orientation (See 2D Example) The first camera's coordinates (X cl , Y c , , Z cl ) are -8.644, -13.338, 13.898 The second camera's coordinates (X c2 , Y c2 , Z_J are 5.797, - 7.466, 9.423 Calculate the projection points The picture coordinates (x, y) of a corner of the box as it appears in the first picture and the real-world coordinates (X, Y, O) of the corner projected onto the control rectangle plane are : (x. y) (x, Y, o) 147,302 -2.109, 4.929, 0 The picture coordinates (x, y) of the same corner of the box as it appears in the second picture and the real-world coordinates (X, Y, O) of the corner projected onto the control rectangle plane are: (x, y) (X, Y, O) 286, 390 -4.061, 5.533, 0 Calculate the X. Y. Z coordinates The projection line for a point in a picture passes through the camera focal point and through that point's projection onto the control plane. If a point in one picture matches a point in another picture taken from a different perspective, then the intersection of their projection lines occurs at the 3D coordinates of the object. Aj = Z cl * (Y c2 -Y ) -Z c2 * (Y cl -Yft (65) A x = 13.898* (-7.466-5.580) -9.423* (-13.338-4.911) = - 8.530 B 2 = Z cl * (Y ? -Yι) (66) Bi = 13.898* (5.580-4.911) = 8.394 Ci = Z t2 * (X cl -Xft -Z cl * (X c2 -X 2 ) (67) C α = 9.423* (-8.644- [-2.124] ) -13.898* (5.797- -4.035) = - -198.586 D 2 = -13.698* (-4.035- [-2.124] ) = 26.129 E, = (Y cl -Yft • (X c2 -X 2 ) - (X cl -Xft • (Y, -Y 2 ) (69) Ej = (-13.338-4.911) • (5.797- [-4.035] ) - (-8.644- [-2.124] ) • (-7.466-5.680) = -265.025 Fi = (Y C . 1 -Y 1 ) * (X -X 1 ) - (X rl -X,) * (Y,-Yft (70) Fi = (-13.338-4.911) • (-4.035- [-2.124] ) - (-8.644- [-2.124] ) • (5.680-4.911) = 39.604 t, = 8.530*8.394- [-198.586] *26.129- [-265.025] *39.604 = 0.145 A 2 = 8.530 B 2 = -9.423* (5.580-4.911) = -5.891 d = 198.588 D 2 = Z c2 * (X 2 -X,) (75) D 2 = 9.423* (-4.035- [-2.124] ) = -18.394 E 2 = 265.025 F 2 = (X 0? -X ? ) * (Y 2 -Y 1 - (Y 02 -Y 2 ) * (X 2 -X 1 ) (77) F 2 = (5.797- [-4.035] ) • (5.580-4.911) - (-7.466-5.580) • (- 4.035- [-2.124] ) = -19.420 t, = -A 2 *B 2 -C 2 *D 2 -E 2 *F 2 (78) t 2 = -8.530*-5.691-198.586*-18.394-285.025*-19.420 = 0.081 X = (Xι+X 2 + t ! * (X c2 -X 2 ) +t 2 * (X cl -X : ) ) /2 (79) X = (-2.124 + -4.035 + 0.145* ( 5.797- [-4.036] ) +0.081* (-8.644- [-2.124] ) )/2 = -2.635 Y = (Yi+Y. + t ! * (Y c2 -Y 2 )+t 2 * (Y cl -Yj) ) /2 (80) Y = (4.911 + 5.580 + 0.145* (-7.466-5.580)+0.081* (-13.338- 4.911) ) /2 = 3.549 Z = (t ! *Z c2 + t 2 *Z cl ) /2 (81) Z = (0.145*9.423+0.081*13.895) /2 = 1.246

Claims: - 32 - What is claimed is : 1. A method for remotely determining the normal view shape of a region of damaged biological tissue such as a wound comprising; a. positioning a target object having known contrasting visual features and dimensions proximate said wound, b. forming an optical image of said target object and said wound using an optical imaging system that has an optical axis that may be inclined at an oblique angle to a normal to the plane of said wound, c. determining a coordinate transformation that maps an oblique image of said target object into a normal view thereof, d. mapping by said coordinate transformation said image of said target object into a normal image thereof, and e. mapping by said coordinate transformation said wound image, thereby obtaining a normal image of said wound. 2. The method of Claim 1 further including the step of determining the size of selected features of said wound. 3. The method of Claim 2 wherein said size determination step includes measuring the length of features of said normal image of said target object, dividing the measured length by the real world length of the corresponding feature of the actual target object to obtain a ratio k, and multiplying the length of selected wound features in said transformed image by k. 4. The method of Claim 1 wherein said target object is further defined as having thereon at least one, first optically imageable pair of intersecting lines. 5. The method of Claim 4 wherein said first pair of intersection line is further defined as defining a first corner of a polygon. 6. The method of Claim 5 wherein said target object is further defined as having thereon a second pair of optically imageable intersecting lines, said second pair of intersecting lines defining a second corner of said polygon. - 33 - 7. The method of Claim 6 wherein said polygon is further defined as being a quadrilateral . 8. The method of Claim 7 wherein said quadrilateral is further defined as being a rectangle. 9. The method of Claim 4 wherein said first pair of optically imageable lines on said target object is further defined as defining on one side thereof an area of one color and on the opposite side thereof an area of a contrasting color. 10. The method of Claim 4 wherein said target object is further defined as having thereon a transparent region. 11. The method of Claim 4 wherein said target object is further defined as having through the thickness dimension thereof a perforation. 12. A method for optically assessing characteristics of a wound comprising; a. positioning a target object comprising a thin sheet having formed thereon optically contrasting features proximate said wound, b. forming an optical image of said target object and said wound using an optical imaging system that has an optical axis that may be inclined at an oblique angle to a normal to the plane of said wound, c. determining a coordinate transformation that maps an oblique image of said target object into a normal view thereof, d. mapping by said coordinate transformation said image of said target object into a normal image thereof, e. mapping by said coordinate transformation said wound image, thereby obtaining a normal image of said wound, and f. observing features of said normal image of said wound. 13. The method of Claim 12 wherein the plane of said target object is positioned approximately coplanar with a tangent plane to said wound. / / / - 34 - 14. The method of Claim 13 wherein said optically contrasting features of target object are further defined as including a rectangular shape . 15 The method of Claim 14 wherein said coordinate transformation is further defined as mapping a quadrilateral image of said rectangle into a rectangle with sides m the same ratio as said rectangular shape of said target object 16. An apparatus for remotely assessing characteristics of a wound comprising; a. a target object having quantifiable visual features for placement proximate said wound, b. means for forming an image of said target object and said wound, c. means for performing a coordinate transformation that maps an oblique image of said target object into a normal view thereof, and d means for mapping by said coordinate transformation said image and said wound image of said target object into normal images thereof, 17. The apparatus of Claim 16 further including means for measuring selected features of said wound image. 18. The apparatus of Claim 16 wherein said target object is further defined as having thereon at least one first optically imageable pair of intersecting lines 19. A target object for imaging simultaneously with an adjacent wound, said target object having thereon at least one, first optically imageable pair of intersecting lines, whereby an obliquely formed image of said target object and said wound may be transformed into normal views of said target object and said wound. 20. The target object of Claim 19 further defined as having thereon a second pair of intersecting lines, said first and second pair of intersecting lines defining a polygon. 21. The target object of Claim 20 wherein the interior of said polygon is transparent.

Inventors: DERKSEN GRANT D, TAYLOR GEOFFREY L, DERKSEN, GRANT, D., TAYLOR, GEOFFREY, L.
Applicants: DERKSEN GRANT D, TAYLOR GEOFFREY L, DERKSEN, GRANT, D., TAYLOR, GEOFFREY, L.
Application year: 1996
Main classification: A61B
IPC codes: A61B5, G01C11