CN107067001A - The determination method and system of salivary gland region of interest - Google Patents
The determination method and system of salivary gland region of interest Download PDFInfo
- Publication number
- CN107067001A CN107067001A CN201710062687.5A CN201710062687A CN107067001A CN 107067001 A CN107067001 A CN 107067001A CN 201710062687 A CN201710062687 A CN 201710062687A CN 107067001 A CN107067001 A CN 107067001A
- Authority
- CN
- China
- Prior art keywords
- interest
- region
- salivary gland
- pixel
- module
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
- G06T2207/10128—Scintigraphy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
This application provides the determination method and system of salivary gland region of interest, it is related to medical diagnosis technical field, wherein, the determination method of the salivary gland region of interest includes：First, the salivary gland dynamic imaging of patient and the template pre-established are subjected to image registration, afterwards, according to the coordinate of salivary gland region of interest in template the first region of interest is delineated in the salivary gland dynamic imaging of patient, wherein, salivary gland region of interest includes the position of bilateral parotid and bilateral salivary gland, so, the quantity of radiological imaging agent in each pixel in the first region of interest was calculated every one minute, whether the quantity calculated by the above method judges the quantity of developer in each pixel into positive growth, if, pixel is determined as the second region of interest, if not, then pixel the second region of interest is not determined as, finally, second region of interest is defined as to the region of interest of salivary gland, the region of interest of salivary gland can be fast and accurately sketched out by this method.
Description
Technical field
The present invention relates to the determination method and system of medical diagnostic field, more particularly to salivary gland region of interest.
Background technology
Salivary gland is the body of gland of people or vertebrate oral cavity endocrine saliva.The function of salivary gland, which is mainly, to salivate,
On the one hand it is used as moistening oral cavity, is on the other hand conducive to swallowing and speaking.Due to, contain amylase in the saliva of human secretory,
Starch that can also be in preliminary exposition food.
With the quickening pace of modern life with the increase of people's life stress, in recent years, the lesion of salivary gland is in rise year by year
Trend, for example, salivary gland inflammation, xerosis, enlargement and other lesions, above-mentioned lesion big city cause salivary gland function disorderly
Disorderly.At present, processing method medically is the salivary gland dynamic imaging for gathering patient, i.e., allow patient to inject developer first, so
Afterwards, patient is allowed to be kept at least 30 minutes with same posture, wherein at the 20th minute, doctor, which gives patient's acidic food, stimulates saliva
Liquid excretion, doctor obtains the counting of the radiological imaging agent in salivary gland in real time during this, finally, by doctor rule of thumb in saliva
The region of interest of salivary gland is sketched out in liquid gland dynamic imaging manually, then salivary gland pair is calculated by manual or calculator etc.
The quantity and sour post-stimulatory discharge percentage of radiological imaging agent intake, the concrete function of salivary gland is evaluated with this.
In addition, during medical treatment, the background around salivary gland is often chosen together as reference, because salivary gland
Color can progressively be deepened in salivary gland dynamic imaging in the presence of radiological imaging agent, and background is in the effect of physiological metabolism
Under in salivary gland dynamic imaging color can progressively shoal.By the comparison with background, can accurate measurement salivary gland it is interested
The quantity of radiological imaging agent in area.
Although using this conventional imageological examination of above-mentioned salivary gland dynamic imaging, the function of salivary gland can be evaluated,
But, rely primarily on the working experience of doctor in the operating process that position determination is carried out to salivary gland region of interest, subjectivity compared with
By force, diagnostic accuracy is poor.
The content of the invention
In view of this, the purpose of the embodiment of the present invention is the provision of the determination method of salivary gland region of interest, by inciting somebody to action
The salivary gland dynamic imaging of patient is matched with the template pre-established so that the position of salivary gland region of interest is determined more
Plus it is accurate.
In a first aspect, the embodiments of the invention provide the determination method of salivary gland region of interest, including：By the saliva of patient
Gland dynamic imaging and the template pre-established carry out image registration；
First to be delineated in the salivary gland dynamic imaging of patient interested according to the coordinate of salivary gland region of interest in template
Area, wherein, salivary gland region of interest includes the position of bilateral parotid and bilateral salivary gland；
The quantity of radiological imaging agent in each pixel in the first region of interest was calculated every one minute；
Judge the quantity of developer in each pixel whether into positive growth；
If it is, pixel is determined as into the second region of interest；
If it is not, then pixel is not determined as into the second region of interest；
Second region of interest is defined as to the region of interest of salivary gland.
With reference in a first aspect, the embodiments of the invention provide the possible embodiment of the first of first aspect, wherein, side
Method also includes：
According to the coordinate of background region of interest in template the 3rd region of interest is delineated in the salivary gland dynamic imaging of patient；
The quantity of radiological imaging agent in each pixel in the 3rd region of interest was calculated every one minute；
Judge the quantity of developer in each pixel whether into negative growth；
If it is, pixel is determined as into the 4th region of interest；
If it is not, then pixel is not determined as into the 4th region of interest；
4th region of interest is defined as to the region of interest of background.
With reference in a first aspect, the embodiments of the invention provide the possible embodiment of second of first aspect, wherein, institute
Stating method also includes：
The salivary gland dynamic imaging of previous many cases patient is chosen by doctor terminal, wherein, salivary gland dynamic imaging is put for injection
The image of the 15th minute after penetrating property developer；
The salivary gland dynamic imaging of previous many cases patient is divided into doctor terminal into head imaging and trunk is imaged；
Doctor terminal judges whether the position of head imaging is placed in the middle；
Doctor terminal judges whether head imaging is complete；
Doctor terminal judges whether head imaging is clear；
Using position, between two parties, completely and clearly a head imaging is used as reference picture to doctor terminal；
The head imaging of remaining and reference picture are carried out image registration by doctor terminal, and make head imaging and reference picture weight
It is folded；
Doctor terminal calculate it is overlapping after head imaging with the salivary gland region of interest in reference picture each pixel it is overlapping
Rate；
The pixel that Duplication exceedes threshold value set in advance is constituted template by doctor terminal.
With reference to the first possible embodiment of first aspect, the embodiments of the invention provide the third of first aspect
Possible embodiment, wherein, method also includes：
The quantity of the radiological imaging agent before syringe injection patient is gathered, and is denoted as c；
The quantity of the radiological imaging agent after syringe injection patient is gathered, and is denoted as d；
When calculating the 15th minute in the second region of interest radiological imaging agent quantity, and be denoted as a；
When calculating the 15th minute in the 4th region of interest radiological imaging agent quantity, and be denoted as b；
The uptake ratio e of salivary gland is calculated, wherein,
With reference to the first possible embodiment of first aspect, the embodiments of the invention provide the 4th of first aspect kind
Possible embodiment, wherein, method also includes：
Calculating acid stimulates the maximum quantity of radiological imaging agent in preceding second region of interest, and is denoted as f；
When radiological imaging agent reaches maximum quantity in preceding second region of interest of calculating acid stimulation, put in the 4th region of interest
The quantity of penetrating property developer, and it is denoted as h；
The minimum number of radiological imaging agent in the second region of interest after acid stimulation is calculated, and is denoted as g；
When radiological imaging agent reaches minimum number in the second region of interest after calculating acid is stimulated, put in the 4th region of interest
The quantity of penetrating property developer, and it is denoted as i；
The discharge percentage k of salivary gland is calculated, wherein,
Second aspect, the embodiments of the invention provide the determination system of salivary gland region of interest, including：
Registration module, for the salivary gland dynamic imaging of patient and the template pre-established to be carried out into image registration；
First region of interest delineates module, for the coordinate according to salivary gland region of interest in template patient salivary gland
The first region of interest is delineated in dynamic imaging, wherein, salivary gland region of interest includes the position of bilateral parotid and bilateral salivary gland；
Quantity statistics module, for calculating in the first region of interest radiological imaging agent in each pixel every one minute
Quantity；
First judge module, for judging the quantity of developer in each pixel whether into positive growth；
Affirmative performing module, for if it is, pixel is determined as into the second region of interest；
Negate performing module, for if it is not, then pixel is not determined as into the second region of interest；
Region of interest determining module, the region of interest for the second region of interest to be defined as to salivary gland.
With reference to second aspect, the embodiments of the invention provide the possible embodiment of the first of second aspect, wherein, also
Including：
3rd region of interest delineates module, and the salivary gland for the coordinate according to background region of interest in template in patient is moved
The 3rd region of interest is delineated in state imaging；
Number statistical module, for calculating in the 3rd region of interest Radionuclide imaging in each pixel every one minute
The quantity of agent；
Second judge module, for judging the quantity of developer in each pixel whether into negative growth；
Affirmative module, for if it is, pixel is determined as into the 4th region of interest；
Negate module, for if it is not, then pixel is not determined as into the 4th region of interest；
Background determining unit, the region of interest for the 4th region of interest to be defined as to background.
With reference to second aspect, the embodiments of the invention provide the possible embodiment of second of second aspect, wherein, also
Including：
Selection module, the salivary gland dynamic imaging for choosing previous many cases patient by doctor terminal, wherein, salivary gland dynamic
Image the image for the 15th minute after injection radiological imaging agent；
Division module, head imaging and trunk are divided into for doctor terminal by the salivary gland dynamic imaging of previous many cases patient
Portion is imaged；
Determination module, judges whether the position of head imaging is placed in the middle, whether complete and whether clear for doctor terminal；
With reference to determining module, for doctor terminal, using position, between two parties, completely and clearly a head imaging is used as reference picture；
Overlay module, carries out image registration, and show head for doctor terminal by the head imaging of remaining and reference picture
As overlapping with reference picture；
Duplication computing module, for doctor terminal calculate it is overlapping after head imaging and salivary gland in reference picture it is interested
The Duplication of each pixel in area；
Template comprising modules, template is constituted for doctor terminal by the pixel that Duplication exceedes threshold value set in advance.
With reference to the first possible embodiment of second aspect, the embodiments of the invention provide the third of second aspect
Possible embodiment, wherein, in addition to：
Former number calculating section, for gathering the quantity of the radiological imaging agent before syringe injection patient, and is denoted as c；
Existing number calculating section, for gathering the quantity of the radiological imaging agent after syringe injection patient, and is denoted as d；
Second region of interest quantity statistics module, radiological imaging agent in the second region of interest during for calculating the 15th minute
Quantity, and be denoted as a；
4th region of interest quantity statistics module, radiological imaging agent in the 4th region of interest during for calculating the 15th minute
Quantity, and be denoted as b；
Uptake ratio computing module, the uptake ratio e for calculating salivary gland, wherein,
With reference to the first possible embodiment of second aspect, the embodiments of the invention provide the 4th of second aspect kind
Possible embodiment, wherein, in addition to：
Quantity statistics module before salivary gland acid is stimulated, radiological imaging agent in preceding second region of interest is stimulated for calculating acid
Maximum quantity, and be denoted as f；
Quantity statistics module before background acid is stimulated, stimulates radiological imaging agent in preceding second region of interest to reach for calculating acid
During to maximum quantity, the quantity of radiological imaging agent in the 4th region of interest, and it is denoted as h；
Quantity statistics module after salivary gland acid is stimulated, for calculating after acid is stimulated radiological imaging agent in the second region of interest
Minimum number, and be denoted as g；
Quantity statistics module after background acid is stimulated, for calculating after acid is stimulated, radiological imaging agent reaches in the second region of interest
During to minimum number, the quantity of radiological imaging agent in the 4th region of interest, and it is denoted as i；Percentage computing module is discharged, is used for
The discharge percentage k of salivary gland is calculated, wherein,
The determination method and system of salivary gland region of interest provided in an embodiment of the present invention, wherein, the salivary gland is interested
The determination method in area includes：First, the salivary gland dynamic imaging of patient and the template pre-established are subjected to image registration, and pressed
The first region of interest is delineated in the salivary gland dynamic imaging of patient according to the coordinate of salivary gland region of interest in template, wherein, saliva
Liquid gland region of interest includes the position of bilateral parotid and bilateral salivary gland, meanwhile, calculated every one minute in the first region of interest
The quantity of radiological imaging agent in each pixel, by the comparison of above-mentioned quantity, to judge the number of developer in each pixel
Whether amount is into positive growth, if it is, pixel is determined as into the second region of interest, if it is not, then pixel is not determined as into the second sense
Region of interest, finally, the second region of interest is defined as the region of interest of salivary gland, and trouble is calculated by what the above method can be quantified
The salivary gland region of interest of person, and then improve the diagnostic accuracy to region of interest.
To enable the above objects, features and advantages of the present invention to become apparent, preferred embodiment cited below particularly, and coordinate
Appended accompanying drawing, is described in detail below.
Brief description of the drawings
, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical scheme of the prior art
The accompanying drawing used required in embodiment or description of the prior art is briefly described, it should be apparent that, in describing below
Accompanying drawing is some embodiments of the present invention, for those of ordinary skill in the art, before creative work is not paid
Put, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 shows the structure connection figure of the determination method for the salivary gland region of interest that the embodiment of the present invention is provided；
Fig. 2 shows the structural representation of the determination method for the salivary gland region of interest that the embodiment of the present invention is provided；
Fig. 3 shows the structural framing figure of the determination system for the salivary gland region of interest that the embodiment of the present invention is provided；
Fig. 4 shows the structure connection figure of the determination system for the salivary gland region of interest that the embodiment of the present invention is provided.
Icon：1- registration modules；The region of interest of 2- first delineates module；3- quantity statistics modules；The judge modules of 4- first；
5- affirms performing module；6- negates performing module；7- region of interest determining modules.
Embodiment
Below in conjunction with accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Ground is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.Generally exist
The component of the embodiment of the present invention described and illustrated in accompanying drawing can be arranged and designed with a variety of configurations herein.Cause
This, the detailed description of the embodiments of the invention to providing in the accompanying drawings is not intended to limit claimed invention below
Scope, but it is merely representative of the selected embodiment of the present invention.Based on embodiments of the invention, those skilled in the art are not doing
The every other embodiment obtained on the premise of going out creative work, belongs to the scope of protection of the invention.
At present, when the salivary gland function of needs assessment patient, medically conventional processing method is to gather the saliva of patient
Liquid gland dynamic imaging, that is, allow patient to inject developer in advance, then, allows patient to be kept at least 30 minutes with same posture, this
In the section time, the counting of the radiological imaging agent in salivary gland is obtained in real time by doctor, finally, by doctor rule of thumb in saliva
The region of interest of salivary gland is sketched out in gland dynamic imaging manually, the function of salivary gland can be evaluated using the above method, still,
The working experience of doctor is relied primarily in the operating process that position determination is carried out to salivary gland region of interest, it is subjective, lead
Cause the precision of diagnosis poor.
In addition, when doctor wants to obtain the intake percentage of patient 15 minutes, it is necessary to salivary gland during manual calculations 15 minutes
Radiocounting before radiocounting, and injection, after injection in syringe, then obtain 15 points by a series of calculating of formula
The intake percentage of clock, calculates trouble, takes.In addition, salivary gland dynamic imaging needs the image of collection 30 minutes, so,
Some patientss are difficult to keep fixed posture in gatherer process, further result in measurement accuracy decline.
Based on this, the embodiments of the invention provide the determination method and system of salivary gland region of interest, below by implementation
Example is described.
Embodiment 1
Referring to Fig. 1, the determination method for the salivary gland region of interest that the present embodiment is proposed specifically includes following steps：
Step S101：The salivary gland dynamic imaging of patient and the template pre-established are subjected to image registration, here, saliva
Gland dynamic imaging is the image that patient injects the 15th minute after radiological imaging agent.For the ease of checking result, diagnosed
It is preceding can patient internal introducing radiopharmaceutical, so, after radiopharmaceutical is introduced in vivo, dense can gather in target organ or tissue
In, the ray that it is launched is detected by imaging instrument, so as to obtain the distributed image of medicine in vivo, on this basis
To diagnose various diseases.Common radiological imaging agent has many kinds, for example, F-DOPA developers, dopamine operating body are imaged
Agent, dopamine receptor imaging agent, hydrocarbon tryptamines transporter developer etc..
Step S102：Coordinate according to salivary gland region of interest in template delineates in the salivary gland dynamic imaging of patient
One region of interest, wherein, salivary gland region of interest includes the position of bilateral parotid and bilateral salivary gland, medically, salivary gland bag
Include the parotid gland, salivary gland, sublingual gland and how each small salivary gland, still, salivary gland around oral cavity it is main with the parotid gland and
Based on salivary gland, therefore, the position that bilateral parotid and bilateral salivary gland are chosen in the method is the first region of interest.Obtaining
Take after the salivary gland dynamic imaging of patient, by being delineated in the salivary gland dynamic imaging and salivary gland region of interest one in template
The coordinate of cause, you can obtain the first region of interest of patient.
Step S103：The quantity of radiological imaging agent in each pixel in the first region of interest was calculated every one minute.
During concrete operations, calculated after patient's injection developer in 1 to 15 minutes every one minute in the first region of interest
The quantity of radiological imaging agent in each pixel, certainly, specific acquisition time interval can also be carried out according to circumstances flexibly
Setting.
Step S104：Judge the quantity of developer in each pixel whether into positive growth.
Because salivary gland color in salivary gland dynamic imaging in the presence of radiological imaging agent can progressively be deepened, root
A time-acttivity curve (X-axis is the time, and Y-axis is radiocounting) is drawn according to the quantity of developer in each pixel, it is right
The time-acttivity curve of each pixel carries out the slope that Linear Quasi merges calculated curve, with this, to judge to show in each pixel
As whether the quantity of agent is into positive growth.
Step S105：, i.e., can be by color in salivary gland dynamic imaging if it is, pixel is determined as into the second region of interest
The pixel progressively deepened in the presence of radiological imaging agent, is determined as the second region of interest.
Step S106：If it is not, then pixel is not determined as into the second region of interest.Color will not met can be in radioactivity
The pixel progressively deepened in the presence of developer, excludes the second region of interest.
Step S107：Second region of interest is defined as to the region of interest of salivary gland.Finally, the second region of interest is determined
For the region of interest of salivary gland.
Due to during the diagnosis of salivary gland, often needing by being compared with background, so, can accurately count
Measure the quantity of the radiological imaging agent in salivary gland region of interest.Therefore, referring to Fig. 2, the determination of the salivary gland region of interest
The region of interest that method also specifically includes background is as follows：
Step S201：According to the coordinate of background region of interest in template the 3rd is delineated in the salivary gland dynamic imaging of patient
Region of interest, specific operation process is with step S102.
Step S202：Radioactivity in each pixel was calculated every one minute in the 3rd region of interest in 1 to 15 minutes
The quantity of developer, it is emerging every the sense of calculating the 3rd in one minute in 1 to 15 minutes after patient's injection developer in concrete operations
In interesting area in each pixel radiological imaging agent quantity, certainly, specific acquisition time interval can also according to circumstances be entered
Row flexibly setting.
Step S203：Whether the quantity of developer in each pixel is judged into negative growth, because background is in Radionuclide imaging
Color can progressively be deepened in salivary gland dynamic imaging in the presence of agent, and one is drawn according to the quantity of developer in each pixel
Time-acttivity curve (X-axis is the time, and Y-axis is radiocounting), is carried out linearly to the time-acttivity curve of each pixel
Be fitted and calculated curve slope, with this, to judge the quantity of developer in each pixel whether into negative growth.
Step S204：If it is, pixel is determined as into the 4th region of interest, when judged result is negative growth, by the picture
Element is determined as the 4th region of interest.
Step S205：If it is not, then pixel is not determined as into the 4th region of interest, i.e., when judged result is negative growth
When, the pixel is excluded into the 4th region of interest.
Step S206：4th region of interest is defined as to the region of interest of background, that is, obtains the region of interest of background.
In the determination method of this salivary gland region of interest, setting up for template is most important, specifically includes：
The salivary gland dynamic imaging of previous many cases patient is chosen by doctor terminal, here, many cases are generally more than 30, wherein,
Salivary gland dynamic imaging is the image of the 15th minute after injection radiological imaging agent.
Concrete processing procedure includes：First, the region of interest for sketching out 2 backgrounds by doctor terminal (is particularly located at patient left
Side and right side head temporo area), the salivary gland dynamic imaging of previous many cases patient is divided into head imaging and trunk by doctor terminal
Imaging.Due to, the head of patient is likely to occur the rotation and twisting of low-angle, therefore, and head is separately handled with trunk, with
The midpoint of bilateral salivary gland lower edge and thyroid gland higher level's line is cut-off rule, and image is divided into head (more than cut-off rule) and trunk
Portion (below cut-off rule).The image for removing trunk only retains head image.By the head image of remaining patient and reference picture
Head image carries out the characteristic image registration based on Affine respectively so that the image of patient passes through most after stretching, rotation, scaling
Possible, holding identical position overlapping with the image of reference picture.
Afterwards, in order to obtain accurate template, judge whether the position of head imaging is placed in the middle, head imaging by doctor terminal
Whether whether the imaging of complete and head is clear.
Using position, between two parties, completely and clearly a head imaging is used as reference picture, i.e. patients head position to occupy to doctor terminal
In, image completely covers head, while bilateral parotid and bilateral salivary gland clearly develop.It has chosen above-mentioned such one
After reference picture, the head imaging of remaining and reference picture are carried out image registration by doctor terminal, and make head imaging and reference picture weight
It is folded.
The region of interest of the salivary gland of all patients after registration is superimposed, it is each in the region of interest after statistical stacking
The number of times that pixel is repeatedly delineated between different patients, the head imaging after doctor terminal calculating is overlapping and the salivary gland in reference picture
The pixel that Duplication exceedes threshold value set in advance is constituted template by the Duplication of each pixel in region of interest, doctor terminal.Tool
When gymnastics is made, selection wherein repetitive rate is more than 95% pixel as the salivary gland region of interest of template, for bilateral parotid and
Bilateral salivary gland generates region of interest respectively.
In addition, the absorption function in order to weigh salivary gland, referring to Fig. 3, the determination method of salivary gland region of interest also includes
The computational methods of the uptake ratio of salivary gland：
Step S301：The quantity of the radiological imaging agent before syringe injection patient is gathered, and is denoted as c.
Step S302：The quantity of the radiological imaging agent after syringe injection patient is gathered, and is denoted as d.Pass through syringe
The quantity and the quantity of the radiological imaging agent before syringe injection patient for injecting the radiological imaging agent after patient make difference
Obtain being injected into the quantity of all radiological imaging agents in patient's body.
Step S303：When calculating the 15th minute in the second region of interest radiological imaging agent quantity, and be denoted as a, that is, count
The quantity of the radiological imaging agent of salivary gland region of interest when calculating the 15th minute.
Step S304：When calculating the 15th minute in the 4th region of interest radiological imaging agent quantity, and be denoted as b, that is, count
The quantity of the radiological imaging agent of background region of interest, passes through the radiological imaging agent of salivary gland region of interest when calculating the 15th minute
Quantity and background region of interest radiological imaging agent quantity make poor, salivary gland region of interest when can accurately obtain 15 minutes
The quantity of interior radiological imaging agent.
Step S305：The uptake ratio e of salivary gland is calculated, wherein,
In addition, it is also an important indicator for investigating salivary gland function that acid, which is stimulated, at the 20th minute, patient is given acid thorn
Swash, for example, chewing and buccal vitamin C, afterwards, calculate the discharge percentage that salivary gland is stimulated acid, comprise the following steps that：
Calculating acid stimulates the maximum quantity of radiological imaging agent in preceding second region of interest, and is denoted as f, and calculates sour stimulation
The quantity of radiological imaging agent is remembered in the 4th region of interest when radiological imaging agent reaches maximum quantity in preceding second region of interest
It is h, also, calculates the minimum number of radiological imaging agent in the second region of interest after acid is stimulated, and be denoted as g, the after acid is stimulated
Radiological imaging agent reaches that the quantity that radiological imaging agent in the 4th region of interest is counted during minimum number is denoted as in two region of interest
I, here, acid carry out being defined at the time of acid is stimulated before stimulating with the sour post-stimulatory moment by patient, typically, often take acid to stimulate the moment
Forward or backward 1 minute.Afterwards, the discharge percentage k of salivary gland is calculated, wherein,
In summary, the determination method for the salivary gland region of interest that the present embodiment is provided includes：First, by the trouble got
The salivary gland dynamic imaging of person and the template pre-established carry out image registration, wherein, template is according to the case pre-entered
And set up, afterwards, first is delineated in the salivary gland dynamic imaging of patient according to the coordinate of salivary gland region of interest in template
Region of interest, wherein, salivary gland region of interest includes the position of bilateral parotid and bilateral salivary gland, i.e., by bilateral parotid and bilateral
Salivary gland all divides salivary gland region of interest into.Because salivary gland is in the presence of radiological imaging agent in salivary gland dynamic imaging
Middle color can progressively deepen, therefore, and radiological imaging agent in each pixel was calculated in the first region of interest every one minute
Quantity, the quantity of developer in each pixel is judged with reference to each above-mentioned quantity whether into positive growth, if into positive growth
, then pixel is determined as the second region of interest, if not into positive growth, then pixel is not determined as that second is interested
Area, finally, the second region of interest is defined as the region of interest of salivary gland, doctor can be eliminated the reliance on by above-mentioned quantization method
Working experience and obtain salivary gland region of interest, so as to improve the diagnostic accuracy to salivary gland region of interest.
Embodiment 2
Referring to Fig. 4, present embodiments providing the determination system of salivary gland region of interest includes：The registration module being sequentially connected
1st, the first region of interest delineates module 2, quantity statistics module 3, certainly the first judge module 4, performing module 5, negative execution mould
Block 6 and region of interest determining module 7, the concrete function of modules are described below：
Registration module 1 is used to the salivary gland dynamic imaging of patient and the template pre-established carrying out image registration, first
Region of interest, which delineates module 2, to be used to hook in the salivary gland dynamic imaging of patient according to the coordinate of salivary gland region of interest in template
The first region of interest is drawn, wherein, salivary gland region of interest includes the position of bilateral parotid and bilateral salivary gland, quantity statistics module
3 quantity for calculating radiological imaging agent in each pixel in the first region of interest every one minute, the first judge module 4
For whether judging the quantity of developer in each pixel into positive growth, if it is determined that be into positive growth, then performing module 5 certainly
Pixel is determined as the second region of interest, then negates that pixel is not determined as then by performing module 6 if it is determined that not being into positive growth
Second region of interest, region of interest determining module 7 is used for the region of interest that the second region of interest is defined as to salivary gland.
In addition, the determination system of salivary gland region of interest also includes with lower module：3rd region of interest, which delineates module, to be used for
According to the coordinate of background region of interest in template the 3rd region of interest, number statistical are delineated in the salivary gland dynamic imaging of patient
Module is used for the number for calculating radiological imaging agent in each pixel in the 3rd region of interest every one minute in 1 to 15 minutes
Whether amount, the second judge module is used for the quantity for judging developer in each pixel into negative growth, if it is, affirming module by picture
Element is determined as the 4th region of interest, if not, pixel is not then determined as the 4th region of interest, background determining unit by negative module
4th region of interest is defined as to the region of interest of background.
In addition, the determination system of salivary gland region of interest also includes：Choosing module is used to choose previous many cases by doctor terminal
The salivary gland dynamic imaging of patient, wherein, salivary gland dynamic imaging is the image of the 15th minute after injection radiological imaging agent, is drawn
Sub-module is used for doctor terminal and the salivary gland dynamic imaging of previous many cases patient is divided into head imaging and trunk imaging, judges
Module is used for doctor terminal and judges whether the position of head imaging is placed in the middle, whether complete and whether clear, is used with reference to determining module
In doctor terminal using position it is placed in the middle, complete and clearly a head imaging as reference picture, overlay module be used for doctor terminal by its
Head imaging and the reference picture of remaining example carry out image registration, and make head imaging overlapping with reference picture, Duplication computing module use
Head imaging and the Duplication of each pixel in the salivary gland region of interest in reference picture after doctor terminal calculating is overlapping, template
Comprising modules are used for doctor terminal and the pixel that Duplication exceedes threshold value set in advance are constituted into template.
In addition, the determination system of salivary gland region of interest also includes：Former number calculating section is used to gather syringe injection
The quantity of radiological imaging agent before patient, and c is denoted as, existing number calculating section is used to gather putting after syringe injection patient
The quantity of penetrating property developer, and d is denoted as, second feels emerging when the second region of interest quantity statistics module 3 is used to calculate the 15th minute
The quantity of radiological imaging agent in interesting area, and be denoted as a, the when the 4th region of interest quantity statistics module is used to calculate the 15th minute
The quantity of radiological imaging agent in four region of interest, and b is denoted as, uptake ratio computing module is used for the uptake ratio e for calculating salivary gland,
Wherein,
Because the post-stimulatory discharge percentage of acid is also the important indicator for weighing salivary gland function, salivary gland region of interest
Determination system also include：Quantity statistics module stimulates Radionuclide imaging in preceding second region of interest for calculating acid before acid is stimulated
The maximum quantity of agent, and f is denoted as, stimulate preceding quantity statistics module to calculate sour stimulate in preceding second region of interest by background acid and radiate
Property developer the quantity of radiological imaging agent is denoted as h, quantity statistics after acid is stimulated in the 4th region of interest when reaching maximum quantity
Module is used for the minimum number for calculating radiological imaging agent in the second region of interest after acid is stimulated, and is denoted as g, is stimulated by background acid
The 4th is interested when radiological imaging agent reaches minimum number in the second region of interest after quantity statistics module calculating acid is stimulated afterwards
The quantity of radiological imaging agent is denoted as i in area, and discharge percentage computing module is used for the discharge percentage k for calculating salivary gland, its
In,
In summary, the determination system for the salivary gland region of interest that the present embodiment is provided includes：What is be sequentially connected matches somebody with somebody quasi-mode
Block 1, the first region of interest delineate module 2, quantity statistics module 3, certainly the first judge module 4, performing module 5, negative execution
Module 6 and region of interest determining module 7, registration module 1 are used for the salivary gland dynamic imaging of patient and the template pre-established
Image registration is carried out, the first region of interest, which delineates module 2, is used for the coordinate according to salivary gland region of interest in template patient's
The first region of interest is delineated in salivary gland dynamic imaging, wherein, salivary gland region of interest includes bilateral parotid and bilateral salivary gland
Position, quantity statistics module 3 is used to calculate in the first region of interest radiological imaging agent in each pixel every one minute
Quantity, whether the first judge module 4 is used to judge the quantity of developer in each pixel into positive growth, if it is, holding certainly
Pixel is determined as the second region of interest by row module 5, if not, pixel is not then determined as that the second sense is emerging by negative performing module 6
Second region of interest, finally, is defined as the region of interest of salivary gland by region of interest determining module 7 by interesting area, by the system,
So that the diagnosis to salivary gland region of interest is more convenient accurate.
Finally it should be noted that：Embodiment described above, is only the embodiment of the present invention, to illustrate the present invention
Technical scheme, rather than its limitations, protection scope of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair
It is bright to be described in detail, it will be understood by those within the art that：Any one skilled in the art
The invention discloses technical scope in, it can still modify to the technical scheme described in previous embodiment or can be light
Change is readily conceivable that, or equivalent substitution is carried out to which part technical characteristic；And these modifications, change or replacement, do not make
The essence of appropriate technical solution departs from the spirit and scope of technical scheme of the embodiment of the present invention, should all cover the protection in the present invention
Within the scope of.Therefore, protection scope of the present invention described should be defined by scope of the claims.
Claims (10)
1. the determination method of salivary gland region of interest, it is characterised in that including：
The salivary gland dynamic imaging of patient and the template pre-established are subjected to image registration；
According to the coordinate of salivary gland region of interest in the template the first sense is delineated in the salivary gland dynamic imaging of the patient
Region of interest, wherein, the salivary gland region of interest includes the position of bilateral parotid and bilateral salivary gland；
The quantity of radiological imaging agent in each pixel in first region of interest was calculated every one minute；
Judge the quantity of developer in each pixel whether into positive growth；
If it is, the pixel is determined as into the second region of interest；
If it is not, then the pixel is not determined as into the second region of interest；
Second region of interest is defined as to the region of interest of salivary gland.
2. the determination method of salivary gland region of interest according to claim 1, it is characterised in that methods described also includes：
It is emerging that coordinate according to background region of interest in the template delineates the 3rd sense in the salivary gland dynamic imaging of the patient
Interesting area；
The quantity of radiological imaging agent in each pixel in the 3rd region of interest was calculated every one minute；
Judge the quantity of developer in each pixel whether into negative growth；
If it is, the pixel is determined as into the 4th region of interest；
If it is not, then the pixel is not determined as into the 4th region of interest；
4th region of interest is defined as to the region of interest of background.
3. the determination method of salivary gland region of interest according to claim 1, it is characterised in that methods described also includes：
The salivary gland dynamic imaging of previous many cases patient is chosen by doctor terminal, wherein, the salivary gland dynamic imaging is injection institute
State the image of the 15th minute after radiological imaging agent；
The salivary gland dynamic imaging of the previous many cases patient is divided into the doctor terminal into head imaging and trunk is imaged；
The doctor terminal judges whether the position of the head imaging is placed in the middle；
The doctor terminal judges whether the head imaging is complete；
The doctor terminal judges whether the head imaging is clear；
Using position, between two parties, completely and clearly a head imaging is used as reference picture to the doctor terminal；
The head imaging of remaining and the reference picture are carried out image registration by the doctor terminal, and image the head
It is overlapping with the reference picture；
The doctor terminal calculate it is described it is overlapping after head imaging with it is every in the salivary gland region of interest in the reference picture
The Duplication of individual pixel；
The pixel that Duplication exceedes threshold value set in advance is constituted the template by the doctor terminal.
4. the determination method of salivary gland region of interest according to claim 2, it is characterised in that methods described also includes：
The quantity for the radiological imaging agent that syringe is injected before the patient is gathered, and is denoted as c；
The quantity for the radiological imaging agent that syringe is injected after the patient is gathered, and is denoted as d；
Described in when calculating the 15th minute in the second region of interest radiological imaging agent quantity, and be denoted as a；
Described in when calculating the 15th minute in the 4th region of interest radiological imaging agent quantity, and be denoted as b；
The uptake ratio e of salivary gland is calculated, wherein,
5. the determination method of salivary gland region of interest according to claim 2, it is characterised in that methods described also includes：
Calculating acid stimulates the maximum quantity of radiological imaging agent in preceding second region of interest, and is denoted as f；
When radiological imaging agent reaches maximum quantity in second region of interest before calculating acid is stimulated, the 4th region of interest
The quantity of interior radiological imaging agent, and it is denoted as h；
The minimum number of radiological imaging agent in second region of interest after acid stimulation is calculated, and is denoted as g；
When radiological imaging agent reaches minimum number in second region of interest after calculating acid is stimulated, the 4th region of interest
The quantity of interior radiological imaging agent, and it is denoted as i；
The discharge percentage k of salivary gland is calculated, wherein,
6. the determination system of salivary gland region of interest, it is characterised in that including：
Registration module, for the salivary gland dynamic imaging of patient and the template pre-established to be carried out into image registration；
First region of interest delineates module, for the coordinate according to salivary gland region of interest in the template the patient saliva
The first region of interest is delineated in liquid gland dynamic imaging, wherein, the salivary gland region of interest is included under bilateral parotid and bilateral jaw
The position of gland；
Quantity statistics module, for calculating in first region of interest radiological imaging agent in each pixel every one minute
Quantity；
First judge module, for judging the quantity of developer in each pixel whether into positive growth；
Affirmative performing module, for if it is, the pixel is determined as into the second region of interest；
Negate performing module, for if it is not, then the pixel is not determined as into the second region of interest；
Region of interest determining module, the region of interest for second region of interest to be defined as to salivary gland.
7. the determination system of salivary gland region of interest according to claim 6, it is characterised in that also include：
3rd region of interest delineates module, for the coordinate according to background region of interest in the template the patient saliva
The 3rd region of interest is delineated in gland dynamic imaging；
Number statistical module, for calculating in the 3rd region of interest radiological imaging agent in each pixel every one minute
Quantity；
Second judge module, for judging the quantity of developer in each pixel whether into negative growth；
Affirmative module, for if it is, the pixel is determined as into the 4th region of interest；
Negate module, for if it is not, then the pixel is not determined as into the 4th region of interest；
Background determining unit, the region of interest for the 4th region of interest to be defined as to background.
8. the determination system of salivary gland region of interest according to claim 7, it is characterised in that also include：
Selection module, the salivary gland dynamic imaging for choosing previous many cases patient by doctor terminal, wherein, the salivary gland dynamic
Image to inject the image of the 15th minute after the radiological imaging agent；
Division module, for the doctor terminal by the salivary gland dynamic imaging of the previous many cases patient be divided into head imaging and
Trunk is imaged；
Determination module, judges whether the position of the head imaging is placed in the middle, whether complete and whether clear for the doctor terminal
It is clear；
With reference to determining module, for the doctor terminal, using position, between two parties, completely and clearly a head imaging is used as ginseng
According to figure；
Overlay module, image registration is carried out for the doctor terminal by the head imaging of remaining and the reference picture, and
Make the head imaging overlapping with the reference picture；
Duplication computing module, for the doctor terminal calculate it is described it is overlapping after the head imaging with the reference picture in
The Duplication of each pixel in salivary gland region of interest；
Template comprising modules, the template is constituted for the doctor terminal by the pixel that Duplication exceedes threshold value set in advance.
9. the determination system of salivary gland region of interest according to claim 7, it is characterised in that also include：
Former number calculating section, for gathering the quantity for the radiological imaging agent that syringe is injected before the patient, and is denoted as c；
Existing number calculating section, for gathering the quantity for the radiological imaging agent that syringe is injected after the patient, and is denoted as d；
Second region of interest quantity statistics module, radiological imaging agent in the second region of interest described in during for calculating the 15th minute
Quantity, and be denoted as a；
4th region of interest quantity statistics module, radiological imaging agent in the 4th region of interest described in during for calculating the 15th minute
Quantity, and be denoted as b；
Uptake ratio computing module, the uptake ratio e for calculating salivary gland, wherein,
10. the determination system of salivary gland region of interest according to claim 7, it is characterised in that also include：
Salivary gland is sour to stimulate preceding quantity statistics module, and radiological imaging agent in preceding second region of interest is stimulated for calculating acid
Maximum quantity, and be denoted as f；
Quantity statistics module before background acid is stimulated, for calculating before acid is stimulated, radiological imaging agent reaches in second region of interest
During to maximum quantity, the quantity of radiological imaging agent in the 4th region of interest, and it is denoted as h；
Quantity statistics module after salivary gland acid is stimulated, for calculating after acid is stimulated radiological imaging agent in second region of interest
Minimum number, and be denoted as g；
Quantity statistics module after background acid is stimulated, for calculating after acid is stimulated, radiological imaging agent reaches in second region of interest
During to minimum number, the quantity of radiological imaging agent in the 4th region of interest, and it is denoted as i；
Discharge percentage computing module, the discharge percentage k for calculating salivary gland, wherein,
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710062687.5A CN107067001A (en) | 2017-01-25 | 2017-01-25 | The determination method and system of salivary gland region of interest |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710062687.5A CN107067001A (en) | 2017-01-25 | 2017-01-25 | The determination method and system of salivary gland region of interest |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107067001A true CN107067001A (en) | 2017-08-18 |
Family
ID=59598939
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710062687.5A Pending CN107067001A (en) | 2017-01-25 | 2017-01-25 | The determination method and system of salivary gland region of interest |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107067001A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108171188A (en) * | 2018-01-04 | 2018-06-15 | 天津大学 | For the MCM methods of Rapid matching identification finger vena |
| CN113712532A (en) * | 2021-09-13 | 2021-11-30 | 核工业总医院 | Evaluation method mainly based on magnetic resonance diffusion weighted imaging after transient taste stimulation |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104414636A (en) * | 2013-08-23 | 2015-03-18 | 北京大学 | Magnetic resonance image based cerebral micro-bleeding computer auxiliary detection system |
| CN104751178A (en) * | 2015-03-31 | 2015-07-01 | 上海理工大学 | Pulmonary nodule detection device and method based on shape template matching and combining classifier |
| CN104819754A (en) * | 2015-05-13 | 2015-08-05 | 山东大学 | Medicine bottle liquid level detection method based on image processing |
-
2017
- 2017-01-25 CN CN201710062687.5A patent/CN107067001A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104414636A (en) * | 2013-08-23 | 2015-03-18 | 北京大学 | Magnetic resonance image based cerebral micro-bleeding computer auxiliary detection system |
| CN104751178A (en) * | 2015-03-31 | 2015-07-01 | 上海理工大学 | Pulmonary nodule detection device and method based on shape template matching and combining classifier |
| CN104819754A (en) * | 2015-05-13 | 2015-08-05 | 山东大学 | Medicine bottle liquid level detection method based on image processing |
Non-Patent Citations (2)
| Title |
|---|
| 何涛等: ""唾液腺显像常见影响因素分析"", 《标记免疫分析与临床》 * |
| 陈松等: ""GE Aladdin语言编程在唾液腺动态显像图像半自动化处理中的应用"", 《中国医科大学学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108171188A (en) * | 2018-01-04 | 2018-06-15 | 天津大学 | For the MCM methods of Rapid matching identification finger vena |
| CN108171188B (en) * | 2018-01-04 | 2021-11-23 | 天津大学 | MCM method for quickly matching and identifying finger veins |
| CN113712532A (en) * | 2021-09-13 | 2021-11-30 | 核工业总医院 | Evaluation method mainly based on magnetic resonance diffusion weighted imaging after transient taste stimulation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| DE102015122604A1 (en) | A periodontal disease diagnosis supportive system, a periodontal disease diagnosis supporting system, a periodontal disease diagnosis supporting program, and a periodontal disease diagnosis supporting method | |
| EP2742865B1 (en) | Method for testing periodontal disease | |
| CN108596861A (en) | A kind of CT image metal artifact minimizing technologies based on the study of depth residual error | |
| CN107067001A (en) | The determination method and system of salivary gland region of interest | |
| DE102012104373B4 (en) | Method for simulation of movement of pines and reconstruction with virtual functional bite registration | |
| CN106562800A (en) | Full oral digital information acquisition and saving method and system | |
| Ristau et al. | Comparison of AudaxCeph®’s fully automated cephalometric tracing technology to a semi-automated approach by human examiners | |
| Mitrirattanakul et al. | Accuracy of the intraoral scanner for detection of tooth wear | |
| Strauss et al. | The use of digital technologies in peri‐implant soft tissue augmentation–A narrative review on planning, measurements, monitoring and aesthetics | |
| Dallazen et al. | Comparison of manual (2D) and digital (3D) methods in the assessment of simulated facial edema | |
| CN110477941B (en) | Method, apparatus, computer device and readable storage medium for correcting intake value | |
| RU2421136C2 (en) | Method of determining maturity of hard dental tissues | |
| US20230215001A1 (en) | Method for facilitating caries detection | |
| CN104507415B (en) | Manufacturing device and diagnosis device for periodontal index, | |
| Seo et al. | The reliability of an easy measuring method for abutment convergence angle with a computer-aided design (CAD) system | |
| Amin et al. | The relationships of the maxillary central incisors and canines to the incisive papilla in Jordanians | |
| Li et al. | Accuracy of two best-fit alignment strategies with different reference areas for wear measurement with an intraoral scanner: an in vitro study. | |
| Nimbalkar | Accuracy of volumetric analysis software packages in assessment of tooth volume using CBCT | |
| Hujoel et al. | Quantification of periodontal attachment at single‐rooted teeth | |
| Şakar et al. | Assessment of the accuracy of orthodontic digital models in dental education | |
| CN106447650A (en) | Human body thyroid weight measuring method based on SPECT (Single Photon Emission Computed Tomography) planar imaging | |
| Abdolali et al. | Fully automated detection of the mandibular canal in cone beam CT images using Lie group based statistical shape models | |
| Park et al. | Impact of marginal defects on the accuracy of automated finish line detection in tooth preparation. | |
| Soesilawati et al. | Comparison of Permanent Teeth Eruption by Chronological Age in Indonesian Children. | |
| WO2011062473A1 (en) | Regression method for modeling widths of teeth |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170818 |