BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention related to computer modeling.

2. Scope Of Related Prior Art

Computer models typically are shown visually to a user in which surfaces are displayed at desired locations and orientations. The user may interact with the model to cutaway, move and otherwise manipulate portions of the model. However, much of the this is done without the aid of tactile feedback.

Typically, when a sculptor is modifying or manipulating a physical model, there is force feedback. The computer model simulation lacks this element.

There has been some limited use of force, or haptic, feedback. In molecular docking studies, a robotic arm was used to supply molecular interaction forces as described in "Project GROPE: Haptic Displays for Scientific Visualization," by F. P. Brooks, P. M. Ouh-Young, J. J Batter, P. J. Kilpatrick, Proceedings of SIGGRAPH '90, pp. 177-186 (August 1990).

In another application, a haptic device was used as a nanomanipulator for a scanning tunneling microscope as described in "The Nanomanipulator: A Virtual-Reality Interface for a Scanning Tunneling Microscope" by R. M. Taylor, W. Robinett, V. L. Chi, F. P. Brooks, W. V. Wright, R. S. Williams, and E. J. Snyder, Proceedings of SIGGRAPH '93, pp. 127-134 (August 1993). This enabled scientists to manipulate individual atoms on a planar surface. This, however, was limited only to this specific use.

A medical planning and training system in "Dynamic Force Feedback in a Virtual Knee Palpation," by G. Burdea, N. Langrana, K. Lange, D. Gomez, and S. Deshpande, Journal Of Artificial Intelligence in Medicine, (1994). has also been developed which simulates knee palpation through the use of visual and haptic feedback. This system was specially designed for, and works only with knee palpation.

Currently there is a need to provide tactile feedback corresponding with a visual representation of a computer model to aid in model manipulation.

SUMMARY OF THE INVENTION

A computer model system starts with a volumetric dataset model. A visual representation is provided to an operator. The operator then employs a haptic device to "examine" the model. The haptic device(s) is monitored to determine its location. This location is equated to locations of the model. A scalar map of values, referred to as density values are employed in determining a normal force to an isosurface, and also viscous forces opposite the velocity of the haptic device. This creates the illusion of physically touching the model.

The operator may then select one of several "virtual tools" which can be used to modify the volumetric model. The model parameters in a small locality of the device are modified to create the illusion of melting, building-up, painting, air brushing, burning, stamping and other modifications to the model. The forces provided to the haptic device are also modified to be consistent with the visual representation and the volumetric model.

OBJECTS OF THE INVENTION

An object of the present invention is to provide a system which provides force feedback sensation to an operator consistent with a visual representation of the computer model.

Another object of the present invention is to provide a means for manipulating a three-dimensional volumetric model in which the user may visualize modifications and "feel" the model as it is being modified.

BRIEF DESCRIPTION OF THE DRAWINGS

While the novel features of the invention are set forth with particularity in the appended claims, the invention, both as to organization and content, will be better understood and appreciated, along with other objects and features thereof, from the following detailed description taken in conjunction with the drawings, in which:

FIG. 1 is a simplified flowchart illustrating the functioning of the present invention.

FIG. 2 is a diagram of the forces being simulated.

FIG. 3 illustrates the use of some model modification tools on a volumetric wall.

FIG. 4 is a simplified block diagram of a computer modeling system according to the present invention.

DETAILED DESCRIPTION OF THE INVENTION

A simplified flowchart of the functioning of the present invention is shown in FIG. 1. The haptic interaction loop begins in step 11 after an initial image of the scene has been computed and displayed. The next step in the interaction loop is to obtain the current position and orientation of a haptic pointer device operated by an operator 1 at step 13. The virtual counterpart of this physical pointer device will be referred to as a "tool", since it will often be used, for example, as a virtual scalpel, chisel, or paintbrush.

An operator may intend to modify the model in step 15. If necessary, the modification computation is performed in step 17, and the volume database is updated in step 19. In addition, the pixel status buffer is updated in step 21 to indicate which pixels of the image have been affected by this modification operation.

Once the optional modification step is complete, the current force is computed in step 23 and supplied to the haptic device.

During the rendering in step 25, some small number of pixels that require updating are rendered by conventional methods, such as a ray casting method.

Finally, in step 27, it is determined if it is time to refresh the physical display device. If so, the current image is copied to the screen and graphics hardware is used to render a geometric object in step 27 that indicates the size, position, and orientation of the current tool.

The data modification operation need not occur for every iteration of the haptic interaction loop. Instead, a timer may be employed indicating elapsed time since the previous modification. The main reason for limiting the rate at which data modification occurs is that in the preferred embodiment, a haptic refresh rate of 1 to 5 kHz is maintained. Therefore, there is only a small amount of computational time left over after the force calculation in each iteration. Increasing the rate of data modification would decrease the amount of time available to update the pixels of the image affected by the modification.

The rate at which the physical display device is refreshed is also limited by an elapsed time threshold. In the case, the limit is imposed since refresh rates much greater than 30 Hz are unnecessary.

