CN107133462A - Data processing method, device and equipment - Google Patents

Data processing method, device and equipment Download PDF

Info

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
Application number
CN201710288367.1A
Other languages
Chinese (zh)
Other versions
CN107133462B (en
Inventor
马润霞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201710288367.1A priority Critical patent/CN107133462B/en
Publication of CN107133462A publication Critical patent/CN107133462A/en
Application granted granted Critical
Publication of CN107133462B publication Critical patent/CN107133462B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject 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

Data processing method, device and equipment
【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.
CN201710288367.1A 2017-04-27 2017-04-27 Data processing method, device and equipment Active CN107133462B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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.