CN107133462A - Data processing method, device and equipment - Google Patents
Data processing method, device and equipment Download PDFInfo
- Publication number
- CN107133462A CN107133462A CN201710288367.1A CN201710288367A CN107133462A CN 107133462 A CN107133462 A CN 107133462A CN 201710288367 A CN201710288367 A CN 201710288367A CN 107133462 A CN107133462 A CN 107133462A
- Authority
- CN
- China
- Prior art keywords
- crystal
- energy peak
- specified
- spectrum curve
- fitting
- Prior art date
- Legal status (The legal status 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 status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
The embodiments of the invention provide a kind of data processing method, device and equipment.The embodiment of the present invention specifies the original gamma-spectrometric data of crystal by obtaining in PET device;According to original gamma-spectrometric data, the first fitting spectrum curve and the first energy peak for specifying crystal are obtained;Set up the relevance of the first fitting spectrum curve and original gamma-spectrometric data, and the crystal energy peak value for specifying crystal is obtained based on the relevance, during the energy peak for specifying crystal is obtained, reduce manual adjustment, so as to reduce the time of data processing, improve data-handling efficiency, and improve the intelligence degree of data handling procedure, solve to a certain extent more about manually being operated in the data handling procedure of energy peak in the prior art, intelligence degree is relatively low, causes the problem of data-handling efficiency is relatively low.
Description
【Technical field】
This programme is related to technical field of image processing, more particularly to a kind of data processing method, device and equipment.
【Background technology】
Currently, PET (Positron Emission Tomography, Positron Emission Computed Tomography) technology exists
Played an important role in medical science detection.PET system detects gamma rays using scintillation crystal detectors, and gamma rays lead to
Cross Compton scattering or photoelectric absorption excites in scintillator crystal materials and sends visible ray.Each crystal dissipates in generation Compton
Penetrate or during photoelectric absorption effect, certain rule is presented in photon energy distribution thereon, and non-scatter photon energy concentrates on peak value
511Kev, mainly scattered photon of the energy in below 420kev.It is right respectively using multiple data processing channels in PET system
The photon numbers of each energy of each crystal are counted, and the original energy of each crystal can be obtained based on these statistics
Modal data.
In PET system, area is come using LLD (Lower energy-Level Discriminator, low energy discriminator)
Divide scattered photon and non-scatter photon, by setting LLD threshold values to exclude scattering events., can be true by some if LLD threshold values are too high
Meet event exclusion, if LLD threshold values are too low, excessive scattering events can be made to be mixed into detection and influence PET system to detect
Accuracy.Therefore, the setting of LLD threshold values is most important for PET system.
LLD threshold values are determined by the energy peak of crystal, and the energy peak of crystal is the original energy by crystal
What modal data was obtained.In the prior art, when gamma-spectrometric data is unimodal (i.e. single peak value), LM is utilized
(Levenberg-Marquardt, Lie Wen-Burger) algorithm obtains the energy peak of crystal, is bimodal (i.e. two for gamma-spectrometric data
Individual peak value) situation, the energy peak of crystal is obtained by manual adjustment mode.So, in the case where bi-modal case is more,
The number of times of manual adjustment increases therewith, causes data-handling efficiency relatively low.
During this programme is realized, inventor has found that at least there are the following problems in the prior art:
More about manually being operated in the data handling procedure of energy peak, intelligence degree is relatively low, causes data processing
It is less efficient.
【The content of the invention】
In view of this, this programme embodiment provides a kind of data processing method, device and equipment, to solve existing skill
More about manually being operated in the data handling procedure of energy peak in art, intelligence degree is relatively low, causes data-handling efficiency
Relatively low the problem of.
In a first aspect, this programme embodiment provides a kind of data processing method, methods described includes:
Obtain the original gamma-spectrometric data that crystal is specified in positron e mission computed tomography equipment;
According to the original gamma-spectrometric data, the first fitting spectrum curve of the specified crystal and first fitting are obtained
The first energy peak in spectrum curve;
The relevance of the first fitting spectrum curve and the original gamma-spectrometric data is set up, and is obtained based on the relevance
The crystal energy peak value of the specified crystal.
Aspect as described above and any possible implementation, it is further provided a kind of implementation, set up described the
The relevance of one fitting spectrum curve and the original gamma-spectrometric data, and based on the crystal of the relevance acquisition specified crystal
Energy peak, including:
Obtain the first corresponding first function of fitting spectrum curve the second letter corresponding with the original gamma-spectrometric data
Several penalties;
The specified operation values of the penalty are compared with specified threshold, comparative result is obtained;
According to the comparative result, the crystal energy peak value of the specified crystal is determined.
Aspect as described above and any possible implementation, it is further provided a kind of implementation, according to the ratio
Relatively result, determines the crystal energy peak value of the specified crystal, including:
In the case where the specified operation values of the penalty exceed the specified threshold, expectation-maximization algorithm is utilized
Obtain the second fitting spectrum curve of the specified crystal;
The second energy peak and photon numbers row that photon numbers ranked first are extracted from the described second fitting spectrum curve
Second the 3rd energy peak;
According to specified principle, choose second energy peak or the 3rd energy peak is used as the specified crystal
Crystal energy peak value.
Aspect as described above and any possible implementation, it is further provided a kind of implementation, according to the ratio
Relatively result, determines the crystal energy peak value of the specified crystal, including:
In the case where the specified operation values of the penalty are not less than the specified threshold, first energy is determined
Peak value is the crystal energy peak value of the specified crystal.
Aspect as described above and any possible implementation, it is further provided a kind of implementation, methods described is also
Including:
Obtain the crystal energy that all crystal in crystal block are specified in the positron e mission computed tomography equipment
Measure peak value;
Pseudo- coloured silkization processing is carried out to the crystal energy peak value of all crystal in the specified crystal block, pseudo- coloured silkization figure is obtained
Picture.
Aspect as described above and any possible implementation, it is further provided a kind of implementation, methods described is also
Including:
Recognize that color-values in the pseudo- coloured silkization image meet the disordered crystal of exceptional condition;
Determine the area-of-interest of the disordered crystal position;
Obtain the median of the color-values of all crystal in the area-of-interest;
The corresponding color-values of the disordered crystal are adjusted to the median;
The crystal energy peak value of the disordered crystal is adjusted according to the corresponding color-values correspondence of the disordered crystal after adjustment.
Aspect as described above and any possible implementation, it is further provided a kind of implementation, methods described is also
Including:
When crystal in the pseudo- coloured silkization image is chosen, display is chosen the original spectrum curve of crystal, the original
Beginning spectrum curve is corresponding with original gamma-spectrometric data;
When energy peak is conditioned in the original spectrum curve of the selected crystal, the correspondence adjustment pseudo- coloured silkization figure
The color-values of crystal are chosen as described in, and correspondence adjusts the crystal energy peak value of the selected crystal.
Aspect as described above and any possible implementation, it is further provided a kind of implementation, according to the original
Beginning gamma-spectrometric data, obtains the first energy in the first fitting spectrum curve and the first fitting spectrum curve of the specified crystal
Measure before peak value, methods described also includes:
Filtration treatment, the original gamma-spectrometric data after being filtered are carried out to the original gamma-spectrometric data;
According to the original gamma-spectrometric data, the first fitting spectrum curve of the specified crystal and first fitting are obtained
The first energy peak in spectrum curve, including:
According to the original gamma-spectrometric data after filtering, the first fitting spectrum curve and described first of the specified crystal is obtained
It is fitted the first energy peak in spectrum curve.
Aspect as described above and any possible implementation, it is further provided a kind of implementation, methods described is also
Including:
According to the crystal energy peak value of the specified crystal, the positron e mission computed tomography equipment is determined
In corresponding low energy discriminator threshold value.
Second aspect, this programme embodiment provides a kind of data processing equipment, and described device includes:
Acquisition module, for obtaining positron e mission computed tomography positron e mission computed tomography
The original gamma-spectrometric data of crystal is specified in equipment;
First fitting module, for according to the original gamma-spectrometric data, obtaining the first fitting power spectrum of the specified crystal
The first energy peak in curve and the first fitting spectrum curve;
Module is set up, the relevance for setting up the first fitting spectrum curve and the original gamma-spectrometric data, and base
The crystal energy peak value of the specified crystal is obtained in the relevance.
The third aspect, this programme embodiment provides a kind of data processing equipment, and the equipment includes:
Processor;
Memory for storing the processor-executable instruction;
The processor is configured as:
Obtain the original gamma-spectrometric data that crystal is specified in positron e mission computed tomography equipment;
According to the original gamma-spectrometric data, the first fitting spectrum curve of the specified crystal and first fitting are obtained
The first energy peak in spectrum curve;
The relevance of the first fitting spectrum curve and the original gamma-spectrometric data is set up, and is obtained based on the relevance
The crystal energy peak value of the specified crystal.
The embodiment of the present invention has the advantages that:
The embodiment of the present invention is obtained according to being fitted by being fitted processing to the original gamma-spectrometric data for specifying crystal
The relevance of first fitting spectrum curve and original gamma-spectrometric data, to obtain the crystal energy peak value of specified crystal so that obtaining
During the energy peak that crystal must be specified, manual adjustment is reduced, so as to reduce the time of data processing, number is improved
According to treatment effeciency, and improve the intelligence degree of data handling procedure.
【Brief description of the drawings】
, below will be attached to what is used required in embodiment in order to illustrate more clearly of the technical scheme of this programme embodiment
Figure is briefly described, it should be apparent that, drawings in the following description are only some embodiments of this programme, for this area
For those of ordinary skill, without having to pay creative labor, it can also be obtained according to these accompanying drawings other attached
Figure.
Fig. 1 is the first pass exemplary plot of data processing method provided in an embodiment of the present invention.
Fig. 2 is the exemplary plot of spectrum curve provided in an embodiment of the present invention.
Fig. 3 is the second procedure exemplary plot of data processing method provided in an embodiment of the present invention.
Fig. 4 is the 3rd flow example figure of data processing method provided in an embodiment of the present invention.
Fig. 5 is the 4th flow example figure of data processing method provided in an embodiment of the present invention.
Fig. 6 is the 5th flow example figure of data processing method provided in an embodiment of the present invention.
Fig. 7 is the 6th flow example figure of data processing method provided in an embodiment of the present invention.
Fig. 8 is the 7th flow example figure of data processing method provided in an embodiment of the present invention.
Fig. 9 is the functional block diagram of data processing equipment provided in an embodiment of the present invention.
Figure 10 is the simplified block diagram of data processing equipment.
【Embodiment】
In order to be better understood from the technical scheme of this programme, this programme embodiment is retouched in detail below in conjunction with the accompanying drawings
State.
It will be appreciated that described embodiment is only a part of embodiment of this programme, rather than whole embodiments.Base
Embodiment in this programme, those of ordinary skill in the art obtained under the premise of creative work is not made it is all its
Its embodiment, belongs to the scope of this programme protection.
The term used in this programme embodiment is the purpose only merely for description specific embodiment, and is not intended to be limiting
This programme." one kind ", " described " and "the" of singulative used in this programme embodiment and appended claims
It is also intended to including most forms, unless context clearly shows that other implications.
It should be appreciated that term "and/or" used herein is only a kind of incidence relation for describing affiliated partner, represent
There may be three kinds of relations, for example, A and/or B, can be represented:Individualism A, while there is A and B, individualism B these three
Situation.In addition, character "/" herein, it is a kind of relation of "or" to typically represent forward-backward correlation object.
Depending on linguistic context, word as used in this " if " can be construed to " ... when " or " when ...
When " or " in response to determining " or " in response to detection ".Similarly, depending on linguistic context, phrase " if it is determined that " or " if detection
(condition or event of statement) " can be construed to " when it is determined that when " or " in response to determine " or " when the detection (condition of statement
Or event) when " or " in response to detection (condition or event of statement) ".
Embodiment one
The embodiments of the invention provide a kind of data processing method.The data processing method can pass through application APP
To realize, the terminal such as computer, medical work station can obtain corresponding data processing function by installing the application program.
Fig. 1 is the first pass exemplary plot of data processing method provided in an embodiment of the present invention.As shown in figure 1, this implementation
In example, data processing method may include steps of:
S101, obtains the original gamma-spectrometric data that crystal is specified in PET device.
S102, according to original gamma-spectrometric data, obtaining the first fitting spectrum curve for specifying crystal and the first fitting can set a song to music
The first energy peak in line.
S103, sets up the relevance of the first fitting spectrum curve and original gamma-spectrometric data, and referred to based on the relevance
Determine the crystal energy peak value of crystal.
, can be by the original gamma-spectrometric data first storage of collection after original gamma-spectrometric data is collected from PET device
In equipment.Then when needing, original gamma-spectrometric data is obtained from storage device.
Wherein, original gamma-spectrometric data is the gamma-spectrometric data of specified crystal in PET device.Original gamma-spectrometric data is one group discrete
The photon numbers that data, such as energy value are a are m, and the photon numbers that energy value is b are n.
Original gamma-spectrometric data can be converted to corresponding original spectrum curve.Original spectrum curve can be by by original energy
The corresponding discrete point of modal data is carried out in smoothing junction and obtained according to the ascending order of energy value.The exemplary plot of spectrum curve
Refer to Fig. 2.In Fig. 2, abscissa represents the corresponding photon energy values of crystal, and ordinate represents the quantity of photon.Wherein power spectrum
The corresponding abscissa of peak of curve is the energy peak of crystal.
Wherein, S102, according to original gamma-spectrometric data, obtains the first fitting spectrum curve for specifying crystal and the first fitting energy
The first energy peak in spectral curve, can include:According to original gamma-spectrometric data, obtained using LM algorithms and specify the first of crystal
It is fitted the first energy peak in spectrum curve and the first fitting spectrum curve.
The LM algorithms in S102 are illustrated below.
It is distributed and is studied according to crystal energy, most PET crystal energies is distributed in the range of 0~511Kev, Mei Geneng
The distribution between counting in amount meets mixed type Gaussian curve model equation below (1):
Go out specified crystal corresponding parameter A, b, μ, α, β, γ match value using LM algorithm the Fitting Calculations, wherein, μ is finger
Determine the corresponding energy peak of crystal, you can determined to specify the corresponding first fitting spectrum curve of crystal by these parameters, this first
The energy peak being fitted in spectrum curve is the first energy peak, i.e. μ.
Assuming that the first corresponding first function of fitting spectrum curve is function f (x), corresponding second letter of original spectrum curve
Number is function t (x), then penalty=| f (x)-t (x) |.
Wherein, specified threshold rule of thumb data can be set.Wherein, the specified operation values of penalty can be punished
The quadratic sum of penalty function.For example, it is assumed that penalty=| f (x)-t (x) |, specified threshold is k, then k can be met | f (x)-
T (x) | quadratic sum be less than k value.
When gamma-spectrometric data is unimodal, the penalty of relatively becoming of usual penalty and specified threshold does not surpass
Specified threshold is crossed, can now be determined to specify the crystal energy peak value of crystal according to the first energy peak obtained in S102.
When gamma-spectrometric data is bimodal, the penalty of relatively becoming of usual penalty and specified threshold exceedes specified threshold
Value, can now determine to specify the crystalline substance of crystal according to EM (Expectation-maximization, expectation maximization) algorithms
Body energy peak.
During a concrete implementation, according to original gamma-spectrometric data, obtain and specify the first fitting of crystal to set a song to music
Before the first energy peak in line and the first fitting spectrum curve, data processing method can also include:To original power spectrum number
According to progress filtration treatment, the original gamma-spectrometric data after being filtered;According to original gamma-spectrometric data, obtain and specify the first of crystal to intend
The first energy peak in spectrum curve and the first fitting spectrum curve is closed, including:According to the original gamma-spectrometric data after filtering, obtain
Fetching determines the first energy peak in the first fitting spectrum curve and the first fitting spectrum curve of crystal.
Wherein, to the purpose that original gamma-spectrometric data carries out filtration treatment removed in gathered data due to mutation or burr
Etc. the mushing error caused, so as to improve the accuracy of data processing.
Embodiment illustrated in fig. 1, by being fitted processing to the original gamma-spectrometric data for specifying crystal, and is obtained according to fitting
The first fitting spectrum curve and original gamma-spectrometric data relevance, to obtain the crystal energy peak value of specified crystal so that
During obtaining the energy peak for specifying crystal, manual adjustment is reduced, so as to reduce the time of data processing, is improved
Data-handling efficiency, and improve the intelligence degree of data handling procedure.
Fig. 3 is the second procedure exemplary plot of data processing method provided in an embodiment of the present invention.As shown in figure 3, this implementation
In example, data processing method may include steps of:
S301, obtains the original gamma-spectrometric data that crystal is specified in PET device.
S302, according to original gamma-spectrometric data, obtaining the first fitting spectrum curve for specifying crystal and the first fitting can set a song to music
The first energy peak in line.
S303, obtains first and is fitted the corresponding first function of spectrum curve second function corresponding with original gamma-spectrometric data
Penalty.
S304, specified operation values and the specified threshold of penalty are compared, comparative result is obtained.
S305, according to comparative result, it is determined that specifying the crystal energy peak value of crystal.
In S305, it is not necessary to manual adjustment, and according to the comparative result of penalty and specified threshold, with regard to that can automatically determine
Specify the crystal energy peak value of crystal.So, the whole process for obtaining the crystal energy peak value for specifying crystal can be according to journey
Sequence is carried out automatically, greatly reduces manual adjustment, so as to can not only improve the intelligence degree of data handling procedure, additionally it is possible to
Processing time is reduced, the efficiency of data processing is improved.
The corresponding spectrum curve bi-modal case of Crystal Spectrum data is more, embodiment illustrated in fig. 3 reduction manual adjustment
Time it is more so that the effect for improving data-handling efficiency is more obvious.
Fig. 4 is the 3rd flow example figure of data processing method provided in an embodiment of the present invention.As shown in figure 4, this implementation
In example, data processing method may include steps of:
S401, obtains the original gamma-spectrometric data that crystal is specified in PET device.
S402, according to original gamma-spectrometric data, obtaining the first fitting spectrum curve for specifying crystal and the first fitting can set a song to music
The first energy peak in line.
S403, obtains first and is fitted the corresponding first function of spectrum curve second function corresponding with original gamma-spectrometric data
Penalty.
S404, specified operation values and the specified threshold of penalty are compared, comparative result is obtained.
S405, judges whether the specified operation values of penalty exceed specified threshold, if the specified computing of penalty
Value exceedes specified threshold, performs S406, if the specified operation values of penalty are not less than specified threshold, performs S409.
S406, the second fitting spectrum curve for specifying crystal is obtained using EM algorithms;
S407, extracts the second energy peak and photon numbers row that photon numbers ranked first from the second fitting spectrum curve
Second the 3rd energy peak.
S408, according to specified principle, chooses the second energy peak or the 3rd energy peak as the crystal energy of specified crystal
Peak value is measured, is terminated.
S409, it is the crystal energy peak value for specifying crystal to determine the first energy peak, is terminated.
The Energy distribution of PET crystal occurs that bimodal (relatively low energy quantitative value peaks, with respect to the quantity of high-energy
Be worth relatively low peak) situation meet mixed Gauss model, now, the spectrum curve of each crystal can use formula below (2) table
Show.
Using EM algorithm the Fitting Calculations go out corresponding parameter A1, b1, the μ 1 of specified crystal, α 1, β 1, γ 1, A2, b2, μ 2, α 2,
β 2, γ 2 match value, wherein, μ 1 and μ 2 are to specify corresponding two energy peaks of crystal, you can is determined by these parameters specified
Two energy peaks in the corresponding second fitting spectrum curve of crystal, the second fitting spectrum curve are respectively μ 1 and μ 2.
In embodiment illustrated in fig. 4, in the case where penalty is not less than specified threshold, will obtain (can be calculated using LM
Method obtain) the first energy peak be defined as specify crystal crystal energy peak value, penalty exceed specified threshold feelings
Under condition, the second energy peak in the second fitting spectrum curve for specifying crystal and the second fitting spectrum curve is obtained using EM algorithms
Value and the 3rd energy peak, and according to specified principle, choose the second energy peak or the 3rd energy peak is used as specified crystal
Crystal energy peak value, so as to occur bimodal situation even in spectrum curve, can also automatically obtain the crystal of crystal according to program
Energy peak, greatly reduces manual adjustment, it is to avoid cause data processing time long because of excessive manual intervention, reduces
The time of data processing, and then improve the efficiency of data processing.
Fig. 5 is the 4th flow example figure of data processing method provided in an embodiment of the present invention.As shown in figure 5, this implementation
In example, data processing method may include steps of:
S501, obtains the original gamma-spectrometric data that crystal is specified in PET device.
S502, according to original gamma-spectrometric data, obtaining the first fitting spectrum curve for specifying crystal and the first fitting can set a song to music
The first energy peak in line.
S503, obtains first and is fitted the corresponding first function of spectrum curve second function corresponding with original gamma-spectrometric data
Penalty.
S504, specified operation values and the specified threshold of penalty are compared, comparative result is obtained.
S505, according to comparative result, it is determined that specifying the crystal energy peak value of crystal.
S506, obtains the crystal energy peak value that all crystal in crystal block are specified in PET device;
S507, carries out pseudo- coloured silkization to the crystal energy peak value for specifying all crystal in crystal block and handles, obtain pseudo- coloured silkization figure
Picture.
In embodiment illustrated in fig. 5, after it is determined that specifying the crystal energy peak value of crystal, own always according in specified crystal block
The crystal energy peak value of crystal, is handled by pseudo- coloured silkization, obtains pseudo- coloured silkization image.In pseudo- coloured silkization image, the color of each crystal
It is worth for rgb value, the value is corresponding with the crystal energy peak value of crystal.By pseudo- coloured silkization image, user can easily pass through
Eyes intuitively identify the crystal of image color value lofty (color-values are abnormal in other words), so as to the crystal energy of these crystal
Amount peak value is suitably adjusted, so that the crystal energy peak value of these crystal more meets actual value.
It should be noted that in data processing, whole crystal of PET device can be carried out into piecemeal first, so
Each piece of crystal is respectively processed afterwards, thus polylith crystal data can parallel be located simultaneously by multiple threads
Reason, so as to improve data processing speed.For example, it is assumed that whole crystal of PET device have 1152*96 (symbol " * " they are multiplication sign),
These crystal are divided into 36 pieces, each piece of crystal number is 64*48, then using multithreading, and 36 crystal blocks are carried out
Parallel processing.For the data volume very complicated data of big and data structure, the efficiency of data processing can be so significantly improved.
Fig. 6 is the 5th flow example figure of data processing method provided in an embodiment of the present invention.As shown in fig. 6, this implementation
In example, data processing method may include steps of:
S601, obtains the original gamma-spectrometric data that crystal is specified in PET device.
S602, according to original gamma-spectrometric data, obtaining the first fitting spectrum curve for specifying crystal and the first fitting can set a song to music
The first energy peak in line.
S603, obtains first and is fitted the corresponding first function of spectrum curve second function corresponding with original gamma-spectrometric data
Penalty.
S604, specified operation values and the specified threshold of penalty are compared, comparative result is obtained.
S605, according to comparative result, it is determined that specifying the crystal energy peak value of crystal.
S606, obtains the crystal energy peak value that all crystal in crystal block are specified in PET device;
S607, carries out pseudo- coloured silkization to the crystal energy peak value for specifying all crystal in crystal block and handles, obtain pseudo- coloured silkization figure
Picture.
Color-values meet the disordered crystal of exceptional condition in S608, the pseudo- coloured silkization image of identification;
S609, determines the area-of-interest of disordered crystal position;
S610, obtains the median of the color-values of all crystal in area-of-interest;
The corresponding color-values of disordered crystal are adjusted to median by S611;
S612, the crystal energy peak value of disordered crystal is adjusted according to the corresponding color-values correspondence of disordered crystal after adjustment.
Embodiment illustrated in fig. 6 provides the mode that a kind of crystal energy peak value to crystal is adjusted.
, will after color-values meet the disordered crystal of exceptional condition in identifying pseudo- coloured silkization image in embodiment illustrated in fig. 6
The corresponding color-values of disordered crystal are adjusted to the median of the color-values of all crystal of related area-of-interest, and according to adjustment
The corresponding color-values correspondence of disordered crystal adjusts the crystal energy peak value of disordered crystal afterwards, so as to the crystalline substance of obvious abnormal crystal
Body energy peak is adjusted, and it is more conformed to actual value, so as to improve the accuracy of data processing.
Fig. 7 is the 6th flow example figure of data processing method provided in an embodiment of the present invention.As shown in fig. 7, this implementation
In example, data processing method may include steps of:
The original gamma-spectrometric data of crystal is specified in S701, collection PET device.
S702, according to original gamma-spectrometric data, obtaining the first fitting spectrum curve for specifying crystal and the first fitting can set a song to music
The first energy peak in line.
S703, obtains first and is fitted the corresponding first function of spectrum curve second function corresponding with original gamma-spectrometric data
Penalty.
S704, specified operation values and the specified threshold of penalty are compared, comparative result is obtained.
S705, according to comparative result, it is determined that specifying the crystal energy peak value of crystal.
S706, obtains the crystal energy peak value that all crystal in crystal block are specified in PET device;
S707, carries out pseudo- coloured silkization to the crystal energy peak value for specifying all crystal in crystal block and handles, obtain pseudo- coloured silkization figure
Picture.
S708, when the crystal in pseudo- coloured silkization image is chosen, display is chosen the original spectrum curve of crystal, wherein,
The original spectrum curve for being chosen crystal is corresponding with the original gamma-spectrometric data of selected crystal.
S709, when energy peak is conditioned in the original spectrum curve for being chosen crystal, the pseudo- coloured silkization image of correspondence adjustment
In be chosen crystal color-values, and correspondence adjustment be chosen crystal crystal energy peak value.
Embodiment illustrated in fig. 7 provides the mode that another crystal energy peak value to crystal is adjusted.
In embodiment illustrated in fig. 7, after obtaining specifying the corresponding pseudo- coloured silkization image of crystal block, the crystalline substance in pseudo- coloured silkization image
When body is chosen, the original spectrum curve of selected crystal can be shown.On this basis, when original of the user to selected crystal
After energy peak in beginning spectrum curve is adjusted, it can be chosen according in the pseudo- coloured silkization image of the regulation of user correspondence adjustment
The color-values of crystal and the crystal energy peak value of crystal.
Fig. 8 is the 7th flow example figure of data processing method provided in an embodiment of the present invention.As shown in figure 8, this implementation
In example, data processing method may include steps of:
S801, obtains the original gamma-spectrometric data that crystal is specified in PET device.
S802, according to original gamma-spectrometric data, obtaining the first fitting spectrum curve for specifying crystal and the first fitting can set a song to music
The first energy peak in line.
S803, obtains first and is fitted the corresponding first function of spectrum curve second function corresponding with original gamma-spectrometric data
Penalty.
S804, specified operation values and the specified threshold of penalty are compared, comparative result is obtained.
S805, according to comparative result, it is determined that specifying the crystal energy peak value of crystal.
S806, according to the crystal energy peak value of specified crystal, determines the corresponding low energy discriminator LLD thresholds in PET device
Value.
In embodiment illustrated in fig. 8, after it is determined that specifying the crystal energy peak value of crystal, determined according to crystal energy peak value
Corresponding low energy discriminator LLD threshold values in PET device so that the setting of LLD threshold values is more reasonable, so as to improve PET
The accuracy of the testing result of equipment.
Data processing method provided in an embodiment of the present invention, by being fitted place to the original gamma-spectrometric data for specifying crystal
Reason, and the first fitting spectrum curve and the specified operation values of the penalty of original spectrum curve that are obtained according to fitting are with specifying
The comparative result of threshold value, to determine the crystal energy peak value for specifying crystal so that obtaining the mistake for the energy peak for specifying crystal
Cheng Zhong, reduces manual adjustment, so as to reduce the time of data processing, improves data-handling efficiency, and improve number
According to the intelligence degree of processing procedure.
Embodiment two
The embodiments of the invention provide a kind of data processing equipment, the data processing equipment can be realized in previous embodiment
Each step of data processing method.
Fig. 9 is the functional block diagram of data processing equipment provided in an embodiment of the present invention.As shown in figure 9, in the present embodiment,
Data processing equipment includes:
Acquisition module 910, the original of crystal is specified for obtaining in positron e mission computed tomography PET device
Gamma-spectrometric data;
First fitting module 920, specified for according to original gamma-spectrometric data, obtaining crystal the first fitting spectrum curve and
The first energy peak in first fitting spectrum curve;
Module 930 is set up, the relevance for setting up the first fitting spectrum curve and original gamma-spectrometric data, and based on the pass
Connection property obtains the crystal energy peak value for specifying crystal.
During a concrete implementation, module 930 is set up for setting up the first fitting spectrum curve and original energy
The relevance of modal data, and when obtaining the crystal energy peak value for specifying crystal based on the relevance, specifically for:First is obtained to intend
Close the penalty of the corresponding first function of spectrum curve second function corresponding with the original gamma-spectrometric data;By penalty
Specified operation values be compared with specified threshold, obtain comparative result;According to comparative result, it is determined that specifying the crystal energy of crystal
Measure peak value.
During a concrete implementation, module 930 is set up for according to comparative result, it is determined that specifying the crystalline substance of crystal
During body energy peak, specifically for:In the case where the specified operation values of penalty exceed specified threshold, obtained using EM algorithms
Fetching determines the second fitting spectrum curve of crystal;The second energy that photon numbers ranked first is extracted from the second fitting spectrum curve
The 3rd energy peak that peak value and photon numbers ranked second;According to specified principle, the second energy peak or the 3rd energy peak are chosen
It is worth as the crystal energy peak value for specifying crystal.
During a concrete implementation, module 930 is set up for according to comparative result, it is determined that specifying the crystalline substance of crystal
During body energy peak, specifically for:In the case where the specified operation values of penalty are not less than specified threshold, the first energy is determined
It is the crystal energy peak value for specifying crystal to measure peak value.
During a concrete implementation, data processing equipment can also include:Peak value acquisition module, for obtaining
The crystal energy peak value of all crystal in crystal block is specified in PET device;Pseudo- coloured silkization processing module, for specifying in crystal block
The crystal energy peak value of all crystal carries out pseudo- coloured silkization processing, obtains pseudo- coloured silkization image.
During a concrete implementation, include peak value acquisition module and pseudo- coloured silkization processing module in data processing equipment
On the basis of, data processing equipment can also include:
Identification module, for recognizing that color-values in pseudo- coloured silkization image meet the disordered crystal of exceptional condition;
Area-of-interest determining module, the area-of-interest for determining the disordered crystal position;
Median acquisition module, the median for obtaining the color-values of all crystal in area-of-interest;
Color-values adjusting module, for the corresponding color-values of disordered crystal to be adjusted into median;
Peak value adjusting module, the crystal for adjusting disordered crystal according to the corresponding color-values correspondence of disordered crystal after adjustment
Energy peak.
During a concrete implementation, include peak value acquisition module and pseudo- coloured silkization processing module in data processing equipment
On the basis of, data processing equipment can also include:
Curve display module, when being chosen for the crystal in pseudo- coloured silkization image, display is chosen the original energy of crystal
Spectral curve, wherein, the original spectrum curve for being chosen crystal is corresponding with the original gamma-spectrometric data of selected crystal;
Color-values and peak value adjusting module, for be chosen crystal original spectrum curve in energy peak be conditioned
When, the color-values of crystal are chosen described in the pseudo- coloured silkization image of correspondence adjustment, and correspondence adjustment is chosen the crystal energy of crystal
Measure peak value.
During a concrete implementation, data processing equipment can also include:Filtering module, for original power spectrum
Data carry out filtration treatment, the original gamma-spectrometric data after being filtered;First fitting module 920 is for according to original power spectrum number
According to specific to use when obtaining the first energy peak in the first fitting spectrum curve for specifying crystal and the first fitting spectrum curve
In:According to the original gamma-spectrometric data after filtering, the first fitting spectrum curve for specifying crystal and the first fitting spectrum curve are obtained
In the first energy peak.
During a concrete implementation, data processing equipment can also include:LLD threshold determination modules, for root
According to the crystal energy peak value of specified crystal, the corresponding low energy discriminator LLD threshold values in PET device are determined.
Because the data processing equipment in the present embodiment is able to carry out the data processing method in previous embodiment one, this reality
The part that example is not described in detail is applied, the related description to data processing method in previous embodiment one is referred to.
Data processing equipment provided in an embodiment of the present invention, by being fitted place to the original gamma-spectrometric data for specifying crystal
Reason, and the relevance of the first fitting spectrum curve and original gamma-spectrometric data obtained according to fitting, to obtain the crystalline substance of specified crystal
Body energy peak so that during the energy peak for specifying crystal is obtained, manual adjustment is reduced, so as to reduce data
The time of processing, data-handling efficiency is improved, and improve the intelligence degree of data handling procedure.
Embodiment three
The embodiment of the present invention provides a kind of data processing equipment, and the equipment includes:Processor;It can be held for storing processor
The memory of row instruction;Processor is configured as:Obtain and crystal is specified in positron e mission computed tomography PET device
Original gamma-spectrometric data;According to original gamma-spectrometric data, the first fitting spectrum curve for specifying crystal and the first fitting power spectrum are obtained
The first energy peak in curve;The relevance of the first fitting spectrum curve and original gamma-spectrometric data is set up, and based on the association
Property obtain specify crystal crystal energy peak value.
Wherein, data processing equipment can be computer.
Figure 10 is the simplified block diagram of data processing equipment.Referring to Figure 10, the data processing equipment 1000 can include and one
The processor 1001 of individual or multiple data storage facility connections, the data storage facility can include storage medium 1006 and internal memory
Unit 1004.Data processing equipment 1000 can also include input interface 1005 and output interface 1007, for another device
Or system is communicated.Internal storage location 1004 or storage medium are storable in by the CPU of processor 1001 program codes performed
In 1006.
Processor 1001 in data processing equipment 1000, which is called, is stored in internal storage location 1004 or storage medium 1006
Program code, performs following each step:
Obtain the original gamma-spectrometric data that crystal is specified in positron e mission computed tomography PET device;
According to original gamma-spectrometric data, obtain in the first fitting spectrum curve for specifying crystal and the first fitting spectrum curve
First energy peak;
The relevance of the first fitting spectrum curve and original gamma-spectrometric data is set up, and specified crystal is obtained based on the relevance
Crystal energy peak value.
In above-described embodiment, storage medium can be read-only storage (Read-Only Memory, ROM), or readable
Write, for example hard disk, flash memory.Internal storage location can be random access memory (Random Access Memory, RAM).Internal memory
Unit can be with processor physical integration or integrated in memory or being configured to single unit.
Processor is the control centre of the said equipment (equipment is above-mentioned server or above-mentioned client), and at offer
Device is managed, for execute instruction, carrying out interrupt operation, there is provided clocking capability and various other functions.Alternatively, processor bag
Include the CPU 0 and CPU 1 shown in one or more CPU (CPU), such as Figure 10.The said equipment includes one
Or multiple processors.Processor can be monokaryon (single CPU) processor or multinuclear (multi -CPU) processor.Unless otherwise stated,
It is described as the part of such as processor or memory for performing task and can realize as universal component, it is temporarily used for given
Time performs task, or is embodied as being manufactured specifically for the particular elements of the execution task.Terminology used herein " processor "
Refer to one or more devices, circuit and/or process cores, for processing data, such as computer program instructions.
It is storable in by the CPU of the processor program codes performed in internal storage location or storage medium.Alternatively, it is stored in
Program code in storage medium can be copied into internal storage location and be performed so as to the CPU of processor.Processor is executable at least
One kernel (such as LINUXTM、UNIXTM、WINDOWSTM、ANDROIDTM、IOSTM), it is well known that the kernel is used for by control
The communication of execution, control and the peripheral unit of other programs or process and the use of control computer device resource is made to control
The operation of the said equipment.
Said elements in the said equipment can be connected to each other by bus, bus such as data/address bus, address bus, control
One of bus, expansion bus and local bus or its any combination.
It is apparent to those skilled in the art that, for convenience and simplicity of description, the system of foregoing description,
The specific work process of device and unit, may be referred to the corresponding process in preceding method embodiment, will not be repeated here.
In several embodiments that this programme is provided, it should be understood that disclosed system, apparatus and method can be with
Realize by another way.For example, device embodiment described above is only schematical, for example, the unit
Divide, only a kind of division of logic function there can be other dividing mode when actually realizing, for example, multiple units or group
Part can combine or be desirably integrated into another system, or some features can be ignored, or not perform.It is another, it is shown
Or the coupling each other discussed or direct-coupling or communication connection can be by some interfaces, device or unit it is indirect
Coupling is communicated to connect, and can be electrical, machinery or other forms.
The unit illustrated as separating component can be or may not be it is physically separate, it is aobvious as unit
The part shown can be or may not be physical location, you can with positioned at a place, or can also be distributed to multiple
On NE.Some or all of unit therein can be selected to realize the mesh of this embodiment scheme according to the actual needs
's.
In addition, each functional unit in this programme each embodiment can be integrated in a processing unit, can also
That unit is individually physically present, can also two or more units it is integrated in a unit.Above-mentioned integrated list
Member can both be realized in the form of hardware, it would however also be possible to employ hardware adds the form of SFU software functional unit to realize.
The preferred embodiment of this programme is the foregoing is only, not to limit this programme, all essences in this programme
God is with principle, and any modification, equivalent substitution and improvements done etc. should be included within the scope of this programme protection.
Claims (10)
1. a kind of data processing method, it is characterised in that methods described includes:
Obtain the original gamma-spectrometric data that crystal is specified in positron e mission computed tomography equipment;
According to the original gamma-spectrometric data, the first fitting spectrum curve and the first fitting power spectrum of the specified crystal are obtained
The first energy peak in curve;
The relevance of the first fitting spectrum curve and the original gamma-spectrometric data is set up, and obtains described based on the relevance
Specify the crystal energy peak value of crystal.
2. according to the method described in claim 1, it is characterised in that described to set up the first fitting spectrum curve and the original
The relevance of beginning gamma-spectrometric data, and included based on the crystal energy peak value that the relevance obtains the specified crystal:
Obtain the first corresponding first function of fitting spectrum curve second function corresponding with the original gamma-spectrometric data
Penalty;
The specified operation values of the penalty are compared with specified threshold, comparative result is obtained;
According to the comparative result, the crystal energy peak value of the specified crystal is determined.
3. method according to claim 2, it is characterised in that according to the comparative result, determine the specified crystal
Crystal energy peak value, including:
In the case where the specified operation values of the penalty exceed the specified threshold, obtained using expectation-maximization algorithm
Second fitting spectrum curve of the specified crystal;
The second energy peak and photon numbers that extraction photon numbers ranked first from the described second fitting spectrum curve ranked second
The 3rd energy peak;
According to specified principle, second energy peak or the 3rd energy peak are chosen as the crystal of the specified crystal
Energy peak.
4. method according to claim 2, it is characterised in that according to the comparative result, determine the specified crystal
Crystal energy peak value, including:
In the case where the specified operation values of the penalty are not less than the specified threshold, first energy peak is determined
For the crystal energy peak value of the specified crystal.
5. method according to claim 2, it is characterised in that methods described also includes:
Obtain the crystal energy peak that all crystal in crystal block are specified in the positron e mission computed tomography equipment
Value;
Pseudo- coloured silkization processing is carried out to the crystal energy peak value of all crystal in the specified crystal block, pseudo- coloured silkization image is obtained.
6. method according to claim 5, it is characterised in that methods described also includes:
Recognize that color-values in the pseudo- coloured silkization image meet the disordered crystal of exceptional condition;
Determine the area-of-interest of the disordered crystal position;
Obtain the median of the color-values of all crystal in the area-of-interest;
The corresponding color-values of the disordered crystal are adjusted to the median;
The crystal energy peak value of the disordered crystal is adjusted according to the corresponding color-values correspondence of the disordered crystal after adjustment.
7. method according to claim 5, it is characterised in that methods described also includes:
When crystal in the pseudo- coloured silkization image is chosen, display is chosen the original spectrum curve of crystal, the original energy
Spectral curve is corresponding with original gamma-spectrometric data;
When energy peak is conditioned in the original spectrum curve of the selected crystal, in the correspondence adjustment pseudo- coloured silkization image
The color-values of the selected crystal, and correspondence adjust the crystal energy peak value of the selected crystal.
8. method according to claim 2, it is characterised in that methods described also includes:
According to the crystal energy peak value of the specified crystal, determine in the positron e mission computed tomography equipment
Corresponding low energy discriminator threshold value.
9. a kind of data processing equipment, it is characterised in that described device includes:
Acquisition module, the original gamma-spectrometric data of crystal is specified for gathering in positron e mission computed tomography equipment;
First fitting module, for according to the original gamma-spectrometric data, obtaining the first fitting spectrum curve of the specified crystal
With the first energy peak in the described first fitting spectrum curve;
Module is set up, for setting up the relevance of the first fitting spectrum curve and the original gamma-spectrometric data, and is based on being somebody's turn to do
Relevance obtains the crystal energy peak value of the specified crystal.
10. a kind of data processing equipment, it is characterised in that the equipment includes:
Processor;
Memory for storing the processor-executable instruction;
The processor is configured as:
Gather the original gamma-spectrometric data that crystal is specified in positron e mission computed tomography equipment;
According to the original gamma-spectrometric data, the first fitting spectrum curve and the first fitting power spectrum of the specified crystal are obtained
The first energy peak in curve;
The relevance of the first fitting spectrum curve and the original gamma-spectrometric data is set up, and obtains described based on the relevance
Specify the crystal energy peak value of crystal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710288367.1A CN107133462B (en) | 2017-04-27 | 2017-04-27 | Data processing method, device and equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710288367.1A CN107133462B (en) | 2017-04-27 | 2017-04-27 | Data processing method, device and equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107133462A true CN107133462A (en) | 2017-09-05 |
CN107133462B CN107133462B (en) | 2020-08-18 |
Family
ID=59716206
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710288367.1A Active CN107133462B (en) | 2017-04-27 | 2017-04-27 | Data processing method, device and equipment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107133462B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109949953A (en) * | 2019-04-02 | 2019-06-28 | 中国科学院高能物理研究所 | Nuclear fuel rod degree of enrichment passive detection system and degree of enrichment reconstructing method |
CN110148991A (en) * | 2019-06-17 | 2019-08-20 | 北京理工大学 | A kind of over-charging of battery method for early warning and system based on big data |
CN110211095A (en) * | 2019-05-06 | 2019-09-06 | 东软医疗***股份有限公司 | Energy diagram and crystal positions search table generating method, device, storage medium |
CN111199081A (en) * | 2018-10-31 | 2020-05-26 | 南京智能仿真技术研究院有限公司 | Energy spectrum data processing device and simulation system and method for simulating energy spectrum curve thereof |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6624422B2 (en) * | 2001-09-25 | 2003-09-23 | Ge Medical Systems Global Technology Company, Llc | Method for dynamic stabilization of PET detector gains |
CN104035122A (en) * | 2014-05-22 | 2014-09-10 | 沈阳东软医疗***有限公司 | Method and device for correcting energy value |
CN104173074A (en) * | 2014-07-21 | 2014-12-03 | 北京辛耕普华医疗科技有限公司 | Geometric calibration method for SPECT (single photon emission computed tomography) system |
CN104599302A (en) * | 2015-01-13 | 2015-05-06 | 上海联影医疗科技有限公司 | Method for obtaining PET crystal energy peak value and setting energy frequency discriminator |
CN105115994A (en) * | 2015-07-22 | 2015-12-02 | 武汉数字派特科技有限公司 | Digital PET energy parameterization calibration method and system |
-
2017
- 2017-04-27 CN CN201710288367.1A patent/CN107133462B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6624422B2 (en) * | 2001-09-25 | 2003-09-23 | Ge Medical Systems Global Technology Company, Llc | Method for dynamic stabilization of PET detector gains |
CN104035122A (en) * | 2014-05-22 | 2014-09-10 | 沈阳东软医疗***有限公司 | Method and device for correcting energy value |
CN104173074A (en) * | 2014-07-21 | 2014-12-03 | 北京辛耕普华医疗科技有限公司 | Geometric calibration method for SPECT (single photon emission computed tomography) system |
CN104599302A (en) * | 2015-01-13 | 2015-05-06 | 上海联影医疗科技有限公司 | Method for obtaining PET crystal energy peak value and setting energy frequency discriminator |
CN105115994A (en) * | 2015-07-22 | 2015-12-02 | 武汉数字派特科技有限公司 | Digital PET energy parameterization calibration method and system |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111199081A (en) * | 2018-10-31 | 2020-05-26 | 南京智能仿真技术研究院有限公司 | Energy spectrum data processing device and simulation system and method for simulating energy spectrum curve thereof |
CN109949953A (en) * | 2019-04-02 | 2019-06-28 | 中国科学院高能物理研究所 | Nuclear fuel rod degree of enrichment passive detection system and degree of enrichment reconstructing method |
CN110211095A (en) * | 2019-05-06 | 2019-09-06 | 东软医疗***股份有限公司 | Energy diagram and crystal positions search table generating method, device, storage medium |
CN110211095B (en) * | 2019-05-06 | 2023-05-30 | 沈阳智核医疗科技有限公司 | Energy map and crystal position lookup table generation method, device and storage medium |
CN110148991A (en) * | 2019-06-17 | 2019-08-20 | 北京理工大学 | A kind of over-charging of battery method for early warning and system based on big data |
CN110148991B (en) * | 2019-06-17 | 2020-12-08 | 北京理工大学 | Big data-based battery overcharge early warning method and system |
Also Published As
Publication number | Publication date |
---|---|
CN107133462B (en) | 2020-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107133462A (en) | Data processing method, device and equipment | |
Niesen et al. | Algorithm 919: A Krylov subspace algorithm for evaluating the ϕ-functions appearing in exponential integrators | |
Zucker et al. | Weighted log rank type statistics for comparing survival curves when there is a time lag in the effectiveness of treatment | |
US20030182597A1 (en) | Method for optimization of memory usage for a computer program | |
Christianson et al. | Automated size‐specific CT dose monitoring program: assessing variability in CT dose | |
WO2019144556A1 (en) | Electronic device, automatic case allocation method for credit card collection cases, and storage medium | |
CN108983279A (en) | A kind of Low background Spectra Unfolding Methods based on Nal detector | |
DE112012005014T5 (en) | Perform arithmetic operations using both large and small floating-point values | |
CN105139282B (en) | A kind of power grid achievement data processing method, device and computing device | |
DE102007009300A1 (en) | Computer system, has register provided for storing connection information, which indicates whether parallel/serial-advanced technology attachment-memory unit is connected with connector | |
CN112017042A (en) | Resource quota determining method and device based on tweed distribution and electronic equipment | |
CN106528755A (en) | Hot topic generation method and device | |
CN115730605B (en) | Data analysis method based on multidimensional information | |
CN108717744A (en) | Identify the method, apparatus and terminal device of the seal serial number on financial documents | |
Gandhi et al. | The XMM large scale structure survey: properties and two-point angular correlations of point-like sources | |
CN110378389A (en) | A kind of Adaboost classifier calculated machine creating device | |
CN110378607A (en) | A kind of automobile financial service system based on algorithm | |
US11075004B2 (en) | Discrimination result providing apparatus, operating method of discrimination result providing apparatus, discrimination result providing program, and discrimination result providing system | |
CN111160394A (en) | Training method and device of classification network, computer equipment and storage medium | |
CN116030950A (en) | Medical data integration management method | |
CN108091398B (en) | Patient grouping method and device | |
Ria et al. | A comparison of COVID-19 and imaging radiation risk in clinical patient populations | |
EP3279697A1 (en) | Technique for normalizing scintigraphy image | |
CN108764301A (en) | A kind of distress in concrete detection method based on reversed rarefaction representation | |
CN101903912A (en) | Rendering using multiple intensity redistribution functions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |