WO2008098284A1 - An edge detection method and a particle counting method - Google Patents
An edge detection method and a particle counting method Download PDFInfo
- Publication number
- WO2008098284A1 WO2008098284A1 PCT/AU2008/000165 AU2008000165W WO2008098284A1 WO 2008098284 A1 WO2008098284 A1 WO 2008098284A1 AU 2008000165 W AU2008000165 W AU 2008000165W WO 2008098284 A1 WO2008098284 A1 WO 2008098284A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- target pixel
- edge
- pixels
- array
- intensity value
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
Definitions
- An aspect of the present invention relates to image processing for edge detection. Embodiments of this aspect of the invention find application, though not exclusively, in the decomposition of complex waveforms.
- Another aspect of the present invention relates to the counting of particles.
- Embodiments of this aspect of the invention find application, though not exclusively, in the medical field for use in the counting of particles in signal images, for example the counting of emboli in signal images generated by the use of a transcranial Doppler machine to analyse a patient's blood flow.
- Transcranial Doppler machines work by transmitting sound above the auditory range into a patient with enough power output to travel through the cranium and into the brain. Signals are reflected within the patient and detected by a probe. It is possible to analyse a specific depth within the patient by assuming a mean speed of sound in tissue and waiting a corresponding time period before listening to echo's reflected back to the probe.
- Doppler effects are used to calculate the velocity relative to the ultrasound beam of the object within the patient that reflected the sound. This relative velocity is proportional to the difference between the transmitted and received frequencies. The actual velocity of the object within the patient that reflected the sound is the relative velocity multiplied by the cosine of the angle of insonation.
- Blood is made up of numerous cells, which typically travel with a range of velocities dependant on the turbulence of the blood flow at the point of interest. Therefore, when analysing blood flow using a TCD machine, the reflected signal received by the probe is typically complex. Known fast Fourier transform techniques are typically performed on the received waveform in order to decompose the signal into the set of velocities from which it was formed. These may be displayed visually with velocity on the y-axis and time on the x-axis. The intensity of the signal is typically shown as different colours; usually with white as the most intense signal.
- High Intensity Transients are short, intense signals that may be generated by particulate material travelling through the Doppler sample gate.
- HITS may be recorded individually, or in large groups that are referred to as “flurries", “showers” or “curtains”. Such events may overload the Doppler spectra, producing a "white-out” of signal.
- the known prior art methods for doing so focus analysis upon the sound signal received at the probe. Using such prior art techniques, it is typically only possible to detect single particles,
- the method further includes the step of defining a respective integral value for the target pixel and each of the neighbouring pixels, wherein each of said respective integral values are dependent upon whether the respective intensity value of the respective pixel exceeds the threshold intensity value.
- a threshold edge detection value is defined.
- the respective integral values are selected such that an aggregate value calculated from a summation of the respective integral values is greater than, or equal to, the threshold edge detection value only if the respective intensity values of the target pixel and no more than six of the neighbouring pixels exceed the threshold intensity value.
- the cellular neural network processing model defines a receptive field having a two dimensional array of three rows and three columns, the receptive field receiving as input the respective intensity values of the target pixel and the eight neighbouring pixels.
- a computer implemented method of processing an image having a plurality of pixels so as to designate at least one edge exceeding a threshold intensity value, each of the pixels having a respective intensity value said method including the steps of: applying a method as described above to a target pixel selected from amongst the plurality of pixels and storing a resultant edge designation for the target pixel; redefining the target pixel from amongst the plurality of pixels and repeating the preceding step until edge designations for substantially all of the pixels have been completed.
- the resultant edge determinations are stored in a two-dimensional Boolean array, which in some embodiments is displayable as an image.
- a computer implemented method of inferring a particle count from a signal image having a plurality of pixels, each of the pixels having a respective intensity value said method including the steps of: processing the signal image using an edge designation algorithm so as to designate edges exceeding a threshold intensity value; formulating an array representative of the edges; processing the array so as to detect peaks defined by said edges; counting the number of peaks; and using a correlation scheme responsive to the number of peaks to infer a particle count.
- the edge designation algorithm is a computer implemented method as described above with reference to the first and/or second aspect of the invention.
- the signal image is derived from an analysis of a blood flow within a patient.
- the signal image may be the output from a transcranial Doppler machine.
- the array is defined within a plurality of memory addresses and/or data storage addresses.
- the array is two dimensional and is displayable as an image.
- the correlation scheme defines a one-to-one relationship between peaks and particles.
- a computer- readable medium containing computer executable code for instructing a computer to perform the method according to any one of the preceding aspects.
- a downloadable or remotely executable file or combination of files containing computer executable code for instructing a computer to perform a method according to any one of the preceding aspects .
- a computing apparatus having a central processing unit, associated memory and storage devices, and input and output devices, said apparatus being configured to perform a method according to any one of the preceding aspects .
- Figure 1 is a schematic depiction of an array of nine pixels
- Figure 2 is a sample of a raw image as produced by a TCD machine
- Figure 3 is a first sample of spectral image data extracted from an image produced by a TCD machine
- Figure 4 is a second sample of spectral image data extracted from an image produced by a TCD machine
- Figure 5 is a sample of the edge designation array image arising from an analysis of the spectral sample of figure 3 ;
- Figure 6 is an example of the edge designation array image arising from an analysis of the spectral sample of figure 4;
- Figure 7 is an example of an edge designation array;
- Figure 8 is an example of a height array
- Figure 9 is a flow chart showing steps performed in an embodiment of the invention from an end -user's perspective; and Figure 10 is a flow chart showing general processing steps performed in an embodiment of the invention.
- a preferred embodiment of the present invention uses a computer to process a target pixel of an image.
- the image is formed from a square or rectangular array of pixels, each of which has an intensity value lying within the range of 0 to 255 inclusive.
- the pixels making up the image typically vary in intensity values from regions having a relatively lower intensity value to regions having a relatively higher intensity value.
- the preferred embodiment aims to designate whether the target pixel lies on an edge at which the intensity starts to exceed a threshold intensity value.
- This processing is performed with reference to not only the target pixel, but also the adjacent eight neighbouring pixels. That is, the processing takes place with reference to a 3 by 3 array as shown schematically in figure 1, in which the pixels are labeled for ease of reference as pixels 1 to 9.
- the target pixel is pixel 5 and its neighbouring eight pixels are pixels 1 to 4 and pixels 6 to 9.
- the processor compares the intensity value of the target pixel to the threshold intensity value. If the intensity value of the target pixel is less than the threshold intensity value, then the pixel is designated as not being an edge pixel. However, if the intensity value of the target pixel is greater than the threshold value, the processor proceeds to compare the respective intensity values of the eight neighbouring pixels to the threshold intensity value. If the respective intensity values of no more than six of the eight neighbouring pixels exceeds the threshold intensity value, then the target pixel is designated as an edge pixel. Alternatively, if the respective intensity values of seven or eight of the eight neighbouring pixels exceeds the threshold intensity value, then the target pixel is designated as not being an edge pixel.
- Table 1 The various possibilities are shown below in table 1.
- the threshold intensity is arbitrarily set to a value of
- the processor compares the intensity value of target pixel 5, which is 160, to the threshold intensity value, which is 150, and thereby determines that the intensity value of target pixel 5 exceeds the threshold intensity value.
- the processor compares the intensity values of the eight neighbouring pixels 1 to 4 and 6 to 9 and determines that six out of eight of those values exceed the threshold of 150 (i.e. pixels 1, 2, 4, 7, 8 and 9).
- the processor designates the target pixel 5 of this example as an edge pixel.
- Another embodiment of the invention implements a cellular neural network processing model.
- This model has a receptive field in the form of a two dimensional, 3 by 3 array of cells, which corresponds to the 3 by 3 array of pixels centered on the target pixel.
- the receptive field receives as input the respective intensity values of the corresponding pixels.
- Each of the cells has an internal state called a "potential", which, in the preferred embodiment, is pre-set to -60V.
- Cells can receive and send "information", in the form of a positive or negative potential to their immediate neighbours. This positive or negative value is added to the target cell's internal potential.
- the potential available to be sent by each of the cells is determined by the processor based upon the position of the cell within the receptive field and upon the intensity values of the corresponding nine pixels.
- the values of the resting potentials, and the amounts by which potentials may be varied are all either positive or negative integral values.
- the central cell starts at the resting potential of -60V, however if the intensity value of the target pixel is greater than the threshold intensity value, then the centre cell can increase its potential by 160V to 100V.
- Each of the outer eight cells has the capacity to send a negative potential of 1 IV to the central cell. This will only occur for any one of the eight outer cells if the intensity value of the corresponding pixel is greater than the threshold intensity value. Once the eight outer cells have sent their "information" to the central cell, an aggregate value is determined, which is the resultant potential of the central cell.
- the processor compares the aggregate value against a threshold edge detection value, which in the preferred embodiment is 30V. This determines whether or not the target pixel should be designated as an edge pixel. If the aggregate value is greater than 30V, then the target pixel is designated as an edge pixel. If the potential of the central cell is less than 30V, then the target pixel is designated as not an edge pixel. It will be appreciated by those skilled in the art that the integral values of -60V, 160V, -11V and 30V, as used in the preferred embodiment, have been selected to ensure that if no more than six surround cells send a negative value, and the centre cell sends a positive value, then the aggregate value will be above the threshold edge detection value. On the other hand, the aggregate value will be less than the threshold edge detection value in the following circumstances:
- this cellular neural network approach makes use of integral mathematics and it will be appreciated by those skilled in the art that computers are adapted to perform integral mathematics at a far higher rate as compared to the processing of floating point arithmetic. Applying this cellular neural network model to the example pixel intensity values shown in figure 1, the centre cell's potential is increased by 160V from a resting potential of -60V due to the target pixel having an intensity value in excess of the threshold intensity value.
- the processor designates the target pixel to be an edge pixel.
- the preferred embodiment may also be used to extract edges that are less than the threshold intensity.
- One technique for doing so is to start by inverting the image (i.e. transposing high intensity pixels into low intensity pixels and vice versa).
- Another technique is to alter the mathematics of the analysis so as to focus upon edges that are less than the threshold.
- the method described above to process an edge designation for a single target pixel may be used iteratively to process edge designations for all, or substantially all, of the pixels within an image.
- the result is stored in a digitally accessible format, such as in the computer's random access memory, or some other digital storage medium.
- the target pixel is then redefined and the process is repeated.
- the re-definition of the target pixel involves indexing along one of the axes of the image, for example the x-axis, and performing an edge designation for each pixel on that axis.
- the processor indexes to the next row of the x-axis and repeats the edge designation process for that row. This process continues along adjacent rows until edge designations for all, or substantially all, of the pixels within the image have been completed.
- no edge designations are processed for the pixels in the far top and bottom rows and the far left and right hand columns, as they are not surrounded by a full compliment of eight neighbouring pixels.
- Each edge determination can have only one of two possible results; either the target pixel is an edge pixel, or it is not.
- the results of all edge determinations for a given image may be stored in a two-dimensional Boolean array having substantially the same dimensions as the original image. If the pixel at position (x,y) of the image is designated as an edge pixel, then a value of 1 is stored in the (x,y) position of the Boolean array. Alternatively, if the pixel at position (x,y) of the image is designated as not an edge pixel, then a value of 0 is stored in the (x,y) position of the Boolean array.
- this type of two dimensional Boolean array is displayable as an image simply by depicting a 0 value in one colour (for example white) and a 1 value in another colour (for example black) . This allows a user to display an image of the edges extracted from the original image.
- another preferred embodiment of the invention is a computer implemented method of inferring a particle count from a signal image.
- the general steps performed by the processor in this method are depicted in figure 10.
- TCD images are derived from an interaction of sound waves with blood flow within a patient.
- the preferred embodiment is particularly adapted for use in analyzing the signal images that are produced by TCD machines so as to infer a count of the number of emboli in the patient's blood that flowed past the TCD sensor during formation of the TCD image.
- the images comprise a plurality of pixels, each of which has a respective intensity value that is representative of the intensity of the signal received by the probe of the TCD machine.
- the processor commences the method by receiving image data 20 from a TCD machine.
- image data 20 An example of such image data 20 is shown in figure 2.
- the processor is programmed to automatically extract the Doppler spectral data portion of this image based upon a known standard positioning of the spectral data within the overall image data.
- the user manually inputs the coordinates of the spectra data, for example by using a mouse to point at the edges of the spectral data sections of the overall image. This step in the processing is shown at step 21 of figure
- the processor then defines a two-dimensional image intensity array within the computer's random access memory or other data storage means.
- the size of the image intensity array corresponds to the size of the extracted spectral image 20 such that each of the dimensions of the image intensity array correspond to the number of pixels in the x and y axes respectively.
- the extracted spectral image is 300 pixels by 200 pixels, then an image intensity array also having 300 columns and 200 rows is defined within the computer's memory.
- step 22 implements step 22 as shown in figure 10 wherein parameters for correlating a particular colour with an intensity value are pre-set based on a known range of colours in the TCD machine output.
- step 23 utilizes step 23 in which parameters for correlating a particular colour with an intensity value are dynamically calculated based upon an analysis of the spectral image data.
- the computer dispenses with the need for an image intensity array and instead reads the value of the various pixel's intensities directly from the image during the subsequent processing steps.
- the computer receives the intensity values of each pixel in the signal image in real time as the signal image is being produced by the TCD machine. In such embodiments the computer typically analyses the received data stream to infer a particle count in real time.
- the next step is the decomposition of the image into edges at step 24 of figure 10.
- This commences with the definition of a threshold intensity value that is to be used in the edge detection analysis.
- One preferred embodiment of the invention provides the operator with an opportunity to select a pixel from the TCD image. This is done with the use of a suitable input device such as a mouse, with which the user points and clicks on the selected image pixel.
- the intensity of the selected pixel is ascertained by the processor and the threshold intensity value is defined as equal to the ascertained value.
- the user should select a pixel having an intensity that is representative of the typical type of High Intensity Transient that the method is ultimately seeking to correlate with a particle count.
- the threshold intensity value is calculated automatically from a mathematical analysis of the intensity values of the pixels of the image. For example, in one such embodiment, the computer calculates an image intensity value that is 7db of the background noise and the threshold intensity value defined as the result.
- the computer firstly determines which type of TCD machine produced the signal image that is to be analysed. This determination may be either automated, based upon an analysis of the image, or it may be responsive to user input.
- the computer then polls a database in which pre -determined threshold values are stored for each of the known TCD machines that are commercially available.
- the database is stored locally on the computer. In yet other embodiments, the database is remotely accessible, for example via the internet.
- the processor applies an edge designation algorithm to the data stored within the image intensity array so as to designate edges exceeding the threshold intensity value.
- the edge designation algorithm is in accordance with the method described above with reference to figure 1.
- Alternative embodiments make use of prior art edge designation algorithms.
- the results of the edge designation algorithm are stored in a two dimensional edge designation array in which the values are representative of the edges.
- the size of this edge designation array substantially corresponds to the size of the original signal image. More particularly, this edge designation array has two less rows and two less columns than the original signal image, since the outer-most rows and columns of the signal image are not processible using the edge designation algorithm described above.
- the values within the edge designation array are defined within a plurality of memory addresses and/or data storage addresses.
- the values stored in this edge designation array form a Boolean data set in which the value 0 is associated with a pixel that is not designated as an edge and the value 1 is associated with a pixel that is designated as an edge.
- the state of this edge designation array may be represented graphically at step 26 by displaying an array of pixels corresponding to the edge designation array, in which edge designation array positions having a 0 value are shown in a light colour and edge designation array positions having a 1 value are shown in a dark colour. Examples of displays arising from such edge designation arrays are shown in figures 5 and 6.
- the display shown in figure 5 corresponds to the edges designated from the signal image shown in figure 3.
- the display shown in figure 6 corresponds to the edges designated from the signal image shown in figure 4.
- the edge designation array is processed so as to detect the peaks defined by the edges. As the peaks are detected, they are counted at step 28.
- a sample edge designation array of six columns by 5 rows is depicted in figure 7.
- the processor commences a transformation of the two dimensional edge designation array into a one dimensional height array by sampling the value that the edge designation array stores at the first row of the first column. This position of the edge designation array corresponds to the edge designation for the top left corner of the signal image.
- the processor increments along y coordinates, keeping x constant, until a value of 1, indicating an edge, is found.
- the example edge designation array shown in figure 7 transforms into the following one dimensional height array, with the values in the array being inversely proportional to the height of the respective edge :
- the one dimensional height array is then processed so as to detect and count the peaks. This commences with the definition of a peak counter variable that is initially set to a zero value.
- the processor indexes its way along the one dimensional height array and for each position 'n' within the one dimensional array, checks whether the value of the one dimensional array at position 'n-1' is greater than the value of the one dimensional array at position 'n' . If this relationship tests positive (i.e. if the height at the position to the left of position 'n' is lower than the height at position 'n'), then the processor checks whether the value of the one dimensional array at position 'n+1 ' is greater than the value of the one dimensional array at position 'n'.
- the example one dimensional array of ⁇ 3,1,3,2,1,3 ⁇ is processed to have two peaks (at the second and fifth positions) and the peak counter variable therefore has a value of two at the conclusion of the processing of the one dimensional height array for this example.
- the peaks that have been detected may be shown in a display at step 29.
- the receptive field is a one dimensional cellular array that is five cells long and the values from five of the positions of the one dimensional height array is provided as input to the receptive field.
- Each cell has an initial resting potential of -30V. From this resting potential, the potential of the centre cell is increased by 95 V to a value of 65V. The potential of the centre cell is then decreased by 13 V for each of the surrounding cells that has a value less than, or equal to, the value of the centre cell. If the resultant potential of the centre cell exceeds a peak detection threshold of 40V, then a peak is detected and the peak counter variable is incremented.
- the processor indexes along the one dimensional height array, applying this cellular neural network approach to each of the positions of the one dimensional height array.
- this implementation of a cellular neural network advantageously uses integral mathematics to speed processing as compared to processing requiring floating point arithmetic.
- This cellular neural network approach also advantageously tends to smoothen ragged peaks and is preferable to simply taking an average since outlying values do not skew the mean.
- the processor uses a correlation scheme responsive to the number of peaks to infer a particle count.
- the correlation scheme is a one-to-one relationship between the number of peaks and the number of particles.
- the number of particles inferred by the preferred embodiment from its analysis of the signal image is equal to the number of peaks counted.
- the correlation scheme involves multiplying the number of peaks by a factor (which may be greater or less than one) to calculate the number of inferred particles. This factor can be used to compensate for the over or under detection of HITS in the signal image.
- Figure 9 illustrates the process flow from an user's perspective.
- the user enters the patient's details at step 40 and these details are stored in the database 41.
- the user identifies the type of TCD machine which prepared the signal images being analysed. This allows the processor to access the appropriate data required to analyse the image data 20.
- the user selects the image data files that are to be analysed and at steps 44 and 45 the processor performs the processing and counting of particles. Once this processing is complete, a pointer to the source image, a pointer to the transformed image and the final particle count are stored in the database 41.
- Preferred embodiments of the invention may be implemented on a range of computing platforms.
- the preferred embodiment utilizes a computing apparatus configured to perform the various processing steps.
- This computing apparatus has a central processing unit (CPU) capable of executing software that is written in at least one of many known programming languages; associated memory, for example RAM and ROM; storage devices such as hard drives, writable CD ROMS and flash memory; input devices such as a keyboard and mouse; output devices, for example a printer; a display in the form of a screen and a communications link in the form of a modem.
- CPU central processing unit
- RAM random access memory
- storage devices such as hard drives, writable CD ROMS and flash memory
- input devices such as a keyboard and mouse
- output devices for example a printer
- the actual platform upon which the invention is implemented will vary depending upon factors such as the amount of processing power required.
- the computing apparatus is a stand alone computer, whilst in other embodiments the computing apparatus is formed from a networked array of inter
- the present invention may be embodied in computer software in the form of executable code for instructing a computer to perform the inventive method.
- the software and its associated data are capable of being stored upon a computer -readable medium in the form of one or more compact disks.
- Alternative embodiments make use of other forms of digital storage media, such as Digital Versatile Discs (DVD's), hard drives, flash memory, Erasable Programmable Read-Only Memory EPROM, and the like.
- DVD's Digital Versatile Discs
- the software and its associated data may be stored as one or more downloadable or remotely executable files that are accessible via a computer communications network such as the internet.
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Image Analysis (AREA)
Abstract
Description
Claims
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU2008215155A AU2008215155A1 (en) | 2007-02-14 | 2008-02-11 | An edge detection method and a particle counting method |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU2007900707A AU2007900707A0 (en) | 2007-02-14 | An edge detection method and a particle counting method | |
AU2007900707 | 2007-02-14 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2008098284A1 true WO2008098284A1 (en) | 2008-08-21 |
Family
ID=39689536
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/AU2008/000165 WO2008098284A1 (en) | 2007-02-14 | 2008-02-11 | An edge detection method and a particle counting method |
Country Status (2)
Country | Link |
---|---|
AU (1) | AU2008215155A1 (en) |
WO (1) | WO2008098284A1 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102171369A (en) * | 2008-10-06 | 2011-08-31 | 东曹株式会社 | Method for measuring survivin mRNA |
US8250021B2 (en) | 2008-07-17 | 2012-08-21 | Luminex Corporation | Methods, storage mediums, and systems for configuring classification regions within a classification matrix of an analysis system and for classifying particles of an assay |
US8717568B2 (en) | 2010-07-16 | 2014-05-06 | Luminex Corporation | Analyzing particle quantity and/or distribution within an assay analysis system |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2200513A (en) * | 1987-01-22 | 1988-08-03 | British Aerospace | Edge detection and enhancement of digitized images |
EP1569170A1 (en) * | 2002-12-05 | 2005-08-31 | Seiko Epson Corporation | Characteristic region extraction device, characteristic region extraction method, and characteristic region extraction program |
US7136515B2 (en) * | 2001-09-13 | 2006-11-14 | Intel Corporation | Method and apparatus for providing a binary fingerprint image |
-
2008
- 2008-02-11 WO PCT/AU2008/000165 patent/WO2008098284A1/en active Application Filing
- 2008-02-11 AU AU2008215155A patent/AU2008215155A1/en not_active Abandoned
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2200513A (en) * | 1987-01-22 | 1988-08-03 | British Aerospace | Edge detection and enhancement of digitized images |
US7136515B2 (en) * | 2001-09-13 | 2006-11-14 | Intel Corporation | Method and apparatus for providing a binary fingerprint image |
EP1569170A1 (en) * | 2002-12-05 | 2005-08-31 | Seiko Epson Corporation | Characteristic region extraction device, characteristic region extraction method, and characteristic region extraction program |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8250021B2 (en) | 2008-07-17 | 2012-08-21 | Luminex Corporation | Methods, storage mediums, and systems for configuring classification regions within a classification matrix of an analysis system and for classifying particles of an assay |
CN102171369A (en) * | 2008-10-06 | 2011-08-31 | 东曹株式会社 | Method for measuring survivin mRNA |
US8717568B2 (en) | 2010-07-16 | 2014-05-06 | Luminex Corporation | Analyzing particle quantity and/or distribution within an assay analysis system |
US9719855B2 (en) | 2010-07-16 | 2017-08-01 | Luminex Corporation | Evaluation of focusing performance in an assay analysis system |
Also Published As
Publication number | Publication date |
---|---|
AU2008215155A1 (en) | 2008-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
White et al. | Completeness in weak-lensing searches for clusters | |
CN101052991A (en) | Feature weighted medical object contouring using distance coordinates | |
DE102016100367A1 (en) | Sparse tracking in sound beam intensity impulse imaging | |
Chainais | Infinitely divisible cascades to model the statistics of natural images | |
CN102930511B (en) | Method for analyzing velocity vector of flow field of heart based on gray scale ultrasound image | |
Adams | Contour mapping and differential systematics of geographic variation | |
WO2008098284A1 (en) | An edge detection method and a particle counting method | |
Kerut et al. | Review of methods for texture analysis of myocardium from echocardiographic images: a means of tissue characterization | |
Roxborough et al. | Tetrahedron based, least squares, progressive volume models with application to freehand ultrasound data | |
JPH08201500A (en) | Underwater acoustic signal detector | |
Preston | Ξ-filters | |
Fakiris et al. | Quantification of regions of interest in swath sonar backscatter images using grey-level and shape geometry descriptors: The TargAn software | |
CN102841341A (en) | Pulse radar moving-target detection method | |
Alhadidi et al. | cDNA microarray genome image processing using fixed spot position | |
Amiri et al. | Segmentation of ultrasound images based on scatterer density using U-Net | |
DE102005051535A1 (en) | System and method for detecting frosted glass objects and nodules | |
CN113963427B (en) | Method and system for rapid in-vivo detection | |
Alipoor et al. | A novel logarithmic edge detection algorithm | |
Bons et al. | Image analysis of microstructures in natural and experimental samples | |
JPH0282947A (en) | Method for detecting and analyzing skin surface conformation | |
Hudaib et al. | New methodology for microarray spot segmentation and gene expression analysis | |
Shai et al. | Spatial bias removal in microarray images | |
Bamber et al. | Fast image processing systems for evaluating the clinical potential of ultrasound speckle suppression and parametric imaging | |
Wayalun et al. | A hybrid technique of noise reduction with periductal fibrosis ultrasound images for periductal fibrosis detection system of cholangiocarcinoma surveillance | |
Xing et al. | Simulated phantom images for optimizing wavelet-based image processing algorithms in mammography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 08706052 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2008215155 Country of ref document: AU |
|
ENP | Entry into the national phase |
Ref document number: 2008215155 Country of ref document: AU Date of ref document: 20080211 Kind code of ref document: A |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 08706052 Country of ref document: EP Kind code of ref document: A1 |