WO2000016123A9 - Systems and processing algorithms for ultrasound time-of-flight digitizer systems - Google Patents

Systems and processing algorithms for ultrasound time-of-flight digitizer systems

Info

Publication number
WO2000016123A9
WO2000016123A9 PCT/US1999/020479 US9920479W WO0016123A9 WO 2000016123 A9 WO2000016123 A9 WO 2000016123A9 US 9920479 W US9920479 W US 9920479W WO 0016123 A9 WO0016123 A9 WO 0016123A9
Authority
WO
WIPO (PCT)
Prior art keywords
time
flight
wave
measurements
pulses
Prior art date
Application number
PCT/US1999/020479
Other languages
French (fr)
Other versions
WO2000016123A1 (en
Inventor
Zahi Zloter
Gideon Shenholz
Ron Serber
Original Assignee
Electronics For Imaging Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from US09/150,251 external-priority patent/US6265676B1/en
Application filed by Electronics For Imaging Inc filed Critical Electronics For Imaging Inc
Priority to EP99946775A priority Critical patent/EP1112513A1/en
Priority to AU59104/99A priority patent/AU751093B2/en
Publication of WO2000016123A1 publication Critical patent/WO2000016123A1/en
Publication of WO2000016123A9 publication Critical patent/WO2000016123A9/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S11/00Systems for determining distance or velocity not using reflection or reradiation
    • G01S11/14Systems for determining distance or velocity not using reflection or reradiation using ultrasonic, sonic, or infrasonic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/06Systems determining the position data of a target
    • G01S15/08Systems for measuring distance only
    • G01S15/10Systems for measuring distance only using transmission of interrupted, pulse-modulated waves
    • G01S15/101Particularities of the measurement of distance
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/87Combinations of sonar systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52004Means for monitoring or calibrating

Definitions

  • the present invention relates to digitizer systems and, in particular, it concerns processing algorithms for ultrasound time-of-flight digitizer systems and systems employing such algorithms.
  • the present invention is a processing algorithm for ultrasound time-of-flight digitizer systems, and a system employing such an algorithm.
  • a method for tracking variations in distance D calculated from time-of-flight measurements of a sequence of pulses of a pressure wave oscillation from a transmitter to a receiver, the pressure wave oscillation having a given wavelength and wave period comprising: (a) identifying a state of synchronous operation by obtaining at least two time-of-flight measurements derived from successive pressure wave pulses which satisfy given synchronicity criteria; (b) monitoring successive time-of-flight measurements to identify a shifted time-of- flight measurement which varies by at least half of the wave period from a predicted time-of-flight value calculated from a number of preceding time-of-flight measurements; (c) identifying a shift factor corresponding to an integer multiple of the wave period by which the shifted time-of-flight measurement must be corrected to obtain a corrected time-of-flight measurement falling within
  • the time-of-flight measurements are made by a technique configured to identify a predetermined point within a cycle, such as by identifying a first zero crossing of a received signal after the signal has exceeded a given threshold value.
  • a shifted time-of-flight measurement for which the shift factor exceeds a predetermined maximum value, typically of less than 3, is discarded.
  • the corrected time-of- flight measurement is employed as a previous time-of-flight measurement for the step of monitoring performed on a subsequent time-of-flight measurement.
  • the state of synchronous operation is identified by obtaining at least three time-of-flight measurements derived from successive pressure wave pulses for which successive time-of-flight measurements vary by less than half of the wave period.
  • the state of synchronous operation is identified by obtaining at least three time-of-flight measurements derived from successive pressure wave pulses which vary substantially linearly.
  • the predicted time-of- flight value is calculated by geometrical extrapolation from at least two previous time- of-flight measurements, and preferably by extrapolation of a second order polynomial fitting the previous three time-of-flight measurements.
  • At least one supplementary shift test is performed, the step of correcting being performed selectively in response to the supplementary shift test.
  • the supplementary shift test may include determining an order in which a positive and a negative signal amplitude threshold are exceeded, or may include: (a) determining at least one peak signal amplitude occurring after a signal amplitude threshold is exceeded; and (b) calculating whether the peak signal amplitude differs from that of a corresponding peak signal amplitude from a previous pulse by more than a predefined ratio.
  • the transmitter is associated with a drawing implement which includes a contact switch for identifying operative contact between the drawing implement and a surface, the sequence of pulses being initiated in response to identification of the operative contact, the method further comprising continuing transmission of the sequence of pulses for a given delay period, typically at least about _ second, after the contact switch has ceased to indicate the operative contact so as to preserve the state of synchronous operation during intermittent contact.
  • a system for processing timing information and a received signal corresponding to a sequence of pulses of a pressure wave oscillation received by a receiver to track variations in a distance D calculated from time-of-flight measurements of the pulses, the pressure wave oscillation having a given wavelength and wave period comprising: (a) a signal processor for processing the received signal to derive an effective time-of-arrival for each pulse; (b) a timing module associated with the signal processor, the timing module being configured to derive a time-of-flight for each pulse from the timing information and the effective time-of-arrival; (c) a synchronous operation module associated with the timing module and configured to analyze the times-of-flight to identify a state of synchronous operation when at least two successive pressure wave pulses satisfy predefined synchronicity criteria; (d) a monitoring module associated with the timing module and configured to monitor successive time-of-flight measurements to identify a shifted time-of-flight measurement which varies by at
  • the predefined synchronicity criteria include that, for first, second and third time-of-flight measurements calculated from a first, a second and a third successive pulse, respectively, a difference between the first and the second time-of-flight measurements and a difference between the second and the third time-of-flight measurements are both less than half of the wave period, or that at least three time-of-flight measurements derived from successive pressure wave pulses vary substantially linearly.
  • the signal processor is configured to identify a predetermined point within a cycle as the effective time-of- arrival. This point typically corresponds to a first zero crossing of a the received signal after the signal has exceeded a given threshold value.
  • the signal processor performs at least one supplementary shift test, the correction module being configured to correct the distance D selectively in response to the supplementary shift test.
  • This test may include determining an order in which a positive and a negative signal amplitude threshold are exceeded, or may include: (a) determining at least one peak signal amplitude occurring after a signal amplitude threshold is exceeded; and (b) calculating whether the peak signal amplitude differs from that of a corresponding peak signal amplitude from a previous pulse by more than a predefined ratio.
  • the monitoring module - calculates the predicted time-of-flight value by geometrical extrapolation from at least two previous time-of-flight measurements, and typically, by extrapolation of a second order polynomial fitting the previous three time-of-flight measurements.
  • a transmitter device for use with a system for digitizing operative strokes of a handheld drawing implement, the drawing implement having a body and an operative tip
  • the transmitter device comprising: (a) a housing with an aperture, the housing being configured for receiving a portion of the body of the drawing implement with its operative tip extending through the aperture; (b) a normally-closed switch deployed so as to be opened by relative movement between the drawing implement and the housing resulting from pressure applied to the operative tip; and (c) a primary spring deployed to bias the drawing implement to a forward position in which the switch is closed.
  • a secondary spring weaker than the primary spring, deployed to act upon the drawing implement in a rearward direction so as to suspend the drawing implement within the housing.
  • a centering element associated with the primary spring and providing an abutment surface configured to align a rear end of the drawing implement centrally within the housing.
  • a system for identifying the position of a transmitter of a pulsed pressure-wave signal in at least two dimensions comprising a plurality of interconnected modular receiver units, each of the modular receiver units including: (a) a pressure- wave receiver configured to convert a received pressure-wave signal into an electrical signal; (b) a signal processor associated with the pressure-wave receiver and configured to process the electrical signal to generate a detection output indicative of reception of a pulse of the pulsed pressure-wave signal; (c) a timing module associated with the signal processor and responsive to the detection output and an externally supplied synchronization input to measure a time-of-flight of the pulse; and (d) an output module associated with the timing module and configured to output data related to the time-of-flight.
  • the plurality of interconnected modular receiver units are interconnected by a common data line, the output module of each of the modular receiver units being configured to provide a cascade trigger function for triggering sequential transmission of the output data by all of the modular receiver units along the common data line.
  • FIG. 1 is a plot of ultrasound signal amplitude against time for the beginning of a pulse in which a threshold is exceeded during a second cycle of the pulse
  • FIG. 2 is a plot of ultrasound signal amplitude against time for the beginning of a pulse in which a threshold is exceeded during a third cycle of the pulse
  • FIG. 3 is a plot of distance against time for successive pulses illustrating a first aspect of an algorithm according to the present invention
  • FIGS. 4A and 4B are plots of ultrasound signal amplitude against time for the beginning of two pulses illustrating the use of a positive and negative threshold technique according to the present invention
  • FIGS. 5A and 5B are plots of ultrasound signal amplitude against time for the beginning of two further pulses illustrating a refinement of the technique of Figures 4 A and 4B employing peak amplitude detection;
  • FIG 6 is a block diagram of a receiver module from a system, constructed and operative according to the teachings of the present invention, for tracking variations in a distance D calculated from time-of-flight measurements of pressure-wave pulses;
  • FIG. 7 is a block diagram of a system including a number of the modules of Figure 6; and FIG. 8 is a schematic cross-sectional view through a transmitter device, constructed and operative according to the teachings of the present invention, for use with a system for digitizing operative strokes of a hand-held drawing implement.
  • the present invention is a processing algorithm for an ultrasound time-of-flight digitizer system, and a corresponding system.
  • the invention also provides a simple technique for encoding information in an electromagnetic synchronization signal.
  • the present invention is applicable to any situation where distance from a transmitter to a receiver is measured by the time-of-flight of a pressure-wave pulse. In a minimal case, this could be a single transmitter/receiver pair. Typically, however, such systems employ at least two receivers for two- dimensional measurement, or at least three receivers for three-dimensional measurement. Alternatively, multiple transmitters may be used with a single receiver.
  • the techniques of the invention relate primarily to processing of signals received by one receiver originating from one transmitter which may be applied to any transmitter/receiver pair within the system.
  • Figures 1 and 2 show the beginnings of two ultrasound pulses as the amplitude grows over several cycles.
  • a fixed amplitude threshold may be exceeded at various different stages of the pulse amplitude envelope.
  • the threshold is exceeded during the second cycle of the pulse in Figure 1 , but only in the third cycle of the pulse in Figure 2. Transitions between cycles such as shown frequently give rise to significant discontinuities or "jumps" in the measured distance.
  • the present invention provides a method for tracking variations in distance D calculated from time-of-flight measurements of a sequence of pulses of a pressure wave oscillation from a transmitter to a receiver.
  • the method includes identifying a state of synchronous operation b y obtaining at least two time-of-flight measurements derived from successive pressure wave pulses satisfying some predefined synchronicity criteria.
  • One preferred set of synchronicity criteria is that at least three time-of-flight measurements derived from successive pressure wave pulses for which successive time-of-flight measurements vary by less than half of the wave period k.
  • the first time- of-flight measurement 20 taken at ⁇ is used to define a range of ⁇ k/2 about the measured value.
  • the subsequent time-of-flight measurement 22 at t 2 is then compared to this range. When two successive measurements each fall within this range from the previous readings, the system operation is defined to be synchronous.
  • the calculation is based on a geometrical extrapolation from at least two previous distance measurements.
  • the extrapolation is based on a second order polynomial fitting the previous three time-of-flight measurements.
  • the method then identifies a shift factor corresponding to an integer multiple of the wave period k by which the shifted time-of-flight measurement must be corrected to obtain a corrected time-of-flight measurement falling within half of the wave period from the preceding time-of-flight measurement. In this case, a shift factor of -1 is required.
  • Distance D calculated from the shifted time-of-flight measurement, is then corrected by the product of the shift factor and the pressure wave wavelength to map shifted time-of-flight measurement 26 to corrected measurement 28.
  • the corrected value 28 is then used in the prediction of the subsequent time-of-flight value.
  • the initial calculation will assume the same shift factor as used for the previous measurement. If, however, as shown here, the subsequent automatically-corrected measurement 32 calculated using the same shift factor lies more than half a wavelength from the predicted value, a reverse shift factor will be used to return the measurement to the predicted value 30.
  • this method provides a particularly simple solution to the problem of discontinuities due to transitions between cycles of detection. Assuming that the period between pulses is sufficiently small that distances moved between - pulses are always less than A 2, a jump in excess of A 2 is a strong indication that a transition between cycles of detection has occurred. In principle, cumulative correction of such jumps to bring each step within the required range from the previous measurement should provide a complete solution to the problem. In practice, certain situations such as extreme acceleration of the transmitter unit may lead to mistaken corrections. Various refinements are preferably provided to deal with such problematic cases, as will now be described.
  • the refinements to the present invention to address problems of this type may be subdivided into two categories. Firstly, steps are taken to standardize and stabilize the detection process, thereby reducing the proportion of anomalous measurements. Secondly, as an alternative, or preferably at the same time, the "correction" of measurements based on the primary algorithm described above may be made conditional upon performance of at least one supplementary shift test designed to verify the correction.
  • time-of-flight measurements according to the present invention are preferably made by a technique configured to identify a predetermined point within a cycle. The simplest technique of this sort is to identify a first zero crossing of the received signal after the signal has exceeded a given threshold value. With respect to supplementary shift tests, a number of possible approaches exist, each of which may be used alone or in combination with others.
  • FIG. 4A shows a signal in which a positive threshold 40 is exceeded in the first cycle at point 42 whereas Figure 4B shows exceedance in the second cycle at point 44.
  • Figure 4A shows a signal in which a positive threshold 40 is exceeded in the first cycle at point 42 whereas Figure 4B shows exceedance in the second cycle at point 44.
  • These two cases are preferably distinguished according to the present invention by adding a corresponding negative threshold 46 and determining an order in which the positive and the negative signal amplitude thresholds are exceeded.
  • positive threshold 40 is exceeded before negative threshold 46 is exceeded at point 43.
  • negative threshold 46 is exceeded first at point 45. This switch-over in order may be taken as indicative that a change n detection cycle has occurred.
  • Figures 5A and 5B show that this last approach alone is not always sufficient. Specifically, Figure 5A shows a negative detection at point 48 followed by a positive detection at point 50. Figure 5B, on the other hand, shows a positive detection at point 52 followed by a negative detection at point 54. In both cases, however, the signal has been detected during the first cycle of the signal.
  • the at least one supplementary shift test preferably includes determining at least one peak signal amplitude occurring after a signal amplitude threshold is exceeded and comparing this peak amplitude with a corresponding amplitude from the previous detection.
  • the peak signal amplitudes measured after detection A 1 and A 2 are relatively close in magnitude.
  • the second peak amplitude A 2 is significantly greater than A r
  • a ratio above a given predefined value between peak amplitude measurements from successive detections may be taken as an indication of change in detection cycle.
  • a twin criteria of positive-negative threshold exceedance reversal and significant variation in peak amplitude after detection, used together, may provide a particularly robust confirmation algorithm for suspected detection cycle transitions.
  • the system of the invention is typically made up from a number of modules, each of which may be implemented in a wide range of forms. These include software modules, which are software programs performing the functions of the system as described below. These software modules should be written in a suitable programming language, such as C or C++, which is compatible with the operating system of the computer on which these modules are being run. The exact selection of such a programming language could easily be made by one of ordinary skill in the art.
  • computer includes, but is not limited to, personal computers (PC) having an operating system such as DOS, WindowsTM, OS/2TM or Linux; MackintoshTM computers; computers having JAVATM-OS as the operating system; and graphical workstations such as the computers of Sun Microsystems TM and Silicon GraphicsTM, and other computers having some version of the UNIX operating system such as AIX or SOLARISTM of Sun MicrosystemsTM; or any other known and available operating system, as well as any other suitable combination of computational hardware and an operating system.
  • PC personal computers
  • an operating system such as DOS, WindowsTM, OS/2TM or Linux
  • MackintoshTM computers computers having JAVATM-OS as the operating system
  • graphical workstations such as the computers of Sun Microsystems TM and Silicon GraphicsTM
  • AIX or SOLARISTM of Sun MicrosystemsTM or any other known and available operating system, as well as any other suitable combination of computational hardware and an operating system.
  • WindowsTM includes but is not limited to Windows95TM, Windows 3.xTM in which "x” is an integer such as "1”, Windows NTTM, Windows98TM, Windows CETM and any new or upgraded versions of these operating systems by Microsoft Inc. (Seattle, Washington, USA), as well as any other operating system from Microsoft Inc.
  • modules could be implemented as hardware alone, or as a combination of hardware and software known as "firmware", which includes software programming instructions burnt onto a ROM (read only memory) chip.
  • firmware which includes software programming instructions burnt onto a ROM (read only memory) chip.
  • Various of the modules may also include volatile or non-volatile memory or other storage devices of any kind as required. The precise implementation could easily be performed b y one of ordinary skill in the art. In any case, the modules are described herein according to function rather than strictly as physically separate entities.
  • a digitizer system is implemented using a number of modular receiver units 100 which can be connected together in a series of any required length to cover a given area. Typically at least two units are used, and preferably at least four. It has also been found advantageous to use an even number of units deployed so that pairs of units are optimally deployed with overlapping "fields of view" to receive from each subsection of the area. More specifically, in the implementation of Figure 6, receiver unit 100 is configured to track variations in a distance D calculated from time-of-flight measurements of pressure-wave pulses according to the techniques of the present invention described above.
  • receiver unit 100 has a signal processor 102 for processing the signal from ultrasound receiver 103 to derive an effective time-of- arrival for each pulse, typically taken as a first zero crossing of a the received signal after the signal has exceeded a given threshold value.
  • a timing module 104 then derives a time-of-flight for each pulse from the effective time-of-arrival and pulse timing information, preferably derived from a timing input 105 typically originating from an infrared receiver.
  • the time-of-flight values are then processed by a synchronous operation module 106, a monitoring module 108, a shift factor module 110 and a correction module 112 to adjust for jumps due to changes in detection cycle.
  • synchronous operation module 106 analyses the times-of-flight to identify a state of synchronous operation, and monitoring module 108 monitors successive time-of-flight measurements to identify a shifted time-of-flight measurements.
  • Shift factor module 110 identifies a shift factor corresponding to an integer multiple of the wave period by which the shifted time-of-flight measurement must be corrected to obtain a corrected time-of-flight measurement falling within half of the wave period from the predicted time-of-flight measurement.
  • shift factor module 110 is configured to designate as erroneous any shifted time-of-flight measurement for which the shift factor exceeds two, or preferably one.
  • Correction module corrects the distance D calculated from the shifted time-of-flight measurement accordingly.
  • the corrected value of D is then preferably fed back to monitoring module 108 to be used as the previous time-of-flight measurement for monitoring a subsequent time-of-flight measurement.
  • these are preferably implemented as part of signal processor 102.
  • a separate supplementary shift test module 114 may be provided.
  • output module 116 preferably has input, output and reset trigger connections for controlling a transmission sequence which enables a large number of units to be connected on common data lines, as will be described.
  • Output module 116 also provides data output and optionally data input connections, which may be implemented with desired connection format such as an RS232 interface.
  • FIG. 7 shows schematically how a number of receiver units 100 are combined in a system 120 with common data lines 122 and 124.
  • the trigger output of each unit is connected to the trigger input of the next, with the trigger input of the first unit ("channel 1") being connected to ground.
  • the output from an IR receiver 126 supplied along IR sync, line 128 to the IR signal timing module 105 of each unit 100 generates a trigger reset.
  • This trigger reset is effective to prime the unit to wait for the next received pressure-wave pulse and to set its trigger output to "high”.
  • each receiver unit 100 is configured to wait until its trigger input goes "low” (i.e., ground) and then to transmit its distance measurement data.
  • each channel transmits its data and then triggers transmission of the next channel.
  • the external device can readily identify which data comes from which receiver unit by the order of transmission which is dictated by the connections of trigger outputs and inputs and always starts from the channel which is permanently connected to ground. Additional identification information may also be included in the data transmission. It will be noted that the structure described provides great structural flexibility and simplicity. Since each receiver unit 100 performs only its own fairly simple signal processing and calculations, the internal structure of the units is simple and low cost to implement. Furthermore, since all receiver units are identical, production costs are minimized.
  • the modular nature of the system allows extension of the system to cover large areas without any structural modification other than addition of receiver units.
  • the corresponding modifications required in triangulation algorithms to be performed b y the external device are trivial.
  • the external device may also be programmed to adapt automatically to the number of receiver units present.
  • Transmitter device 60 for use with a system for digitizing operative strokes of a hand-held drawing implement 62.
  • Transmitter device 60 includes a housing 64 configured for receiving at least a portion of the body of the drawing implement 62 with its operative tip extending through an aperture 66.
  • a particular feature of transmitter device 60 is that operative contact of the drawing implement tip against a surface is identified by use of a normally-closed switch 68.
  • the term "normally-closed switch” is used to refer to a switch structure in which movement is detected by the breaking of a circuit normally completed by the switch.
  • normally-closed switch 68 is deployed so as to be opened by relative movement between drawing implement 62 and housing 64 resulting from pressure applied to the operative tip.
  • a primary spring 70 is deployed to bias drawing implement 62 to a forward position in which switch 68 is closed.
  • the use of a normally-closed switch provides immediate detection of contact with a surface by detecting the onset of relative movement between the drawing implement and the housing.
  • a secondary spring 72 weaker than primary spring 70, is deployed to act upon drawing implement 62 in a rearward direction, i.e., tending to retract the operative point.
  • the opposition and alignment of springs 70 and 72 serves to suspend the drawing implement properly aligned within housing 64 and allows the housing to accommodate drawing implements with a range of lengths.
  • the stronger primary spring 70 ensures that switch 68 returns to its closed state whenever contact force is not applied to the operative tip of drawing implement 62.
  • normally-closed switch 68 is deployed in a removable cover portion 64a of housing 64.
  • switch 68 is located at other positions, for example, near aperture 66.
  • a further preferred feature of transmitter device 60 is a centering element 74, typically associated with primary spring 70, which provides an abutment surface 76 shaped to align a rear end of drawing implement 62 centrally within the cylindrical hollow of housing 64. This ensures alignment of drawing implements within housing 64 independent of their outer dimensions and avoids frictional contact between the drawing implement and the walls of housing 64 which could otherwise interfere with operation of switch 68.
  • abutment surface 76 is formed with a tapered projection 78 which, under the action of springs 70 and 72, is self centering in the axial recess common to the rear end of most marker type drawing implements.
  • the circuitry of the transmitter unit this is typically relatively simple, requiring only synchronous pulsed transmission of the pressure-wave and electromagnetic transmitters.
  • the electromagnetic signal can carry additional information such as pen color and the states of one or more switches.
  • a specific problem in implementation of the transmitter relates to maintaining synchronicity during brief operative strokes.
  • drawing implements such as pens or erasers typically use a contact switch for identifying operative contact between the drawing implement and a surface. The sequence of ultrasound pulses is then initiated in when operative contact is formed. Since effective functioning of the present invention requires establishing "synchronous operation", the repeated stopping and starting of the ultrasound transmitter normally caused by short contact strokes leads to repeated loss of synchronicity.
  • transmission of the sequence of pulses is preferably continued for a given delay period, typically at least about _ second, after the contact switch has ceased to indicate the operative contact.
  • a given delay period typically at least about _ second
  • the status of the contact switch is then conveyed to the base unit separately, either through a hard- wired connection when present, or preferably via an electromagnetic link from the transmitter unit to a corresponding electromagnetic receiver. In the latter case, transmission both of the pressure-wave signals and of the electromagnetic signal are continued for the delay period.

Abstract

A method and corresponding system for tracking variations in distance D calculated from time-of-flight measurements of a sequence of pulses of a pressure wave oscillation from a transmitter to a receiver identifies a state of synchronous operation by obtaining at least two time-of-flight measurements derived from successive pressure wave pulses which satisfy given synchronicity criteria. Successive time-of-flight measurements are then monitored to identify a shifted time-of-flight measurement which varies by at least half of the wave period from a predicted time-of-flight value calculated from a number of preceding time-of-flight measurements. A shift factor is then identified corresponding to an integer multiple of the wave period by which the shifted time-of-flight measurement must be corrected to obtain a corrected time-of-flight measurement falling within half of the wave period from the predicted time-of-flight value. The distance D calculated from the shifted time-of-flight measurement is then corrected by the product of the shift factor and the pressure wave wavelength.

Description

SYSTEMS AND PROCESSING ALGORITHMS FOR ULTRASOUND
TIME-OF-FLIGHT DIGITIZER SYSTEMS
This application is a continuation-in-part of serial no. 09/030,825 which was filed on 26 February, 1998.
FIELD AND BACKGROUND OF THE INVENTION
The present invention relates to digitizer systems and, in particular, it concerns processing algorithms for ultrasound time-of-flight digitizer systems and systems employing such algorithms.
It is known to employ ultrasound time-of-flight measurements to measure distances in a two- or three-dimensional digitizer system. Such systems employ one transmitter which transmits ultrasound pulses which are received by receivers at a number of positions. Alternative configurations employ multiple transmitters with a single receiver. In either case, by measuring the time-of-flight of the pulses from the transmitter to the receiver, the distance between them can be calculated. Synchronization of the transmitter and receiver may be achieved either by a hard wired connection or by a wireless electromagnetic link.
Identification of the beginning of each ultrasound pulse received is non trivial. The reaction characteristic of transducer elements generally used to generate the pulses, together with dispersion of the signal during transit, cause the pulse to rise over a number of cycles before reaching its maximum amplitude. As a result, a threshold level set to reliably distinguish the pulse signals over background noise may be triggered at different stages of the pulse. Depending on the transmitter- receiver distance and various environmental conditions, the threshold may be exceeded sometimes during the first cycle, sometimes during the second, and sometimes during the third, leading to considerable imprecision (see Figures 1 and 2). For a typical operating frequency of about 40 kHz and taking the speed of sound to be 330 m/s, each cycle corresponds to a distance of about 8 mm. Such a range of error is unacceptable for typical applications such as digitizers for writing implements, computer mice and the like.
Conventional ultrasound digitizer systems are also somewhat inflexible in their hardware configurations. Typically, a predefined receiver arrangement is produced for each given application. The electronic components are then centralized in a control box. Such arrangements leave little or no flexibility to adapt systems to applications with larger dimensions.
There is therefore a need for systems and methods for accurately tracking variations in distance calculated from time-of-flight measurements of a sequence of pulses of a pressure wave from a transmitter to a receiver. It would also be highly advantageous to provide a modular receiver system in which the number of receiver units may be increased to cover any desired area.
SUMMARY OF THE INVENTION
The present invention is a processing algorithm for ultrasound time-of-flight digitizer systems, and a system employing such an algorithm. According to the teachings of the present invention there is provided, a method for tracking variations in distance D calculated from time-of-flight measurements of a sequence of pulses of a pressure wave oscillation from a transmitter to a receiver, the pressure wave oscillation having a given wavelength and wave period, the method comprising: (a) identifying a state of synchronous operation by obtaining at least two time-of-flight measurements derived from successive pressure wave pulses which satisfy given synchronicity criteria; (b) monitoring successive time-of-flight measurements to identify a shifted time-of- flight measurement which varies by at least half of the wave period from a predicted time-of-flight value calculated from a number of preceding time-of-flight measurements; (c) identifying a shift factor corresponding to an integer multiple of the wave period by which the shifted time-of-flight measurement must be corrected to obtain a corrected time-of-flight measurement falling within half of the wave period from the predicted time-of-flight value; and (d) correcting the distance D calculated from the shifted time-of-flight measurement by the product of the shift factor and the pressure wave wavelength.
According to a further feature of the present invention, the time-of-flight measurements are made by a technique configured to identify a predetermined point within a cycle, such as by identifying a first zero crossing of a received signal after the signal has exceeded a given threshold value.
According to a further feature of the present invention, a shifted time-of-flight measurement for which the shift factor exceeds a predetermined maximum value, typically of less than 3, is discarded.
According to a further feature of the present invention, the corrected time-of- flight measurement is employed as a previous time-of-flight measurement for the step of monitoring performed on a subsequent time-of-flight measurement.
According to a further feature of the present invention, the state of synchronous operation is identified by obtaining at least three time-of-flight measurements derived from successive pressure wave pulses for which successive time-of-flight measurements vary by less than half of the wave period.
According to a further feature of the present invention, the state of synchronous operation is identified by obtaining at least three time-of-flight measurements derived from successive pressure wave pulses which vary substantially linearly.
According to a further feature of the present invention, the predicted time-of- flight value is calculated by geometrical extrapolation from at least two previous time- of-flight measurements, and preferably by extrapolation of a second order polynomial fitting the previous three time-of-flight measurements.
According to a further feature of the present invention, at least one supplementary shift test is performed, the step of correcting being performed selectively in response to the supplementary shift test. The supplementary shift test may include determining an order in which a positive and a negative signal amplitude threshold are exceeded, or may include: (a) determining at least one peak signal amplitude occurring after a signal amplitude threshold is exceeded; and (b) calculating whether the peak signal amplitude differs from that of a corresponding peak signal amplitude from a previous pulse by more than a predefined ratio.
According to a further feature of the present invention, the transmitter is associated with a drawing implement which includes a contact switch for identifying operative contact between the drawing implement and a surface, the sequence of pulses being initiated in response to identification of the operative contact, the method further comprising continuing transmission of the sequence of pulses for a given delay period, typically at least about _ second, after the contact switch has ceased to indicate the operative contact so as to preserve the state of synchronous operation during intermittent contact.
There is also provided according to the teachings of the present invention, a system for processing timing information and a received signal corresponding to a sequence of pulses of a pressure wave oscillation received by a receiver to track variations in a distance D calculated from time-of-flight measurements of the pulses, the pressure wave oscillation having a given wavelength and wave period, the system comprising: (a) a signal processor for processing the received signal to derive an effective time-of-arrival for each pulse; (b) a timing module associated with the signal processor, the timing module being configured to derive a time-of-flight for each pulse from the timing information and the effective time-of-arrival; (c) a synchronous operation module associated with the timing module and configured to analyze the times-of-flight to identify a state of synchronous operation when at least two successive pressure wave pulses satisfy predefined synchronicity criteria; (d) a monitoring module associated with the timing module and configured to monitor successive time-of-flight measurements to identify a shifted time-of-flight measurement which varies by at least half of the wave period from a predicted time- of-flight value calculated from a number of preceding time-of-flight measurements; (e) a shift factor module associated with the monitoring module and configured to identify a shift factor corresponding to an integer multiple of the wave period b y which the shifted time-of-flight measurement must be corrected to obtain a corrected time-of-flight measurement falling within half of the wave period from the predicted time-of-flight value; and (f) a correction module associated with the timing module and configured to correct the distance D calculated from the shifted time-of-flight measurement by the product of the shift factor and the pressure wave wavelength.
According to a further feature of the present invention, the predefined synchronicity criteria include that, for first, second and third time-of-flight measurements calculated from a first, a second and a third successive pulse, respectively, a difference between the first and the second time-of-flight measurements and a difference between the second and the third time-of-flight measurements are both less than half of the wave period, or that at least three time-of-flight measurements derived from successive pressure wave pulses vary substantially linearly. According to a further feature of the present invention, the signal processor is configured to identify a predetermined point within a cycle as the effective time-of- arrival. This point typically corresponds to a first zero crossing of a the received signal after the signal has exceeded a given threshold value.
According to a further feature of the present invention, the signal processor performs at least one supplementary shift test, the correction module being configured to correct the distance D selectively in response to the supplementary shift test. This test may include determining an order in which a positive and a negative signal amplitude threshold are exceeded, or may include: (a) determining at least one peak signal amplitude occurring after a signal amplitude threshold is exceeded; and (b) calculating whether the peak signal amplitude differs from that of a corresponding peak signal amplitude from a previous pulse by more than a predefined ratio. According to a further feature of the present invention, the monitoring module - calculates the predicted time-of-flight value by geometrical extrapolation from at least two previous time-of-flight measurements, and typically, by extrapolation of a second order polynomial fitting the previous three time-of-flight measurements.
There is also provided according to the teachings of the present invention, a transmitter device for use with a system for digitizing operative strokes of a handheld drawing implement, the drawing implement having a body and an operative tip, the transmitter device comprising: (a) a housing with an aperture, the housing being configured for receiving a portion of the body of the drawing implement with its operative tip extending through the aperture; (b) a normally-closed switch deployed so as to be opened by relative movement between the drawing implement and the housing resulting from pressure applied to the operative tip; and (c) a primary spring deployed to bias the drawing implement to a forward position in which the switch is closed.
According to a further feature of the present invention, there is also provided a secondary spring, weaker than the primary spring, deployed to act upon the drawing implement in a rearward direction so as to suspend the drawing implement within the housing.
According to a further feature of the present invention, there is also provided a centering element associated with the primary spring and providing an abutment surface configured to align a rear end of the drawing implement centrally within the housing.
There is also provided according to the teachings of the present invention, a system for identifying the position of a transmitter of a pulsed pressure-wave signal in at least two dimensions, the system comprising a plurality of interconnected modular receiver units, each of the modular receiver units including: (a) a pressure- wave receiver configured to convert a received pressure-wave signal into an electrical signal; (b) a signal processor associated with the pressure-wave receiver and configured to process the electrical signal to generate a detection output indicative of reception of a pulse of the pulsed pressure-wave signal; (c) a timing module associated with the signal processor and responsive to the detection output and an externally supplied synchronization input to measure a time-of-flight of the pulse; and (d) an output module associated with the timing module and configured to output data related to the time-of-flight. According to a further feature of the present invention, the plurality of interconnected modular receiver units are interconnected by a common data line, the output module of each of the modular receiver units being configured to provide a cascade trigger function for triggering sequential transmission of the output data by all of the modular receiver units along the common data line.
BRIEF DESCRIPTION OF THE DRAWINGS
The invention is herein described, by way of example only, with reference to the accompanying drawings, wherein:
FIG. 1 is a plot of ultrasound signal amplitude against time for the beginning of a pulse in which a threshold is exceeded during a second cycle of the pulse; FIG. 2 is a plot of ultrasound signal amplitude against time for the beginning of a pulse in which a threshold is exceeded during a third cycle of the pulse;
FIG. 3 is a plot of distance against time for successive pulses illustrating a first aspect of an algorithm according to the present invention;
FIGS. 4A and 4B are plots of ultrasound signal amplitude against time for the beginning of two pulses illustrating the use of a positive and negative threshold technique according to the present invention;
FIGS. 5A and 5B are plots of ultrasound signal amplitude against time for the beginning of two further pulses illustrating a refinement of the technique of Figures 4 A and 4B employing peak amplitude detection; FIG 6 is a block diagram of a receiver module from a system, constructed and operative according to the teachings of the present invention, for tracking variations in a distance D calculated from time-of-flight measurements of pressure-wave pulses;
FIG. 7 is a block diagram of a system including a number of the modules of Figure 6; and FIG. 8 is a schematic cross-sectional view through a transmitter device, constructed and operative according to the teachings of the present invention, for use with a system for digitizing operative strokes of a hand-held drawing implement.
DESCRIPTION OF THE PREFERRED EMBODIMENTS
The present invention is a processing algorithm for an ultrasound time-of-flight digitizer system, and a corresponding system. The invention also provides a simple technique for encoding information in an electromagnetic synchronization signal.
The principles and operation of the systems, algorithms and techniques according to the present invention may be better understood with reference to the drawings and the accompanying description. Before turning to the techniques and systems of the present invention, it should be noted that the present invention is applicable to any situation where distance from a transmitter to a receiver is measured by the time-of-flight of a pressure-wave pulse. In a minimal case, this could be a single transmitter/receiver pair. Typically, however, such systems employ at least two receivers for two- dimensional measurement, or at least three receivers for three-dimensional measurement. Alternatively, multiple transmitters may be used with a single receiver. In any case, the techniques of the invention relate primarily to processing of signals received by one receiver originating from one transmitter which may be applied to any transmitter/receiver pair within the system. Referring now to the drawings, Figures 1 and 2 show the beginnings of two ultrasound pulses as the amplitude grows over several cycles. As explained earlier, a fixed amplitude threshold may be exceeded at various different stages of the pulse amplitude envelope. In the example shown here, the threshold is exceeded during the second cycle of the pulse in Figure 1 , but only in the third cycle of the pulse in Figure 2. Transitions between cycles such as shown frequently give rise to significant discontinuities or "jumps" in the measured distance. To address this problem, the present invention provides a method for tracking variations in distance D calculated from time-of-flight measurements of a sequence of pulses of a pressure wave oscillation from a transmitter to a receiver. In general terms, the method includes identifying a state of synchronous operation b y obtaining at least two time-of-flight measurements derived from successive pressure wave pulses satisfying some predefined synchronicity criteria. One preferred set of synchronicity criteria is that at least three time-of-flight measurements derived from successive pressure wave pulses for which successive time-of-flight measurements vary by less than half of the wave period k. Thus, referring to Figure 3, the first time- of-flight measurement 20 taken at ^ is used to define a range of ±k/2 about the measured value. The subsequent time-of-flight measurement 22 at t2 is then compared to this range. When two successive measurements each fall within this range from the previous readings, the system operation is defined to be synchronous.
Parenthetically, it should be noted that these criteria for synchronous operation are neither absolutely necessary nor exclusive. Thus, for example, it may be desired to supplement these criteria by allowing alternative criteria also defined as "synchronous operation". In one example, linear variation of D over three successive measurements may be considered synchronous operation even though the change in D may exceed ±k/2 between measurements. After entering synchronous operation, successive time-of-flight measurements are monitored to identify the first time-of-flight measurement to lie outside the A 2 range from a predicted time-of-flight value calculated from a number of preceding time-of-flight measurements. In this case, the first so-called "shifted time-of-flight measurement" 26 occurs at t4. The predicted time-of-flight value may be calculated in a number of ways. Preferably, the calculation is based on a geometrical extrapolation from at least two previous distance measurements. In one particularly effective implementation, the extrapolation is based on a second order polynomial fitting the previous three time-of-flight measurements. The method then identifies a shift factor corresponding to an integer multiple of the wave period k by which the shifted time-of-flight measurement must be corrected to obtain a corrected time-of-flight measurement falling within half of the wave period from the preceding time-of-flight measurement. In this case, a shift factor of -1 is required. Distance D, calculated from the shifted time-of-flight measurement, is then corrected by the product of the shift factor and the pressure wave wavelength to map shifted time-of-flight measurement 26 to corrected measurement 28. An assumption is normally made that subsequent measurements should also be corrected by the same shift factor until a new mismatch to the predicted value is encountered. Clearly, the shift factor used in this method should be limited to some reasonable value. It is much more likely that a sudden unpredicted change in measured distance of multiple wavelengths originates from sound interference or a momentary malfunction than from multiple shifts in the cycle of detection. Typically, a variation of only one, or at most, two cycles between detections would be expected.
It should be noted that the corrected value 28 is then used in the prediction of the subsequent time-of-flight value. The initial calculation will assume the same shift factor as used for the previous measurement. If, however, as shown here, the subsequent automatically-corrected measurement 32 calculated using the same shift factor lies more than half a wavelength from the predicted value, a reverse shift factor will be used to return the measurement to the predicted value 30.
Conceptually, this method provides a particularly simple solution to the problem of discontinuities due to transitions between cycles of detection. Assuming that the period between pulses is sufficiently small that distances moved between - pulses are always less than A 2, a jump in excess of A 2 is a strong indication that a transition between cycles of detection has occurred. In principle, cumulative correction of such jumps to bring each step within the required range from the previous measurement should provide a complete solution to the problem. In practice, certain situations such as extreme acceleration of the transmitter unit may lead to mistaken corrections. Various refinements are preferably provided to deal with such problematic cases, as will now be described.
The refinements to the present invention to address problems of this type may be subdivided into two categories. Firstly, steps are taken to standardize and stabilize the detection process, thereby reducing the proportion of anomalous measurements. Secondly, as an alternative, or preferably at the same time, the "correction" of measurements based on the primary algorithm described above may be made conditional upon performance of at least one supplementary shift test designed to verify the correction. With respect to standardization and stabilization of the detection process, time-of-flight measurements according to the present invention are preferably made by a technique configured to identify a predetermined point within a cycle. The simplest technique of this sort is to identify a first zero crossing of the received signal after the signal has exceeded a given threshold value. With respect to supplementary shift tests, a number of possible approaches exist, each of which may be used alone or in combination with others.
A first supplemental shift test is illustrated in Figures 4A and 4B. Figure 4A shows a signal in which a positive threshold 40 is exceeded in the first cycle at point 42 whereas Figure 4B shows exceedance in the second cycle at point 44. These two cases are preferably distinguished according to the present invention by adding a corresponding negative threshold 46 and determining an order in which the positive and the negative signal amplitude thresholds are exceeded. Thus, in Figure 4A, positive threshold 40 is exceeded before negative threshold 46 is exceeded at point 43. In Figure 4B, on the other hand, negative threshold 46 is exceeded first at point 45. This switch-over in order may be taken as indicative that a change n detection cycle has occurred.
Figures 5A and 5B, on the other hand, show that this last approach alone is not always sufficient. Specifically, Figure 5A shows a negative detection at point 48 followed by a positive detection at point 50. Figure 5B, on the other hand, shows a positive detection at point 52 followed by a negative detection at point 54. In both cases, however, the signal has been detected during the first cycle of the signal.
To address this shortcoming, as an alternative to, or further refinement of, the order-based algorithm of Figures 4A and 4B, the at least one supplementary shift test preferably includes determining at least one peak signal amplitude occurring after a signal amplitude threshold is exceeded and comparing this peak amplitude with a corresponding amplitude from the previous detection. Thus, taking Figures 5A and 5B as detections of successive pulses, the peak signal amplitudes measured after detection A1 and A2, respectively, are relatively close in magnitude. In the successive pulses of Figures 4A and 4B, on the other hand, the second peak amplitude A2 is significantly greater than Ar Thus, a ratio above a given predefined value between peak amplitude measurements from successive detections may be taken as an indication of change in detection cycle. A twin criteria of positive-negative threshold exceedance reversal and significant variation in peak amplitude after detection, used together, may provide a particularly robust confirmation algorithm for suspected detection cycle transitions.
Turning now to Figure 6, it will be appreciated that the algorithms described above may be implemented in a wide range of ways. The system of the invention is typically made up from a number of modules, each of which may be implemented in a wide range of forms. These include software modules, which are software programs performing the functions of the system as described below. These software modules should be written in a suitable programming language, such as C or C++, which is compatible with the operating system of the computer on which these modules are being run. The exact selection of such a programming language could easily be made by one of ordinary skill in the art. Hereinafter, the term "computer" includes, but is not limited to, personal computers (PC) having an operating system such as DOS, Windows™, OS/2™ or Linux; Mackintosh™ computers; computers having JAVA™-OS as the operating system; and graphical workstations such as the computers of Sun Microsystems ™ and Silicon Graphics™, and other computers having some version of the UNIX operating system such as AIX or SOLARIS™ of Sun Microsystems™; or any other known and available operating system, as well as any other suitable combination of computational hardware and an operating system. Hereinafter, the term "Windows™" includes but is not limited to Windows95™, Windows 3.x™ in which "x" is an integer such as "1", Windows NT™, Windows98™, Windows CE™ and any new or upgraded versions of these operating systems by Microsoft Inc. (Seattle, Washington, USA), as well as any other operating system from Microsoft Inc.
In addition, the modules could be implemented as hardware alone, or as a combination of hardware and software known as "firmware", which includes software programming instructions burnt onto a ROM (read only memory) chip. Various of the modules may also include volatile or non-volatile memory or other storage devices of any kind as required. The precise implementation could easily be performed b y one of ordinary skill in the art. In any case, the modules are described herein according to function rather than strictly as physically separate entities.
In one particularly advantageous implementation which will be illustrated in Figures 6 and 7, a digitizer system is implemented using a number of modular receiver units 100 which can be connected together in a series of any required length to cover a given area. Typically at least two units are used, and preferably at least four. It has also been found advantageous to use an even number of units deployed so that pairs of units are optimally deployed with overlapping "fields of view" to receive from each subsection of the area. More specifically, in the implementation of Figure 6, receiver unit 100 is configured to track variations in a distance D calculated from time-of-flight measurements of pressure-wave pulses according to the techniques of the present invention described above. Thus, receiver unit 100 has a signal processor 102 for processing the signal from ultrasound receiver 103 to derive an effective time-of- arrival for each pulse, typically taken as a first zero crossing of a the received signal after the signal has exceeded a given threshold value. A timing module 104 then derives a time-of-flight for each pulse from the effective time-of-arrival and pulse timing information, preferably derived from a timing input 105 typically originating from an infrared receiver. The time-of-flight values are then processed by a synchronous operation module 106, a monitoring module 108, a shift factor module 110 and a correction module 112 to adjust for jumps due to changes in detection cycle. Specifically, synchronous operation module 106 analyses the times-of-flight to identify a state of synchronous operation, and monitoring module 108 monitors successive time-of-flight measurements to identify a shifted time-of-flight measurements. Shift factor module 110 then identifies a shift factor corresponding to an integer multiple of the wave period by which the shifted time-of-flight measurement must be corrected to obtain a corrected time-of-flight measurement falling within half of the wave period from the predicted time-of-flight measurement. Typically, shift factor module 110 is configured to designate as erroneous any shifted time-of-flight measurement for which the shift factor exceeds two, or preferably one. Correction module corrects the distance D calculated from the shifted time-of-flight measurement accordingly. The corrected value of D is then preferably fed back to monitoring module 108 to be used as the previous time-of-flight measurement for monitoring a subsequent time-of-flight measurement. With regard to the supplementary shift tests described earlier, these are preferably implemented as part of signal processor 102. Alternatively, or additionally, a separate supplementary shift test module 114 may be provided. To facilitate the modular connection of receiver units 100 as mentioned above, communication of distance measurements to an external device, such as a computer, is preferably managed by an output module 116. Output module 116 preferably has input, output and reset trigger connections for controlling a transmission sequence which enables a large number of units to be connected on common data lines, as will be described. Output module 116 also provides data output and optionally data input connections, which may be implemented with desired connection format such as an RS232 interface.
Figure 7 shows schematically how a number of receiver units 100 are combined in a system 120 with common data lines 122 and 124. The trigger output of each unit is connected to the trigger input of the next, with the trigger input of the first unit ("channel 1") being connected to ground. The output from an IR receiver 126 supplied along IR sync, line 128 to the IR signal timing module 105 of each unit 100 generates a trigger reset. This trigger reset is effective to prime the unit to wait for the next received pressure-wave pulse and to set its trigger output to "high". After receipt of a pressure-wave signal and the calculation of a corresponding distance measurement, each receiver unit 100 is configured to wait until its trigger input goes "low" (i.e., ground) and then to transmit its distance measurement data. Immediately after transmission, it sets its trigger output to low. As a result, a chain response is generated in which each channel transmits its data and then triggers transmission of the next channel. The external device can readily identify which data comes from which receiver unit by the order of transmission which is dictated by the connections of trigger outputs and inputs and always starts from the channel which is permanently connected to ground. Additional identification information may also be included in the data transmission. It will be noted that the structure described provides great structural flexibility and simplicity. Since each receiver unit 100 performs only its own fairly simple signal processing and calculations, the internal structure of the units is simple and low cost to implement. Furthermore, since all receiver units are identical, production costs are minimized.
The modular nature of the system allows extension of the system to cover large areas without any structural modification other than addition of receiver units. The corresponding modifications required in triangulation algorithms to be performed b y the external device are trivial. If desired, the external device may also be programmed to adapt automatically to the number of receiver units present.
Turning finally to Figure 8, this shows a preferred structure of transmitter device 60 for use with a system for digitizing operative strokes of a hand-held drawing implement 62. Transmitter device 60 includes a housing 64 configured for receiving at least a portion of the body of the drawing implement 62 with its operative tip extending through an aperture 66. A particular feature of transmitter device 60 is that operative contact of the drawing implement tip against a surface is identified by use of a normally-closed switch 68. The term "normally-closed switch" is used to refer to a switch structure in which movement is detected by the breaking of a circuit normally completed by the switch. In this case, normally-closed switch 68 is deployed so as to be opened by relative movement between drawing implement 62 and housing 64 resulting from pressure applied to the operative tip. A primary spring 70 is deployed to bias drawing implement 62 to a forward position in which switch 68 is closed. In contrast to conventional microswitches which switch after a predefined distance of travel, the use of a normally-closed switch provides immediate detection of contact with a surface by detecting the onset of relative movement between the drawing implement and the housing.
Preferably, a secondary spring 72, weaker than primary spring 70, is deployed to act upon drawing implement 62 in a rearward direction, i.e., tending to retract the operative point. The opposition and alignment of springs 70 and 72 serves to suspend the drawing implement properly aligned within housing 64 and allows the housing to accommodate drawing implements with a range of lengths. At the same time, the stronger primary spring 70 ensures that switch 68 returns to its closed state whenever contact force is not applied to the operative tip of drawing implement 62.
In the example shown here, normally-closed switch 68 is deployed in a removable cover portion 64a of housing 64. However, one ordinarily skilled in the art will appreciate that alternative implementations may also be used in which switch 68 is located at other positions, for example, near aperture 66.
A further preferred feature of transmitter device 60 is a centering element 74, typically associated with primary spring 70, which provides an abutment surface 76 shaped to align a rear end of drawing implement 62 centrally within the cylindrical hollow of housing 64. This ensures alignment of drawing implements within housing 64 independent of their outer dimensions and avoids frictional contact between the drawing implement and the walls of housing 64 which could otherwise interfere with operation of switch 68. Preferably, abutment surface 76 is formed with a tapered projection 78 which, under the action of springs 70 and 72, is self centering in the axial recess common to the rear end of most marker type drawing implements.
With regard to the circuitry of the transmitter unit, this is typically relatively simple, requiring only synchronous pulsed transmission of the pressure-wave and electromagnetic transmitters. Optionally, the electromagnetic signal can carry additional information such as pen color and the states of one or more switches. A specific problem in implementation of the transmitter relates to maintaining synchronicity during brief operative strokes. In order to avoid wasteful continuous transmission, drawing implements such as pens or erasers typically use a contact switch for identifying operative contact between the drawing implement and a surface. The sequence of ultrasound pulses is then initiated in when operative contact is formed. Since effective functioning of the present invention requires establishing "synchronous operation", the repeated stopping and starting of the ultrasound transmitter normally caused by short contact strokes leads to repeated loss of synchronicity. To address this problem, transmission of the sequence of pulses is preferably continued for a given delay period, typically at least about _ second, after the contact switch has ceased to indicate the operative contact. In this way the state of synchronous operation is preserved during intermittent contact. The status of the contact switch is then conveyed to the base unit separately, either through a hard- wired connection when present, or preferably via an electromagnetic link from the transmitter unit to a corresponding electromagnetic receiver. In the latter case, transmission both of the pressure-wave signals and of the electromagnetic signal are continued for the delay period.
It will be appreciated that the above descriptions are intended only to serve as examples, and that many other embodiments are possible within the spirit and the scope of the present invention.

Claims

WHAT IS CLAIMED IS:
1. A method for tracking variations in distance D calculated from time-of- flight measurements of a sequence of pulses of a pressure wave oscillation from a transmitter to a receiver, the pressure wave oscillation having a given wavelength and wave period, the method comprising: (One) identifying a state of synchronous operation by obtaining at least two time-of-flight measurements derived from successive pressure wave pulses which satisfy given synchronicity criteria; (Two) monitoring successive time-of-flight measurements to identify a shifted time-of-flight measurement which varies by at least half of the wave period from a predicted time-of-flight value calculated from a number of preceding time-of-flight measurements; (Three) identifying a shift factor corresponding to an integer multiple of the wave period by which said shifted time-of-flight measurement must be corrected to obtain a corrected time-of-flight measurement falling within half of the wave period from the predicted time-of-flight value; and (Four) correcting the distance D calculated from said shifted time-of-flight measurement by the product of said shift factor and the pressure wave wavelength.
2. The method of claim 1 , wherein the time-of-flight measurements are made by a technique configured to identify a predetermined point within a cycle.
3. The method of claim 2, wherein the time-of-flight measurements are made by identifying a first zero crossing of a received signal after the signal has exceeded a given threshold value.
4. The method of claim 1 , further comprising disregarding a shifted time- of-flight measurement for which said shift factor exceeds a predetermined maximum value.
5. The method of claim 4, wherein said predetermined maximum value is less than 3.
6. The method of claim 1 , wherein said corrected time-of-flight measurement is employed as a previous time-of-flight measurement for said step of monitoring performed on a subsequent time-of-flight measurement.
7. The method of claim 1 , wherein state of synchronous operation is identified by obtaining at least three time-of-flight measurements derived from successive pressure wave pulses for which successive time-of-flight measurements vary by less than half of the wave period.
8. The method of claim 1 , wherein state of synchronous operation is identified by obtaining at least three time-of-flight measurements derived from successive pressure wave pulses which vary substantially linearly.
9. The method of claim 1 , wherein said predicted time-of-flight value is calculated by geometrical extrapolation from at least two previous time-of-flight measurements.
10. The method of claim 1 , wherein said predicted time-of-flight value is calculated by extrapolation of a second order polynomial fitting the previous three time-of-flight measurements.
11. The method of claim 1 , further comprising performing at least one supplementary shift test, said step of correcting being performed selectively response to said supplementary shift test.
12. The method of claim 11 , wherein said at least one supplementary shift test includes determining an order in which a positive and a negative signal amplitude threshold are exceeded.
13. The method of claim 11 , wherein said at least one supplementary shift test includes:
(One) determining at least one peak signal amplitude occurring after a signal amplitude threshold is exceeded; and (Two) calculating whether said peak signal amplitude differs from that of a corresponding peak signal amplitude from a previous pulse by more than a predefined ratio.
14. The method of claim 1 , wherein the transmitter is associated with a drawing implement which includes a contact switch for identifying operative contact between the drawing implement and a surface, the sequence of pulses being initiated in response to identification of said operative contact, the method further comprising continuing transmission of the sequence of pulses for a given delay period after said contact switch has ceased to indicate said operative contact so as to preserve said state of synchronous operation during intermittent contact.
15. The method of claim 14, wherein said given delay period is at least about _ second.
16. A system for processing timing information and a received signal corresponding to a sequence of pulses of a pressure wave oscillation received by a receiver to track variations in a distance D calculated from time-of-flight measurements of the pulses, the pressure wave oscillation having a given wavelength and wave period, the system comprising:
(One) a signal processor for processing the received signal to derive an effective time-of-arrival for each pulse; (Two) a timing module associated with said signal processor, said timing module being configured to derive a time-of-flight for each pulse from the timing information and said effective time-of-arrival; (Three)a synchronous operation module associated with said timing module and configured to analyze said times-of-flight to identify a state of synchronous operation when at least two successive pressure wave pulses satisfy predefined synchronicity criteria; (Four) a monitoring module associated with said timing module and configured to monitor successive time-of-flight measurements to identify a shifted time-of-flight measurement which varies by at least half of the wave period from a predicted time-of-flight value calculated from a number of preceding time-of-flight measurements; (Five) a shift factor module associated with said monitoring module and configured to identify a shift factor corresponding to an integer multiple of the wave period by which said shifted time-of-flight measurement must be corrected to obtain a corrected time-of-flight measurement falling within half of the wave period from the predicted time-of-flight value; and (Six) a correction module associated with said timing module and configured to correct the distance D calculated from said shifted time-of-flight measurement by the product of said shift factor and the pressure wave wavelength.
17. The system of claim 16, wherein said predefined synchronicity criteria include that, for first, second and third time-of-flight measurements calculated from a first, a second and a third successive pulse, respectively, a difference between said first and said second time-of-flight measurements and a difference between said second and said third time-of-flight measurements are both less than half of the wave period.
18. The system of claim 16, wherein said predefined synchronicity criteria include that at least three time-of-flight measurements derived from successive pressure wave pulses vary substantially linearly.
19. The system of claim 16, wherein said signal processor is configured to identify a predetermined point within a cycle as said effective time-of-arrival.
20. The system of claim 19, wherein said predetermined point corresponds to a first zero crossing of a the received signal after the signal has exceeded a given threshold value.
21. The system of claim 16, wherein said shift factor module is configured to designate as erroneous any shifted time-of-flight measurement for which said shift factor exceeds a predetermined maximum value.
22. The system of claim 21 , wherein said predetermined maximum value is less than 3.
23. The system of claim 16, wherein said monitoring module is configured to employ said corrected time-of-flight measurement as the previous time-of-flight measurement for monitoring a subsequent time-of-flight measurement.
24. The system of claim 16, wherein said signal processor is configured to perform at least one supplementary shift test, said correction module being configured to correct the distance D selectively in response to said supplementary shift test.
25. The system of claim 24, wherein said at least one supplementary shift test includes determining an order in which a positive and a negative signal amplitude threshold are exceeded.
26. The system of claim 24, wherein said at least one supplementary shift test includes:
(One) determining at least one peak signal amplitude occurring after a signal amplitude threshold is exceeded; and
(Two) calculating whether said peak signal amplitude differs from that of a corresponding peak signal amplitude from a previous pulse by more than a predefined ratio.
27. The system of claim 16, wherein said monitoring module calculates said predicted time-of-flight value by geometrical extrapolation from at least two previous time-of-flight measurements.
28. The system of claim 16, wherein said monitoring module calculates said predicted time-of-flight value by extrapolation of a second order polynomial fitting the previous three time-of-flight measurements.
29. A transmitter device for use with a system for digitizing operative strokes of a hand-held drawing implement, the drawing implement having a body and an operative tip, the transmitter device comprising: (One) a housing with an aperture, the housing being configured for receiving a portion of the body of the drawing implement with its operative tip extending through said aperture;
(Two) a normally-closed switch deployed so as to be opened by relative movement between the drawing implement and said housing resulting from pressure applied to the operative tip; and
(Three)a primary spring deployed to bias the drawing implement to a forward position in which said switch is closed.
30. The transmitter device of claim 29, further comprising a secondary spring, weaker than said primary spring, deployed to act upon the drawing implement in a rearward direction so as to suspend the drawing implement within said housing.
31. The transmitter device of claim 29, further comprising a centering element associated with said primary spring and providing an abutment surface configured to align a rear end of the drawing implement centrally within said housing.
32. A system for identifying the position of a transmitter of a pulsed pressure-wave signal in at least two dimensions, the system comprising a plurality of interconnected modular receiver units, each of said modular receiver units including: (One) a pressure-wave receiver configured to convert a received pressure- wave signal into an electrical signal; (Two) a signal processor associated with said pressure-wave receiver and configured to process said electrical signal to generate a detection output indicative of reception of a pulse of the pulsed pressure-wave signal; (Three)a timing module associated with said signal processor and responsive to said detection output and an externally supplied synchronization input to measure a time-of-flight of said pulse; and (Four) an output module associated with said timing module and configured to output data related to said time-of-flight.
33. The system of claim 32, wherein said plurality of interconnected modular receiver units are interconnected by a common data line, said output module of each of said modular receiver units being configured to provide a cascade trigger function for triggering sequential transmission of said output data by all of said modular receiver units along said common data line.
PCT/US1999/020479 1998-09-10 1999-09-09 Systems and processing algorithms for ultrasound time-of-flight digitizer systems WO2000016123A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
EP99946775A EP1112513A1 (en) 1998-09-10 1999-09-09 Systems and processing algorithms for ultrasound time-of-flight digitizer systems
AU59104/99A AU751093B2 (en) 1998-09-10 1999-09-09 Systems and processing algorithms for ultrasound time-of-flight digitizer systems

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US09/150,251 US6265676B1 (en) 1997-03-10 1998-09-10 Systems and processing algorithms for ultrasound time-of-flight digitizer systems
US09/150,251 1998-09-10
US09/288,058 US6151014A (en) 1998-02-26 1999-04-07 Systems and processing algorithms for ultrasound time-of-flight digitizer systems

Publications (2)

Publication Number Publication Date
WO2000016123A1 WO2000016123A1 (en) 2000-03-23
WO2000016123A9 true WO2000016123A9 (en) 2000-08-03

Family

ID=26847466

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US1999/020479 WO2000016123A1 (en) 1998-09-10 1999-09-09 Systems and processing algorithms for ultrasound time-of-flight digitizer systems

Country Status (3)

Country Link
US (1) US6151014A (en)
AU (1) AU751093B2 (en)
WO (1) WO2000016123A1 (en)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6335723B1 (en) * 1998-10-02 2002-01-01 Tidenet, Inc. Transmitter pen location system
IL129450A (en) * 1999-04-14 2004-09-27 Pegasus Technologies Ltd Presentation board digitizers
US6925881B1 (en) 2002-01-17 2005-08-09 Southwest Research Institute Time shift data analysis for long-range guided wave inspection
CA2481643C (en) 2002-04-15 2012-09-11 Epos Technologies Limited Method and system for obtaining positioning data
EP1759268A2 (en) 2004-05-17 2007-03-07 Epos Technologies Limited Acoustic robust synchronization signaling for acoustic positioning system
US20110098554A1 (en) 2004-09-29 2011-04-28 Tel Hashomer Medical Research Infrastructure And Services Ltd. Monitoring of convection enhanced drug delivery
US7367944B2 (en) 2004-12-13 2008-05-06 Tel Hashomer Medical Research Infrastructure And Services Ltd. Method and system for monitoring ablation of tissues
ES2376588T3 (en) 2005-03-23 2012-03-15 Epos Development Ltd. PROCEDURE AND SYSTEM FOR A DIGITAL FEATHER ASSEMBLY.
CN103399675B (en) * 2007-03-14 2016-12-28 高通股份有限公司 MEMS microphone
US9181555B2 (en) * 2007-07-23 2015-11-10 Ramot At Tel-Aviv University Ltd. Photocatalytic hydrogen production and polypeptides capable of same
US8431841B2 (en) * 2009-01-29 2013-04-30 Pegasus Technologies Ltd. Rangefinder pivotable between obverse and reverse positions
US8064290B2 (en) * 2009-04-28 2011-11-22 Luidia, Inc. Digital transcription system utilizing small aperture acoustical sensors
US8120994B2 (en) * 2009-04-28 2012-02-21 Luidia, Inc. Digital transcription system utilizing acoustical detectors having apertures with a vertical orientation relative to the work surface
EP2671093B1 (en) * 2011-01-31 2019-01-16 Sunnybrook Health Sciences Centre Ultrasonic probe with ultrasonic transducers addressable on common electrical channel
US8907931B2 (en) 2011-06-20 2014-12-09 Sony Corporation Electronic terminal, input correction method, and program
CN102566797A (en) 2011-12-07 2012-07-11 深圳易方数码科技股份有限公司 Electronic pen and tracker
WO2016040218A1 (en) * 2014-09-08 2016-03-17 Luidia, Inc. Pen-location-determining and transcription method and apparatus with automatic page-flip detection
US11197962B2 (en) * 2018-02-26 2021-12-14 Verily Life Sciences Llc Waveform reconstruction for ultrasound time of flight measurements

Family Cites Families (54)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US33936A (en) * 1861-12-17 Improvement in pumps
US3134099A (en) * 1962-12-21 1964-05-19 Ibm Ultrasonic data converter
US3396366A (en) * 1965-05-18 1968-08-06 Lab For Electronics Inc Ultrasonic wave directive assembly
US3382480A (en) * 1966-07-06 1968-05-07 Thomas K. Tsao Method of detecting traffic
US3613066A (en) * 1968-10-22 1971-10-12 Cii Computer input equipment
US3838212A (en) * 1969-07-11 1974-09-24 Amperex Electronic Corp Graphical data device
US3706850A (en) * 1971-04-23 1972-12-19 Bell Telephone Labor Inc Telewriting system
US3914754A (en) * 1971-05-18 1975-10-21 Martonair Ltd Machine element position detection system
US4103278A (en) * 1976-12-06 1978-07-25 Kanje Satake Obstacle detecting apparatus using ultrasonic waves
NL7803764A (en) * 1978-04-10 1979-10-12 Philips Nv ACOUSTIC WRITING COMBINATION CONTAINING A WRITING INSTRUMENT WITH AN ASSOCIATED WRITING TABLET.
JPS5828554B2 (en) * 1978-11-22 1983-06-16 松下電器産業株式会社 ultrasonic distance meter
US4309098A (en) * 1978-12-01 1982-01-05 Polaroid Corporation Autofocus cine camera having automatic variable speed focusing
ES487980A1 (en) * 1979-03-27 1980-10-01 Contraves Ag Additional device for a particle analyzer.
US4336987A (en) * 1980-12-22 1982-06-29 Polaroid Corporation Camera employing zone focusing
US4439846A (en) * 1981-11-18 1984-03-27 Polaroid Corporation Sonar rangefinder system
US4506354A (en) * 1982-09-30 1985-03-19 Position Orientation Systems, Ltd. Ultrasonic position detecting system
US4488000A (en) * 1982-09-30 1984-12-11 New York Institute Of Technology Apparatus for determining position and writing pressure
US4440482A (en) * 1982-09-30 1984-04-03 Polaroid Corporation Sonic autofocus camera having variable sonic beamwidth
US4564928A (en) * 1982-09-30 1986-01-14 New York Institute Of Technology Graphical data apparatus
US4490814A (en) * 1982-09-30 1984-12-25 Polaroid Corporation Sonic autofocus camera having variable sonic beamwidth
US4480485A (en) * 1982-10-01 1984-11-06 Panametrics, Inc. Acoustic flowmeter with envelope midpoint tracking
JPS59107297U (en) * 1983-01-08 1984-07-19 富士通株式会社 Eraser for electronic blackboard
AU1371784A (en) * 1983-03-28 1984-10-25 Hitachi Construction Machinery Numerical value measuring instrument using oscilloscope device
US4578674A (en) * 1983-04-20 1986-03-25 International Business Machines Corporation Method and apparatus for wireless cursor position control
US4570033A (en) * 1983-09-09 1986-02-11 Numonics Corporation Polyphase digitizer
EP0136697A3 (en) * 1983-10-03 1987-03-25 Nissan Motor Co., Ltd. Ultrasonic vehicle rangefinder
US4552991A (en) * 1983-11-03 1985-11-12 Numonics Corporation Absolute position coordinate determining device employing a single phase difference measurement to determine cursor position
US4772764A (en) * 1983-12-01 1988-09-20 Xecutek Corporation Soundwave position digitizer
US4739860A (en) * 1984-05-29 1988-04-26 Nissan Motor Co., Ltd. Ultrasonic rangefinder
US4682159A (en) * 1984-06-20 1987-07-21 Personics Corporation Apparatus and method for controlling a cursor on a computer display
DE3582968D1 (en) * 1984-07-25 1991-07-04 Hitachi Ltd TABLET-COORDINATE INPUT DEVICE USING ELASTIC SHAFT.
JPS61167889A (en) * 1985-01-18 1986-07-29 Nippon Soken Inc Range measuring device
EP0229637A3 (en) * 1986-01-09 1988-10-05 Wacom Co., Ltd. Electronic blackboard apparatus
CH669676A5 (en) * 1986-01-23 1989-03-31 Zellweger Uster Ag
US4777329A (en) * 1987-08-24 1988-10-11 Microfield Graphics, Inc. Graphic input system
CN1018590B (en) * 1987-09-14 1992-10-07 株式会社瓦科姆 Electronic blackboard device and writer and eraser thereof
US4853715A (en) * 1988-06-17 1989-08-01 Numonics Corporation Plotter head control device
JPH01320521A (en) * 1988-06-22 1989-12-26 Wacom Co Ltd Electronic blackboard device and its writing tool or the like
DE3905956A1 (en) * 1989-02-25 1990-09-13 Fraunhofer Ges Forschung DEVICE FOR MEASURING ULTRASONIC TIME
US4963703A (en) * 1989-07-18 1990-10-16 Numonics Corporation Coordinate determining device using spatial filters
JPH06502509A (en) * 1990-10-22 1994-03-17 ロジテク,インコーポレイテッド Ultrasonic positioning method and device
GB2250093A (en) * 1990-11-20 1992-05-27 Graeme William Younger Ultrasonic position measuring system
CA2058219C (en) * 1991-10-21 2002-04-02 Smart Technologies Inc. Interactive display system
US5380959A (en) * 1992-06-15 1995-01-10 Carroll Touch, Inc. Controller for an acoustic touch panel
WO1994002911A1 (en) * 1992-07-24 1994-02-03 Toda Koji Ultrasonic touch system
US5248856A (en) * 1992-10-07 1993-09-28 Microfield Graphics, Inc. Code-based, electromagnetic-field-responsive graphic data-acquisition system
FR2698191A1 (en) * 1992-11-17 1994-05-20 Lectra Systemes Sa Acquisition and processing of graphic data
US5205826A (en) * 1992-11-30 1993-04-27 Chen Long Hsiung Automatically blocked safety syringe
US5434370A (en) * 1993-11-05 1995-07-18 Microfield Graphics, Inc. Marking system with pen-up/pen-down tracking
DE69423296T2 (en) * 1993-11-30 2000-11-30 Hewlett Packard Co Arrangement for data entry
JPH07306746A (en) * 1994-03-15 1995-11-21 Fujitsu Ltd Ultrasonic coordinate input device
JP3007933B2 (en) * 1994-12-15 2000-02-14 富士通株式会社 Ultrasonic coordinate input device
DE19602810A1 (en) * 1996-01-26 1997-07-31 Siemens Ag Method and device for measuring the transit time of an electrical, electromagnetic or acoustic signal
IL120417A (en) * 1997-03-10 2000-09-28 Electronics For Imaging Inc Presentation board digitizer systems

Also Published As

Publication number Publication date
AU751093B2 (en) 2002-08-08
WO2000016123A1 (en) 2000-03-23
AU5910499A (en) 2000-04-03
US6151014A (en) 2000-11-21

Similar Documents

Publication Publication Date Title
US6265676B1 (en) Systems and processing algorithms for ultrasound time-of-flight digitizer systems
US6151014A (en) Systems and processing algorithms for ultrasound time-of-flight digitizer systems
US6414673B1 (en) Transmitter pen location system
EP0491132B1 (en) Coordinate input apparatus
US20110261654A1 (en) Optimum pseudo random sequence determining method, position detection system, position detection method, transmission device and reception device
US9146307B2 (en) Ultrasonic detection device
JP2007500348A (en) Distance measuring method and apparatus using ultrasonic waves
CA2316208A1 (en) Improved pen positioning system
US8583381B2 (en) Ultrasonic propagation time measurement system
US8498839B2 (en) Position determination system, transmission device and reception device, and position determination method
US9157984B2 (en) Position detection system, transmission device and reception device, and position detection method
GB2318232A (en) Adaptive threshold circuit
Elmer et al. High resolution ultrasonic distance measurement in air using coded signals
US5046053A (en) Acoustic signal detection circuit
JP2003139848A (en) Two-frequency type microwave sensor
JP2006145403A (en) Ultrasonic measurement circuit and liquid-level detection system using the same
CN107576964B (en) Echo time measuring method of linear frequency conversion signal
EP1112513A1 (en) Systems and processing algorithms for ultrasound time-of-flight digitizer systems
KR102261686B1 (en) Ultrasonic Wave Distance Meter
JP2007179507A (en) System for preventing erroneous ultrasonic pen coordinate input
US8707786B2 (en) Object detection apparatus
CN113330686A (en) Transmission of values by means of pulse-width-modulated signals
JPS6144378A (en) Ultrasonic digitizer
JPS61116680A (en) Distance detector
Papageorgiou et al. Accurate Displacement Measurement based on the Frequency Variation Monitoring of Ultrasonic Signals

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AL AM AT AU AZ BA BB BG BR BY CA CH CN CU CZ DE DK EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT UA UG UZ VN YU ZA ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW SD SL SZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
AK Designated states

Kind code of ref document: C2

Designated state(s): AL AM AT AU AZ BA BB BG BR BY CA CH CN CU CZ DE DK EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT UA UG UZ VN YU ZA ZW

AL Designated countries for regional patents

Kind code of ref document: C2

Designated state(s): GH GM KE LS MW SD SL SZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG

COP Corrected version of pamphlet

Free format text: PAGES 2/5 AND 4/5, DRAWINGS, REPLACED BY NEW PAGES 2/5 AND 4/5; DUE TO LATE TRANSMITTAL BY THE RECEIVING OFFICE

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
WWE Wipo information: entry into national phase

Ref document number: 59104/99

Country of ref document: AU

WWE Wipo information: entry into national phase

Ref document number: 1999946775

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 1999946775

Country of ref document: EP

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

WWG Wipo information: grant in national office

Ref document number: 59104/99

Country of ref document: AU

WWW Wipo information: withdrawn in national office

Ref document number: 1999946775

Country of ref document: EP