Each voxel of the volumetric data set represents some properties of the data in the small three-dimensional region surrounding that voxel. Typically, volumetric data consists of scalar values representing material density, from which a gradient direction and magnitude can be computed. Shading properties such as color and ambient, diffuse, and specular lighting coefficients may also be stored in each voxel. A voxel may contain additional material properties such as stiffness and viscosity.

Example data which could be used to represent a voxel of the computer model are shown in Table 1. In this example, a single index value is used to represent several material properties through the use of a look-up table. This reduces the amount of computer memory required to store the model, but requires all properties derived from the look-up table to be modified simultaneously.

              TABLE 1______________________________________Property      Type          Size (bytes)______________________________________Density       Scalar        1Gradient Direction         Encoded Unit Vector                       2Gradient Magnitude         Scalar        1Color         R,G,B         3Material Properties         Look-up Table Index                       1______________________________________

The force equations are based on two principal requirements. First, the interaction forces must be calculated fast enough to be used within an interactive system. Second, the forces imparted to the user should be consistent with the rendering of the volumetric object.

In order to meet the speed requirement and since the haptic device used in the present implementation could only handle translation forces, the force calculation is simplified to a point contact. This has been shown to be a reasonable simplification for many tasks. The general equation used for feeling an object using a point contact model is:

F=A+R(V)+S(N)                                              (1)

and is illustrated in FIG. 2. The force F supplied to operator 1 located at position P and moving in direction V is equal to the vector sum of an ambient force A, a motion retarding force R(V), and a stiffness force S(N) normal to object D₂ which is inside object D₁. The ambient force A is the sum of all global forces acting on the tool, such as a gravitational or buoyant force. The motion retarding force is proportional to velocity and can be used to represent a viscous force. The last term captures the stiffness of the object and is always in the direction of the gradient. When simulating interaction on rigid surfaces, which are generally governed by Hooke's law, this term can then be set to a linear force function in the direction of the surface normal and proportional to the penetration distance of point P.

This general equation for force feedback is the basis for calculating forces which are consistent with different types of rendering methods.

The display of volume data may require segmentation of a volume into its parts. In a similar manner, a segmentation step which produces force feedback properties to the model may also be employed.

In Table 1, please notice that the material properties take into account both viscous and stiffness forces.

While it is possible to precompute distance to an isosurface for every voxel in the volume, it was decided to investigate techniques for approximating stiffness and retarding forces based only on the density field. There are two reasons for this. First, the system allows for the interactive modification of the volume. Creating a new distance map for the volume would be prohibitive. Second, for small penetration distances, the density field itself can give a reasonable approximation of the distance to an isosurface. Similar to volume rendering, the retarding and stiffness force functions used to feel an isosurface become dependent on transfer functions: ##EQU1##

Here the density d is used as an indicator of penetration distance between the isosurface density values d_i and d_j, where d_i <d_j. The function f_r (d) maps density values into retarding force magnitudes while f_s (d) maps density value into stiffness force magnitudes. For example, these transfer functions may be: ##EQU2##

The retarding force f_r (d) is set to a linear function proportional to the difference in density above d_i. Similar to haptic volume rendering, the coefficients E, F, and G specify a linear mapping from density values to force magnitudes.

The stiffness force f_s (d) varies from zero to G depending linearly on where the value d lies between d_i and d_j. This can be viewed as a penetrable shell model with viscous internal properties. A nice property of this model is that it allows the user to feel subsurface structure since the density and normal vector may change below the surface. Many other operator-selected transfer functions may be used.

To aid the user during volume modification, synthetic guide forces can be applied. For example, a virtual plane perpendicular to a surface can be used as a guide when attempting to cut a straight line.

The data modification component of this haptic visualization method is an optional component that can be used to modify any local property of the data stored in the volume database. A local property is one that affects only a small region of the volume, and therefore will cause only a small region of the image to require updating. All properties represented by a value stored in each voxel are local properties of the model. Values which are derived through the use of a look-up table represent dependent properties that cannot be modified independently. For example, if opacity is derived from a look-up table indexed by material density, it cannot be modified using this method since any change to the look-up table could potentially result in wide spread changes in the model. The index value itself could be modified, thereby modifying all properties that are derived from the look-up table. If opacity is instead stored in each voxel, then this is a local property that can be modified with this method.

Local data properties can be modified by setting them to a constant value, or by applying a filter to the current values. With the constant value method, the current color value in each voxel that the tool encounters could be replaced by some constant color value in order to "paint" the object.

With the filter method, an object could be "melted" by updating density values in a small region around the tool location according to d_i =(1-α)d_i-1, where d_i is the new density value, d_i-1 is the current density value, and α is obtained by sampling a spherical filter with a Gaussian distribution centered at the tool location.

In contrast to melting, an object can be "constructed" using, d_i =αD+(1-α)d_i-1 where D is the density of the material that we are adding to the volume. Note that melting is just a special case of constructing with D=0. In fact, constructing will appear like melting whenever the opacity of d_i is less than the opacity of the d_i-1 that it is replacing.

The two modification methods described above can lead to a wide variety of tools by varying the constant value or filter type, and the data properties that are modified. In the present system, a modification operation is defined by providing the modification method, the extent of modification, the properties affected, and any necessary filter and constant values. A tool is composed of one or more modification operations, with the condition that no property is affected by more than one modification operation in a tool.

FIG. 3 illustrates the use of some model modification tools on a volumetric wall. These tools are described briefly in Table 2, where the first column indicates a tool name, the second column lists the modification method used for each operation, the third column defines the properties that are modified for each operation, and the last column describes the modification process. Specific constants for the instance of the tools shown in FIG. 2 are given in parenthesis in the last column.

              TABLE 2______________________________________                Property  DescriptionTool Name  Method    Affected  (example)______________________________________Melt       Filter    Density   remove densityConstruct  Filter    Density   add density                          (63% dense)Burn       Filter    Density   remove density      Filter    Color     blend in blackSquirt     Filter    Density   add density                          (63% dense)      Filter    Color     add colorStamp      Filter    Density   add shape (cross-                          hair filter)      Filter    Color     add color (green)      Constant  Index     set material (shiny)Paint      Constant            set color (yellow)      Constant            set material (shiny)Airbrush   Filter    Color     blend color (purple)      Constant  Index     set material (dull)______________________________________

In FIG. 4 a computer modeling system according to the present invention is shown. Operator 1 interacts with a haptic device 33, a control panel 31 and a display device 35. The display device may be any apparatus which can display two-dimensional images and computer models to operator 1. The control panel may be any conventional input device used for providing information to a computer system, or may be the haptic device itself.

The haptic device may be an input/output type device such as a "pencil-like" segment which operator 1 holds, connected to a series of arm segments in which the segments are connected to each other to make a continuous mechanical arm. Each segment meets at a joint which has both a sensor for determining the angular position of the joint, and a motor for causing a feedback force to be applied between the arm segments, and ultimately to the operator. A tip of a "pencil-like" arm segment of a haptic device 33 may be used as the functional point at which forces are computed and modifications are made.

A computer model is created by a model creation device 3 and provided to model memory 43. Optionally, model memory may receive a computer model in its final form and have it prestored. A tracking device 37 (which may alternatively be provided by haptic device 33) monitors sensors of haptic device 33 to determine a location of a desired point of haptic device 33. Typically the point monitored would be the end of the "pencil-like" arm segments of haptic device 33.

A force computation device 39 is connected to tracking device 37 and model memory 43. Force computation device 39 receives the information from tracking device 37 to determine the location of the functional point of the haptic device. When this point has a location which intersects the model stored in model memory 43, force computation device provides a force value to haptic device 33 causing haptic device 33 to exert that amount of force in the specified direction against operator 1. This gives the illusion to operator 1 that the tip is hitting the model.

Model memory has stored densities which represent stiffness and viscosity coefficients for different locations of the model to simulate different materials. For example, a very high stiffness force and a low viscosity force would simulate ice and conversely a very low stiffness force and a high viscosity would simulate a putty- like object.

A renderer 45 reads the graphics information from model memory 43 and provides this information to a display device 35. Display device 35 then displays model surfaces, and optionally the haptic device in correct relation to those model surfaces which is correctly associated with the forces experienced by operator 1 from haptic device 33. This provides the illusion through both sight and touch of actually interacting with the virtual model.

Not only is it possible to interact with, feel, and visualize the model, but when operator 1 interacts with control panel 31, a mode may be selected in which the model may be modified and forces computed and images rendered according to the modified model. The modes determined by a control panel 31 may allow operator 1, using haptic device 33 to cut away, burn, melt, build up, pull, or push portions of the model in the vicinity of the tip of pencil-like arm of haptic device 33.

Based on the mode selected by operator 1 through control panel 31, a model modification device 41 modifies parameters of the model consistent with the mode and location of haptic device 33. Examples of possible modifications which may be made to the model are described later.

The present invention provides force feedback based upon the characteristics of the model itself, however, the present invention does not have to be limited to only non-moving objects. The present invention, with proper force calculations, may also provide force feedback for interaction with moving objects. For instance, haptic device 33 may be used to collide with, or "hit" an object causing the object to move to a different location. This may provide feedback such as a "feel" of hitting an object which may be possibly moving when it is hit. The material of the object may also be taken into account along with the motion and momentum of the object in order to have the object crush or deform as it is hit and/or moved in a direction. In the instance of a collision, an inertial force is also computed and added to the other existing forces to produce the resultant force.

The more processing and the more complex the calculations, the more processing power is required which may become restrictive. Theoretically, however, it is possible to determine, and provide this force and visual feedback of moving bodies.

While specific embodiments of the invention have been illustrated and described herein, it is realized that modifications and changes will occur to those skilled in the art. It is therefore to be understood that the appended claims are intended to cover all such modification and changes as fall within the true spirit and scope of the invention.