WO2006107947A2 - Hyperspectral imaging in diabetes and peripheral vascular disease - Google Patents

Hyperspectral imaging in diabetes and peripheral vascular disease Download PDF

Info

Publication number
WO2006107947A2
WO2006107947A2 PCT/US2006/012461 US2006012461W WO2006107947A2 WO 2006107947 A2 WO2006107947 A2 WO 2006107947A2 US 2006012461 W US2006012461 W US 2006012461W WO 2006107947 A2 WO2006107947 A2 WO 2006107947A2
Authority
WO
WIPO (PCT)
Prior art keywords
image
tissue
color
saturation
interest
Prior art date
Application number
PCT/US2006/012461
Other languages
French (fr)
Other versions
WO2006107947A3 (en
Inventor
Jenny E. Freeman
Svetlana V. Panasyuk
Michael Hopmeier
Kevin Schomacker
Derek Brand
Original Assignee
Hypermed, 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
Application filed by Hypermed, Inc. filed Critical Hypermed, Inc.
Priority to CA2604829A priority Critical patent/CA2604829C/en
Publication of WO2006107947A2 publication Critical patent/WO2006107947A2/en
Publication of WO2006107947A3 publication Critical patent/WO2006107947A3/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/742Details of notification to user or communication with user or patient ; user input means using visual displays
    • A61B5/743Displaying an image simultaneously with additional graphical information, e.g. symbols, charts, function plots
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
    • A61B5/14552Details of sensors specially adapted therefor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/41Detecting, measuring or recording for evaluating the immune or lymphatic systems
    • A61B5/412Detecting or monitoring sepsis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/44Detecting, measuring or recording for evaluating the integumentary system, e.g. skin, hair or nails
    • A61B5/441Skin evaluation, e.g. for skin disorder diagnosis
    • A61B5/444Evaluating skin marks, e.g. mole, nevi, tumour, scar
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/44Detecting, measuring or recording for evaluating the integumentary system, e.g. skin, hair or nails
    • A61B5/441Skin evaluation, e.g. for skin disorder diagnosis
    • A61B5/445Evaluating skin irritation or skin trauma, e.g. rash, eczema, wound, bed sore
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7275Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/0205Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
    • G01J3/0224Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using polarising or depolarising elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/0272Handheld
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/0283Details using a charging unit
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/10Arrangements of light sources specially adapted for spectrometry or colorimetry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/12Generating the spectrum; Monochromators
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/2823Imaging spectrometer
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/40ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Definitions

  • the invention is directed to methods and systems of hyperspectral and multispectral imaging of medical tissues.
  • the invention is directed to new devices, tools and processes for the detection and evaluation of diseases and disorders such as diabetes and peripheral vascular disease that incorporate hyperspectral/multispectral imaging.
  • Type 1 diabetes Diabetes afflicts an estimated 194 million people worldwide, affecting 7.9% of Americans (over 21 million people) and 7.8% of Europeans. Between 85% and 95% of all diabetics suffer from Type 2 diabetes, although nearly 5 million people worldwide suffer from Type 1 diabetes, affecting an estimated 1.27 million people in Europe and another 1.04 million people in the United States 1 . Both Type 1 and Type 2 diabetic patients are at higher risk for a wide array of complications including heart disease, kidney disease (e.g. nephropathy), ocular diseases (e.g. glaucoma), and neuropathy and nerve damages to name a few 2 . The feet of diabetic patients are at risk for a wide array of complications, which are discussed below. Problems with the foot that affect the ambulatory nature of the patient are not only important from the standpoint of physical risk, but also convey an emotional risk as well, as these problems disrupt the fundamental independence of the patient by limiting his or her ability to walk.
  • kidney disease e.g. nephropathy
  • ocular diseases e.
  • Peripheral arterial disease affects primarily people older than 55. There are currently 59.3 million Americans older than 55, and over 12 million of them have symptomatic peripheral vascular disease. It is estimated that only 20% of all patients with PAD have been diagnosed at this time. This represents a dramatically underpenetrated market. Although pharmacologic treatments for PAD have traditionally been poor, 2.1 million nevertheless receive pharmacologic treatment for the symptoms of PAD, and current diagnostic tests are not considered to be very sensitive indicators of disease progression or response to therapy. Additionally, 443,000 patients undergo vascular procedures such as peripheral arterial bypass surgery (100,000) or peripheral angioplasty (343,000) annually and are candidates for pre and post surgical testing. One difficulty in diagnosing PAD is that in the general population, only about 10% of persons with PAD experience classic symptoms of intermittent claudication. About 40% of patients do not complain of leg pain, while the remaining 50% have leg symptoms which differ from classic claudication.
  • Venous occlusive disease due to incompetent valves in veins designed to prevent backfiow and deep vein thrombosis, results in venous congestion and eventually stasis ulcers. Approximately 70% of leg ulcers are due to venous occlusion. Many of these ulcers are found at the medial malleolus. The foot is generally swollen and the skin near the ulcer site is brownish in appearance.
  • Diabetic feet are at risk for a wide range of pathologies, including microcirculatory changes, peripheral vascular disease, ulceration, infection, deep tissue destruction and metabolic complications.
  • the development of an ulcer in the diabetic foot is commonly a result of a break in the barrier between the dermis of the skin and the subcutaneous fat that cushions the foot during ambulation. This, in turn, can lead to increased pressure on the dermis, resulting in tissue ischemia and eventual death, and ultimately result in an ulcer. 4
  • Musculoskeletal Abnormalities Include altered foot mechanics, limited joint mobility, and bony deformities, and can lead to harmful changes in biomechanics and gait. This increases pressures associated with various regions of the foot. Alteration or atrophy of fat pads from increased pressure can lead to skin loss or callus, both of which increase the risk of ulceration by two orders of magnitude.
  • Peripheral Vascular Disease Caused by atherosclerotic obstruction of large vessels resulting in arterial insufficiency 10 is common in the elderly populations and is yet more common and severe in diabetics. 11 Diabetics may develop atherosclerotic disease of large-sized and medium-sized arteries, however, significant atherosclerotic disease of the infrapopliteal segments is particularly common. The reason for this is thought to result from a number of metabolic abnormalities in diabetics, including high LDL and VLDL levels, elevated plasma von Willebrand factor, inhibition of prostacyclin synthesis, elevated plasma fibrinogen levels, and increased platelet adhesiveness.
  • Venous occlusions are typically observed in the elderly who typically presented with swollen lower extremities and foot ulcers typically at the medial malleolus.
  • Previous studies have shown that a foot ulcer precedes roughly 85% of all lower extremity amputations in diabetic patients ' and that 15% of all diabetic patients will develop a foot ulcer during the course of their lifetimes. 14 More than 88,000 amputations performed annually on diabetics, 15 and roughly an additional 30,000 amputations are performed on non- diabetics, mostly related to peripheral vascular disease.
  • Foot pathology is major source of morbidity among diabetics and is a leading cause of hospitalization.
  • the infected and/or ischemic diabetic foot ulcer accounts for about 25% of all hospital days among people with diabetes, and the costs of foot disorder diagnosis and management are estimated at several billion dollars annually. 16 ' 17
  • the first step in the assessment of the diabetic foot is the clinical examination 18 ' 19 . All patients with diabetes require a thorough pedal examination at least once a year, even without signs of neuropathy. Evaluation of the diabetic patient with peripheral vascular disease should include a thorough medical history, vascular history, physical examination, neurologic evaluation for neuropathy and a thorough vascular examination. 20
  • the present invention overcomes the problems and disadvantages associated with current strategies and designs, and provides new tools and methods for detecting tissue at risk of developing into an ulcer, for detecting problems with diabetic foot disease, and for evaluating the potential for wounds to heal.
  • One embodiment of the invention is directed to a medical instrument comprising a first stage optic responsive to illumination of tissue, a spectral separator, one or more polarizers, an imaging sensor, a diagnostic processor, a filter control interface, a general purpose operating module to assess the state of tissue in diabetic subjects following a set of instructions, and a calibrator.
  • the instrument further comprises a second stage optic responsive to illumination of tissue.
  • the set of instructions comprises preprocessing hyperspectral information, building a visual image, defining a region of interest in tissue, converting the visual image into units of optical density by taking a negative logarithm of each decimal base, decomposing a spectra for each pixel into several independent components, determining three planes for an RGB pseudo-color image, determining a sharpness factor plane, converting the RGB pseudo-color image to a hue-saturation-value/intensity image having a plane, adjusting the hue- saturation-value/intensity image plane with the sharpness factor plane, converting the hue- saturation-value/intensity image back to the RGB pseudo-color image, removing outliers beyond a standard deviation and stretching image between 0 and 1, displaying the region of interest in pseudo-colors; and characterizing a metabolic state of the tissue of interest.
  • the region of interest is one of a pixel, a specified region or an entire field of view.
  • determining three planes for an RGB pseudo-color image comprises one or more characteristic features of the spectra
  • determining a sharpness factor plane comprises a combination of the images at different wavelengths
  • removing outliers beyond a standard deviation comprises three standard deviations
  • displaying the region of interest in pseudo-colors comprises one of performing one in combination with a color photoimage of a subject, in addition to a color photo image of a subject, and projecting onto the tissue of interest.
  • defining the color intensity plane as apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation, defining the color intensity plane as reflectance in blue-green-orange region, adjusting the hue saturation comprises adjusting a color resolution of the pseudo-color image according to quality of apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation, adjusting the hue saturation further comprises one or a combination of reducing resolution of hue and saturation color planes by binning the image, resizing the image, and smoothing the image through filtering higher frequency components out, and further interpolating the smoothed color planes on a grid of higher resolution intensity plane.
  • Another embodiment of the invention is directed to a method for assessing the state of tissue of a diabetic subject comprising, preprocessing hyperspectral information, building a visual image, defining a region of interest in tissue, converting the visual image into units of optical density by taking a negative logarithm of each decimal base, decomposing a spectra for each pixel into several independent components, determining three planes for an RGB pseudocolor image, determining a sharpness factor plane, converting the RGB pseudo-color image to a hue-saturation- value/intensity image having a plane, adjusting the hue- saturation- value/intensity image plane with the sharpness factor plane, converting the hue-saturation-value/intensity image back to the RGB pseudo-color image, removing outliers beyond a standard deviation and stretching image between 0 and 1, displaying the region of interest in pseudo-colors; and characterizing a metabolic state of the tissue of interest.
  • the region of interest is one of a pixel, a specified region or an entire field of view.
  • determining three planes for an RGB pseudo-color image comprises one or more characteristic features of the spectra
  • determining a sharpness factor plane comprises a combination of the images at different wavelengths
  • removing outliers beyond a standard deviation comprises three standard deviations
  • displaying the region of interest in pseudo-colors comprises one of performing one in combination with a color photoimage of a subject, in addition to a color photo image of a subject, and projecting onto the tissue of interest.
  • defining the color intensity plane as apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation, defining the color intensity plane as reflectance in blue-green-orange region, adjusting the hue saturation comprises adjusting a color resolution of the pseudo-color image according to quality of apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation, adjusting the hue saturation further comprises one or a combination of reducing resolution of hue and saturation color planes by binning the image, resizing the image, and smoothing the image through filtering higher frequency components out, and further interpolating the smoothed color planes on a grid of higher resolution intensity plane.
  • FIG 1 Block diagram depicting a portable hyperspectral imaging apparatus.
  • FIG 2 Basic specifications of the MHSI system.
  • FIG 3 OxyHb and DeoxyHb HSV/I color chart. Schematic representation of the MHSI display (left) showing the interplay between the oxyHb and deoxyHb coefficients and describing some of the potential physiological consequences of values of the MHSI.
  • FIG 4 Representative data from dorsal surface of foot showing individual oxyHb and deoxyHb values and how they can be used to evaluate regions of the tissue.
  • FIG 5 Representative data from tissue showing sensitivity of MHSI to drug-induced changes in the vasculature, (left to right) Visible image of foot surface post iontophoresis (IP), representative spectra pre and post iontophoresis with Acetylcholine (IP) showing greater oxyHb levels after IP. Images of increased oxyHb coefficient ring where IP occurred, image of deoxyHb, showing little change post IP.
  • IP foot surface post iontophoresis
  • IP Acetylcholine
  • FIG 6 Representative data from an ulcer located on the sole (ulcer 1) and dorsal surface (ulcer 2) of the foot.
  • FIG 7 MHSI information from the soles and dorsal surfaces of four patients. Each row of images represents data from one patient. The two columns on the left represent data from the soles of the feet, while the columns on the right represent data from the dorsal surfaces of the feet.
  • FIG 8 MHSI image of diabetic foot ulcer with 200 segment radial profile.
  • FIG 9 MHSI of wounds during healing.
  • FIG 10 color image (a) and pseudo-color image (b) of a rabbit ear Description of the Invention
  • HSI or hyperspectral imaging is a novel method of "imaging spectroscopy” that generates a "gradient map" of a region of interest based on local chemical composition.
  • HSI has been used in satellite investigation of suspected chemical weapons production areas 22 , geological features 23 , and the condition of agricultural fields 24 and has recently been applied to the investigation of physiologic and pathologic changes in living tissue in animal and human studies to provide information as to the health or disease of tissue that is otherwise unavailable.
  • 25 MHSI for medical applications has been shown to accurately predict viability and survival of tissue deprived of adequate perfusion, and to differentiate diseased (e.g. tumor) and ischemic tissue from normal tissue.
  • Spectroscopy is used in medicine to monitor metabolic status in a variety of tissues.
  • One of the most common spectroscopic applications is in pulse oximetry, which utilize the different oxyhemoglobin (oxyHb) and deoxyhemoglobin (deoxyHb) absorption bands to estimate arterial hemoglobin oxygen saturation.
  • oxyHb oxyhemoglobin
  • deoxyHb deoxyhemoglobin
  • HSI spatial information for the individual components and the ratio
  • imaging spectroscopy where the multi- dimensional (spatial & spectral) data are represented in what is called a “hypercube.”
  • the spectrum of reflected light is acquired for each pixel in a region, and each such spectrum is subjected to standard analysis. This allows the creation of an image based on the metabolic state of the region of interest (ROI).
  • ROI region of interest
  • Non-invasive measurements of oxygen or blood flow have been demonstrated previously, with investigators using thermometry, 42 point diffuse reflectance spectroscopy, 43 ' 44 and laser Doppler imaging. 45 Sheffield et al, have also reviewed laser Doppler and TcPO 2 measurements and their specific applications to wound healing. 46 While other techniques have been utilized in both the research lab and the clinic and have the advantage of a longer experience base, MHSI is superior to other technologies and can provide predictive information on the onset and outcomes of diabetic foot ulcers, venous stasis ulcers and peripheral vascular disease.
  • MHSI has the ability to show anatomically relevant information that is useful in the assessment of local, regional and systemic disease. This is important in the assessment of people with diabetes and/or peripheral vascular disease. MHSI shows the oxygen delivery and oxygen extraction of each pixel in the image collected. These images with pixels ranging from 20 microns to 120 microns have been useful in several ways.
  • MHSI shows the effects on the microcirculation of systemic diabetes, smoking, a variety of medications such as all of the classes of antihypertensives (ACE inhibitors, ARBs, Beta blockers, Peripheral arterial and arteriolar dilators), vasodilators (such as nitroglycerine, quinine, morphine), vasoconstrictors (including coffee, tobacco, pseudephedrine, Ritalin, epinephrine, levophedrine, neosynepherine), state of hydration, state of cardiac function (baseline, exercise, congestive heart failure), systemic infection or sepsis as well as other viral or bacterial infections and parasitic diseases.
  • ACE inhibitors ACE inhibitors, ARBs, Beta blockers
  • vasodilators such as nitroglycerine, quinine, morphine
  • vasoconstrictors including coffee, tobacco, pseudephedrine, Rital
  • the size of the pixels used is important in that it is smaller than the spacing of the perforating arterioles ( ⁇ 0.8 mm) 47 of the dermis and therefore permits the visualization of the distribution of mottling or other patterns associated with the anatomy of the microcirculation and its responses.
  • the image delivers information about the oxygen delivery and oxygen extraction for a particular region as it is influenced by blood flow through the larger vessels of that region of the body. For example an image of the top of the foot reflects both the systemic microvascular status and the status of the large (macrovascular) vessels supplying the leg.
  • MHSI shows the actual effect of the combination of systemic, regional and local effects on small pieces of tissue. This combines the effects of systemic and regional effects described above with the effects of local influences on the tissue including pressure, neuropathy, localized small vessel occlusion, localized trauma or wounding, pressure sore, inflammation, and wound healing. Angiogenesis during wound healing is readily monitored with MHSI.
  • Wounds other than on the foot can be similarly assessed, such as sacral decubiti, other areas of pressure necrosis, prosthesis stumps, skin flap tissue before, after or during surgery, areas of tissue breakdown after surgery, and burn injuries.
  • Current optical methods for evaluating tissues for the conditions described above include:
  • TcPO 2 Transcutaneous PO2
  • TcPO 2 measurements appeared cumbersome, lengthy (-20-30 minutes), highly operator dependant, and carried data only from skin directly under the probe (with little ability to distinguish the spatial characteristics of the ischemic area). While TcPO 2 has been shown to carry statistically significant information in terms of quantifying tissue at risk for ulceration, 48 TcPO 2 was not encouraging as a useful clinical device.
  • Non-imaging techniques such as near-infrared absorption spectroscopy (NIRS) or TcPO 2 , rely on measurements at a single point in tissue which may not accurately reflect overall tissue condition or provide anatomically relevant data, and probe placement on the skin can alter blood flow and cannot deliver accurate information in the area of an ulcer or directly surrounding it. Because MHSI is truly remote sensing, data are acquired at a distance, eliminating probe placement errors and allowing the investigation of the wound itself, which some techniques can not accomplish due to infection risk.
  • NIRS near-infrared absorption spectroscopy
  • TcPO 2 TcPO 2
  • Oxyhemoglobin in the arms and feet of ulcerated subjects is lower than oxyhemoglobin in diabetics without the ulceration.
  • the strong signal suggests ability to distinguish diabetics at lower and high risk.
  • Oxygen saturation level in the skin of arms and feet of diabetics is lower than oxygen saturation in the skin of controls. This is at a statistically significant level allowing separation between diabetics and controls.
  • MHSI quantitatively assesses different areas of tissue metabolism on both dorsal and plantar foot surfaces of any curvature.
  • MHSI evaluates state of tissue as a function of distance away from ulcer to assess the viability of surrounding tissue, and evaluate the degree of risk of further ulceration.
  • MHSI can be classified with a 4-quadrant system to determine the metabolic state of tissue using oxygen delivery and oxygen extraction: low/low, low/high, high/high, and high/low. This metric is used in distinguishing healthy tissue from ulcerated, or from a tissue at risk of ulceration.
  • MHSI is a unique visualization method that produces an image that combines spatial information from three independent parameters characterizing tissue: oxygenated and deoxygenated hemoglobin concentrations and light absorption.
  • MHSI evaluates skin metabolism at high resolution of 20-120 microns per image pixel.
  • Areas of decreased MHSI indicate tissue at risk for non-healing, ulcer extension, or primary ulceration.
  • MHSI differentiates between regions of tissue associated with a present foot ulcer on the basis of biomarkers such as oxyHb and deoxyHb coefficients.
  • MHSI evaluates temporal changes in oxygen delivery and extraction to particular areas, both, on local and systemic scale. The trend in the change of oxyHb and deoxyHb are used to predict healing status of a wound/ulcer as well as progression of diabetic complications.
  • MHSI examines tissue for gross features that may be indicative of global risks of complications, such as poor perfusion or the inability of the microcirculation to react and compensate in tissue.
  • MHSI has potential in diagnosing global microcirculatory insufficiencies and impacting on other complications of diabetes associated with the microvasculature besides foot ulcers.
  • MHSI is superior to other modalities for assessing the healing potential of tissue adjacent to ulcers.
  • MHSI provides more direct measurements of oxyHb and deoxyHb activities of the affected tissue. Hence, the discrimination is not markedly improved by adding iontophoresis results to refine prediction as is required for Laser Doppler to do so.
  • MHSI has significant advantages over laser Doppler and TcPO 2 measurements. Whereas MHSI is able to deliver spatially relevant data with high spatial resolution, TcPO 2 delivers only single point data. Laser Doppler data has poor spatial resolution and is frequently reported as a single mean numerical value across the region of interest.
  • the major clinical advantage of hyperspectral imaging is the delivery of metabolic information derived from the tissue's spectral properties in an easily interpretable image format with high spatial resolution.
  • This 2 -D information allows gradients in biomarker levels to be assessed spatially. Multiple images taken over time allow the gradiant to be measured temporally. This adds new dimensions to the assessment of ulceration risk and tissue healing in that it will allow the physician to target therapy and care to specific at risk areas much earlier than previously possible.
  • the reporting of biomarkers such as oxyHb and deoxyHb levels in tissue individually and in an image format where spatial distributions can be assessed has not been done before. Typically the two numbers are combined in a ratio and reported as percent hemoglobin oxygen saturation (O 2 Sat).
  • MHSI has the clear potential to be developed into a cost effective, easy to use, turn-key camera-based metabolic sensor given the availability and relatively low price of components.
  • MHSI information can be used to predict the onset of foot ulcers before there are clinical indications, and provides early detection, diagnosis, and quantification of progression of microcirculatory complications such as neuropathy in diabetic patients.
  • MHSI technology can evaluate the ulcer and surrounding area to predict whether that will heal or require surgical intervention.
  • the present invention also provides MHSI that is useful in the prediction and monitoring of peripheral venous disease including venous ulcers.
  • MHSI provides anatomically relevant spectral information
  • its use of spectral data of reflected electro-magnetic radiation provides detailed tissue information. Since different types of tissue reflect, absorb, and scatter light differently, in theory the hyperspectral cubes contain enough information to differentiate between tissue types and conditions.
  • MHSI is more robust than conventional analyses since it is based on a few general properties of the spectral profiles (slope, offset, water, oxyHb, deoxyHb, and its ratio) and is therefore flexible with respect to spectral coverage and not sensitive to a particular light wavelength.
  • MHSI is faster than conventional analyses because it uses fast image processing • techniques that allow superposition of absorbance, scattering, and oxygenation information in one pseudo-color image. Visible MHSI is useful because it clearly depicts oxyHb and deoxyHb which are important, physiologically relevant biomarkers in a spatially relevant fashion. Similarly, NIR shows water, oxyHb and deoxyHb.
  • MHSI main purposes include 1) expand human capabilities beyond the ordinary array of senses; 2) expand the human brain capabilities by pre-analyzing the spectral characteristics of the observable subject; 3) perform these tasks with real or near-real time data acquisition.
  • the aim of MHSI is to facilitate the diagnosis and assessment of the metabolic state of tissue.
  • MHSI delivers results in an intuitive form by pairing MHSI pseudo-color image with a high quality color picture composed from the same hyperspectral data. Identification and assessment of a region of interest (ROI) is easily achieved by flipping between color and MHSI images, and zooming onto the ROI. The images can be seen on a computer screen or projector, and/or stored and transported as any other digital information, and/or printed out.
  • the MHSI image preserves the high resolution of the hyperspectral imager thereby allowing further improvement with upgraded hardware.
  • MHSI transcribes vast 3D spectral information sets into one image preserving biological complexity via millions of color shades.
  • the particular color and distinct shape of features in the pseudo-color image allow discrimination between tissue types such as ulcers, callus, intact skin, hematoma, and superficial blood vessels.
  • the algorithm presents oxyHb, deoxyHb and SH S IO 2 to the user to conclude characteristics of the tissue including, but not limited to, discerning whether the tissue is healing or whether it is at a high risk of ulceration.
  • a particular color code contains adequate information for diagnosis and is presented as such.
  • MHSI by itself is not a definite decision making algorithm; it is a tool that a medical professional can use in order to give a confident diagnosis.
  • MHSI contains a decision making algorithm that provides the physician with a diagnosis.
  • MHSI provides currently unavailable information to the doctor, preferably to be used in conjunction with other clinical assessments to provide an accurate diagnosis.
  • MHSI provides images for further analysis by the user. As more information is gathered, a spectral library is preferably compiled to allow MHSI to be a true diagnostic device.
  • MHSI is preferably used to quantify medical therapies in order to measure the effectiveness of new therapeutic " agents or procedures.
  • a typical subject population can be broken down into one of three groups: those that will heal independent of therapy, those that will not heal independent of therapy, and the borderline cases that may benefit from the therapy.
  • MHSI preferably is used to select borderline subjects for these studies where the treatment if effective most likely benefits the subject.
  • MHSI is used to quantify wound progression or prevention in order to identify new therapeutic agents and to develop individual therapeutic regiments depending on subject response.
  • One embodiment of the invention is directed to a medical instrument comprising a first- stage optic responsive to illumination of a tissue, a spectral separator, one or more polarizers, an imaging sensor, a diagnostic processor, a filter control interface, and a general-purpose operating module (FIG 1).
  • a medical instrument comprising a first- stage optic responsive to illumination of a tissue, a spectral separator, one or more polarizers, an imaging sensor, a diagnostic processor, a filter control interface, and a general-purpose operating module (FIG 1).
  • the spectral separator is optically responsive to the first-stage optic and has a control input
  • the polarizer filters a plurality of light beams into a plane of polarization before entering the imaging sensor
  • the imaging sensor is optically responsive to the spectral separator and has an image data output
  • the diagnostic processor comprises an image acquisition interface with an input responsive to the imaging sensor and one or more diagnostic protocol modules wherein each diagnostic protocol module contains a set of instructions for operating the spectral separator and for operating the filter control interface
  • the filter control interface comprises a control output provided to the control input of the spectral separator, which directs the spectral separator independently of the illumination to receive one or more wavelengths of the illumination to provide multispectral or hyperspectral information as determined by the set of instructions provided by the one or more diagnostic protocol module
  • the general-purpose operating module performs filtering and acquiring steps one or more times depending on the set of instructions provided by the one or more diagnostic protocol modules.
  • the instrument may also comprise a second-stage optic responsive to illumination of the tissue.
  • the one or more wavelengths of illumination are one or a combination of UV, visible, NIR, and IR.
  • the multispectral or hyperspectral information determines one or more of the metabolic state of tissue to assess areas at high risk of developing into a foot ulcer or other wounded tissue to assess the potential of an ulcer or the tissue to heal.
  • Preferred embodiments include multispectral or hyperspectral information gathered remotely and noninvasively.
  • an imaging system could be affixed to a wounded area to track its progress over time. Such a system could be attached to or embedded in a dressing, skin covering or a device used to impact wound healing or maintain tissue integrity such as a vacuum suction system or a bed upon which a patient is lying or a shoe, boot or offloading device.
  • Another embodiment is directed to the set of instructions comprising: preprocessing the hyperspectral information, building a visual image, defining a region of interest of the tissue, converting all hyperspectral image intensities into units of optical density by taking a negative logarithm of each decimal base, decomposing a spectra for each pixel into several independent components, determining three planes for an RGB pseudo-color image, determining a sharpness factor plane, converting the RGB pseudo-color image to a hue-saturation-value/intensity (HSV/I) image having a plane, scaling the hue-saturation-value/intensity image plane with the sharpness factor plane, converting the hue-saturation-value/intensity image back to the RGB pseudo-color image, removing outliers beyond a standard deviation and stretching image between 0 and 1, displaying the region of interest in pseudo-colors; and characterizing a metabolic state of the tissue of interest.
  • HSV/I hue-saturation-value/intensity
  • the region of interest may be a pixel, a group of pixels in a prespecified region of a prespecified shape or a handoutlined shape or an entire field of view.
  • determining the three planes for an RGB pseudo-color image comprises one or more characteristic features of the spectra.
  • determining a sharpness factor plane comprises a combination of the images at different wavelengths, preferably by taking a ratio of a yellow plane in the range of about 550-580 nm to a green plane in the range of about 495-525 nm, or by taking a combination of oxyHb and deoxyHb spectral components, or by taking a ratio between a wavelength in the red region in the range 615-710 nm and a wavelength in the yellow region in the range of about 550-580 nm or in the orange region in the range of about 580-615 nm.
  • outliers are removed beyond a standard deviation, preferably three standard deviations.
  • the region of interest is displayed in pseudo-colors, performed with one of in combination with a color photo image of a subject, or in addition to a color photo image of a subject, or by projecting the pseudo-color image onto the observed surface.
  • Another embodiment of the invention is directed to a method for evaluating DFU or area of tissue at risk comprising preprocessing the hyperspectral information, building a visual image, defining a region of interest of the tissue, converting all hyperspectral image intensities into units of optical density by taking a negative logarithm of each decimal base, decomposing a spectra for each pixel into several independent components, determining three planes for an RGB pseudo-color image, determining a sharpness factor plane, converting the RGB pseudo-color image to a hue-saturation-value/intensity (HSV/I) image having a plane, scaling the hue- saturation-value/intensity image plane with the sharpness factor plane, converting the hue- saturation-value/intensity image back to the RGB pseudo-color image, removing outliers beyond a standard deviation and stretching image between 0 and 1, displaying the region of interest in pseudo-colors, and characterizing a metabolic state of the tissue of interest.
  • HSV/I hue-saturation-value/in
  • Another embodiment is directed to a medical instrument comprising an image projector, an illumination source, a remote control device and a real-time data processing package.
  • a medical instrument comprising an image projector, an illumination source, a remote control device and a real-time data processing package.
  • Such a system could project the colorized or other kind of image with relevant information back onto the tissue from which it was taken to assist the physician in diagnosis and treatment such as wound debridement.
  • information can be transmitted to the physician using multiple means, one of such can be a heads-up display.
  • Another embodiment is intended to help tell the doctor level of amputation, safety of debriding tissue, likelihood for tissue to heal, selection and monitoring of specific therapy including topical pharmaceuticals, skin-like coverings, vacuum suction apparatus, systemic pharmaceuticals, adequacy of surgical, stenting or atherectomy procedure, extension of infection vs inflammation of tissue to assist in therapy, identification of organism responsible for local or systemic infection.
  • Yet another embodiment can give information about tissue hydration and potentially information about oxyHb and deoxyHb from deeper tissue using NIR wavelengths. These can be used as a stand alone device or as paired with the more standard Visible wavelength MHSI device as shown in figure 2.
  • Yet another embodiment can derive and present information from changes seen radiating from an area of wounded, ulcerated or otherwise abnormal tissue or from any change in tissue characteristics over a distance .
  • a "gradient map" thus produced can be used to generate a diagnosis, predict the capability of the tissue to heal, define a level for amputation, define the infection vs inflammation, define areas of ischemia, define areas of tissue at risk for ulceration etc.
  • Another embodiment can involve dividing the region of interest into radial segments, pie like segments or a combination of the two or into squares or other geometric shapes and using these segments to compare and contrast different regions of tissue in the same field of view or as compared to a similar field of view on the contralateral extremity or on another part of the body (such as the forearm, the upper leg, etc.).
  • the radial segments can also be compared to similar locations at different time points to demonstrate change over time in response to different therapeutic interventions, changes in tissue physiology , either local, regional or systemic due either to progression or remission of disease or of the effects of topical or systemic medications or therapies.
  • Such measurements can be used to evaluate wound healing, tissue regeneration, angiogenesis, vasculogenesis, arteriogenesis, infection, inflammation, microvascular disease or alterations, or other changes in tissue characteristics or physiology associated with the implementation of negative pressure (vacuum suction applied to the wound), hyperbaric therapy, grafting of autologuous, heterograft, xenograft or biological or synthetic skin substituetes, administration of topical agents including antibiotics, cleansers, growth factors, surgical intervention, angioplasty, stenting, atherectomy, laser therapy, vasodilator therapy, offloading, compression, effects of pressure due to orthotic or prosthetic, effects of electromagnetic, acupuncture, massage, infrared, vibration or other therapies. .
  • Such measurements can be considered as biomarkers representing tissue oxygen delivery and oxygen extraction, tissue oxygenation, tissue perfusion, tissue metabolism or other characteristics correlated with MHSI measurements.
  • Such measurements can be used in association with the implementatiton of hyperbaric therapy delivered to assist in the healing of ulceration in diabetic or other foot ulceration, or other wounds in other parts of the body.
  • the tissue can be monitored before and at specified intervals during therapy or continuously during therapy to determine when the tissue has been adequately modified (oxygenated) by the therapy or that there has been sufficient change in tissue metabolism as described by the MHSI measurements of oxyHb, deoxy Hb or other measured parameters or whether no benefit is being delivered.
  • MHSI can be used to determine the appropriate duration of HBO therapy during a given session and as to whether sufficient benefit has been delivered from a course of therapy that it can safely be discontinued and that the wound will then be likely to heal with more standard methods.
  • MHSI can be used to determine the capability of tissue to heal after debridement and hence the relative safety of pursuing such an approach. Similarly, MHSI can be used to help determine the lowest level of amputation that can be performed with successful healing. Similarly MHSI can be used to determine whether elective surgery to the foot, lower extremity or other body part where evaluation and or quantitation of perfusion, oxygenation, or tissue metabolism would assist in determination of the safety of undertaking such a procedure or the location in which to direct such a procedure. MHSI can be utilized before debridement, amputation or other surgery to make this determination or during debridement, amputatioriv or other surgery to better assess tissue to improve surgical outcomes.
  • Such measurements can be used for the determination of which patients or which wounds are likely to improve with any of the above mentioned therapies, which patients or wounds or portions of wounds are healing or worsening, when a given therapy is sufficient (this could be during or immedialty after application of a therapy such as hyperbaric therapy or a debridement or a particular cleansing or pharmaceutical regimen or after a longer course of several days of therapy such as a vacuum therapy .
  • MHSI criteria can be used to determine when a tissue will accept a skingraft or benf ⁇ t from an allograft or other skin replacement.
  • One embodiment uses a single system that employs light wavelengths derived from the UV, visible, the near infrared, short wave infrared, mid infrared or far infrared portion of the electromagnetic spectrum. Another embodiment uses a system that uses one or more wavelengths from more than one of these wavelength regimes. One such system using wavelengths from more than one of these wavelength groupings is shown in figure two. In other embodiments, a single sensor could be used to collect light from more than one wavelength regime.
  • FIG. 1 A portable hyperspectral imaging apparatus according to an embodiment of the invention is depicted in FIG. 1.
  • Portable apparatus 10 weighs less than 100 pounds, preferably less than 25 pounds, and more preferably less than 10 pounds.
  • the portable apparatus may be battery operated, have some other form of portable power source or more preferably, may have a connector adapted to connect to an existing power source.
  • Portable apparatus 10 comprises an optical acquisition system 36 and a diagnostic processor 38.
  • Optical acquisition system 36 comprises means to acquire broadband data, visible data, ultraviolet data, infrared data, hyperspectral data, or any combination thereof.
  • optical acquiring means comprises a first-stage imaging optic 40, a spectral separator 42, a second-stage optic 44, and an imaging sensor 46.
  • optical acquiring means may be any acquisition system suited for acquiring broadband data, visible data, ultraviolet data, infrared data, hyperspectral data, or any combination thereof.
  • one or more polarizers 41, 43 are included in the acquisition system to compile the light into a plane of polarization before entering the imaging sensor.
  • a calibrator is also included in the system.
  • the first polarizer 43 is placed anywhere in the optical path, preferably in front of the receiving camera 46.
  • the second polarizer 41 is placed in front of illuminating lights 20 such that the incident light polarization is controlled.
  • the incident light is crossed polarized with the light recorded by the camera 46 to reduce specular reflection or polarization at different angles to vary intensity of the reflected light recorded by the camera.
  • the illumination is provided by the remote light(s) 20, various sources tailored or adapted to the need of the instrument, preferably positioned around the light receiving opening of the system, or otherwise placed to afford optimal performance.
  • the light can be a circular array of focused LED lights that emit light at the particular wavelengths (or ranges) that are used in the processing algorithm, or in the ranges of wavelengths (e.g., visible and/or near-infrared).
  • the circular arrangement of the light sources provides even illumination that reduces shadowing.
  • the light wavelength selectivity reduces effect of the observation on the observing subject.
  • the configuration may also vary depending on the particular needs and operation of the system.
  • an optical head is mounted to the wall of the examination room, more preferably, an overhead light structure is located in the operating room, or more preferably, the system has a portable table with an observational window overlooking the operating site.
  • the preferred embodiment may also be used as part of another instrument.
  • another instrument For example, as an adjunct to an endoscope.
  • the first-stage optic receives light collected from a tissue sample through a polarizer and focuses the light onto the surface of the spectral separator.
  • the spectral separator is a liquid crystal tunable filter (LCTF).
  • LCTF 42 is a programmable filter that sequentially provides light from selected wavelength bands with small (for example, 7-10 nm) bandwidth from the light collected from the sample.
  • Second-stage optic 44 receives the narrow band of light passing through the spectral separator and focuses the light onto the image sensor 46.
  • the image sensor is preferably, although not necessarily, a two-dimensional array sensor, such as a charge-coupled device array (CCD) or CMOS, which delivers an image signal to the diagnostic processor 38.
  • CCD charge-coupled device array
  • CMOS complementary metal-coupled device array
  • Diagnostic processor 38 includes an image acquisition interface 50, that has an input responsive to an output of the image sensor 46 and an output provided to a general-purpose operating module 54.
  • the general-purpose operating module includes routines that perform image processing, and that operates and controls the various parts of the system.
  • the general-purpose operating module also controls the light source(s) (e.g. LED array) allowing for switching on and off during measurement as required by the algorithm.
  • the general-purpose operating module has control output provided to a filter control interface 52, which in turn has an output provided to the spectral separator 42.
  • the general-purpose operating module also interacts with a number of diagnostic protocol modules 56A, 56B, . . . 54N, and has an output provided to a video display.
  • the diagnostic process includes special purpose hardware, general-purpose hardware with special-purpose software, or a combination of the two.
  • the diagnostic processor also includes an input device 58, which is operatively connected to the general-purpose operating module.
  • a storage device 60 and printer 62 also are operatively connected to the general- purpose operating module.
  • a portable or semi-portable apparatus is employed within line of site (or with optical access) of the object or area of interest, e.g., diabetic foot with or without an ulcer, or general area of interest.
  • An operator begins by selecting a diagnostic protocol module using the input device.
  • Each diagnostic protocol module is adapted to detect particular tissue characteristics of the target.
  • the diagnostic module could be specific for diabetes, for peripheral vascular disease, for venous stasis disease or for a combination of these disease states.
  • the apparatus may contain only one diagnostic module adapted for general medical diagnosis.
  • Diagnostic processor 38 responds to the operator's input by obtaining a series of transfer functions and an image processing protocol and an image processing protocol from the selected diagnostic protocol module 56.
  • the diagnostic processor provides the filtering transfer functions to the spectral separator 42 via its filter control interface 52 and then instructs the image acquisition interface 50 to acquire and store the resulting filtered image from the image sensor 46.
  • the general-purpose operating module 54 repeats these filtering and acquiring steps one or more times, depending on the number of filter transfer functions stored in the selected diagnostic protocol module.
  • the filtering transfer functions can represent bandpass, multiple bandpass, or other filter characteristics and can include wavelengths in preferably the UV, preferably the visible, preferably the NIR and preferably, the IR electromagnetic spectrum.
  • the light source delivering light to the target of interest can be filtered as opposed to the returned light collected by the detector.
  • a tunable source delivers the information.
  • both a tunable source and a tunable detector may be utilized.
  • tuning takes the form of LCTF, acousto-optical tunable filter (AOTF), filter wheels, matched filters, diffraction gratings or other spectral separators.
  • the light source may be a tungsten halogen or xenon lamp, but is preferably a light emitting diode (LED).
  • the unique cool illumination provided by the LED prevents overheating of skin which may result in poor imaging resolution.
  • the LED provides sufficient light while producing minimal or no increase in skin temperature.
  • This lighting system in combination with the polarizer allows adequate illumination while preventing surface glare from internal organs and overheating of skin.
  • illumination can arise from more passive sources such as room lights or from sunlight.
  • the image acquisition interface 50 begins processing these image planes based on the image processing protocol from the selected diagnostic protocol module 56N. Processing operations can include general image processing of combined images, such as comparing the relative amplitude of the collected light at different wavelengths, adding amplitudes of the collected light at different wavelengths, or computing other combinations of signals corresponding to the acquired planes.
  • the computed image is displayed on the display 12. Other preferred embodiments include storing the computed image in the storage device 60 or printing the computed image out on printer 62.
  • a calibrator is included in the system.
  • Calibrator has an area colored with a pattern of two (or more) colors.
  • colors are chosen that have a distinct absorption band in the wavelength range similar to oxyHb and deoxyHb - preferably in the range 500-600 nm.
  • the colors are placed into a pattern, preferably, a checker-board pattern, where 1 out of 4 squares has color 1, and 3 out of 4 squares have color2. Thus, approximately 25% of the squares are colorl and 75% of the squares are color2.
  • the system takes a hypercube being slightly out of focus - that provides blurring of colors into each pixel.
  • the system simulates the way the biological mixture (oxyHb + deoxyHb) is observed by using "optical" mixture via combination of pattern (with known spatial concentrations) and analog blurring (defocusing - for speed. Defocusing can also be done in the software through the use of computational filters).
  • diagnostic protocol modules 56, printer 62, display 12, or any combination thereof may be omitted from portable device 10.
  • acquired images are stored in storage device 60 during the medical procedure. At a later time, these images are transferred via a communications link to a second device or computer located at a remote location, for example, hospital medical records, for backup or reviewing at a later time.
  • This second device can have the omitted diagnostic protocol modules, printer, display, or any combination thereof.
  • the stored images are transferred from portable device 10, located in the clinic, via a communications link to a remote second device in real time.
  • the system has facility to project real-time hyperspectral data onto the operation field, region of interest, or viewing window positioned above the operating site. Heads Up Display. Also, display completely separately for remote guidance (i.e. on a wall screen for a group of people to review in real time, or post procedure).
  • the projected information has precise one-to-one mapping to the illuminated surface (e.g. wound, operating surface, tissue) and provides the surgeon or podiatrist with necessary information in efficient and non-distractive way.
  • the images can be zoomed in/out to provide variable magnification.
  • This subsystem consists of the following elements: 1) image projector with field-of-view precisely co-aligned with the field-of-view of the hyperspectral imager, 2) miniature remote control device which allows the surgeon or podiatrist to switch projected image on and off without turning from the site of debridement and change highlight structure and/or translucency on the projected image to improve visibility of the features of interest as well as projected image brightness and intensity, 3) real-time data processing package which constructs projected image based on hyperspectral data and operator/surgeon input, 4) optional viewing window positioned above the operating site that is translucent for real observation or opaque for projecting pseudo-color solution or higher resolution images.
  • the MHSI system consists of three functional modules - a Spectral Imager (SI), supporting Controller and Power Module (CPM) and Control and Data Acquisition Computer (CDAC).
  • the MHSI also includes a thermometer that remotely measures the temperature at the tissue surface.
  • the Spectral Imager is mounted on suspension arm which neutralizes device weight and allows for easy positioning and focusing of the instrument.
  • the suspension arm is attached to wheeled cart which supports CPM and CDAC as well. This configuration is very mobile and permits wide range of device spatial and directional motions.
  • FIG 2 shows the preferred system specifications along with a diagram of our focusing methodology and the optical design of the Spectral Imager.
  • a liquid crystal tunable filters (LCTF 's) was used as the wavelength selector and are coupled to complementary metal oxide semiconductor (CMOS) imaging sensors.
  • CMOS complementary metal oxide semiconductor
  • a major issue in the collection of hyperspectral imaging data is the position and focusing of the instrument. While our Spectral Imaging Module is positioned on a ball joint that allows free rotation and virtually any angle of incidence to the patient, it is imperative that there be a system in place for targeting the image to a particular spot on the tissue and ensuring that the instrument will be at the proper distance from the tissue to achieve optimal focus. Positioning and focusing with our system are facilitated by two mirrored collimated light beams or lasers that cross precisely at the instrument focal plane (FIG 2), and so bringing the spectral imaging module into position where the two light spots overlap on the tissue to insure optimal focus. To achieve precisely calibrated images, the system may use a specially designed calibration pad placed at the focal plane of the system and measured prior to each patient measurement.
  • the calibration pad includes a diffusely reflective surface to quantify the intensity of the illumination at each wavelength and color bars to validate wavelength accuracy of the system. Calibration data measured at a preset time such as during maintenance calibrations can be stored and compared to with each use to decide whether the system is within specifications and should proceed to patient measurements.
  • the system may use a fiducial label or target placed in the field of view which the image registration module can perform a self-alignment procedure before or during the operation as necessary.
  • Devices of the present invention allow for the creation and unique identification of patterns in data that highlight the information of interest.
  • the data sets in this case may be discrete images, each tightly bounded in spectra that can then be analyzed. This is analogous to looking at a scene through various colored lenses, each filtering out all but a particular color, and then a recombining of these images into something new.
  • Such techniques as false color analysis (assigning new colors to an image that don't represent the true color but are an artifact designed to improve the image analysis by a human) are also applicable.
  • optics can be modified to provide a zoom function, or to transition from a micro environment to a macro environment and a macro environment to a micro environment. Further, commercially available features can be added to provide real-time or near real-time functioning. Data analysis can be enhanced by triangulation with two or more optical acquisition systems. Polarization may be used as desired to enhance signatures for various targets.
  • the present invention also encompasses the ability to combine the data in various manners including vision fusion, summation, subtraction and other, more complex processes whereby certain unique signatures for information of interest can be defined so that background data and imagery can be removed, thereby highlighting features or information of interest. This can also be combined with automated ways of noting or highlighting items, areas or information of interest in the display of the information.
  • the hyperspectrally resolved image in the present invention is comprised of a plurality of spectral bands. Each spectral band is adjacent to another forming a continuous set.
  • each spectral band having a bandwidth of less than 50 nm, more preferably less than 30 nm, more preferably less than 20 nm, more preferably, from about 20 - 40 nm, more preferably, from about 20 - 30 nm, more preferably, from about 10 - 20 nm, more preferably from about 10 - 15 nm, and more preferably from about 5 - 12 nm.
  • MHSI medical hyperspectral imager
  • Preprocess the MHSI data Preferably, by removing background radiation by subtracting the calibrated background radiation from each newly acquired image while accounting for uneven light distribution by dividing each image by the reflectance calibrator image and registering images across a hyperspectral cube.
  • FIG 2 shows examples of spectra taking from single pixels at different tissue sites within an image. Tissue sites include connective tissues, oxygenated tissues, muscle, tumor, and blood.
  • the color hue plane as apparent concentration of oxygenated Hb, or deoxygenated Hb, or their mathematical combination, e.g. total Hb, oxygen saturation, etc.
  • the color saturation plane as apparent concentration of oxygenated Hb, or deoxygenated Hb, or their mathematical combination, e.g. total Hb, oxygen saturation, etc.
  • the color intensity (value) plane as reflectance in blue-green-orange region (preferably in the range of light at about 450-580 nm).
  • reduce resolution of hue and saturation color planes by binning the image (e.g. by 2, 3, 4, etc. pixels), or/and by resizing the image, or/and by smoothing the image through filtering higher frequency components out. Interpolate the smoothed color planes on the grid of higher resolution intensity (value) plane.
  • HSV/I hue-saturation- value/intensity
  • Characterize the metabolic state of the tissue of interest e.g. risk for ulceration, potential to heal.
  • the saturation and/or intensity of the assigned color Preferably, by using the saturation and/or intensity of the assigned color and provide a qualitative color scale bar.
  • the set of instructions includes only the steps of preprocessing the hyperspectral information, building a visual image, using the entire field of view, converting all hyperspectral image intensities into units of optical density by taking a negative logarithm of each decimal base, and characterizing a metabolic state of the tissue of interest. More preferably, the set of instructions comprises preprocessing the hyperspectral information, defining a region of interest of the tissue, and characterizing a state of the tissue of interest.
  • Another preferred embodiment entails reducing the hyperspectral data in the spectral dimension into a small set of physiologic parameters involves resolving the spectral images into several linearly independent images (e.g. oxyhemoglobin, deoxyhemoglobin, an offset coefficient encompassing multiple scattering (MS) properties and a slope coefficient) in the visible regime.
  • Another embodiment determines four images (e.g. oxyhemoglobin, deoxyhemoglobin, offset/ scattering coefficient, and water absorption) in the near infrared region of the spectrum.
  • linear regression fit coefficients cl, c2, c3 and c4 will be calculated for reference oxy-Hb, deoxy-Hb, and MS spectra, respectively, for each spectrum (Sij) in an image cube:
  • FIG 3 summarizes the display of the MHSI, showing a schematic diagram explaining the scenarios of low and high oxy and deoxy hemoglobin coefficients as well as a color scale that indicates a color plat that shows the vertical color scale and the horizontal brightness scale.
  • tissue By measuring an MHSI where the oxyhemoglobin component is high and the deoxyhemoglobin component is low (upper left hand corner of FIG 3), it could be concluded that that particular area of tissue has adequate perfusion and oxygenation, and is able to satisfy its metabolic needs with the oxygen that is being delivered. That this tissue has the lowest level of risk for ulceration and the highest probability of healing. If tissue demonstrates a low oxyhemoglobin level in addition to a low deoxyhemoglobin level (lower left corner of FIG 3), this would imply that the tissue was receiving low total volume of blood.
  • tissue demonstrates a low oxyhemoglobin level in addition to a high deoxyhemoglobin level (lower right corner of FIG 3) this would imply that the tissue has metabolic requirements exceeding available oxygen delivery. In both of these regions there is expected to be a higher risk of ulceration or difficulties with wound healing. If the tissue has a high oxyhemoglobin coefficient and also has a high deoxyhemoglobin content, (lower left corner of FIG 3) this tissue was receiving a larger total volume of blood, and that the oxygen extracted from the blood stream was adequate to support tissue metabolism. This could be indicative of inflammation. Our technique will uniquely permit discrimination between each of these disparate physiologic conditions. For example, if the value is faded purple (upper left hand quadrant) the tissue has very high oxygenation, as discussed above, and is very likely to heal. The color map (right) gives an indication of how the MHSI would be represented in an image format.
  • hyperspectral imaging for use in the peripheral vascular and diabetes clinic, designed both to be mobile for ease of use and to facilitate the most accurate data collection possible for this project.
  • This system provides fast and precise measurement of reflectance spectra, and is characterized by high spatial and spectral resolution, the ability to process spectral data in real time. It has been equipped with a turn-key software interface for the user.
  • Proprietary image registration software insures image stability when measuring spectra of animated objects.
  • the system does not rely on external illumination, rather it contains very efficient internal visible (and NIR in certain versions) light sources, which allow to achieve high signal to noise ratios in measured data without putting noticeable heat load on a biological subject (variations in skin temperature during acquisition are on the order of 0.1 C).
  • the values for x and y can be taken from a single pixel or from a ROI defined by the user.
  • the hemoglobin oxygen saturation (S HSI O 2 ), x/(x+y) can be presented for a pixel or a ROI.
  • FIG 4 shows examples of tissue with low and high oxyHb and deoxyHb values corresponding to tissue at risk of ulceration and a wound that is likely to heal, respectively.
  • Another way of presenting linearly independent variables is through their sum and difference: (x+y) and (x-y). The first would be THb, and the second would "hint" on oxygen extraction - which indicates the kind of Hb that is predominant at the site.
  • oxyHb, deoxyHb, and S HSI O 2 represent typical oxyHb, deoxyHb, and S HSI O 2 values for two body positions, forearm and foot, and for various stages of diabetes: nondiabetics, diabetics without peripheral neuropathy, and diabetics with peripheral neuropathy.
  • the value for oxyHb and S HSI O 2 are lower in the feet of diabetic subjects with neuropathy compared to the other two groups, a group at high risk for developing foot ulcers.
  • the values for oxyHb, deoxyHb, and S HSI O 2 depend on body location, that once calibrated can be accounted for by the diagnostic module.
  • MHSI oximetry values at baseline prior to iontophoresis of acetylcholine
  • MHSI provides relevant physiological information at the systemic, regional and local levels.
  • Forearm data measures systemic microvasculature changes since the forearm is not affected by macrovasculature or somatic neuropathy as found in the lower extremities.
  • Dorsal foot measurements are indicative of microvascular and macro vascular effects including atherosclerotic changes occurring in large vessels exacerbated by diabetes.
  • MHSI data from the right and left lower extremity can be compared to help differentiate the stages of the damage.
  • MHSI can be used to find local information that can be associated to the risk of developing a foot ulcer of the progression of disease by examining the area around an ulcer.
  • MHSI can not only be used for determining risk of foot ulceration, but also for determining systemic progression of diabetic microvascular disease. In one embodiment this is determined by mean oxyHb, deoxyHb and/or other values for a region of interest. In another embodiment this is determined by heterogeneity of oxyHb, deoxyHb and/or other values for a region of interest. In another embodiment this is determined by the patterning of oxyHb, deoxyHb and/or other values for a region of interest.
  • this is determined by changes over a given time period within a measurement session (between 1 picosecond and one hour, preferably between 100 microseconds and 10 minutes and more preferably between 100 microseconds and 15 seconds) and in the mean values or patterns of oxyHb, deoxyHb and/or other values for a region of interest.
  • These measurements can be used to determine a diabetes progression index (DPI).
  • a DPI can be calculated by comparing a MHSI value or set of values from a single point in time with another point in time (preferably 1 month to two years, more preferably 2 months to one year and most preferably 3-6 months)
  • HMI Hyperspectral Microvascular Index
  • MSHI images have the ability to differentiate between regions of tissue associated with a present foot ulcer on the basis of biomarkers such as the oxyHb and deoxyHb coefficients.
  • FIG 6 shows an ulcer on the sole of the foot of a type 1 diabetic patient (ulcer 1). From the visible image on the left, little distinguishes one area of the ulcer from another. However when looking at the image with the MHSI, there is obvious discriminatory power between the state of tissue seen in the purple oval, which is likely to heal, and that surrounded by the black oval, which is tissue at risk for further ulceration. It is important to note that the skin on the sole of this patient's feet is highly calloused, with a thick stratum corneum, but one is still able to differentiate tissue based on its spectral signatures.
  • the device works on all naturally or surgically exposed tissue or tissue otherwise visualized with laparoscopy, endoscopy, retinoscopy or other visualization techniques. Ulcer 2 was located on the dorsal surface of the foot, on the patient's big toe (FIG 6). These images further show the ability to differentiate between tissue at risk and tissue likely to heal. Additionally, tissue surrounding a fungal infection on the patient's middle toe (bottom right-hand corner of the image) has an MHSI that can demonstrate inflamed or infected tissue.
  • tissue can be examined for gross features indicative of global risks of complications, such as poor perfusion or the inability of the microcirculation to react and compensate in tissue.
  • iontophoretic application of the vasodilator acetylcholine (ACH) or nitroprusside was used to stimulate the vasodilation of the microvasculature on the dorsal surface of the foot and on the forearm of the patients and measured the reaction with MHSI (FIG 7).
  • ulcer images were divided into 25 concentric circles lmm apart and 8 pie segments forming 200 sectors per ulcer (FIG 8).
  • a radial profile analysis was undertaken where the ulcer center was defined at the first visit, and registered images from subsequent visits to this.
  • OxyHb, deoxyHb, total-hemoglobin and C ⁇ Sat were calculated for each sector.
  • MHSI microvascular abnormalities in the diabetic foot and provides early information assist in managing foot ulceration and predict outcomes in patients with diabetes.
  • MHSI is used to monitor angiogenesis during wound healing.
  • An example of wound healing in a diabetic rabbit wound model shows that during the healing process, images of oxyHb and deoxyHb show patterns that change in shape, area, and amplitude with time. Similar patterns were noted in experimental models of shock, but the changes observed for shock occurred on a shorter time scale; minutes rather than several days as the wound heals. The rabbit's ears were observed at days 1, 2, 5, and 10.
  • MHSI is ideally suited for characterization of the local heterogeneity in oxyHb and deoxyHb and their spatial changes with time. For example, the zone of hyperemia surrounding a wound as measured by the oxyHb coefficient decreases with time in wounds that heal (FIG 9).
  • the black rings denote location of a future wound.
  • the color hue represents apparent oxyHb concentrations
  • color saturation represents apparent deoxyHb concentrations.
  • Both, oxyHb and deoxyHb vary predominantly between 40 and 90 MHSI units (color bar to the right).
  • the remaining images to the right show change in a region of interest 17-by-17 mm (black box in (a) and (b)) over 10 days.
  • the color image was reconstructed from MHSI data, showing a party of the observed area 50-by-40 mm, recorded at the baseline on day 1.
  • the pseudo-image (b) was obtained as a result of hyperspectral processing, showing a distribution of the oxygenated (oxy) and deoxygenated (deoxy) hemoglobin in the underlying tissue at the same time.
  • Panasyuk SV Freeman JE, Cooke WI, Hopmeier MJ, Convertino VA.
  • MHSI medical hyperspectral imaging

Abstract

The invention is directed to methods and systems of hyperspectral and multispectral imaging of medical tissues. In particular, the invention is directed to new devices, tools and processes for the detection and evaluation of diseases and disorders such as, but not limited to diabetes and peripheral vascular disease, that incorporate hyperspectral or multispectral imaging.

Description

HYPERSPECTRAL IMAGING IN DIABETES AND PERIPHERAL VASCULAR DISEASE
Reference to Related Applications
The present application claims priority to U.S. Provisional Patent Application Serial Number 60/667,677 entitled Hyperspectral Imaging in Diabetes, filed April 4, 2005, and U.S. Provisional Patent Application Serial Number 60/xxx,xxx entitled Hyperspectral Imaging of Angiogenesis, filed March 27, 2006which is hereby incorporated by reference.
Background
1. Field of the Invention
The invention is directed to methods and systems of hyperspectral and multispectral imaging of medical tissues. In particular, the invention is directed to new devices, tools and processes for the detection and evaluation of diseases and disorders such as diabetes and peripheral vascular disease that incorporate hyperspectral/multispectral imaging.
2. Background of the Invention
Diabetes afflicts an estimated 194 million people worldwide, affecting 7.9% of Americans (over 21 million people) and 7.8% of Europeans. Between 85% and 95% of all diabetics suffer from Type 2 diabetes, although nearly 5 million people worldwide suffer from Type 1 diabetes, affecting an estimated 1.27 million people in Europe and another 1.04 million people in the United States1. Both Type 1 and Type 2 diabetic patients are at higher risk for a wide array of complications including heart disease, kidney disease (e.g. nephropathy), ocular diseases (e.g. glaucoma), and neuropathy and nerve damages to name a few2. The feet of diabetic patients are at risk for a wide array of complications, which are discussed below. Problems with the foot that affect the ambulatory nature of the patient are not only important from the standpoint of physical risk, but also convey an emotional risk as well, as these problems disrupt the fundamental independence of the patient by limiting his or her ability to walk.
Peripheral arterial disease (PAD) affects primarily people older than 55. There are currently 59.3 million Americans older than 55, and over 12 million of them have symptomatic peripheral vascular disease. It is estimated that only 20% of all patients with PAD have been diagnosed at this time. This represents a dramatically underpenetrated market. Although pharmacologic treatments for PAD have traditionally been poor, 2.1 million nevertheless receive pharmacologic treatment for the symptoms of PAD, and current diagnostic tests are not considered to be very sensitive indicators of disease progression or response to therapy. Additionally, 443,000 patients undergo vascular procedures such as peripheral arterial bypass surgery (100,000) or peripheral angioplasty (343,000) annually and are candidates for pre and post surgical testing. One difficulty in diagnosing PAD is that in the general population, only about 10% of persons with PAD experience classic symptoms of intermittent claudication. About 40% of patients do not complain of leg pain, while the remaining 50% have leg symptoms which differ from classic claudication.
Relying on medial history and physical examination alone is unsatisfactory. In one study, 44 percent of PAD diagnoses were false positive and 19 percent were false negative when medical history and physical examination alone were used.3 For this reason, physicians have looked for other means to help in providing diagnosis. As in the case of diabetic foot disease, current technologies have fallen short. Nonetheless, patients are frequently sent to peripheral vascular laboratories for non-invasive studies. While the test results are known to be inaccurate, these results do provide some additional information to physicians for assistance in diagnosis or treatment decisions.
Another problem face by physicians is disease of the peripheral veins. Venous occlusive disease due to incompetent valves in veins designed to prevent backfiow and deep vein thrombosis, results in venous congestion and eventually stasis ulcers. Approximately 70% of leg ulcers are due to venous occlusion. Many of these ulcers are found at the medial malleolus. The foot is generally swollen and the skin near the ulcer site is brownish in appearance.
Pathology
Diabetic feet are at risk for a wide range of pathologies, including microcirculatory changes, peripheral vascular disease, ulceration, infection, deep tissue destruction and metabolic complications. The development of an ulcer in the diabetic foot is commonly a result of a break in the barrier between the dermis of the skin and the subcutaneous fat that cushions the foot during ambulation. This, in turn, can lead to increased pressure on the dermis, resulting in tissue ischemia and eventual death, and ultimately result in an ulcer.4 There are a number of factors that weigh heavily in the process of ulceration5 - affecting different aspects of the foot - that lead to a combination of effects that greatly increase the risk of ulceration:6
Neuropathy - Results in a loss of protective sensation in the foot, exposing patients to undue, sudden or repetitive stress. Can cause a lack of awareness of damage to the foot as it be occurs and physical defects and deformities7 which lead to even greater physical stresses on the foot. It can also lead to increased risk of cracking and the development of fissures in calluses, creating a potential entry for bacteria and increased risk of infection.8
Microcirculatory Changes - Often seen in association with hyperglycemic damage.9 Functional abnormalities occur at several levels, including hyaline basement membrane thickening and capillary leakage. On a histologic level, it is well known that diabetes causes a thickening of the endothelial basement membrane which in turn may lead to impaired endothelial cell function.
Musculoskeletal Abnormalities - Include altered foot mechanics, limited joint mobility, and bony deformities, and can lead to harmful changes in biomechanics and gait. This increases pressures associated with various regions of the foot. Alteration or atrophy of fat pads from increased pressure can lead to skin loss or callus, both of which increase the risk of ulceration by two orders of magnitude.
Peripheral Vascular Disease - Caused by atherosclerotic obstruction of large vessels resulting in arterial insufficiency10 is common in the elderly populations and is yet more common and severe in diabetics.11 Diabetics may develop atherosclerotic disease of large-sized and medium-sized arteries, however, significant atherosclerotic disease of the infrapopliteal segments is particularly common. The reason for this is thought to result from a number of metabolic abnormalities in diabetics, including high LDL and VLDL levels, elevated plasma von Willebrand factor, inhibition of prostacyclin synthesis, elevated plasma fibrinogen levels, and increased platelet adhesiveness.
Venous Disease - Caused by incompetent valves controlling backfiow between the deep veins and the more superficial veins or thrombosis of the deep veins. Venous occlusions are typically observed in the elderly who typically presented with swollen lower extremities and foot ulcers typically at the medial malleolus. Previous studies have shown that a foot ulcer precedes roughly 85% of all lower extremity amputations in diabetic patients ' and that 15% of all diabetic patients will develop a foot ulcer during the course of their lifetimes.14 More than 88,000 amputations performed annually on diabetics, 15 and roughly an additional 30,000 amputations are performed on non- diabetics, mostly related to peripheral vascular disease. Estimations have shown that between 2- 6% of diabetic patients will develop a foot ulcer every year13' 16 and that the attributable cost for an adult male between 40 and 65 years old is over $27,000 (1995 US dollars) for the two years after diagnosis of the foot ulcer.16 In conjunction with the increased total costs of care, Ramsey et al showed that diabetic patients incurred more visits to the emergency room (more than twice as many as control patients), more outpatient hospital visits (between 2X and 3X as many as control subjects) and more inpatient hospital days (between 3X and 4X as many as control patients) during the course of an average year.
Foot pathology is major source of morbidity among diabetics and is a leading cause of hospitalization. The infected and/or ischemic diabetic foot ulcer accounts for about 25% of all hospital days among people with diabetes, and the costs of foot disorder diagnosis and management are estimated at several billion dollars annually.16' 17
Current Diagnostic Procedures
The first step in the assessment of the diabetic foot is the clinical examination18' 19. All patients with diabetes require a thorough pedal examination at least once a year, even without signs of neuropathy. Evaluation of the diabetic patient with peripheral vascular disease should include a thorough medical history, vascular history, physical examination, neurologic evaluation for neuropathy and a thorough vascular examination.20
The next step in the work up of a patient with significant peripheral vascular or diabetic foot disease is non-invasive testing.21 Current clinical practice can include ankle brachial index (ABI), transcutaneous oxygen measurements (TcPO2), pulse volume recordings (PVR) and laser Doppler flowmetry. All of these clinical assessments are highly subjective with significant inter- and intra-observer variability especially in longitudinal studies. None of these methods are discriminatory for feet at risk, and none of them provide any information about the spatial variability across the foot. Doppler ultrasound with B-mode realtime imaging is typically used to diagnose deep vein thrombosis while photo and air plethysmography are used to measure volume refill rates as a means of locating and diagnosing valvular insufficiency. Currently there is no method to accurately assess the predisposition to serious foot complications, to define the real extent of disease or to track the efficacy of therapeutics over time.
Summary of Invention
The present invention overcomes the problems and disadvantages associated with current strategies and designs, and provides new tools and methods for detecting tissue at risk of developing into an ulcer, for detecting problems with diabetic foot disease, and for evaluating the potential for wounds to heal.
One embodiment of the invention is directed to a medical instrument comprising a first stage optic responsive to illumination of tissue, a spectral separator, one or more polarizers, an imaging sensor, a diagnostic processor, a filter control interface, a general purpose operating module to assess the state of tissue in diabetic subjects following a set of instructions, and a calibrator. Preferably, the instrument further comprises a second stage optic responsive to illumination of tissue. Preferably, the set of instructions comprises preprocessing hyperspectral information, building a visual image, defining a region of interest in tissue, converting the visual image into units of optical density by taking a negative logarithm of each decimal base, decomposing a spectra for each pixel into several independent components, determining three planes for an RGB pseudo-color image, determining a sharpness factor plane, converting the RGB pseudo-color image to a hue-saturation-value/intensity image having a plane, adjusting the hue- saturation-value/intensity image plane with the sharpness factor plane, converting the hue- saturation-value/intensity image back to the RGB pseudo-color image, removing outliers beyond a standard deviation and stretching image between 0 and 1, displaying the region of interest in pseudo-colors; and characterizing a metabolic state of the tissue of interest.
Preferably, the region of interest is one of a pixel, a specified region or an entire field of view. Preferably, determining three planes for an RGB pseudo-color image comprises one or more characteristic features of the spectra, determining a sharpness factor plane comprises a combination of the images at different wavelengths, removing outliers beyond a standard deviation comprises three standard deviations, displaying the region of interest in pseudo-colors comprises one of performing one in combination with a color photoimage of a subject, in addition to a color photo image of a subject, and projecting onto the tissue of interest. Preferably, defining the color intensity plane as apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation, defining the color intensity plane as reflectance in blue-green-orange region, adjusting the hue saturation comprises adjusting a color resolution of the pseudo-color image according to quality of apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation, adjusting the hue saturation further comprises one or a combination of reducing resolution of hue and saturation color planes by binning the image, resizing the image, and smoothing the image through filtering higher frequency components out, and further interpolating the smoothed color planes on a grid of higher resolution intensity plane.
Another embodiment of the invention is directed to a method for assessing the state of tissue of a diabetic subject comprising, preprocessing hyperspectral information, building a visual image, defining a region of interest in tissue, converting the visual image into units of optical density by taking a negative logarithm of each decimal base, decomposing a spectra for each pixel into several independent components, determining three planes for an RGB pseudocolor image, determining a sharpness factor plane, converting the RGB pseudo-color image to a hue-saturation- value/intensity image having a plane, adjusting the hue- saturation- value/intensity image plane with the sharpness factor plane, converting the hue-saturation-value/intensity image back to the RGB pseudo-color image, removing outliers beyond a standard deviation and stretching image between 0 and 1, displaying the region of interest in pseudo-colors; and characterizing a metabolic state of the tissue of interest.
Preferably, the region of interest is one of a pixel, a specified region or an entire field of view. Preferably, determining three planes for an RGB pseudo-color image comprises one or more characteristic features of the spectra, determining a sharpness factor plane comprises a combination of the images at different wavelengths, removing outliers beyond a standard deviation comprises three standard deviations, displaying the region of interest in pseudo-colors comprises one of performing one in combination with a color photoimage of a subject, in addition to a color photo image of a subject, and projecting onto the tissue of interest.
Preferably, defining the color intensity plane as apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation, defining the color intensity plane as reflectance in blue-green-orange region, adjusting the hue saturation comprises adjusting a color resolution of the pseudo-color image according to quality of apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation, adjusting the hue saturation further comprises one or a combination of reducing resolution of hue and saturation color planes by binning the image, resizing the image, and smoothing the image through filtering higher frequency components out, and further interpolating the smoothed color planes on a grid of higher resolution intensity plane.
Other embodiments and advantages of the invention are set forth in part in the description, which follows, and in part, may be obvious from this description, or may be learned from the practice of the invention.
Brief Description of Drawings:
FIG 1: Block diagram depicting a portable hyperspectral imaging apparatus. FIG 2: Basic specifications of the MHSI system.
FIG 3: OxyHb and DeoxyHb HSV/I color chart. Schematic representation of the MHSI display (left) showing the interplay between the oxyHb and deoxyHb coefficients and describing some of the potential physiological consequences of values of the MHSI.
FIG 4: Representative data from dorsal surface of foot showing individual oxyHb and deoxyHb values and how they can be used to evaluate regions of the tissue.
FIG 5: Representative data from tissue showing sensitivity of MHSI to drug-induced changes in the vasculature, (left to right) Visible image of foot surface post iontophoresis (IP), representative spectra pre and post iontophoresis with Acetylcholine (IP) showing greater oxyHb levels after IP. Images of increased oxyHb coefficient ring where IP occurred, image of deoxyHb, showing little change post IP.
FIG 6: Representative data from an ulcer located on the sole (ulcer 1) and dorsal surface (ulcer 2) of the foot.
FIG 7: MHSI information from the soles and dorsal surfaces of four patients. Each row of images represents data from one patient. The two columns on the left represent data from the soles of the feet, while the columns on the right represent data from the dorsal surfaces of the feet.
FIG 8: MHSI image of diabetic foot ulcer with 200 segment radial profile. FIG 9: MHSI of wounds during healing.
FIG 10: color image (a) and pseudo-color image (b) of a rabbit ear Description of the Invention
Background of Hyperspectral Imaging
HSI or hyperspectral imaging is a novel method of "imaging spectroscopy" that generates a "gradient map" of a region of interest based on local chemical composition. HSI has been used in satellite investigation of suspected chemical weapons production areas22, geological features23, and the condition of agricultural fields24 and has recently been applied to the investigation of physiologic and pathologic changes in living tissue in animal and human studies to provide information as to the health or disease of tissue that is otherwise unavailable.25 MHSI for medical applications (MHSI) has been shown to accurately predict viability and survival of tissue deprived of adequate perfusion, and to differentiate diseased (e.g. tumor) and ischemic tissue from normal tissue.
Spectroscopy is used in medicine to monitor metabolic status in a variety of tissues. One of the most common spectroscopic applications is in pulse oximetry, which utilize the different oxyhemoglobin (oxyHb) and deoxyhemoglobin (deoxyHb) absorption bands to estimate arterial hemoglobin oxygen saturation.28 One of the drawbacks of these systems is that they provide no information about the spatial distribution or heterogeneity of the data. In addition, these systems report the ratio of oxyHb and deoxyHb together losing diagnostic information that can be garnered by evaluating the state of the individual components. Such spatial information for the individual components and the ratio is provided by HSI, which is considered a method of "imaging spectroscopy", where the multi- dimensional (spatial & spectral) data are represented in what is called a "hypercube." The spectrum of reflected light is acquired for each pixel in a region, and each such spectrum is subjected to standard analysis. This allows the creation of an image based on the metabolic state of the region of interest (ROI). In vivo, MHSI has been used to demonstrate otherwise unobserved changes in pathophysiology. Specific studies have evaluated the macroscopic distribution of skin oxygen saturation,29 the in-situ detection of tumor during breast cancer resection in the rat,27 the determination of tissue viability following plastic surgery & burns,30' 3I claudication and foot ulcers in diabetic patients,32"37 and applications to shock and lower body negative pressure (LBNP) in pigs and humans, respectively.38"40 In a skin pedicle flap model in the rat, tissue that has insufficient oxygenation to remain viable is readily apparent from local oxygen saturation maps calculated from hyperspectral images acquired immediately following surgery; by contrast, clinical signs of impending necrosis do not become apparent for 12 hours after surgery.41
Non-invasive measurements of oxygen or blood flow have been demonstrated previously, with investigators using thermometry,42 point diffuse reflectance spectroscopy,43' 44 and laser Doppler imaging.45 Sheffield et al, have also reviewed laser Doppler and TcPO2 measurements and their specific applications to wound healing.46 While other techniques have been utilized in both the research lab and the clinic and have the advantage of a longer experience base, MHSI is superior to other technologies and can provide predictive information on the onset and outcomes of diabetic foot ulcers, venous stasis ulcers and peripheral vascular disease.
Because MHSI has the ability to show anatomically relevant information that is useful in the assessment of local, regional and systemic disease. This is important in the assessment of people with diabetes and/or peripheral vascular disease. MHSI shows the oxygen delivery and oxygen extraction of each pixel in the image collected. These images with pixels ranging from 20 microns to 120 microns have been useful in several ways. In the case of systemic disease, MHSI shows the effects on the microcirculation of systemic diabetes, smoking, a variety of medications such as all of the classes of antihypertensives (ACE inhibitors, ARBs, Beta blockers, Peripheral arterial and arteriolar dilators), vasodilators (such as nitroglycerine, quinine, morphine), vasoconstrictors (including coffee, tobacco, pseudephedrine, Ritalin, epinephrine, levophedrine, neosynepherine), state of hydration, state of cardiac function (baseline, exercise, congestive heart failure), systemic infection or sepsis as well as other viral or bacterial infections and parasitic diseases. The size of the pixels used is important in that it is smaller than the spacing of the perforating arterioles (~0.8 mm)47 of the dermis and therefore permits the visualization of the distribution of mottling or other patterns associated with the anatomy of the microcirculation and its responses. In the case of the use of MHSI for regional assessment, in addition to the above systemic effects at play, the image delivers information about the oxygen delivery and oxygen extraction for a particular region as it is influenced by blood flow through the larger vessels of that region of the body. For example an image of the top of the foot reflects both the systemic microvascular status and the status of the large (macrovascular) vessels supplying the leg. This can reflect atherosclerotic or other blockage of the vessel, potential injury to the vessel with narrowing, or spasm of some of the smaller vessels. It can also reflect other regionalized processes such as neuropathy or venous occlusion or compromise or stasis. In the case of local disease MHSI shows the actual effect of the combination of systemic, regional and local effects on small pieces of tissue. This combines the effects of systemic and regional effects described above with the effects of local influences on the tissue including pressure, neuropathy, localized small vessel occlusion, localized trauma or wounding, pressure sore, inflammation, and wound healing. Angiogenesis during wound healing is readily monitored with MHSI.
Wounds other than on the foot can be similarly assessed, such as sacral decubiti, other areas of pressure necrosis, prosthesis stumps, skin flap tissue before, after or during surgery, areas of tissue breakdown after surgery, and burn injuries. Current optical methods for evaluating tissues for the conditions described above include:
Laser Doppler (LD) - In early iontophoresis experiments as well as recent efforts both LD and MHSI data were collected, and some changes in our images (total hemoglobin) are primarily a consequence of changes in perfusion which was roughly correlated to LD. However, important other changes in MHSI images that report specifically O2 extraction and tissue metabolism (02Sat) are not related to perfusion or LD readings per-se. Superior spatial resolution with MHSI, and O2 extraction information adds highly important clinical information.
Transcutaneous PO2 (TCPO2) - TcPO2 data collected in subjects with peripheral vascular disease and ischemia study as well as in patients with diabetes both with and without foot ulcers. TcPO2 measurements appeared cumbersome, lengthy (-20-30 minutes), highly operator dependant, and carried data only from skin directly under the probe (with little ability to distinguish the spatial characteristics of the ischemic area). While TcPO2 has been shown to carry statistically significant information in terms of quantifying tissue at risk for ulceration,48 TcPO2 was not encouraging as a useful clinical device.
Non-imaging techniques - Techniques such as near-infrared absorption spectroscopy (NIRS) or TcPO2, rely on measurements at a single point in tissue which may not accurately reflect overall tissue condition or provide anatomically relevant data, and probe placement on the skin can alter blood flow and cannot deliver accurate information in the area of an ulcer or directly surrounding it. Because MHSI is truly remote sensing, data are acquired at a distance, eliminating probe placement errors and allowing the investigation of the wound itself, which some techniques can not accomplish due to infection risk.
In short, analysis of the present invention supports the following conclusions:
1. Level of oxygenated hemoglobin in the tissue of arms and feet of diabetic subjects is lower than the level of oxygenated hemoglobin in the skin of control subjects. This is statistically significant result with separation between diabetics and controls.36
2. Oxyhemoglobin in the arms and feet of ulcerated subjects is lower than oxyhemoglobin in diabetics without the ulceration. The strong signal suggests ability to distinguish diabetics at lower and high risk.
3. Oxygen saturation level in the skin of arms and feet of diabetics is lower than oxygen saturation in the skin of controls. This is at a statistically significant level allowing separation between diabetics and controls.
4. MHSI quantitatively assesses different areas of tissue metabolism on both dorsal and plantar foot surfaces of any curvature.
5. MHSI evaluates state of tissue as a function of distance away from ulcer to assess the viability of surrounding tissue, and evaluate the degree of risk of further ulceration.
6. MHSI can be classified with a 4-quadrant system to determine the metabolic state of tissue using oxygen delivery and oxygen extraction: low/low, low/high, high/high, and high/low. This metric is used in distinguishing healthy tissue from ulcerated, or from a tissue at risk of ulceration.
I l 7. MHSI is a unique visualization method that produces an image that combines spatial information from three independent parameters characterizing tissue: oxygenated and deoxygenated hemoglobin concentrations and light absorption.
8. MHSI evaluates skin metabolism at high resolution of 20-120 microns per image pixel.
9. Specific MHSI regions associated with the margins of the ulcer correlate to inflammation (and/or infection).
10. Areas of decreased MHSI indicate tissue at risk for non-healing, ulcer extension, or primary ulceration.
11. MHSI differentiates between regions of tissue associated with a present foot ulcer on the basis of biomarkers such as oxyHb and deoxyHb coefficients.
12. MHSI evaluates temporal changes in oxygen delivery and extraction to particular areas, both, on local and systemic scale. The trend in the change of oxyHb and deoxyHb are used to predict healing status of a wound/ulcer as well as progression of diabetic complications.
13. Specific results from MHSI are indicative of inflamed tissue.
14. MHSI examines tissue for gross features that may be indicative of global risks of complications, such as poor perfusion or the inability of the microcirculation to react and compensate in tissue.
15. MHSI has potential in diagnosing global microcirculatory insufficiencies and impacting on other complications of diabetes associated with the microvasculature besides foot ulcers.
MHSI is superior to other modalities for assessing the healing potential of tissue adjacent to ulcers. MHSI provides more direct measurements of oxyHb and deoxyHb activities of the affected tissue. Hence, the discrimination is not markedly improved by adding iontophoresis results to refine prediction as is required for Laser Doppler to do so. MHSI has significant advantages over laser Doppler and TcPO2 measurements. Whereas MHSI is able to deliver spatially relevant data with high spatial resolution, TcPO2 delivers only single point data. Laser Doppler data has poor spatial resolution and is frequently reported as a single mean numerical value across the region of interest. The major clinical advantage of hyperspectral imaging is the delivery of metabolic information derived from the tissue's spectral properties in an easily interpretable image format with high spatial resolution. This 2 -D information allows gradients in biomarker levels to be assessed spatially. Multiple images taken over time allow the gradiant to be measured temporally. This adds new dimensions to the assessment of ulceration risk and tissue healing in that it will allow the physician to target therapy and care to specific at risk areas much earlier than previously possible. The reporting of biomarkers such as oxyHb and deoxyHb levels in tissue individually and in an image format where spatial distributions can be assessed has not been done before. Typically the two numbers are combined in a ratio and reported as percent hemoglobin oxygen saturation (O2Sat). MHSI has the clear potential to be developed into a cost effective, easy to use, turn-key camera-based metabolic sensor given the availability and relatively low price of components.
Surprisingly, MHSI information according to this invention can be used to predict the onset of foot ulcers before there are clinical indications, and provides early detection, diagnosis, and quantification of progression of microcirculatory complications such as neuropathy in diabetic patients. For patients with foot ulcers, MHSI technology can evaluate the ulcer and surrounding area to predict whether that will heal or require surgical intervention. The present invention also provides MHSI that is useful in the prediction and monitoring of peripheral venous disease including venous ulcers.
There are many advantages to using MHSI. Not only does MHSI provide anatomically relevant spectral information, its use of spectral data of reflected electro-magnetic radiation (ultraviolet - UV, visible, near infrared - NIR, and infrared - IR) provides detailed tissue information. Since different types of tissue reflect, absorb, and scatter light differently, in theory the hyperspectral cubes contain enough information to differentiate between tissue types and conditions. MHSI is more robust than conventional analyses since it is based on a few general properties of the spectral profiles (slope, offset, water, oxyHb, deoxyHb, and its ratio) and is therefore flexible with respect to spectral coverage and not sensitive to a particular light wavelength. MHSI is faster than conventional analyses because it uses fast image processing techniques that allow superposition of absorbance, scattering, and oxygenation information in one pseudo-color image. Visible MHSI is useful because it clearly depicts oxyHb and deoxyHb which are important, physiologically relevant biomarkers in a spatially relevant fashion. Similarly, NIR shows water, oxyHb and deoxyHb.
The simplicity of the presented false color images representing distribution of various chemical species, either singly or in combination (such as ratioed), or in other more sophisticated image processing techniques allow for the display of results in real to near-real time. Another advantage of MHSI is easy interpretation. Color changes show the different tissue types or condition, but the distinction is not a yes/no type. MHSI color scheme allows the surgeon or podiatrist to differentiate between different tissue types and states. In addition, the color and the shape of structures depict different composition and level of viability of the tissue. The data is then represented in a developed MHSI standard format. OxyHb and deoxyHb are presented in a format similar to a blood pressure reading that is easy for physicians to understand. Additionally, a tissue oxygen saturation value denoted as SHSIO2 is also provided.
MHSI main purposes include 1) expand human capabilities beyond the ordinary array of senses; 2) expand the human brain capabilities by pre-analyzing the spectral characteristics of the observable subject; 3) perform these tasks with real or near-real time data acquisition. In summary, the aim of MHSI is to facilitate the diagnosis and assessment of the metabolic state of tissue.
Results of analysis have to be presented in an easily accessible and interpretable form. MHSI delivers results in an intuitive form by pairing MHSI pseudo-color image with a high quality color picture composed from the same hyperspectral data. Identification and assessment of a region of interest (ROI) is easily achieved by flipping between color and MHSI images, and zooming onto the ROI. The images can be seen on a computer screen or projector, and/or stored and transported as any other digital information, and/or printed out. The MHSI image preserves the high resolution of the hyperspectral imager thereby allowing further improvement with upgraded hardware.
Additionally, MHSI transcribes vast 3D spectral information sets into one image preserving biological complexity via millions of color shades. The particular color and distinct shape of features in the pseudo-color image allow discrimination between tissue types such as ulcers, callus, intact skin, hematoma, and superficial blood vessels. Initially, the algorithm presents oxyHb, deoxyHb and SHSIO2 to the user to conclude characteristics of the tissue including, but not limited to, discerning whether the tissue is healing or whether it is at a high risk of ulceration. In another embodiment, a particular color code contains adequate information for diagnosis and is presented as such. In one iteration, MHSI by itself is not a definite decision making algorithm; it is a tool that a medical professional can use in order to give a confident diagnosis. In another iteration, MHSI contains a decision making algorithm that provides the physician with a diagnosis.
Due to the complexity of the biological system, medical personnel desire as much information as possible in order to make the most-reliable diagnosis. MHSI provides currently unavailable information to the doctor, preferably to be used in conjunction with other clinical assessments to provide an accurate diagnosis. MHSI provides images for further analysis by the user. As more information is gathered, a spectral library is preferably compiled to allow MHSI to be a true diagnostic device.
MHSI is preferably used to quantify medical therapies in order to measure the effectiveness of new therapeutic" agents or procedures. For example, in wound healing studies, a typical subject population can be broken down into one of three groups: those that will heal independent of therapy, those that will not heal independent of therapy, and the borderline cases that may benefit from the therapy. MHSI preferably is used to select borderline subjects for these studies where the treatment if effective most likely benefits the subject. MHSI is used to quantify wound progression or prevention in order to identify new therapeutic agents and to develop individual therapeutic regiments depending on subject response.
One embodiment of the invention is directed to a medical instrument comprising a first- stage optic responsive to illumination of a tissue, a spectral separator, one or more polarizers, an imaging sensor, a diagnostic processor, a filter control interface, and a general-purpose operating module (FIG 1). Preferably, the spectral separator is optically responsive to the first-stage optic and has a control input, the polarizer filters a plurality of light beams into a plane of polarization before entering the imaging sensor, the imaging sensor is optically responsive to the spectral separator and has an image data output, the diagnostic processor comprises an image acquisition interface with an input responsive to the imaging sensor and one or more diagnostic protocol modules wherein each diagnostic protocol module contains a set of instructions for operating the spectral separator and for operating the filter control interface, the filter control interface comprises a control output provided to the control input of the spectral separator, which directs the spectral separator independently of the illumination to receive one or more wavelengths of the illumination to provide multispectral or hyperspectral information as determined by the set of instructions provided by the one or more diagnostic protocol module, and the general-purpose operating module performs filtering and acquiring steps one or more times depending on the set of instructions provided by the one or more diagnostic protocol modules.
The instrument may also comprise a second-stage optic responsive to illumination of the tissue. Preferably, the one or more wavelengths of illumination are one or a combination of UV, visible, NIR, and IR. In preferred embodiments, the multispectral or hyperspectral information determines one or more of the metabolic state of tissue to assess areas at high risk of developing into a foot ulcer or other wounded tissue to assess the potential of an ulcer or the tissue to heal. Preferred embodiments include multispectral or hyperspectral information gathered remotely and noninvasively. Alternatively, an imaging system could be affixed to a wounded area to track its progress over time. Such a system could be attached to or embedded in a dressing, skin covering or a device used to impact wound healing or maintain tissue integrity such as a vacuum suction system or a bed upon which a patient is lying or a shoe, boot or offloading device.
Another embodiment is directed to the set of instructions comprising: preprocessing the hyperspectral information, building a visual image, defining a region of interest of the tissue, converting all hyperspectral image intensities into units of optical density by taking a negative logarithm of each decimal base, decomposing a spectra for each pixel into several independent components, determining three planes for an RGB pseudo-color image, determining a sharpness factor plane, converting the RGB pseudo-color image to a hue-saturation-value/intensity (HSV/I) image having a plane, scaling the hue-saturation-value/intensity image plane with the sharpness factor plane, converting the hue-saturation-value/intensity image back to the RGB pseudo-color image, removing outliers beyond a standard deviation and stretching image between 0 and 1, displaying the region of interest in pseudo-colors; and characterizing a metabolic state of the tissue of interest.
The region of interest may be a pixel, a group of pixels in a prespecified region of a prespecified shape or a handoutlined shape or an entire field of view. Preferably, determining the three planes for an RGB pseudo-color image comprises one or more characteristic features of the spectra. Preferably, determining a sharpness factor plane comprises a combination of the images at different wavelengths, preferably by taking a ratio of a yellow plane in the range of about 550-580 nm to a green plane in the range of about 495-525 nm, or by taking a combination of oxyHb and deoxyHb spectral components, or by taking a ratio between a wavelength in the red region in the range 615-710 nm and a wavelength in the yellow region in the range of about 550-580 nm or in the orange region in the range of about 580-615 nm. Preferably, outliers are removed beyond a standard deviation, preferably three standard deviations. The region of interest is displayed in pseudo-colors, performed with one of in combination with a color photo image of a subject, or in addition to a color photo image of a subject, or by projecting the pseudo-color image onto the observed surface.
Another embodiment of the invention is directed to a method for evaluating DFU or area of tissue at risk comprising preprocessing the hyperspectral information, building a visual image, defining a region of interest of the tissue, converting all hyperspectral image intensities into units of optical density by taking a negative logarithm of each decimal base, decomposing a spectra for each pixel into several independent components, determining three planes for an RGB pseudo-color image, determining a sharpness factor plane, converting the RGB pseudo-color image to a hue-saturation-value/intensity (HSV/I) image having a plane, scaling the hue- saturation-value/intensity image plane with the sharpness factor plane, converting the hue- saturation-value/intensity image back to the RGB pseudo-color image, removing outliers beyond a standard deviation and stretching image between 0 and 1, displaying the region of interest in pseudo-colors, and characterizing a metabolic state of the tissue of interest.
Another embodiment is directed to a medical instrument comprising an image projector, an illumination source, a remote control device and a real-time data processing package. Such a system could project the colorized or other kind of image with relevant information back onto the tissue from which it was taken to assist the physician in diagnosis and treatment such as wound debridement. Alternatively, information can be transmitted to the physician using multiple means, one of such can be a heads-up display.
Another embodiment is intended to help tell the doctor level of amputation, safety of debriding tissue, likelihood for tissue to heal, selection and monitoring of specific therapy including topical pharmaceuticals, skin-like coverings, vacuum suction apparatus, systemic pharmaceuticals, adequacy of surgical, stenting or atherectomy procedure, extension of infection vs inflammation of tissue to assist in therapy, identification of organism responsible for local or systemic infection.
Yet another embodiment can give information about tissue hydration and potentially information about oxyHb and deoxyHb from deeper tissue using NIR wavelengths. These can be used as a stand alone device or as paired with the more standard Visible wavelength MHSI device as shown in figure 2.
Yet another embodiment can derive and present information from changes seen radiating from an area of wounded, ulcerated or otherwise abnormal tissue or from any change in tissue characteristics over a distance . A "gradient map" thus produced can be used to generate a diagnosis, predict the capability of the tissue to heal, define a level for amputation, define the infection vs inflammation, define areas of ischemia, define areas of tissue at risk for ulceration etc.
Another embodiment can involve dividing the region of interest into radial segments, pie like segments or a combination of the two or into squares or other geometric shapes and using these segments to compare and contrast different regions of tissue in the same field of view or as compared to a similar field of view on the contralateral extremity or on another part of the body (such as the forearm, the upper leg, etc.). The radial segments can also be compared to similar locations at different time points to demonstrate change over time in response to different therapeutic interventions, changes in tissue physiology , either local, regional or systemic due either to progression or remission of disease or of the effects of topical or systemic medications or therapies.
Such measurements can be used to evaluate wound healing, tissue regeneration, angiogenesis, vasculogenesis, arteriogenesis, infection, inflammation, microvascular disease or alterations, or other changes in tissue characteristics or physiology associated with the implementation of negative pressure (vacuum suction applied to the wound), hyperbaric therapy, grafting of autologuous, heterograft, xenograft or biological or synthetic skin substituetes, administration of topical agents including antibiotics, cleansers, growth factors, surgical intervention, angioplasty, stenting, atherectomy, laser therapy, vasodilator therapy, offloading, compression, effects of pressure due to orthotic or prosthetic, effects of electromagnetic, acupuncture, massage, infrared, vibration or other therapies. .
Such measurements can be considered as biomarkers representing tissue oxygen delivery and oxygen extraction, tissue oxygenation, tissue perfusion, tissue metabolism or other characteristics correlated with MHSI measurements.
Such measurements can be used in association with the implementatiton of hyperbaric therapy delivered to assist in the healing of ulceration in diabetic or other foot ulceration, or other wounds in other parts of the body. In the case of hyperbaric oxygen therapy, the tissue can be monitored before and at specified intervals during therapy or continuously during therapy to determine when the tissue has been adequately modified (oxygenated) by the therapy or that there has been sufficient change in tissue metabolism as described by the MHSI measurements of oxyHb, deoxy Hb or other measured parameters or whether no benefit is being delivered. MHSI can be used to determine the appropriate duration of HBO therapy during a given session and as to whether sufficient benefit has been delivered from a course of therapy that it can safely be discontinued and that the wound will then be likely to heal with more standard methods.
MHSI. can be used to determine the capability of tissue to heal after debridement and hence the relative safety of pursuing such an approach. Similarly, MHSI can be used to help determine the lowest level of amputation that can be performed with successful healing. Similarly MHSI can be used to determine whether elective surgery to the foot, lower extremity or other body part where evaluation and or quantitation of perfusion, oxygenation, or tissue metabolism would assist in determination of the safety of undertaking such a procedure or the location in which to direct such a procedure. MHSI can be utilized before debridement, amputation or other surgery to make this determination or during debridement, amputatioriv or other surgery to better assess tissue to improve surgical outcomes.
Such measurements can be used for the determination of which patients or which wounds are likely to improve with any of the above mentioned therapies, which patients or wounds or portions of wounds are healing or worsening, when a given therapy is sufficient (this could be during or immedialty after application of a therapy such as hyperbaric therapy or a debridement or a particular cleansing or pharmaceutical regimen or after a longer course of several days of therapy such as a vacuum therapy . MHSI criteria can be used to determine when a tissue will accept a skingraft or benfϊt from an allograft or other skin replacement.
Systemic or regional disease
One embodiment uses a single system that employs light wavelengths derived from the UV, visible, the near infrared, short wave infrared, mid infrared or far infrared portion of the electromagnetic spectrum. Another embodiment uses a system that uses one or more wavelengths from more than one of these wavelength regimes. One such system using wavelengths from more than one of these wavelength groupings is shown in figure two. In other embodiments, a single sensor could be used to collect light from more than one wavelength regime.
A portable hyperspectral imaging apparatus according to an embodiment of the invention is depicted in FIG. 1. Portable apparatus 10 weighs less than 100 pounds, preferably less than 25 pounds, and more preferably less than 10 pounds. Preferably, the portable apparatus may be battery operated, have some other form of portable power source or more preferably, may have a connector adapted to connect to an existing power source.
Portable apparatus 10 comprises an optical acquisition system 36 and a diagnostic processor 38. Optical acquisition system 36 comprises means to acquire broadband data, visible data, ultraviolet data, infrared data, hyperspectral data, or any combination thereof. In a preferred embodiment, optical acquiring means comprises a first-stage imaging optic 40, a spectral separator 42, a second-stage optic 44, and an imaging sensor 46. Alternatively, optical acquiring means may be any acquisition system suited for acquiring broadband data, visible data, ultraviolet data, infrared data, hyperspectral data, or any combination thereof. Preferably, one or more polarizers 41, 43 are included in the acquisition system to compile the light into a plane of polarization before entering the imaging sensor. Preferably, a calibrator is also included in the system.
If the spectral separator 42 does not internally polarizes the light, the first polarizer 43 is placed anywhere in the optical path, preferably in front of the receiving camera 46. The second polarizer 41 is placed in front of illuminating lights 20 such that the incident light polarization is controlled. The incident light is crossed polarized with the light recorded by the camera 46 to reduce specular reflection or polarization at different angles to vary intensity of the reflected light recorded by the camera.
The illumination is provided by the remote light(s) 20, various sources tailored or adapted to the need of the instrument, preferably positioned around the light receiving opening of the system, or otherwise placed to afford optimal performance. The light can be a circular array of focused LED lights that emit light at the particular wavelengths (or ranges) that are used in the processing algorithm, or in the ranges of wavelengths (e.g., visible and/or near-infrared). The circular arrangement of the light sources provides even illumination that reduces shadowing. The light wavelength selectivity reduces effect of the observation on the observing subject. The configuration may also vary depending on the particular needs and operation of the system.
Although the preferred embodiment describes the system as portable, a non-portable system may also be utilized. Preferably, an optical head is mounted to the wall of the examination room, more preferably, an overhead light structure is located in the operating room, or more preferably, the system has a portable table with an observational window overlooking the operating site.
The preferred embodiment may also be used as part of another instrument. For example, as an adjunct to an endoscope.
The first-stage optic receives light collected from a tissue sample through a polarizer and focuses the light onto the surface of the spectral separator. Preferably, the spectral separator is a liquid crystal tunable filter (LCTF). LCTF 42 is a programmable filter that sequentially provides light from selected wavelength bands with small (for example, 7-10 nm) bandwidth from the light collected from the sample. Second-stage optic 44 receives the narrow band of light passing through the spectral separator and focuses the light onto the image sensor 46. The image sensor is preferably, although not necessarily, a two-dimensional array sensor, such as a charge-coupled device array (CCD) or CMOS, which delivers an image signal to the diagnostic processor 38.
Diagnostic processor 38 includes an image acquisition interface 50, that has an input responsive to an output of the image sensor 46 and an output provided to a general-purpose operating module 54. The general-purpose operating module includes routines that perform image processing, and that operates and controls the various parts of the system. The general- purpose operating module also controls the light source(s) (e.g. LED array) allowing for switching on and off during measurement as required by the algorithm. The general-purpose operating module has control output provided to a filter control interface 52, which in turn has an output provided to the spectral separator 42. The general-purpose operating module also interacts with a number of diagnostic protocol modules 56A, 56B, . . . 54N, and has an output provided to a video display. The diagnostic process includes special purpose hardware, general-purpose hardware with special-purpose software, or a combination of the two. The diagnostic processor also includes an input device 58, which is operatively connected to the general-purpose operating module. A storage device 60 and printer 62 also are operatively connected to the general- purpose operating module.
In operation, a portable or semi-portable apparatus is employed within line of site (or with optical access) of the object or area of interest, e.g., diabetic foot with or without an ulcer, or general area of interest. An operator begins by selecting a diagnostic protocol module using the input device. Each diagnostic protocol module is adapted to detect particular tissue characteristics of the target. The diagnostic module could be specific for diabetes, for peripheral vascular disease, for venous stasis disease or for a combination of these disease states. As another example, a screening protocol for feet without ulcers or a potential for healing protocol for feet with ulcers. In an alternative embodiment, the apparatus may contain only one diagnostic module adapted for general medical diagnosis.
Diagnostic processor 38 responds to the operator's input by obtaining a series of transfer functions and an image processing protocol and an image processing protocol from the selected diagnostic protocol module 56. The diagnostic processor provides the filtering transfer functions to the spectral separator 42 via its filter control interface 52 and then instructs the image acquisition interface 50 to acquire and store the resulting filtered image from the image sensor 46. The general-purpose operating module 54 repeats these filtering and acquiring steps one or more times, depending on the number of filter transfer functions stored in the selected diagnostic protocol module. The filtering transfer functions can represent bandpass, multiple bandpass, or other filter characteristics and can include wavelengths in preferably the UV, preferably the visible, preferably the NIR and preferably, the IR electromagnetic spectrum.
In a preferred embodiment, the light source delivering light to the target of interest can be filtered as opposed to the returned light collected by the detector. Thus, a tunable source delivers the information. Alternatively, both a tunable source and a tunable detector may be utilized. Such tuning takes the form of LCTF, acousto-optical tunable filter (AOTF), filter wheels, matched filters, diffraction gratings or other spectral separators. The light source may be a tungsten halogen or xenon lamp, but is preferably a light emitting diode (LED).
The unique cool illumination provided by the LED prevents overheating of skin which may result in poor imaging resolution. Preferably, the LED provides sufficient light while producing minimal or no increase in skin temperature. This lighting system in combination with the polarizer allows adequate illumination while preventing surface glare from internal organs and overheating of skin. In certain embodiments, illumination can arise from more passive sources such as room lights or from sunlight.
Once the image acquisition interface 50 has stored images for all of the image planes specified by the diagnostic protocol chosen by the operator, the image acquisition interface begins processing these image planes based on the image processing protocol from the selected diagnostic protocol module 56N. Processing operations can include general image processing of combined images, such as comparing the relative amplitude of the collected light at different wavelengths, adding amplitudes of the collected light at different wavelengths, or computing other combinations of signals corresponding to the acquired planes. The computed image is displayed on the display 12. Other preferred embodiments include storing the computed image in the storage device 60 or printing the computed image out on printer 62.
In a preferred embodiment, a calibrator is included in the system. Calibrator has an area colored with a pattern of two (or more) colors. To optimize use of the calibrator for this particular application where oxyHb and deoxyHb are important components of the solution, colors are chosen that have a distinct absorption band in the wavelength range similar to oxyHb and deoxyHb - preferably in the range 500-600 nm. The colors are placed into a pattern, preferably, a checker-board pattern, where 1 out of 4 squares has color 1, and 3 out of 4 squares have color2. Thus, approximately 25% of the squares are colorl and 75% of the squares are color2. The system takes a hypercube being slightly out of focus - that provides blurring of colors into each pixel. From the spectra for each pixel, a linear composition of two spectra: one from colorl and another from color2 are observed. The recorded spectra are decomposed in a manner similar to a system that decomposes skin spectra into oxyHb & deoxyHb components. However, in this instance it takes pure color 1 and color2 spectra from library instead of oxyHb & deoxyHb. Valid calibration reports concentrations of 75% for color2 and 25% for colorl. Results are similar to skin analysis, where the output is approximately 90% of oxyHb and 10% of deoxyHb. Other embodiments include but are not limited to, changes to the pattern, the color concentration & intensity, and the number of colors.
In summary, the system simulates the way the biological mixture (oxyHb + deoxyHb) is observed by using "optical" mixture via combination of pattern (with known spatial concentrations) and analog blurring (defocusing - for speed. Defocusing can also be done in the software through the use of computational filters).
If the correct result is obtained, confirmation of the lighting distribution and collection throughput, and the wavelength accuracy of the system given confidence in the spectra (wavelengths and intensity) that are being collected are provided. This provides additional assurance that the data recorded off the patient is acceptable.
In another preferred embodiment, diagnostic protocol modules 56, printer 62, display 12, or any combination thereof, may be omitted from portable device 10. In this embodiment, acquired images are stored in storage device 60 during the medical procedure. At a later time, these images are transferred via a communications link to a second device or computer located at a remote location, for example, hospital medical records, for backup or reviewing at a later time. This second device can have the omitted diagnostic protocol modules, printer, display, or any combination thereof. In another embodiment, the stored images are transferred from portable device 10, located in the clinic, via a communications link to a remote second device in real time.
In a preferred embodiment the system has facility to project real-time hyperspectral data onto the operation field, region of interest, or viewing window positioned above the operating site. Heads Up Display. Also, display completely separately for remote guidance (i.e. on a wall screen for a group of people to review in real time, or post procedure). The projected information has precise one-to-one mapping to the illuminated surface (e.g. wound, operating surface, tissue) and provides the surgeon or podiatrist with necessary information in efficient and non-distractive way. When projected onto an overhang viewing window, the images (real-color and/or pseudo-color) can be zoomed in/out to provide variable magnification. This subsystem consists of the following elements: 1) image projector with field-of-view precisely co-aligned with the field-of-view of the hyperspectral imager, 2) miniature remote control device which allows the surgeon or podiatrist to switch projected image on and off without turning from the site of debridement and change highlight structure and/or translucency on the projected image to improve visibility of the features of interest as well as projected image brightness and intensity, 3) real-time data processing package which constructs projected image based on hyperspectral data and operator/surgeon input, 4) optional viewing window positioned above the operating site that is translucent for real observation or opaque for projecting pseudo-color solution or higher resolution images.
The MHSI system consists of three functional modules - a Spectral Imager (SI), supporting Controller and Power Module (CPM) and Control and Data Acquisition Computer (CDAC). The MHSI also includes a thermometer that remotely measures the temperature at the tissue surface. The Spectral Imager is mounted on suspension arm which neutralizes device weight and allows for easy positioning and focusing of the instrument. The suspension arm is attached to wheeled cart which supports CPM and CDAC as well. This configuration is very mobile and permits wide range of device spatial and directional motions.
FIG 2 shows the preferred system specifications along with a diagram of our focusing methodology and the optical design of the Spectral Imager. In this embodiment, a liquid crystal tunable filters (LCTF 's) was used as the wavelength selector and are coupled to complementary metal oxide semiconductor (CMOS) imaging sensors. Fitted with macro lenses and the positional light focusing system described below, the system has a preferred working focal length of roughly 1 to 2 feet.
A major issue in the collection of hyperspectral imaging data is the position and focusing of the instrument. While our Spectral Imaging Module is positioned on a ball joint that allows free rotation and virtually any angle of incidence to the patient, it is imperative that there be a system in place for targeting the image to a particular spot on the tissue and ensuring that the instrument will be at the proper distance from the tissue to achieve optimal focus. Positioning and focusing with our system are facilitated by two mirrored collimated light beams or lasers that cross precisely at the instrument focal plane (FIG 2), and so bringing the spectral imaging module into position where the two light spots overlap on the tissue to insure optimal focus. To achieve precisely calibrated images, the system may use a specially designed calibration pad placed at the focal plane of the system and measured prior to each patient measurement. The calibration pad includes a diffusely reflective surface to quantify the intensity of the illumination at each wavelength and color bars to validate wavelength accuracy of the system. Calibration data measured at a preset time such as during maintenance calibrations can be stored and compared to with each use to decide whether the system is within specifications and should proceed to patient measurements.
To achieve precise co-registration between hyperspectral image and operating surface, the system may use a fiducial label or target placed in the field of view which the image registration module can perform a self-alignment procedure before or during the operation as necessary.
Devices of the present invention allow for the creation and unique identification of patterns in data that highlight the information of interest. The data sets in this case may be discrete images, each tightly bounded in spectra that can then be analyzed. This is analogous to looking at a scene through various colored lenses, each filtering out all but a particular color, and then a recombining of these images into something new. Such techniques as false color analysis (assigning new colors to an image that don't represent the true color but are an artifact designed to improve the image analysis by a human) are also applicable. Optionally, optics can be modified to provide a zoom function, or to transition from a micro environment to a macro environment and a macro environment to a micro environment. Further, commercially available features can be added to provide real-time or near real-time functioning. Data analysis can be enhanced by triangulation with two or more optical acquisition systems. Polarization may be used as desired to enhance signatures for various targets.
In addition to having the ability to gather data, the present invention also encompasses the ability to combine the data in various manners including vision fusion, summation, subtraction and other, more complex processes whereby certain unique signatures for information of interest can be defined so that background data and imagery can be removed, thereby highlighting features or information of interest. This can also be combined with automated ways of noting or highlighting items, areas or information of interest in the display of the information. The hyperspectrally resolved image in the present invention is comprised of a plurality of spectral bands. Each spectral band is adjacent to another forming a continuous set. Preferably, each spectral band having a bandwidth of less than 50 nm, more preferably less than 30 nm, more preferably less than 20 nm, more preferably, from about 20 - 40 nm, more preferably, from about 20 - 30 nm, more preferably, from about 10 - 20 nm, more preferably from about 10 - 15 nm, and more preferably from about 5 - 12 nm.
It is clear to one skilled in the art that there are many uses for a medical hyperspectral imager (MHSI) according to the invention. The MHSI offers the advantages of performing the functions for such uses faster, more economically, and with less equipment and infrastructure/logistics tail than other conventional techniques. Many similar examples can be ascertained by one of ordinary skill in the art from this disclosure for circumstances where medical personal relies on their visual analysis of the biological system. The MHSI acts like "magic glasses" to help human to see inside and beyond.
Algorithm description
The embodiment of diabetes algorithm involves the following steps:
1. Preprocess the MHSI data. Preferably, by removing background radiation by subtracting the calibrated background radiation from each newly acquired image while accounting for uneven light distribution by dividing each image by the reflectance calibrator image and registering images across a hyperspectral cube.
2. Build a color-photo-quality visual image. Preferably, by concatenating three planes from the hyperspectral cube at the wavelengths that approximately correspond to red (preferably in the range of about 580-800 nm, more preferably in the range of about 600- 700 nm, more preferably in the range of about 625-675 nm and more preferably at about 650 nm), green (preferably in the range of about 480-580 nm, more preferably in the range of about 500-550 nm, more preferably in the range of about 505-515 nm, and more preferably at about 510 nm), and blue (preferably in the range of about 350-490 nm, more preferably in the range of about 400-480 nm, more preferably in the range of about 450- 475 nm, and more preferably at about 470 nm) color along the third dimension to be scaled for RGB image. 3. Define a region of interest (ROI), preferably, where the solution is to be calculated unless the entire field of view to be analyzed.
4. Convert all hyperspectral image intensities into units of optical density. Preferably, by taking the negative logarithm of the decimal base. FIG 2 shows examples of spectra taking from single pixels at different tissue sites within an image. Tissue sites include connective tissues, oxygenated tissues, muscle, tumor, and blood.
5. Decompose the spectra for each pixel (or ROI averaged across several pixels). Preferably, decompose into several independent components, more preferably, two of which are oxyhemoglobin and deoxyhemoglobin.
6. Determine three planes for pseudo-color image. Preferably, define the color hue plane as apparent concentration of oxygenated Hb, or deoxygenated Hb, or their mathematical combination, e.g. total Hb, oxygen saturation, etc. Preferably, define the color saturation plane as apparent concentration of oxygenated Hb, or deoxygenated Hb, or their mathematical combination, e.g. total Hb, oxygen saturation, etc. Preferably, define the color intensity (value) plane as reflectance in blue-green-orange region (preferably in the range of light at about 450-580 nm).
7. Adjust the color resolution of the pseudo-color image according to quality of apparent concentration of oxygenated Hb, or deoxygenated Hb, or their mathematical combination, e.g. total Hb, oxygen saturation, etc. Preferably, reduce resolution of hue and saturation color planes by binning the image (e.g. by 2, 3, 4, etc. pixels), or/and by resizing the image, or/and by smoothing the image through filtering higher frequency components out. Interpolate the smoothed color planes on the grid of higher resolution intensity (value) plane.
8. Convert hue-saturation- value/intensity (HSV/I) image to red-green-blue (RGB) image.
9. Remove outliers in the resulting image, defining an outlier as color intensity deviating from a typical range beyond certain number of standard deviations, preferably three. Stretch the resulting image to fill entire color intensity range, e.g. between 0 and 1 for a double precision image. 10. Display ROI in pseudo-colors, preferably, in combination with the color photo image of the subject, or preferably, in addition to the color photo image of the subject, or more preferably, by projecting the pseudo-color image onto the observed surface. Additional information can be conveyed through images portraying the individual coefficients from oxyHb, deoxyHb, slope and offset coefficients, or any linear or nonlinear combination such as the oxyhemoglobin to deoxyhemoglobin ratio.
11. Characterize the metabolic state of the tissue of interest (e.g. risk for ulceration, potential to heal). Preferably, by using the saturation and/or intensity of the assigned color and provide a qualitative color scale bar.
As is clear to a person of ordinary skill in the art, one or more of the above steps in the algorithm can be performed in a different order or eliminated entirely and still produce adequate and desired results. Preferably, the set of instructions includes only the steps of preprocessing the hyperspectral information, building a visual image, using the entire field of view, converting all hyperspectral image intensities into units of optical density by taking a negative logarithm of each decimal base, and characterizing a metabolic state of the tissue of interest. More preferably, the set of instructions comprises preprocessing the hyperspectral information, defining a region of interest of the tissue, and characterizing a state of the tissue of interest.
Another preferred embodiment entails reducing the hyperspectral data in the spectral dimension into a small set of physiologic parameters involves resolving the spectral images into several linearly independent images (e.g. oxyhemoglobin, deoxyhemoglobin, an offset coefficient encompassing multiple scattering (MS) properties and a slope coefficient) in the visible regime. Another embodiment determines four images (e.g. oxyhemoglobin, deoxyhemoglobin, offset/ scattering coefficient, and water absorption) in the near infrared region of the spectrum. As an example for the visible region of the spectrum, linear regression fit coefficients cl, c2, c3 and c4 will be calculated for reference oxy-Hb, deoxy-Hb, and MS spectra, respectively, for each spectrum (Sij) in an image cube:
Sy = Cx OxyHb + C2 DeoxyHb + C3 Offset + c^Slope
Individual images of the oxyhemoglobin and deoxyhemoglobin components, the slope and offset or any combination, linear or nonlinear, of these terms, for example the oxy- to deoxyhemoglobin ratio, can be presented in addition to producing the pseudo-colored image to the user.
In order to present the MHSI effectively, a display method was developed that has the potential to convey a 2-dimensional index, and convey the values for both the oxy and deoxyhemoglobin coefficients independently. The method of displaying our index uses a color scale, in one iteration this ranges from purple values (high) to brown values (low) to indicate the concentration of oxyhemoglobin in the tissue, and a brightness scale, ranging from very bright (high) to faded (low) associated with the tissue concentration of deoxyhemoglobin. FIG 3 summarizes the display of the MHSI, showing a schematic diagram explaining the scenarios of low and high oxy and deoxy hemoglobin coefficients as well as a color scale that indicates a color plat that shows the vertical color scale and the horizontal brightness scale. By measuring an MHSI where the oxyhemoglobin component is high and the deoxyhemoglobin component is low (upper left hand corner of FIG 3), it could be concluded that that particular area of tissue has adequate perfusion and oxygenation, and is able to satisfy its metabolic needs with the oxygen that is being delivered. That this tissue has the lowest level of risk for ulceration and the highest probability of healing. If tissue demonstrates a low oxyhemoglobin level in addition to a low deoxyhemoglobin level (lower left corner of FIG 3), this would imply that the tissue was receiving low total volume of blood. If tissue demonstrates a low oxyhemoglobin level in addition to a high deoxyhemoglobin level (lower right corner of FIG 3) this would imply that the tissue has metabolic requirements exceeding available oxygen delivery. In both of these regions there is expected to be a higher risk of ulceration or difficulties with wound healing. If the tissue has a high oxyhemoglobin coefficient and also has a high deoxyhemoglobin content, (lower left corner of FIG 3) this tissue was receiving a larger total volume of blood, and that the oxygen extracted from the blood stream was adequate to support tissue metabolism. This could be indicative of inflammation. Our technique will uniquely permit discrimination between each of these disparate physiologic conditions. For example, if the value is faded purple (upper left hand quadrant) the tissue has very high oxygenation, as discussed above, and is very likely to heal. The color map (right) gives an indication of how the MHSI would be represented in an image format.
Describe here is hyperspectral imaging for use in the peripheral vascular and diabetes clinic, designed both to be mobile for ease of use and to facilitate the most accurate data collection possible for this project. This system provides fast and precise measurement of reflectance spectra, and is characterized by high spatial and spectral resolution, the ability to process spectral data in real time. It has been equipped with a turn-key software interface for the user. Proprietary image registration software insures image stability when measuring spectra of animated objects. The system does not rely on external illumination, rather it contains very efficient internal visible (and NIR in certain versions) light sources, which allow to achieve high signal to noise ratios in measured data without putting noticeable heat load on a biological subject (variations in skin temperature during acquisition are on the order of 0.1 C).
All MHSI data were corrected for background and uneven illumination, and normalized by the integration time. The data were ratioed to the reflectance of a calibration standard, and negative decimal logarithm was taken to obtain the absorption data. Images at all wavelengths were co-registered using proprietary software developed by HYPERMED to ensure that each pixel represents the same point on the skin throughout all wavelengths. The spectra were then deconvolved into four linearly-independent spectral components with coefficients representing the amount of hemoglobin (both oxyHb and deoxyHb) in the observed skin. Typically two numbers are presented x/y wherein x represents oxyHb and y is deoxyHb. The values for x and y can be taken from a single pixel or from a ROI defined by the user. In addition the hemoglobin oxygen saturation (SHSIO2), x/(x+y), can be presented for a pixel or a ROI. FIG 4 shows examples of tissue with low and high oxyHb and deoxyHb values corresponding to tissue at risk of ulceration and a wound that is likely to heal, respectively. Another way of presenting linearly independent variables is through their sum and difference: (x+y) and (x-y). The first would be THb, and the second would "hint" on oxygen extraction - which indicates the kind of Hb that is predominant at the site.
Data in the following table represent typical oxyHb, deoxyHb, and SHSIO2 values for two body positions, forearm and foot, and for various stages of diabetes: nondiabetics, diabetics without peripheral neuropathy, and diabetics with peripheral neuropathy. In general, the value for oxyHb and SHSIO2 are lower in the feet of diabetic subjects with neuropathy compared to the other two groups, a group at high risk for developing foot ulcers. In addition, the values for oxyHb, deoxyHb, and SHSIO2 depend on body location, that once calibrated can be accounted for by the diagnostic module. MHSI oximetry values at baseline (prior to iontophoresis of acetylcholine)
Site Group (N) Oxy Deoxy SHSlO2 (%)
Forearm Control (21) 29 ±7* 41 ± 16 42 ± 17**
Diabetic Non-
20 ± 5 44 + 10 32 ±8** Neuropathic (36)
Diabetic
19 ± 7 49 + 10 28 ± 8** Neuropathic (51)
Dorsum of foot Control (21) 25 ± 13 44 ± 18 38 ± 22
Diabetic Non-
24 ± 9 41 ± 11 37 ± 12 Neuropathic (36)
Diabetic
19 ± 9*** 45 ± 13 30 ± 12**** Neuropathic (51)
* p < 0.0001 compared to diabetics with and without neuropathy ** p < 0.0001 for all three groups
*** p < 0.025 when compared to control and nonneoropathic **** p < 0.027 when compared to control and nonneoropathic
In summarizing these data, MHSI provides relevant physiological information at the systemic, regional and local levels. Forearm data measures systemic microvasculature changes since the forearm is not affected by macrovasculature or somatic neuropathy as found in the lower extremities. Dorsal foot measurements are indicative of microvascular and macro vascular effects including atherosclerotic changes occurring in large vessels exacerbated by diabetes. MHSI data from the right and left lower extremity can be compared to help differentiate the stages of the damage. Finally, MHSI can be used to find local information that can be associated to the risk of developing a foot ulcer of the progression of disease by examining the area around an ulcer.
MHSI can not only be used for determining risk of foot ulceration, but also for determining systemic progression of diabetic microvascular disease. In one embodiment this is determined by mean oxyHb, deoxyHb and/or other values for a region of interest. In another embodiment this is determined by heterogeneity of oxyHb, deoxyHb and/or other values for a region of interest. In another embodiment this is determined by the patterning of oxyHb, deoxyHb and/or other values for a region of interest. In other embodiments this is determined by changes over a given time period within a measurement session (between 1 picosecond and one hour, preferably between 100 microseconds and 10 minutes and more preferably between 100 microseconds and 15 seconds) and in the mean values or patterns of oxyHb, deoxyHb and/or other values for a region of interest. These measurements can be used to determine a diabetes progression index (DPI). Alternatively, a DPI can be calculated by comparing a MHSI value or set of values from a single point in time with another point in time (preferably 1 month to two years, more preferably 2 months to one year and most preferably 3-6 months)
Using the active stimulus of acetylcholine as a vasodilator, and the known effects of this on LD measurements data was derived in which changes in MHSI could be observed under known alterations in physiology and compare these with baseline images and with LD data (FIG 5). Using data collected from diabetic patients, as well as previous data from human shock studies and iontophoresis studies, an algorithm was derived that clearly discriminates regions vasodilated by iontophoresis and also discriminates ulcer from non ulcer with a proprietary formula that includes terms for oxyHb and deoxyHb. This was further developed as a Hyperspectral Microvascular Index (HMI), which is a metric of tissue physiology and have explored the use of the MHSI in evaluating tissue of the foot. In circumstances when an ulcer has been present, tissue was examined within the ulcer, directly adjacent to the ulcer, surrounding the ulcer and at various other regions of the foot.
With the aid of MHSI, a quantitative metric is demonstrated with superb separation between ulcerated or wounded and non-ulcerated or wounded tissue. Areas of different tissue metabolism can be seen with 60 micron (20-120) spatial resolution. Regions with an increased MHSI associated with the margins of the ulcer can be seen which correlate to inflammation (and/or infection). Areas of decreased MHSI can be seen in other areas which from previous work in ischemia is considered to be tissue at risk for non-healing, ulcer extension, or primary ulceration. These data validate the capability of our measurement system to have the resolution and appropriate range to quantitatively assess different areas of tissue metabolism on both dorsal and plantar foot surfaces as well as skin on other body areas or other tissues visible through endoscopic techniques or at the time of open surgery of the foot, leg, arm or any other body part including internal organs at laparoscopy or the retina at retinoscopy . The invention provides the capability to perform this quantitative assessment on tissue that demonstrates no visible differences on clinical examination to the skilled examiner.
MSHI images have the ability to differentiate between regions of tissue associated with a present foot ulcer on the basis of biomarkers such as the oxyHb and deoxyHb coefficients. FIG 6 shows an ulcer on the sole of the foot of a type 1 diabetic patient (ulcer 1). From the visible image on the left, little distinguishes one area of the ulcer from another. However when looking at the image with the MHSI, there is obvious discriminatory power between the state of tissue seen in the purple oval, which is likely to heal, and that surrounded by the black oval, which is tissue at risk for further ulceration. It is important to note that the skin on the sole of this patient's feet is highly calloused, with a thick stratum corneum, but one is still able to differentiate tissue based on its spectral signatures. Given that the sole of the foot is often the site of the thickest stratum corneum on the body, the device works on all naturally or surgically exposed tissue or tissue otherwise visualized with laparoscopy, endoscopy, retinoscopy or other visualization techniques. Ulcer 2 was located on the dorsal surface of the foot, on the patient's big toe (FIG 6). These images further show the ability to differentiate between tissue at risk and tissue likely to heal. Additionally, tissue surrounding a fungal infection on the patient's middle toe (bottom right-hand corner of the image) has an MHSI that can demonstrate inflamed or infected tissue.
In addition to the differentiation of local tissue, tissue can be examined for gross features indicative of global risks of complications, such as poor perfusion or the inability of the microcirculation to react and compensate in tissue. In another embodiment, iontophoretic application of the vasodilator acetylcholine (ACH) or nitroprusside was used to stimulate the vasodilation of the microvasculature on the dorsal surface of the foot and on the forearm of the patients and measured the reaction with MHSI (FIG 7).
The potential of hyperspectral imaging in diagnosing global microcirculatory insufficiencies and impacting on other complications of diabetes associated with the microvasculature besides foot ulcers. In FIG 7, hyperspectral measurements from the feet of four patients, with the first two columns of images showing the MHSI of the soles of both feet, and the second two columns showing images of the dorsal surfaces of both feet after the application of ACH via iontophoresis. In the first three patients, an MHSI is seen that is much healthier than that of the fourth patient. Consequently, the fourth patient had a foot ulcer at the time of this study and has a previous history of ulceration. While the contrast between the data from the soles in these patients is striking, there is complementary information in the data from the microvascular response shown in the two columns on the right. Note that the first three patients all have MHSI scores that contain purple information in response to vasodilation, while the fourth patient shows what would be considered an MHSI that was indicative of tissue that was at risk. Microcirculatory changes associated with the progression of diabetes can also be modified by different treatment and therapeutic regimens and with the overlay of other systemic diseases (such as congestive heart failure or hypertension) or treatments or therapies for systemic diseases.
To analyze the ulcer data further, ulcer images were divided into 25 concentric circles lmm apart and 8 pie segments forming 200 sectors per ulcer (FIG 8). A radial profile analysis was undertaken where the ulcer center was defined at the first visit, and registered images from subsequent visits to this. OxyHb, deoxyHb, total-hemoglobin and C^Sat were calculated for each sector.
Each radial pie segment was evaluated for signs of healing, nonhealing or progression in subsequent visits. MHSI measurements and clinical healing results were compared. MHSI algorithms were developed to identify changes associated with ulcer healing, nonhealing and progression. A primary endpoint evaluated the specific sectors of tissue around an ulcer that would heal, not heal or progress. The group estimates for oxyHb, deoxyHb, and O2Sat are given in the following table using a linear mixed effects regression model. Significant differences were seen for healing for the oxyHb and deoxyHb values. Patients who did not heal also demonstrated increased heterogeneity in distant foot and in arm measurements. For the 21 ulcers studied, the algorithm predicted 6 of 7 ulcers that did not heal and 10 of 14 ulcers that healed. Conclusion: MHSI identifies microvascular abnormalities in the diabetic foot and provides early information assist in managing foot ulceration and predict outcomes in patients with diabetes.
Figure imgf000036_0001
Figure imgf000037_0001
MHSI is used to monitor angiogenesis during wound healing. An example of wound healing in a diabetic rabbit wound model shows that during the healing process, images of oxyHb and deoxyHb show patterns that change in shape, area, and amplitude with time. Similar patterns were noted in experimental models of shock, but the changes observed for shock occurred on a shorter time scale; minutes rather than several days as the wound heals. The rabbit's ears were observed at days 1, 2, 5, and 10. MHSI is ideally suited for characterization of the local heterogeneity in oxyHb and deoxyHb and their spatial changes with time. For example, the zone of hyperemia surrounding a wound as measured by the oxyHb coefficient decreases with time in wounds that heal (FIG 9). The color image (a), reconstructed from MHSI data, shows a part of the observed area 50-by-40 mm, recorded at the baseline on day 1. The black rings denote location of a future wound. The pseudocolor image (b), obtained as a result of hyperspectral processing, shows distribution of the oxyHb and deoxyHb in the underlying tissue at the same time. The color hue represents apparent oxyHb concentrations, whereas color saturation (from fade to bright) represents apparent deoxyHb concentrations. Both, oxyHb and deoxyHb vary predominantly between 40 and 90 MHSI units (color bar to the right). The remaining images to the right show change in a region of interest 17-by-17 mm (black box in (a) and (b)) over 10 days. At day 2, the oxy concentrations increased significantly in the area as far as 10 mm away from the wound border. By day 5, the increase in oxygenation became more local (purple area, shrunken to about 5 mm) and new microvasculature formed to feed the area in need (red fork- like vessels in the right top corners appearing in days 5 and 10 images). By the 10th day, the area of increased oxyHb has not changed much, but the peak in oxy amplitude decreased, suggesting a period of steady healing. As depicted in Fig. 10, 50-micron resolution images of a rabbit's ear were taken with MHSI over a ten day period. In Fig. 10(a), the color image was reconstructed from MHSI data, showing a party of the observed area 50-by-40 mm, recorded at the baseline on day 1. The pseudo-image (b) was obtained as a result of hyperspectral processing, showing a distribution of the oxygenated (oxy) and deoxygenated (deoxy) hemoglobin in the underlying tissue at the same time.
Other embodiments and uses of the invention will be apparent to those skilled in the art from consideration of the specification and practice of the invention disclosed herein. All references cited herein, including all publications, U.S. and foreign patents, and patent and provisional applications, and all publications and documents cited herein for any reason, are specifically and entirely incorporated by reference. It is intended that the specification and examples be considered exemplary only. Literature Cited
1. International Diabetes Federation. Information on Global Diabetes Prevalence. In: http ://www. eatlas . idf . or g/Prevalence/; 2004.
2. American Diabetes Association, http ://www. diabetes . or g/tvpe- 1 - diabetes/complications.jsp. 2005.
3. Brearley S, Shearman CP, Simms MH. Peripheral pulse palpation: an unreliable physical sign. Ann R Coll Surg Engl 1992;74(3): 169-71.
4. Frykberg RG, Lavery LA, Pham H, Harvey C, Harkless L, Veves A. Role of neuropathy and high foot pressures in diabetic foot ulceration. Diabetes Care 1998;21(10): 1714-9.
5. Sumpio BE. Foot ulcers. N Engl J Med 2000;343(l l):787-93.
6. Young MJ, Breddy JL, Veves A, Boulton AJ. The prediction of diabetic neuropathic foot ulceration using vibration perception thresholds. A prospective study. Diabetes Care 1994;17(6):557-60.
7. Cavanagh P, Young M, Adams J, al. e. Correlates of structure and function in neuropathic diabetic feet. Diabetologia 1991;34(Suppl 2):A39 (abstract).
8. Lavery LA, Armstrong DG, Vela SA, Quebedeaux TL, Fleischli JG. Practical criteria for screening patients at high risk for diabetic foot ulceration. Arch Intern Med 1998;158(2): 157-62.
9. McMillan DE. Development of vascular complications in diabetes. Vase Med 1997;2(2): 132-42.
10. Novo S. Classification, epidemiology, risk factors, and natural history of peripheral arterial disease. Diabetes Obes Metab 2002;4 Suppl 2:Sl-6. 11. Hittel N, Donnelly R. Treating peripheral arterial disease in patients with diabetes. Diabetes Obes Metab 2002;4 Suppl 2:S26-31.
12. Pecoraro RE, Reiber GE, Burgess EM. Pathways to diabetic limb amputation. Basis for prevention. Diabetes Care 1990;13(5):513-21.
13. Reiber GE, Boyko EJ, Smith DC. Lower extremity foot ulcers and amputations in diabetes. In: Harris, Cowie, Stern, Boyko EJ, Reiber GE, Bennet, editors. Diabetes in America. 2nd ed. Washington, DC: US Government Printing Office; 1995. p. 402-428.
14. Palumbo PJ, Melton LJ. Peripheral vascular disease and diabetes. In: Harris, Hamman, editors. Diabetes in America. 1st ed. Washington, DC: US Government Printing Office; 1985.
15. Frykberg RG, Armstrong DG, Giurinϊ J, Edwards A, Kravette M, Kravitz S, et al. Diabetic foot disorders: a clinical practice guideline. American College of Foot and Ankle Surgeons. J Foot Ankle Surg 2000;39(5 Suρpl):Sl-60.
16. Ramsey SD, Newton K, Blough D, McCulloch DK, Sandhu N, Reiber GE, et al. Incidence, outcomes, and cost of foot ulcers in patients with diabetes. Diabetes Care 1999;22(3):382-7.
17. Harrington C, Zagari MJ, Corea J, Klitenic J. A cost analysis of diabetic lower-extremity ulcers. Diabetes Care 2000;23(9): 1333-8.
18. Caputo GM, Cavanagh PR, Ulbrecht JS, Gibbons GW, Karchmer AW. Assessment and management of foot disease in patients with diabetes. N Engl J Med 1994;331(13):854- 60.
19. Lavery L, Gazewood JD. Assessing the feet of patients with diabetes. J Fam Pract 2000;49(l l Suppl):S9-16.
20. Frykberg RG. Diabetic foot ulcers: pathogenesis and management. Am Fam Physician 2002;66(9): 1655-62.
21. Sykes MT, Godsey JB. Vascular evaluation of the problem diabetic foot. Clin Podiatr Med Surg 1998; 15(l):49-83.
22. Treado PJ, Morris MD. Infrared and Raman spectroscopic imaging. Appl Spectrosc Rev 1994;29:l-38.
23. Riaza A, Strobl P, Muller A, Beisl U, Hausold A. Spectral mapping of rock weathering degrees on granite using hyperspectral DAIS 7915 Spectrometer Data. Internl J Applied Earth Observation and Geoinformation Special issue: Applications of imaging spectroscopy 2001;3-4:345-354.
24. Thenkabail PS, Smith RB, De Pauw E. Hyperspectral vegetation indices and their relationships with agricultural crop characteristics. Remote Sens Environ 2000;71(REMOTE SENS ENVIRON): 158-182.
25. Colarusso P, Kidder LH, Levin IW, et al. Infrared spectroscopic imaging: from planetary to cellular systems. Appl Spectrosc 1998;52:106A-120A.
26. HyperMed, Inc., Issued US patents 6,937,885; 6,810,279, 6,741,884; 6,640,130; and 6,640,132; and others pending. 27. Freeman JE, Shi Y, Panasyuk SV, Rogers AE. Medical hyperspectral imaging (MHSI) of l,2-dimethylbenz(a)-anthracene (DMBA)-induced breast tumors in rats. Poster #1001. In: 27th Annual San Antonio Breast Cancer Symposium; 2004; San Antonio, Texas: Breast Cancer Research and Treatment; 2004. p. S51.
28. Carlson KA, Jahr JS. A historical overview and update on pulse oximetry. Anesthesiol Rev 1993;20: 173-181.
29. Zuzak KJ, Schaeberle MD, Gladwin MT, et al. Noninvasive determination of spatially resolved and time-resolved tissue perfusion in humans during nitric oxide inhibition and inhalation by use of a visible-reflectance hyperspectral imaging technique. Circulation 2001;104:2905-2910.
30. Afromowitz MA, Callis JB, Heimbach DM, DeSoto LA, Norton MK. Multispectral imaging of burn wounds: a new clinical instrument for evaluating burn depth. IEEE Trans Biomed Eng 1988;35(10):842-50.
31. Sowa MG, Leonardi L, Payette JR, Fish JS, Mantsch HH. Near infrared spectroscopic assessment of hemodynamic changes in the early post-burn period. Burns 2001; 27:241- 9. Burns 2001;27:241-249.
32. Dinh T, Panasyuk S, Freeman J, Panasyuk AA, Lew R, Brand D, et al. Medical hyperspectral imaging (MHSI) evaluation of microcirculatory changes to predict clinical outcomes: Application to diabetic foot ulcers. Society of Vascular Medicine and Biology 17th Annual Scientific Session 2006:(Abstract submitted).
33. Dinh T, Panasyuk SV, Jiang C, Freeman J, Panasyuk AA, Nerney M, et al. The use of medical hyperspectral imaging (MHSI) to evaluate microcirculatory changes in diabetic foot ulcers and predict clinical outcomes. American Diabetes Association 66th Scientific Session 2006:(abstract accepted).
34. Dinh T, Panasyuk SV, Tracey BH, Khaodhiar L, Lyons TE, Rosenblum BL, et al. The use of medical hyperspectral imaging (MHSI) to identify patients at risk for developing diabetic foot ulcers. Diabetes 2005;54(Sl):A270.
35. Dinh T, Panasyuk SV, Tracey BH, Khaodhiar L, Lyons TE, Rosenblum BL, et al. The use of medical hyperspectral imaging (MHSI) to identify patients at risk for developing diabetic foot ulcers. American Diabetes Association 65th annual session 2005:Poster #1106-P, June.
36. Greenman RI, Panasyuk S, Wang X, Lyons TE, Dinh T, Longorio L, et al. Early changes in the skin microcirculation and muscle metabolism of the diabetic foot. Lancet 2005; 366: 1711-1718. Lancet 2005;366: 1711-1718.
37. Mansfield JR, Lew RA, Lewis EN, Freeman JE, Kellicut D, Sidawy A. Hyperspectral imaging of diabetic foot conditions. In: Eastern Analytical Society Meeting; 2003; Somerset, NJ; 2003.
38. Freeman JE, Panasyuk SV, Hopmeier MJ, Lew RA, Batchinski AI, Cancio LC. The evaluation of new methods of hyperspectral image analysis for the diagnosis of hemorrhagic shock, Poster #17. In: American Association for the Surgery of Trauma; 2005; 2005. 39. Gillies R, Freeman JE, Cancio LC, Brand D, Hopmeier M, Mansfield JR. Systemic effects of shock and resuscitation monitored by visible hyperspectral imaging. Diabetes Technol Therapeut 2003;5(5):847-855.
40. Panasyuk SV, Freeman JE, Cooke WI, Hopmeier MJ, Convertino VA. Initial demonstration in human subjects of medical hyperspectral imaging (MHSI) as a novel stand-off non-invasive method for diagnosing and measuring hemodynamic compromise, Poster #18. In: American Association for the Surgery of Trauma; 2005; 2005.
41. Payette JR, Sowa MG, Germscheid SL, Stranc MF, Abdulrauf B, Mantsch HH. Noninvasive diagnostics: predicting flap viability with near-IR spectroscopy and imaging. Am Clinical Laboratory 1999;18:4-6.
42. Armstrong DG, Lavery LA. Predicting neuropathic ulceration with infrared dermal thermometry. J Am Podiatr Med Assoc 1997;87(7):336-7.
43. Beckert S, Witte MB, Konigsrainer A, Coerper S. The Impact of the Micro-Lightguide O2C for the Quantification of Tissue Ischemia in Diabetic Foot Ulcers. Diabetes Care 2004;27(12):2863-2867.
44. Rajbhandari SM, Harris ND, Tesfaye S, Ward JD. Early identification of diabetic foot ulcers that may require intervention using the micro lightguide spectrophotometer. Diabetes Care 1999;22(8): 1292-5.
45. Khan F, Newton DJ. Laser Doppler imaging in the investigation of lower limb wounds. Int J Low Extrem Wounds 2003;2(2):74-86.
46. Sheffield PJ, Fife CE, Smith APS. Laser Doppler Flowmetry. In: Wound Care Practice. Flagstaff, AZ: Best Publishing Company; 2004.
47. Wardlaw JM, Dennis MS, Warlow CP, Sandercock PA. Imaging appearance of the symptomatic perforating artery in patients with lacunar infarction: occlusion or other vascular pathology? Ann Neurol 2001;50(2):208-15.
48. Zimny S, Dessel F, Ehren M, Pfohl M, Schatz H. Early detection of microcirculatory impairment in diabetic patients with foot at risk. Diabetes Care 2001;24(10):1810-4.

Claims

What is Claimed:
1. A medical instrument comprising: a first stage optic responsive to illumination of tissue; a spectral separator; one or more polarizers; an imaging sensor; a diagnostic processor; a filter control interface; a general purpose operating module to assess the state of tissue in diabetic subjects following a set of instructions; and a calibrator.
2. The medical instrument of claim 1 further comprising a second stage optic responsive to illumination of tissue.
3. The medical instrument of claim 1 wherein the set of instructions comprises: preprocessing hyperspectral information; building a visual image; defining a region of interest in tissue; converting the visual image into units of optical density by taking a negative logarithm of each decimal base; decomposing a spectra for each pixel into several independent components; determining three planes for an RGB pseudo-color image; determining a sharpness factor plane; converting the RGB pseudo-color image to a hue-saturation- value/intensity image having a plane; adjusting the hue- saturation- value/intensity image plane with the sharpness factor plane; converting the hue-saturation-value/intensity image back to the RGB pseudo-color image; removing outliers beyond a standard deviation and stretching image between 0 and 1; displaying the region of interest in pseudo-colors; and characterizing a metabolic state of the tissue of interest.
4. The medical instrument of claim 3 wherein the region of interest is one of a pixel, a specified region or an entire field of view.
5. The medical instrument of claim 3 wherein determining three planes for an RGB pseudocolor image comprises one or more characteristic features of the spectra.
6. The medical instrument of claim 3 wherein determining a sharpness factor plane comprises a combination of the images at different wavelengths.
7. The medical instrument of claim 3 wherein removing outliers beyond a standard deviation comprises three standard deviations.
8. The medical instrument of claim 3 wherein displaying the region of interest in pseudocolors comprises one of performing one in combination with a color photoimage of a subject, in addition to a color photo image of a subject, and projecting onto the tissue of interest.
9. The medical instrument of claim 3 wherein defining the color intensity plane as apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation.
10. The medical instrument of claim 3 wherein defining the color intensity plane as reflectance in blue-green-orange region.
11. The medical instrument of claim 3 wherein adjusting the hue saturation comprises adjusting a color resolution of the pseudo-color image according to quality of apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation.
12. The medical instrument of claim 11 wherein adjusting the hue saturation further comprises one or a combination of reducing resolution of hue and saturation color planes by binning the image, resizing the image, and smoothing the image through filtering higher frequency components out.
13. The medical instrument of claim 12 further interpolating the smoothed color planes on a grid of higher resolution intensity plane.
14. A method for assessing the state of tissue of a diabetic subject comprising: preprocessing hyperspectral information; building a visual image; defining a region of interest in tissue; converting the visual image into units of optical density by taking a negative logarithm of each decimal base; decomposing a spectra for each pixel into several independent components; determining three planes for an RGB pseudo-color image; determining a sharpness factor plane; converting the RGB pseudo-color image to a hue-saturation-value/intensity image having a plane; adjusting the hue- saturation-value/intensity image plane with the sharpness factor plane; converting the hue-saturation-value/intensity image back to the RGB pseudo-color image; removing outliers beyond a standard deviation and stretching image between 0 and 1; displaying the region of interest in pseudo-colors; and characterizing a metabolic state of the tissue of interest.
15. The method of claim 14 wherein the region of interest is one of a pixel, a specified region or an entire field of view.
16. The method of claim 14 wherein determining three planes for an RGB pseudo-color image comprises one or more characteristic features of the spectra.
17. The method of claim 14 wherein determining a sharpness factor plane comprises a combination of the images at different wavelengths.
18. The method of claim 14 wherein removing outliers beyond a standard deviation comprises three standard deviations.
19. The method of claim 14 wherein displaying the region of interest in pseudo-colors comprises one of performing one in combination with a color photoimage of a subject, in addition to a color photo image of a subject, and projecting onto the tissue of interest.
20. The method of claim 14 wherein defining the color intensity plane as apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation.
21. The method of claim 14 wherein defining the color intensity plane as reflectance in blue- green-orange region.
22. The method of claim 14 wherein adjusting the hue saturation comprises adjusting a color resolution of the pseudo-color image according to quality of apparent concentration of one or a mathematical combination of oxygenated Hb, deoxygenated Hb, and total Hb, oxygen saturation.
23. The method of claim 14 wherein adjusting the hue saturation further comprises one or a combination of reducing resolution of hue and saturation color planes by binning the image, resizing the image, and smoothing the image through filtering higher frequency components out.
24. The method of claim 14 further interpolating the smoothed color planes on a grid of higher resolution intensity plane.
PCT/US2006/012461 2005-04-04 2006-04-04 Hyperspectral imaging in diabetes and peripheral vascular disease WO2006107947A2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CA2604829A CA2604829C (en) 2005-04-04 2006-04-04 Hyperspectral imaging in diabetes and peripheral vascular disease

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US66767705P 2005-04-04 2005-04-04
US60/667,677 2005-04-04
US78597706P 2006-03-27 2006-03-27
US60/785,977 2006-03-27

Publications (2)

Publication Number Publication Date
WO2006107947A2 true WO2006107947A2 (en) 2006-10-12
WO2006107947A3 WO2006107947A3 (en) 2007-11-29

Family

ID=37074032

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2006/012461 WO2006107947A2 (en) 2005-04-04 2006-04-04 Hyperspectral imaging in diabetes and peripheral vascular disease

Country Status (3)

Country Link
US (3) US8374682B2 (en)
CA (2) CA2604829C (en)
WO (1) WO2006107947A2 (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008103918A1 (en) * 2007-02-22 2008-08-28 Wisconsin Alumni Research Foundation Hyperspectral imaging spectrometer for early detection of skin cancer
WO2008112135A2 (en) * 2007-03-09 2008-09-18 Nellcor Puritan Bennett Llc Method and aparatus for spectroscopic tissue analyte measurement
EP2181334A1 (en) * 2007-08-22 2010-05-05 Phadia AB Read-out method and apparatus
JP2013240401A (en) * 2012-05-18 2013-12-05 Hoya Corp Electronic endoscope apparatus
US9274046B2 (en) 2010-04-30 2016-03-01 Chemimage Corporation System and method for gross anatomic pathology using hyperspectral imaging
CN106108857A (en) * 2016-08-23 2016-11-16 王伏声 Diabetic foot classification of severity evaluates system
JP2017000836A (en) * 2016-09-27 2017-01-05 Hoya株式会社 Electronic endoscope apparatus
CN108903962A (en) * 2018-04-27 2018-11-30 沈阳东软医疗系统有限公司 A kind of data processing method and device
EP4194821A1 (en) * 2018-02-15 2023-06-14 Viavi Solutions Inc. Sensor device and methods of use

Families Citing this family (171)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6036924A (en) * 1997-12-04 2000-03-14 Hewlett-Packard Company Cassette of lancet cartridges for sampling blood
US6391005B1 (en) 1998-03-30 2002-05-21 Agilent Technologies, Inc. Apparatus and method for penetration with shaft having a sensor for sensing penetration depth
US8641644B2 (en) 2000-11-21 2014-02-04 Sanofi-Aventis Deutschland Gmbh Blood testing apparatus having a rotatable cartridge with multiple lancing elements and testing means
US8337419B2 (en) 2002-04-19 2012-12-25 Sanofi-Aventis Deutschland Gmbh Tissue penetration device
US7981056B2 (en) 2002-04-19 2011-07-19 Pelikan Technologies, Inc. Methods and apparatus for lancet actuation
US9427532B2 (en) 2001-06-12 2016-08-30 Sanofi-Aventis Deutschland Gmbh Tissue penetration device
US9795747B2 (en) 2010-06-02 2017-10-24 Sanofi-Aventis Deutschland Gmbh Methods and apparatus for lancet actuation
ATE497731T1 (en) * 2001-06-12 2011-02-15 Pelikan Technologies Inc DEVICE FOR INCREASING THE SUCCESS RATE OF BLOOD YIELD OBTAINED BY A FINGER PICK
ATE485766T1 (en) 2001-06-12 2010-11-15 Pelikan Technologies Inc ELECTRICAL ACTUATING ELEMENT FOR A LANCET
WO2002100254A2 (en) 2001-06-12 2002-12-19 Pelikan Technologies, Inc. Method and apparatus for lancet launching device integrated onto a blood-sampling cartridge
US7344507B2 (en) 2002-04-19 2008-03-18 Pelikan Technologies, Inc. Method and apparatus for lancet actuation
US7025774B2 (en) * 2001-06-12 2006-04-11 Pelikan Technologies, Inc. Tissue penetration device
US7316700B2 (en) 2001-06-12 2008-01-08 Pelikan Technologies, Inc. Self optimizing lancing device with adaptation means to temporal variations in cutaneous properties
CA2448905C (en) * 2001-06-12 2010-09-07 Pelikan Technologies, Inc. Blood sampling apparatus and method
US9226699B2 (en) 2002-04-19 2016-01-05 Sanofi-Aventis Deutschland Gmbh Body fluid sampling module with a continuous compression tissue interface surface
US7331931B2 (en) * 2002-04-19 2008-02-19 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US8221334B2 (en) 2002-04-19 2012-07-17 Sanofi-Aventis Deutschland Gmbh Method and apparatus for penetrating tissue
US7232451B2 (en) * 2002-04-19 2007-06-19 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7674232B2 (en) * 2002-04-19 2010-03-09 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US8784335B2 (en) 2002-04-19 2014-07-22 Sanofi-Aventis Deutschland Gmbh Body fluid sampling device with a capacitive sensor
US7226461B2 (en) 2002-04-19 2007-06-05 Pelikan Technologies, Inc. Method and apparatus for a multi-use body fluid sampling device with sterility barrier release
US7297122B2 (en) * 2002-04-19 2007-11-20 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7481776B2 (en) * 2002-04-19 2009-01-27 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7901362B2 (en) 2002-04-19 2011-03-08 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US8579831B2 (en) 2002-04-19 2013-11-12 Sanofi-Aventis Deutschland Gmbh Method and apparatus for penetrating tissue
US7909778B2 (en) 2002-04-19 2011-03-22 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US9248267B2 (en) 2002-04-19 2016-02-02 Sanofi-Aventis Deustchland Gmbh Tissue penetration device
US7976476B2 (en) * 2002-04-19 2011-07-12 Pelikan Technologies, Inc. Device and method for variable speed lancet
US7485128B2 (en) * 2002-04-19 2009-02-03 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7371247B2 (en) * 2002-04-19 2008-05-13 Pelikan Technologies, Inc Method and apparatus for penetrating tissue
US7648468B2 (en) 2002-04-19 2010-01-19 Pelikon Technologies, Inc. Method and apparatus for penetrating tissue
US9795334B2 (en) 2002-04-19 2017-10-24 Sanofi-Aventis Deutschland Gmbh Method and apparatus for penetrating tissue
US7491178B2 (en) * 2002-04-19 2009-02-17 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7547287B2 (en) 2002-04-19 2009-06-16 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US9314194B2 (en) 2002-04-19 2016-04-19 Sanofi-Aventis Deutschland Gmbh Tissue penetration device
US7410468B2 (en) * 2002-04-19 2008-08-12 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7141058B2 (en) * 2002-04-19 2006-11-28 Pelikan Technologies, Inc. Method and apparatus for a body fluid sampling device using illumination
US7717863B2 (en) 2002-04-19 2010-05-18 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7229458B2 (en) 2002-04-19 2007-06-12 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7892183B2 (en) 2002-04-19 2011-02-22 Pelikan Technologies, Inc. Method and apparatus for body fluid sampling and analyte sensing
US7563232B2 (en) * 2002-04-19 2009-07-21 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7291117B2 (en) 2002-04-19 2007-11-06 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US8702624B2 (en) 2006-09-29 2014-04-22 Sanofi-Aventis Deutschland Gmbh Analyte measurement device with a single shot actuator
US8267870B2 (en) 2002-04-19 2012-09-18 Sanofi-Aventis Deutschland Gmbh Method and apparatus for body fluid sampling with hybrid actuation
US8574895B2 (en) 2002-12-30 2013-11-05 Sanofi-Aventis Deutschland Gmbh Method and apparatus using optical techniques to measure analyte levels
AU2003300154A1 (en) * 2002-12-31 2004-07-29 Pelikan Technologies Inc. Method and apparatus for loading penetrating members
US7582060B2 (en) * 2003-05-05 2009-09-01 Vioptix Inc. Diagnosing peripheral vascular disease by monitoring oxygen saturation changes during a hyperemia phase
US7569017B2 (en) * 2003-05-05 2009-08-04 Vioptix Inc. Diagnosing peripheral vascular disease by monitoring oxygen saturation changes during an ischemia phase
US7553285B2 (en) * 2003-05-05 2009-06-30 Vioptix Inc. Diagnosing peripheral vascular disease by monitoring oxygen saturation changes during an accumulation phase
US7247142B1 (en) 2003-05-05 2007-07-24 Vioptix, Inc. Diagnosis of peripheral vascular disease using oxygen saturation
WO2004107964A2 (en) 2003-06-06 2004-12-16 Pelikan Technologies, Inc. Blood harvesting device with electronic control
US20060241666A1 (en) * 2003-06-11 2006-10-26 Briggs Barry D Method and apparatus for body fluid sampling and analyte sensing
WO2006001797A1 (en) 2004-06-14 2006-01-05 Pelikan Technologies, Inc. Low pain penetrating
US8133177B2 (en) * 2003-09-23 2012-03-13 Vasamed, Inc. System and method for assessing capillary vitality
EP1671096A4 (en) 2003-09-29 2009-09-16 Pelikan Technologies Inc Method and apparatus for an improved sample capture device
EP1680014A4 (en) * 2003-10-14 2009-01-21 Pelikan Technologies Inc Method and apparatus for a variable user interface
EP1706026B1 (en) 2003-12-31 2017-03-01 Sanofi-Aventis Deutschland GmbH Method and apparatus for improving fluidic flow and sample capture
US7822454B1 (en) 2005-01-03 2010-10-26 Pelikan Technologies, Inc. Fluid sampling device with improved analyte detecting member configuration
US8828203B2 (en) 2004-05-20 2014-09-09 Sanofi-Aventis Deutschland Gmbh Printable hydrogels for biosensors
EP1765194A4 (en) 2004-06-03 2010-09-29 Pelikan Technologies Inc Method and apparatus for a fluid sampling device
US8548570B2 (en) * 2004-11-29 2013-10-01 Hypermed Imaging, Inc. Hyperspectral imaging of angiogenesis
CA2631564A1 (en) * 2004-11-29 2006-06-01 Hypermed, Inc. Medical hyperspectral imaging for evaluation of tissue and tumor
US8652831B2 (en) 2004-12-30 2014-02-18 Sanofi-Aventis Deutschland Gmbh Method and apparatus for analyte measurement test time
US7355688B2 (en) * 2005-09-08 2008-04-08 Vioptix, Inc. Optical probe for optical imaging system
WO2007044599A2 (en) * 2005-10-06 2007-04-19 Hamilton Scott E Pod connected data monitoring system
US8478386B2 (en) * 2006-01-10 2013-07-02 Accuvein Inc. Practitioner-mounted micro vein enhancer
US9854977B2 (en) 2006-01-10 2018-01-02 Accuvein, Inc. Scanned laser vein contrast enhancer using a single laser, and modulation circuitry
US8255040B2 (en) 2006-06-29 2012-08-28 Accuvein, Llc Micro vein enhancer
US10813588B2 (en) 2006-01-10 2020-10-27 Accuvein, Inc. Micro vein enhancer
US8838210B2 (en) 2006-06-29 2014-09-16 AccuView, Inc. Scanned laser vein contrast enhancer using a single laser
US11278240B2 (en) 2006-01-10 2022-03-22 Accuvein, Inc. Trigger-actuated laser vein contrast enhancer
US8489178B2 (en) 2006-06-29 2013-07-16 Accuvein Inc. Enhanced laser vein contrast enhancer with projection of analyzed vein data
US9492117B2 (en) 2006-01-10 2016-11-15 Accuvein, Inc. Practitioner-mounted micro vein enhancer
US11253198B2 (en) 2006-01-10 2022-02-22 Accuvein, Inc. Stand-mounted scanned laser vein contrast enhancer
US20070224694A1 (en) * 2006-02-10 2007-09-27 Puchalski Daniel M Method and system for hyperspectral detection of animal diseases
US20120145906A1 (en) * 2006-03-03 2012-06-14 Chemlmage Corporation Portable system for detecting explosives and a method of use thereof
WO2009098696A2 (en) * 2008-02-10 2009-08-13 Empire Bio-Medical Devices Inc. Device and method for enhanced wound treatment
US8665507B2 (en) * 2006-06-29 2014-03-04 Accuvein, Inc. Module mounting mirror endoscopy
US8244333B2 (en) * 2006-06-29 2012-08-14 Accuvein, Llc Scanned laser vein contrast enhancer
US8730321B2 (en) 2007-06-28 2014-05-20 Accuvein, Inc. Automatic alignment of a contrast enhancement system
US8463364B2 (en) * 2009-07-22 2013-06-11 Accuvein Inc. Vein scanner
US8594770B2 (en) 2006-06-29 2013-11-26 Accuvein, Inc. Multispectral detection and presentation of an object's characteristics
US10238294B2 (en) 2006-06-29 2019-03-26 Accuvein, Inc. Scanned laser vein contrast enhancer using one laser
US9254220B1 (en) 2006-08-29 2016-02-09 Vasamed, Inc. Method and system for assessing severity and stage of peripheral arterial disease and lower extremity wounds using angiosome mapping
FR2916339A1 (en) * 2007-05-25 2008-11-28 Commissariat Energie Atomique BI-SPECTRAL PER-OPERATIVE OPTICAL PROBE.
US20080306356A1 (en) * 2007-06-05 2008-12-11 Kenneth Darryl Kemp Vascular status monitoring system
GB0717365D0 (en) * 2007-09-07 2007-10-17 Univ Loughborough Imaging
US20090287109A1 (en) * 2008-05-14 2009-11-19 Searete Llc, A Limited Liability Corporation Of The State Of Delaware Circulatory monitoring systems and methods
US20090287120A1 (en) * 2007-12-18 2009-11-19 Searete Llc, A Limited Liability Corporation Of The State Of Delaware Circulatory monitoring systems and methods
US9717896B2 (en) * 2007-12-18 2017-08-01 Gearbox, Llc Treatment indications informed by a priori implant information
US8636670B2 (en) 2008-05-13 2014-01-28 The Invention Science Fund I, Llc Circulatory monitoring systems and methods
US20090287101A1 (en) * 2008-05-13 2009-11-19 Searete Llc, A Limited Liability Corporation Of The State Of Delaware Circulatory monitoring systems and methods
US20090292212A1 (en) * 2008-05-20 2009-11-26 Searete Llc, A Limited Corporation Of The State Of Delaware Circulatory monitoring systems and methods
US20090287094A1 (en) * 2008-05-15 2009-11-19 Seacrete Llc, A Limited Liability Corporation Of The State Of Delaware Circulatory monitoring systems and methods
US20090287191A1 (en) * 2007-12-18 2009-11-19 Searete Llc, A Limited Liability Corporation Of The State Of Delaware Circulatory monitoring systems and methods
US20090163785A1 (en) * 2007-12-20 2009-06-25 Nelson Riley H Virtual non-invasive blood analysis device workstation and associated methods
EP2265324B1 (en) 2008-04-11 2015-01-28 Sanofi-Aventis Deutschland GmbH Integrated analyte measurement system
US8606344B2 (en) * 2008-04-30 2013-12-10 Board Of Regents, The University Of Texas System Integrated patient bed system
AU2009246917A1 (en) 2008-05-13 2009-11-19 Spectral Image, Inc. Systems and methods for hyperspectral medical imaging using real-time projection of spectral information
WO2009142758A1 (en) * 2008-05-23 2009-11-26 Spectral Image, Inc. Systems and methods for hyperspectral medical imaging
US9117133B2 (en) * 2008-06-18 2015-08-25 Spectral Image, Inc. Systems and methods for hyperspectral imaging
WO2010004365A1 (en) 2008-07-10 2010-01-14 Ecole Polytechnique Federale De Lausanne (Epfl) Functional optical coherent imaging
US8406859B2 (en) * 2008-08-10 2013-03-26 Board Of Regents, The University Of Texas System Digital light processing hyperspectral imaging apparatus
US9375169B2 (en) 2009-01-30 2016-06-28 Sanofi-Aventis Deutschland Gmbh Cam drive for managing disposable penetrating member actions with a single motor and motor and control system
US20100312076A1 (en) 2009-06-03 2010-12-09 Bly Deborah C Pop box methods and systems for assessing risk of pressure related skin injury and determining a recommended treatment thereof
US9061109B2 (en) 2009-07-22 2015-06-23 Accuvein, Inc. Vein scanner with user interface
US8290558B1 (en) 2009-11-23 2012-10-16 Vioptix, Inc. Tissue oximeter intraoperative sensor
US9726487B2 (en) * 2009-12-18 2017-08-08 Vito Nv Geometric referencing of multi-spectral data
US20110242533A1 (en) * 2010-04-05 2011-10-06 Chemimage Corporation System and Method for Detecting Hazardous Agents Including Explosives
US8965476B2 (en) 2010-04-16 2015-02-24 Sanofi-Aventis Deutschland Gmbh Tissue penetration device
US7884933B1 (en) 2010-05-05 2011-02-08 Revolutionary Business Concepts, Inc. Apparatus and method for determining analyte concentrations
US9245336B2 (en) 2010-12-15 2016-01-26 Koninklijke Philips N.V. Contour guided deformable image registration
EP2757933B1 (en) 2011-09-22 2019-02-06 The George Washington University Systems for visualizing ablated tissue
KR20140102643A (en) 2011-09-22 2014-08-22 더 조지 워싱턴 유니버시티 Systems and methods for visualizing ablated tissue
LV14444B (en) * 2011-09-30 2012-02-20 Latvijas Universitāte Method and device for the determination of efficiency of human anesthesia and for registration the same, using contactless photoplethysmography
US8655091B2 (en) * 2012-02-24 2014-02-18 Raytheon Company Basis vector spectral image compression
JP6088770B2 (en) * 2012-04-20 2017-03-01 国立大学法人広島大学 Grain component analysis apparatus and grain component analysis method
EP2840957B1 (en) * 2012-04-27 2024-01-03 Stryker European Operations Limited Optical coherent imaging medical device and method
US8897522B2 (en) * 2012-05-30 2014-11-25 Xerox Corporation Processing a video for vascular pattern detection and cardiac function analysis
US9329086B2 (en) 2012-05-30 2016-05-03 Chemimage Technologies Llc System and method for assessing tissue oxygenation using a conformal filter
RU2616653C2 (en) * 2012-06-05 2017-04-18 Хайпермед Имэджинг, Инк. Methods and device for coaxial image forming with multiple wavelengths
CA2914780C (en) 2012-07-10 2020-02-25 Aimago S.A. Perfusion assessment multi-modality optical medical device
US9072426B2 (en) 2012-08-02 2015-07-07 AccuVein, Inc Device for detecting and illuminating vasculature using an FPGA
US9955900B2 (en) * 2012-10-31 2018-05-01 Quaerimus, Inc. System and method for continuous monitoring of a human foot
US9901298B2 (en) 2012-11-01 2018-02-27 Quaerimus Medical Incorporated System and method for prevention of diabetic foot ulcers using total internal reflection imaging
US10517483B2 (en) 2012-12-05 2019-12-31 Accuvein, Inc. System for detecting fluorescence and projecting a representative image
JP6356141B2 (en) 2012-12-10 2018-07-11 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Medical device or system for measuring hemoglobin levels during an accident using a camera-projector system
JP6290248B2 (en) 2012-12-21 2018-03-07 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. System and method for extracting physiological information from remotely detected electromagnetic radiation
US9107567B2 (en) 2012-12-27 2015-08-18 Christie Digital Systems Usa, Inc. Spectral imaging with a color wheel
US9544582B2 (en) 2013-03-01 2017-01-10 Boston Scientific Scimed, Inc. Image sensor calibration
CA2894133C (en) * 2013-03-15 2016-11-01 Synaptive Medical (Barbados) Inc. Surgical imaging systems
JP6737705B2 (en) 2013-11-14 2020-08-12 ザ・ジョージ・ワシントン・ユニバーシティThe George Washingtonuniversity Method of operating system for determining depth of injury site and system for generating images of heart tissue
JP2017500550A (en) 2013-11-20 2017-01-05 ザ・ジョージ・ワシントン・ユニバーシティThe George Washingtonuniversity System and method for hyperspectral analysis of cardiac tissue
US9841322B1 (en) 2014-06-02 2017-12-12 Kemeny Associates LLC Spectral imaging with multiple illumination sources
CA2954625C (en) 2014-06-18 2022-12-13 Innopix, Inc. Spectral imaging system for remote and noninvasive detection of target substances using spectral filter arrays and image capture arrays
JPWO2015199067A1 (en) * 2014-06-24 2017-06-15 株式会社ニコン Image analysis apparatus, imaging system, surgery support system, image analysis method, and image analysis program
US9968285B2 (en) * 2014-07-25 2018-05-15 Christie Digital Systems Usa, Inc. Multispectral medical imaging devices and methods thereof
WO2016048911A1 (en) * 2014-09-22 2016-03-31 Invuity, Inc. Hyperspectral imager
CN107205624B (en) 2014-10-29 2019-08-06 光谱Md公司 Reflective multispectral time discrimination optics imaging method and equipment for tissue typing
WO2016073476A1 (en) 2014-11-03 2016-05-12 The George Washington University Systems and methods for lesion assessment
CN107427213B (en) 2014-11-03 2021-04-16 460医学股份有限公司 System and method for evaluation of contact quality
US10070796B2 (en) * 2015-02-04 2018-09-11 General Electric Company Systems and methods for quantitative microcirculation state monitoring
JP6346576B2 (en) * 2015-02-27 2018-06-20 Hoya株式会社 Image processing device
US10506961B1 (en) 2015-04-15 2019-12-17 Eden Medical, Inc. Diagnostic transducer and method
CN107735015B (en) 2015-05-07 2021-09-10 史赛克欧洲运营有限公司 Method and system for laser speckle imaging of tissue using a color image sensor
US10779904B2 (en) 2015-07-19 2020-09-22 460Medical, Inc. Systems and methods for lesion formation and assessment
US10244987B2 (en) * 2015-08-13 2019-04-02 Pixart Imaging Inc. Physiological detection system with adjustable signal source and operating method thereof
US9949821B2 (en) * 2015-12-22 2018-04-24 Biosense Webster (Israel) Ltd. Colored silicone for implant safety
US10362998B2 (en) 2016-02-25 2019-07-30 Samsung Electronics Co., Ltd. Sensor-based detection of changes in health and ventilation threshold
US10172517B2 (en) 2016-02-25 2019-01-08 Samsung Electronics Co., Ltd Image-analysis for assessing heart failure
US10420514B2 (en) 2016-02-25 2019-09-24 Samsung Electronics Co., Ltd. Detection of chronotropic incompetence
US11164596B2 (en) 2016-02-25 2021-11-02 Samsung Electronics Co., Ltd. Sensor assisted evaluation of health and rehabilitation
CN109313929A (en) * 2016-04-15 2019-02-05 皇家飞利浦有限公司 Annotation applies associated data point with clinical decision support
EP3275361A1 (en) * 2016-07-29 2018-01-31 Leica Instruments (Singapore) Pte. Ltd. Medical imaging apparatus and method for the imaging of a light-sensitive object, such as biological tissue
EP3284396B1 (en) * 2016-08-16 2020-02-12 Leica Instruments (Singapore) Pte. Ltd. Observation apparatus and method for visual enhancement of an observed object
CA3034761A1 (en) 2016-08-24 2018-03-01 Mimosa Diagnostics Inc. Multispectral mobile tissue assessment
WO2018160963A1 (en) 2017-03-02 2018-09-07 Spectral Md, Inc. Machine learning systems and techniques for multispectral amputation site analysis
SG11201908706QA (en) 2017-03-22 2019-10-30 Modulated Imaging Inc Method and apparatus for assessing diabetic circulatory complications
WO2018209090A1 (en) * 2017-05-11 2018-11-15 Kent State University Microcirculation assessment device
WO2019060733A2 (en) * 2017-09-22 2019-03-28 Intuitive Surgical Operations, Inc. Enhancing visible differences between different tissues in computer-assisted tele-operated surgery
US10740884B2 (en) 2018-12-14 2020-08-11 Spectral Md, Inc. System and method for high precision multi-aperture spectral imaging
US10783632B2 (en) * 2018-12-14 2020-09-22 Spectral Md, Inc. Machine learning systems and method for assessment, healing prediction, and treatment of wounds
BR112021011113A2 (en) 2018-12-14 2021-08-31 Spectral Md, Inc. SYSTEM AND METHOD FOR FORMING A SPECTRAL IMAGE OF MULTIPLE HIGH PRECISION APERTURES
KR102632479B1 (en) * 2019-06-17 2024-01-31 삼성전자주식회사 Ppg sensor and method for operating ppg sensor
US10878568B1 (en) 2019-08-08 2020-12-29 Neuroptica, Llc Systems and methods for imaging disease biomarkers
JP2023502951A (en) * 2019-11-15 2023-01-26 ブルーイン、バイオメトリクス、リミテッド、ライアビリティー、カンパニー spectral imaging
WO2021245374A1 (en) 2020-06-03 2021-12-09 King's College London Method and system for joint demosaicking and spectral signature estimation
WO2021262844A1 (en) * 2020-06-23 2021-12-30 Chemimage Corporation Methods for determining blood oxygenation and tissue perfusion levels and devices thereof
US20210400211A1 (en) * 2020-06-23 2021-12-23 Woundtech Multi-modal mobile thermal imaging system
WO2023055025A1 (en) * 2021-09-29 2023-04-06 Samsung Electronics Co., Ltd. Method and electronic device for determining skin information using hyper spectral reconstruction
CN113796852B (en) * 2021-09-30 2023-09-08 太原理工大学 Diabetes foot prediction method based on gradient lifting decision tree model algorithm

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6556853B1 (en) * 1995-12-12 2003-04-29 Applied Spectral Imaging Ltd. Spectral bio-imaging of the eye
US20040209237A1 (en) * 2003-04-18 2004-10-21 Medispectra, Inc. Methods and apparatus for characterization of tissue samples

Family Cites Families (83)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL69528A (en) 1983-08-19 1988-06-30 Brig Res Ltd Apparatus for the prevention of pressure sores
US4647918A (en) 1985-01-16 1987-03-03 Goforth William P Multi-event notification system for monitoring critical pressure points on persons with diminished sensation of the feet
US4805016A (en) * 1986-08-25 1989-02-14 Kabushiki Kaisha Toshiba Endoscopic system for converting primary color images into hue, saturation and intensity images
US4871252A (en) * 1986-10-30 1989-10-03 The Regents Of The University Of California Method and apparatus for object positioning
DE3854259D1 (en) 1987-05-21 1995-09-07 Peter Seitz METHOD AND DEVICE FOR PRODUCING INSERTS OR THE LIKE.
US4947850A (en) 1988-03-11 1990-08-14 Trustees Of The University Of Pennsylvania Method and apparatus for imaging an internal body portion of a host animal
US5438989A (en) 1990-08-10 1995-08-08 Hochman; Darryl Solid tumor, cortical function, and nerve tissue imaging methods and device
US5198977A (en) * 1990-11-27 1993-03-30 Jesse Salb System and method for localization of functional activity in the human brain
JPH04106748U (en) * 1991-02-28 1992-09-14 株式会社島津製作所 Optical biomeasuring device
US5645081A (en) 1991-11-14 1997-07-08 Wake Forest University Method of treating tissue damage and apparatus for same
US5277181A (en) * 1991-12-12 1994-01-11 Vivascan Corporation Noninvasive measurement of hematocrit and hemoglobin content by differential optical analysis
US5642096A (en) 1992-03-20 1997-06-24 Paromed Medizintechnik Gmbh Device for prevention of ulcers in the feet of diabetes patients
EP0767361B1 (en) * 1993-07-22 2000-02-23 Applied Spectral Imaging Ltd. Method and apparatus for spectral imaging
US5349954A (en) 1993-07-23 1994-09-27 General Electric Company Tumor tissue characterization apparatus and method
US5995645A (en) * 1993-08-18 1999-11-30 Applied Spectral Imaging Ltd. Method of cancer cell detection
EP0649117A3 (en) * 1993-10-15 1996-01-31 George S Allen Method for providing medical images.
US5782770A (en) 1994-05-12 1998-07-21 Science Applications International Corporation Hyperspectral imaging methods and apparatus for non-invasive diagnosis of tissue for cancer
US5566473A (en) 1994-08-16 1996-10-22 Salminen; Reijo K. Processing roll apparatus and method for web drying
US5722287A (en) 1995-05-31 1998-03-03 Forstein; Micah Aaron Video pedobarograph system
US20040127800A1 (en) 1995-07-06 2004-07-01 Kimball Victor E. Device for assessing perfusion failure in a patient by measurement of blood flow
JP3568280B2 (en) 1995-07-12 2004-09-22 富士写真フイルム株式会社 Surgical operation support system
US6573063B2 (en) 1995-10-04 2003-06-03 Cytoscan Sciences, Llc Methods and systems for assessing biological materials using optical and spectroscopic detection techniques
US6240309B1 (en) * 1995-10-06 2001-05-29 Hitachi, Ltd. Optical measurement instrument for living body
CA2235772C (en) 1995-10-23 2002-12-31 Cytometrics, Inc. Method and apparatus for reflected imaging analysis
US6553327B2 (en) * 1998-09-16 2003-04-22 Yeda Research & Development Co., Ltd. Apparatus for monitoring a system with time in space and method therefor
US6078681A (en) * 1996-03-18 2000-06-20 Marine Biological Laboratory Analytical imaging system and process
US6544193B2 (en) 1996-09-04 2003-04-08 Marcio Marc Abreu Noninvasive measurement of chemical substances
US6208749B1 (en) * 1997-02-28 2001-03-27 Electro-Optical Sciences, Inc. Systems and methods for the multispectral imaging and characterization of skin tissue
WO1998044839A1 (en) 1997-04-03 1998-10-15 National Research Council Of Canada Method of assessing tissue viability using near-infrared spectroscopy
WO1999022640A2 (en) 1997-10-30 1999-05-14 Hypermed Imaging, Inc. Multispectral/hyperspectral medical instrument
US6937885B1 (en) 1997-10-30 2005-08-30 Hypermed, Inc. Multispectral/hyperspectral medical instrument
US6070093A (en) 1997-12-02 2000-05-30 Abbott Laboratories Multiplex sensor and method of use
JPH11274818A (en) * 1998-03-20 1999-10-08 Fujitsu Ltd High frequency circuit device
US6487440B2 (en) 1998-07-08 2002-11-26 Lifespex, Inc. Optical probe having and methods for difuse and uniform light irradiation
US6083158A (en) 1998-09-01 2000-07-04 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Real-time visualization of tissue ischemia
US6741884B1 (en) 1998-09-03 2004-05-25 Hypermed, Inc. Infrared endoscopic balloon probes
US8024027B2 (en) 1998-09-03 2011-09-20 Hyperspectral Imaging, Inc. Infrared endoscopic balloon probes
US6438396B1 (en) 1998-11-05 2002-08-20 Cytometrics, Inc. Method and apparatus for providing high contrast imaging
US6859275B2 (en) 1999-04-09 2005-02-22 Plain Sight Systems, Inc. System and method for encoded spatio-spectral information processing
US20030026762A1 (en) 1999-05-05 2003-02-06 Malmros Mark K. Bio-spectral imaging system and methods for diagnosing cell disease state
EP1196081B1 (en) 1999-07-02 2013-08-21 HyperMed Imaging, Inc. Integrated imaging apparatus
US6587711B1 (en) 1999-07-22 2003-07-01 The Research Foundation Of Cuny Spectral polarizing tomographic dermatoscope
EP1211977B1 (en) 1999-07-26 2011-05-18 Cardiosense Ltd. Apparatus for the detection of physiological distress in a patient
US7483733B2 (en) 1999-07-26 2009-01-27 Cardiosense, Ltd. Non-invasive method and apparatus to detect and monitor early medical shock, and related conditions
US6750964B2 (en) 1999-08-06 2004-06-15 Cambridge Research And Instrumentation, Inc. Spectral imaging methods and systems
US7840257B2 (en) 2003-01-04 2010-11-23 Non Invasive Technology, Inc. Examination of biological tissue using non-contact optical probes
US6122846A (en) 1999-08-30 2000-09-26 Frank B. Gray Force monitoring shoe
US6640132B1 (en) * 1999-11-17 2003-10-28 Hypermed, Inc. Forensic hyperspectral apparatus and method
AU1547101A (en) * 1999-11-26 2001-06-04 Applied Spectral Imaging Ltd. System and method for functional brain mapping and an oxygen saturation difference map algorithm for effecting same
US6665557B1 (en) * 1999-12-01 2003-12-16 The Research Foundation Of City College Of New York Sprectroscopic and time-resolved optical methods and apparatus for imaging objects in turbed media
US6918883B2 (en) 1999-12-15 2005-07-19 Cannon Research Institute Of Carolinas Medical Center Sock for detection of pressure points on feet
US7239909B2 (en) * 2000-01-19 2007-07-03 Luminetx Technologies Corp. Imaging system using diffuse infrared light
US6556858B1 (en) 2000-01-19 2003-04-29 Herbert D. Zeman Diffuse infrared light imaging system
US6834122B2 (en) 2000-01-22 2004-12-21 Kairos Scientific, Inc. Visualization and processing of multidimensional data using prefiltering and sorting criteria
EP1267709B1 (en) 2000-03-28 2009-04-29 Board of Regents, The University of Texas System Method and apparatus for diagnostic multispectral digital imaging
CA2401234A1 (en) 2000-04-13 2001-10-25 Mark D. Hewko Tissue viability/health monitor utilizing near infrared spectroscopy
DE60130170T2 (en) * 2000-06-06 2008-05-21 Fujifilm Corp. Fluorescent light Image display method and apparatus therefor
ATE417543T1 (en) 2000-07-13 2009-01-15 Univ Virginia Commonwealth USE OF ULTRAVIOLET, NEAR ULTRAVIOLET AND NEAR INFRARED RESONANCE RMAN SPECTROSCOPY AND FLUORESENSENCE SPECTOSCOPY FOR EXAMINING TISSUE FOR SHOCK, CRITICAL ILLNESSES OR OTHER MEDICAL CONDITIONS
US6678398B2 (en) * 2000-09-18 2004-01-13 Sti Medical Systems, Inc. Dual mode real-time screening and rapid full-area, selective-spectral, remote imaging and analysis device and process
AU2002253784A1 (en) * 2000-11-07 2002-08-28 Hypermed, Inc. Hyperspectral imaging calibration device
JP2002158893A (en) * 2000-11-22 2002-05-31 Minolta Co Ltd Device and method for correcting picture and recording medium
US6958076B2 (en) 2001-04-16 2005-10-25 Biomedical Research Associates Inc. Implantable venous valve
US6728561B2 (en) 2001-08-14 2004-04-27 University Of Alabama In Huntsville Multispectral image processing method, apparatus and computer program product for determining the blood oxygen saturation in a vessel
US7049597B2 (en) 2001-12-21 2006-05-23 Andrew Bodkin Multi-mode optical imager
JP2003334163A (en) * 2002-03-14 2003-11-25 Olympus Optical Co Ltd Endoscopic image processor
US6826336B2 (en) 2002-05-22 2004-11-30 The Boeing Company Fiber optic LED illuminator
AU2004227886A1 (en) 2003-04-04 2004-10-21 Lumidigm, Inc. Multispectral biometric sensor
JP4406226B2 (en) * 2003-07-02 2010-01-27 株式会社東芝 Biological information video device
WO2005008198A2 (en) 2003-07-18 2005-01-27 Chemimage Corporation Method and apparatus for compact birefringent interference imaging spectrometer
US20050049467A1 (en) 2003-08-28 2005-03-03 Georgios Stamatas Method for assessing pigmented skin
US7920908B2 (en) 2003-10-16 2011-04-05 David Hattery Multispectral imaging for quantitative contrast of functional and structural features of layers inside optically dense media such as tissue
AU2005226667B2 (en) * 2004-03-19 2010-04-01 Board Of Supervisors Of Louisiana State University And Agricultural And Mechanical College A method for evaluating relative oxygen saturation in body tissues
US7166852B2 (en) * 2004-04-06 2007-01-23 Accuray, Inc. Treatment target positioning system
US20090163809A1 (en) 2004-06-03 2009-06-25 Kane Scott D Medical method and associated apparatus utilizable in accessing internal organs through skin surface
EP1818016B1 (en) * 2004-07-20 2019-02-20 Toshinori Kato Biofunction diagnosis device, biofunction diagnosis method, bioprobe, bioprobe wearing tool, bioprobe support tool, and bioprobe wearing assisting tool
US8548570B2 (en) * 2004-11-29 2013-10-01 Hypermed Imaging, Inc. Hyperspectral imaging of angiogenesis
CA2631564A1 (en) 2004-11-29 2006-06-01 Hypermed, Inc. Medical hyperspectral imaging for evaluation of tissue and tumor
CA2590227A1 (en) 2004-12-07 2006-06-07 Clean Earth Technologies, Llc Method and apparatus for standoff detection of liveness
WO2006086085A2 (en) 2004-12-28 2006-08-17 Hypermed, Inc. Hyperspectral/multispectral imaging in determination, assessment and monitoring of systemic physiology and shock
US20060241495A1 (en) 2005-03-23 2006-10-26 Eastman Kodak Company Wound healing monitoring and treatment
US8971984B2 (en) * 2005-04-04 2015-03-03 Hypermed Imaging, Inc. Hyperspectral technology for assessing and treating diabetic foot and tissue disease
DE602006020618D1 (en) 2005-12-22 2011-04-21 Visen Medical Inc COMBINED X-RAY AND OPTICAL TOMOGRAPHY IMAGING SYSTEM
US9044179B2 (en) 2006-01-04 2015-06-02 The Trustees Of The University Of Pennslyvania Oxygen sensor for internal monitoring of tissue oxygen in vivo

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6556853B1 (en) * 1995-12-12 2003-04-29 Applied Spectral Imaging Ltd. Spectral bio-imaging of the eye
US20040209237A1 (en) * 2003-04-18 2004-10-21 Medispectra, Inc. Methods and apparatus for characterization of tissue samples

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008103918A1 (en) * 2007-02-22 2008-08-28 Wisconsin Alumni Research Foundation Hyperspectral imaging spectrometer for early detection of skin cancer
US8315692B2 (en) 2007-02-22 2012-11-20 Sheinis Andrew I Multi-spectral imaging spectrometer for early detection of skin cancer
WO2008112135A2 (en) * 2007-03-09 2008-09-18 Nellcor Puritan Bennett Llc Method and aparatus for spectroscopic tissue analyte measurement
WO2008112135A3 (en) * 2007-03-09 2008-11-27 Nellcor Puritan Bennett Llc Method and aparatus for spectroscopic tissue analyte measurement
EP2181334A1 (en) * 2007-08-22 2010-05-05 Phadia AB Read-out method and apparatus
EP2181334A4 (en) * 2007-08-22 2010-10-27 Phadia Ab Read-out method and apparatus
US9274046B2 (en) 2010-04-30 2016-03-01 Chemimage Corporation System and method for gross anatomic pathology using hyperspectral imaging
EP2850993A4 (en) * 2012-05-18 2016-01-27 Hoya Corp Electronic endoscope device
CN104486983A (en) * 2012-05-18 2015-04-01 Hoya株式会社 Electronic endoscope device
JP2013240401A (en) * 2012-05-18 2013-12-05 Hoya Corp Electronic endoscope apparatus
CN105996965A (en) * 2012-05-18 2016-10-12 Hoya株式会社 Electronic endoscope device
CN106108857A (en) * 2016-08-23 2016-11-16 王伏声 Diabetic foot classification of severity evaluates system
JP2017000836A (en) * 2016-09-27 2017-01-05 Hoya株式会社 Electronic endoscope apparatus
EP4194821A1 (en) * 2018-02-15 2023-06-14 Viavi Solutions Inc. Sensor device and methods of use
US11796389B2 (en) 2018-02-15 2023-10-24 Viavi Solutions Inc. Sensor device and methods of use
CN108903962A (en) * 2018-04-27 2018-11-30 沈阳东软医疗系统有限公司 A kind of data processing method and device
CN108903962B (en) * 2018-04-27 2022-02-11 东软医疗系统股份有限公司 Data processing method and device

Also Published As

Publication number Publication date
US8655433B2 (en) 2014-02-18
US20070016079A1 (en) 2007-01-18
US20130245455A1 (en) 2013-09-19
US20140112559A1 (en) 2014-04-24
WO2006107947A3 (en) 2007-11-29
CA2947613C (en) 2019-11-05
US8374682B2 (en) 2013-02-12
CA2604829A1 (en) 2006-10-12
CA2604829C (en) 2018-05-15
CA2947613A1 (en) 2006-10-12
US10779773B2 (en) 2020-09-22

Similar Documents

Publication Publication Date Title
US10779773B2 (en) Hyperspectral imaging in diabetes and peripheral vascular disease
US10321869B2 (en) Systems and methods for combining hyperspectral images with color images
US8224425B2 (en) Hyperspectral imaging in diabetes and peripheral vascular disease
CA2592691C (en) Hyperspectral/multispectral imaging in determination, assessment and monitoring of systemic physiology and shock
US11026582B2 (en) OxyVu-1 hyperspectral tissue oxygenation (HTO) measurement system
US10117582B2 (en) Medical hyperspectral imaging for evaluation of tissue and tumor
US10470694B2 (en) Systems and methods for measuring tissue oxygenation
US20150265195A1 (en) Systems and methods for measuring tissue oxygenation
US10013811B2 (en) Hyperspectral image visualization in patients with medical conditions
US9326715B1 (en) OxyVu-1 hyperspectral tissue oxygenation (HTO) measurement system
CA2979384C (en) Systems and methods for measuring tissue oxygenation
WO2007022508A2 (en) Hyperspectral technology for assessing and treating diabetic foot and tissue disease
AU2013202796B2 (en) Hyperspectral/multispectral imaging in determination, assessment and monitoring of systemic physiology and shock
CA2656016A1 (en) Hyperspectral technology for assessing and treating diabetic foot and tissue disease

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application
ENP Entry into the national phase

Ref document number: 2604829

Country of ref document: CA

NENP Non-entry into the national phase

Ref country code: DE

NENP Non-entry into the national phase

Ref country code: RU

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: COMMUNICATION PURSUANT TO RULE 112(1) EPC SENT ON 13.12.07

122 Ep: pct application non-entry in european phase

Ref document number: 06749224

Country of ref document: EP

Kind code of ref document: A2