CN106667513A - CT measurement method for detecting renal filtration function - Google Patents
CT measurement method for detecting renal filtration function Download PDFInfo
- Publication number
- CN106667513A CN106667513A CN201611159034.0A CN201611159034A CN106667513A CN 106667513 A CN106667513 A CN 106667513A CN 201611159034 A CN201611159034 A CN 201611159034A CN 106667513 A CN106667513 A CN 106667513A
- Authority
- CN
- China
- Prior art keywords
- kidney
- contrast medium
- injection
- voxel
- values
- 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
- 238000000691 measurement method Methods 0.000 title abstract 3
- 230000013878 renal filtration Effects 0.000 title abstract 3
- 210000003734 kidney Anatomy 0.000 claims abstract description 128
- 239000002872 contrast media Substances 0.000 claims abstract description 98
- 238000002347 injection Methods 0.000 claims abstract description 46
- 239000007924 injection Substances 0.000 claims abstract description 46
- 238000009825 accumulation Methods 0.000 claims abstract description 39
- 238000002591 computed tomography Methods 0.000 claims abstract description 35
- 230000024924 glomerular filtration Effects 0.000 claims abstract description 14
- 238000012417 linear regression Methods 0.000 claims abstract description 13
- 238000000034 method Methods 0.000 claims description 41
- 238000001914 filtration Methods 0.000 claims description 33
- 238000006243 chemical reaction Methods 0.000 claims description 20
- 238000001514 detection method Methods 0.000 claims description 20
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 claims description 19
- 229910052740 iodine Inorganic materials 0.000 claims description 19
- 239000011630 iodine Substances 0.000 claims description 19
- 230000009977 dual effect Effects 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 8
- 210000000709 aorta Anatomy 0.000 claims description 5
- PNDPGZBMCMUPRI-UHFFFAOYSA-N iodine Chemical group II PNDPGZBMCMUPRI-UHFFFAOYSA-N 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 18
- 239000000193 iodinated contrast media Substances 0.000 description 17
- 238000005259 measurement Methods 0.000 description 12
- 230000003907 kidney function Effects 0.000 description 10
- 210000002700 urine Anatomy 0.000 description 8
- 210000005239 tubule Anatomy 0.000 description 6
- 239000003795 chemical substances by application Substances 0.000 description 4
- 239000003814 drug Substances 0.000 description 4
- 210000000244 kidney pelvis Anatomy 0.000 description 4
- 238000002360 preparation method Methods 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- UVKNREAAHQRBKA-IEOVAKBOSA-I 2-[bis[2-[bis(carboxylatomethyl)amino]ethyl]amino]acetate;technetium-99(4+) Chemical compound [99Tc+4].[O-]C(=O)CN(CC([O-])=O)CCN(CC(=O)[O-])CCN(CC([O-])=O)CC([O-])=O UVKNREAAHQRBKA-IEOVAKBOSA-I 0.000 description 3
- 241000483399 Ipimorpha retusa Species 0.000 description 3
- 208000000223 Solitary Kidney Diseases 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 230000002503 metabolic effect Effects 0.000 description 3
- 210000002381 plasma Anatomy 0.000 description 3
- 238000011158 quantitative evaluation Methods 0.000 description 3
- 239000002002 slurry Substances 0.000 description 3
- 239000002699 waste material Substances 0.000 description 3
- 238000010989 Bland-Altman Methods 0.000 description 2
- 229920001202 Inulin Polymers 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 238000010219 correlation analysis Methods 0.000 description 2
- 239000003792 electrolyte Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000029142 excretion Effects 0.000 description 2
- 238000010253 intravenous injection Methods 0.000 description 2
- JYJIGFIDKWBXDU-MNNPPOADSA-N inulin Chemical compound O[C@H]1[C@H](O)[C@@H](CO)O[C@@]1(CO)OC[C@]1(OC[C@]2(OC[C@]3(OC[C@]4(OC[C@]5(OC[C@]6(OC[C@]7(OC[C@]8(OC[C@]9(OC[C@]%10(OC[C@]%11(OC[C@]%12(OC[C@]%13(OC[C@]%14(OC[C@]%15(OC[C@]%16(OC[C@]%17(OC[C@]%18(OC[C@]%19(OC[C@]%20(OC[C@]%21(OC[C@]%22(OC[C@]%23(OC[C@]%24(OC[C@]%25(OC[C@]%26(OC[C@]%27(OC[C@]%28(OC[C@]%29(OC[C@]%30(OC[C@]%31(OC[C@]%32(OC[C@]%33(OC[C@]%34(OC[C@]%35(OC[C@]%36(O[C@@H]%37[C@@H]([C@@H](O)[C@H](O)[C@@H](CO)O%37)O)[C@H]([C@H](O)[C@@H](CO)O%36)O)[C@H]([C@H](O)[C@@H](CO)O%35)O)[C@H]([C@H](O)[C@@H](CO)O%34)O)[C@H]([C@H](O)[C@@H](CO)O%33)O)[C@H]([C@H](O)[C@@H](CO)O%32)O)[C@H]([C@H](O)[C@@H](CO)O%31)O)[C@H]([C@H](O)[C@@H](CO)O%30)O)[C@H]([C@H](O)[C@@H](CO)O%29)O)[C@H]([C@H](O)[C@@H](CO)O%28)O)[C@H]([C@H](O)[C@@H](CO)O%27)O)[C@H]([C@H](O)[C@@H](CO)O%26)O)[C@H]([C@H](O)[C@@H](CO)O%25)O)[C@H]([C@H](O)[C@@H](CO)O%24)O)[C@H]([C@H](O)[C@@H](CO)O%23)O)[C@H]([C@H](O)[C@@H](CO)O%22)O)[C@H]([C@H](O)[C@@H](CO)O%21)O)[C@H]([C@H](O)[C@@H](CO)O%20)O)[C@H]([C@H](O)[C@@H](CO)O%19)O)[C@H]([C@H](O)[C@@H](CO)O%18)O)[C@H]([C@H](O)[C@@H](CO)O%17)O)[C@H]([C@H](O)[C@@H](CO)O%16)O)[C@H]([C@H](O)[C@@H](CO)O%15)O)[C@H]([C@H](O)[C@@H](CO)O%14)O)[C@H]([C@H](O)[C@@H](CO)O%13)O)[C@H]([C@H](O)[C@@H](CO)O%12)O)[C@H]([C@H](O)[C@@H](CO)O%11)O)[C@H]([C@H](O)[C@@H](CO)O%10)O)[C@H]([C@H](O)[C@@H](CO)O9)O)[C@H]([C@H](O)[C@@H](CO)O8)O)[C@H]([C@H](O)[C@@H](CO)O7)O)[C@H]([C@H](O)[C@@H](CO)O6)O)[C@H]([C@H](O)[C@@H](CO)O5)O)[C@H]([C@H](O)[C@@H](CO)O4)O)[C@H]([C@H](O)[C@@H](CO)O3)O)[C@H]([C@H](O)[C@@H](CO)O2)O)[C@@H](O)[C@H](O)[C@@H](CO)O1 JYJIGFIDKWBXDU-MNNPPOADSA-N 0.000 description 2
- 229940029339 inulin Drugs 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000028327 secretion Effects 0.000 description 2
- 238000002603 single-photon emission computed tomography Methods 0.000 description 2
- 206010062237 Renal impairment Diseases 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000013211 curve analysis Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- -1 freely filter Chemical compound 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000001990 intravenous administration Methods 0.000 description 1
- NTHXOOBQLCIOLC-UHFFFAOYSA-N iohexol Chemical compound OCC(O)CN(C(=O)C)C1=C(I)C(C(=O)NCC(O)CO)=C(I)C(C(=O)NCC(O)CO)=C1I NTHXOOBQLCIOLC-UHFFFAOYSA-N 0.000 description 1
- 229960001025 iohexol Drugs 0.000 description 1
- 208000017169 kidney disease Diseases 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036470 plasma concentration Effects 0.000 description 1
- 231100000857 poor renal function Toxicity 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000012217 radiopharmaceutical Substances 0.000 description 1
- 229940121896 radiopharmaceutical Drugs 0.000 description 1
- 230000002799 radiopharmaceutical effect Effects 0.000 description 1
- 230000009103 reabsorption Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000000626 ureter Anatomy 0.000 description 1
- 230000002485 urinary effect Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Optics & Photonics (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention provides a CT measurement method for detecting a unilateral renal filtration function. The measurement method comprises the steps that CT scanning is performed at a certain moment after injection of a contrast agent, the accumulation amount of the contrast agent in the kidney is measured, the total contrast agent injection amount is divided by the accumulation amount to obtain an accumulation fraction FRA of the kidney, the FRA can directly serve as a pointer for assessing the GFR of the renal filtration function, or the mathematic relation of the FRA and the GFR is determined through linear-regression analysis, and a formula is established to calculate glomerular filtration rate (GFR) based on the FRA.
Description
Technical field
The present invention relates to detect the CT measuring methods of kidney filtering function.
Background technology
Blood Jing glomerular filtrations generate crude urine, then Jing renal tubules and collecting tubule absorb part moisture content therein and electrolyte
Urine discharge is condensed into, kidney filtering function is to removing internal metabolic waste, maintaining the balance of water and electrolyte most important.
Usually said kidney filtering function is the ability that kidney filters blood plasma in the unit interval, because filtration process occurs
In glomerule, so referred to as glomerular filtration rate (GFR), the GFR of adult normal is 120ml/min, that is to say, that meeting per minute
The blood plasma Jing glomerular filtrations for having 120ml form crude urine (being about 60ml per side kidney), and through renal tubules, collecting duct system weight
The new most moisture (about 98%, that is, concentrate) absorbed in crude urine ultimately form urine, and the metabolic waste in blood plasma is in company with urine
Liquid is excreted, and such as impaired renal function, GFR decline, and metabolic waste can not effectively be removed, will accumulate to form uremic in vivo
Disease.
GFR is the pointer of quantitative evaluation renal function.If a certain medicine can freely be filtered from glomerule, and in renal tubules
It is interior to be reuptaked, secretion is there will not be, the kidney clearance rate of this medicine is numerically equal to GFR, at this moment GFR=filters
Excessive/(the plasma concentration * time).The characteristics of inulin has such in physiological test thus be used to measure GFR, but
Because injection, drawing blood, staying the inconveniences such as urine, without actual application value, radiopharmaceutical 99mTc-DTPA in clinical practice
Also there are similar characteristics, while its filtration in kidney, accumulation can in vitro be detected by measurement apparatus such as SPECT, because
This is used to measure GFR (Gates methods).
Most of CT contrast medium such as Iopromides, Iohexol, etc. all have the characteristic of similar inulin, i.e., freely filter,
Finally will not be discharged with urine by reabsorption and secretion, therefore be also good GFR measurement markers things.And CT scan is surveyed
Medicine (iodine) concentration in CT values and tissue for obtaining is linear, therefore the CT scan (being commonly called as enhanced CT) of injection of contrast agent exists
May also be used for measuring GFR while observation pathological change form, i.e., for evaluating the filtering function of kidney.
However, Patlak patterns CT measurement glomerular filtration rate (renal function) of existing reported in literature are because a variety of limit to simultaneously
May not apply to clinic.
The content of the invention
In view of this, it is an object of the invention to provide a kind of CT measurements for quickly and easily detecting kidney filtering function
Method.
To reach above-mentioned purpose, the present invention adopts following technical proposals：
Certain time t Jing after injection contrast medium1Carry out CT scan and measure the accumulation amount of contrast medium in kidney, it is poly- with this
Accumulated amount obtain divided by the injection rate of total contrast medium kidney accumulation fraction (fractional renal accumulation,
FRA), the FRA can be as the index of assessment renal function GFR height；The mathematical formulae that simplicity can further be set up is utilized
FRA is calculating GFR.
The CT measuring methods of above-mentioned detection kidney filtering function specifically include following steps：
Injection of contrast agent will be started as time zero, through the moment t of stipulated time1, CT scan is carried out to kidney and is obtained
To t1The enhanced CT data of moment kidney, the data include constituting all voxels of kidney 3-D view, can therefrom obtain every
The CT values of individual voxel, that is, represent contrast medium concentrations, and all voxel volume sums, that is, represent the volume (V) of kidney；
T is obtained according to the enhanced CT data of the kidney and the volume (V) of kidney1The contrast medium accumulation amount of moment kidney
(Q)；
By the accumulation amount (Q) of contrast medium in the t1 moment kidney divided by total contrast medium injection rate, the poly- of kidney is obtained
Fraction (FRA)；
The accumulation fraction (FRA) of the kidney can be directly as the pointer of assessment kidney filtering function height.
Further, methods described also includes setting up linear regression model (LRM) to calculate based on the accumulation fraction (FRA) of kidney
Glomerular filtration rate (GFR).
Preferably, according to a large amount of early-stage Studies of applicant, based on the accumulation fraction of kidney set up linear regression model (LRM) come
The model for calculating glomerular filtration rate GFR is as shown in Equation 1：
GFR=1416.42 × FRA+2.58 ... (formula 1).
Preferably, the moment t1Span be injection of contrast agent after 2-3 minutes.
It is highly preferred that the moment t1Span be injection of contrast agent after 2.5 minutes.
Moment t1Span selection principle when being CT scan gathered data accumulation quantity of the contrast medium in kidney it is larger
And the amount of outflow ureter and bladder is less, now carries out CT scan gathered data most accurate.
Preferably, before injection of contrast agent, also include carrying out kidney the unenhanced steps of CT, to obtain the unenhanced of kidney
CT data, the data include constituting all voxels of kidney 3-D view, can therefrom obtain the CT values of each voxel, that is, contrast
The background of CT values when agent concentration is zero.
The t1Contrast medium accumulation amount (Q) of moment kidney is by the t1Own in the enhanced CT data of moment kidney
Mean CT-number (the HU of voxel1) deduct the mean CT-number (HU of all voxels in the plain CT data0), then the body for being multiplied by kidney
What product (V) was obtained,
Or by the t1CT value (the HU of each voxel in the enhanced CT data of moment kidneyi) and its volume
(Vvoxel) the sum of products, deduct the CT value (HU of each voxel in the plain CT datak) and its volume (Vvoxel) product
Sum is obtained.
It should be noted that about the CT values of voxel and the statistics of volume or calculating being all that third party is soft in CT measurement data
What part was completed, or the poster processing soft on CT machines completes, and can be completed by business software or CT machine manufacturer's software kits.
Preferably, when the CT scan is dual intensity CT scan or power spectrum CT scan, the t1The contrast medium of moment kidney
Accumulation amount (Q) has another computational methods, subtracts shadow by dual energy subtraction or power spectrum and obtains, specially：
When dual intensity CT scan or power spectrum CT scan, it is not necessary to injection contrast medium before kidney is carried out it is unenhanced, directly by t1
The dual energy subtraction or power spectrum at moment subtracts the incrementss (HU that shadow obtains each voxel CT values in kidneyy), it is multiplied by the volume of each voxel
(Vvoxel) and summation obtains contrast medium accumulation amount (Q) of kidney,
Or by the incrementss (HU of each voxel CT valuesy) obtain the average increase of all voxel CT values in kidney
(HUm), then be multiplied by the volume (V) of kidney and obtain contrast medium accumulation amount (Q) of kidney.
Preferably, total contrast medium injection rate is to be multiplied by contrast medium concentrations by the volume of contrast medium to be multiplied by contrast again
What the conversion coefficient (k) of agent and CT values was obtained.
Preferably, the contrast medium that measuring method is used is iodine contrast medium.
Conversion coefficient k (the units of iodine concentration and CT values：HU*ml/mgI), (mgI is milligram iodine), and the bulb of CT scan
Magnitude of voltage kV is relevant, and bulb magnitude of voltage is higher, and the value of the conversion coefficient is less, when bulb voltage is a certain fixed value, this turn
Coefficient relative constancy is changed, can be determined by the scanning of CT body mould.
For the tube voltage and CT scanner that determine, the conversion coefficient is more stable, i.e.,：
When the bulb voltage of CT scan is 120KV, the conversion coefficient of iodine contrast medium and CT values is 25, and unit is HU*ml/
MgI, i.e., 1mgI is equivalent to 25HU*ml；Can be understood as 1 bold and unconstrained gram iodine can make the CT values that 1 person of outstanding talent rises tissue lift 25HU.
When the bulb voltage of CT scan is 100KV, the conversion coefficient of iodine contrast medium and CT values is 30, and unit is HU*ml/
MgI, i.e., 1mgI is equivalent to 30HU*ml；Can be understood as 1 bold and unconstrained gram iodine can make the CT values that 1 person of outstanding talent rises tissue lift 30HU.
On the other hand, total contrast medium injection rate can also be come by the product of the volume of contrast medium and contrast medium concentrations
Represent, calculate in order to unify dimension when building up fraction (FRA), contrast medium accumulation amount (Q) of the kidney should be divided by conversion system
Number (k).
The dimension of now described total iodine contrast medium injection rate is mgI (bold and unconstrained gram iodine), the amount of above-mentioned contrast medium accumulation amount (Q)
Guiding principle (HU*ml) should be by being scaled mgI (bold and unconstrained gram iodine) divided by conversion coefficient k (dimension of k is HU*ml/mgI)
Preferably, the time used by contrast medium injection should control in a rational scope, contrast medium injection should be 20
Complete in second, it is highly preferred that completing contrast medium injection in 15 seconds.
When further, in order to avoid occurring in endaortic from beginning injection of contrast agent to contrast medium between Different Individual
Between impact of the difference to measurement result, to the time t1Time zero do regulation further below：
After injection of contrast agent, will start contrast medium occur in aorta as time zero, elapsed time t1To kidney
CT scan is carried out, t is obtained1The enhanced CT data of moment kidney, remaining calculates the measuring method of the accumulation fraction (FRA) of kidney
And step with will start injection of contrast agent it is identical as the measuring method and step of time zero.
Do so not only can avoid between Different Individual from start injection of contrast agent occur in contrast medium it is endaortic
Impact of the time difference to measurement result, and the stability and reliability of result can be strengthened.But need after injection contrast medium
Carry out CT monolayers dynamic scan to follow the trail of the change of CT values in aorta, the moment that contrast medium occurs is determined, to a certain degree
On increased the complexity of this technology.
Because the present invention is using the CT scan measuring method of iconography, the CT scan data of every side kidney can be obtained, therefore
Can be used to calculate the FRA of unilateral kidney, and then for assessing the renal function GFR height of one side.
According to the present invention it is proposed that a kind of CT measuring methods of new detection kidney filtering function, are strengthened using conventional CT and are swept
The data information assessment renal function retouched.The method simple, intuitive, the conventional CT that is expert at being capable of quantitative evaluation kidney while diagnosis
Function, and will not additionally increase the x beta radiation amounts of patient, with higher clinical value.
Description of the drawings
Fig. 1 is the image of the kidney CT unenhanced phases of the CT measuring methods based on the detection kidney filtering function of the present invention.
Fig. 2 is CT measuring methods t based on the detection kidney filtering function of the present invention1The figure of moment kidney Enhanced CT
Picture.
Fig. 3 is CT measuring methods t based on the detection kidney filtering function of the present invention1The double kidneys of moment Enhanced CT are three-dimensional
Reconstruction image, shows double kidney volumes and respective mean CT-number.
Fig. 4 is that the linear relationship chart between the iodine contrast medium concentrations and CT values for drawing, CT scan pipe are tested based on phantom study
Voltage is 100kV, and both linear relationships are y=30.01x+45.14, can obtain accordingly contrast medium iodine concentration and CT values it
Between conversion coefficient k=30 (HU*ml/mgI).
Fig. 5 is based on the linear correlation/regression curve analysis chart between FRA and single kidney Gates-GFR.
Fig. 6 is based on the correlation analysiss between list kidney CT-GFR and single kidney Gates-GFR.
Fig. 7 is the renal function (Gates- that renal function (CT-GFR) and ING measured by the inventive method are measured
GFR result consistency analysis figure (Bland-Altman figures), can draw from figure between), and both concordance are very high, and measurement is missed
Difference is little.
Specific embodiment
【Embodiment】
In the embodiment of the CT measuring methods of the detection kidney filtering function of the present invention, before injection of contrast agent,
It is that CT is unenhanced to subject's kidney, excretion phase (2.5 minutes) then after injection of contrast agent, CT scan is carried out again, this
It is referred to as t in application1Moment kidney Enhanced CT.By comparing t1The total injection of moment kidney CT values incrementss and iodine contrast medium
Measure to calculate glomerular filtration rate GFR.
The CT scan of present embodiment is conventional CT scan, comprising unenhanced and enhancing.
Referring to Fig. 1,2,3, the first embodiment of the present invention is described in detail.
The specific implementation step that person under inspection implements is carried out in the following order.
(unenhanced)
To person under inspection implement intravenous injection contrast medium before, kidney (tissue) is carried out first with CT scanner it is unenhanced,
And obtain unenhanced phase CT data.
Fig. 1 shows the unenhanced images of kidney CT, and the contour line in figure is the edge of kidney.
(injection of contrast agent)
After the unenhanced phase CT values of kidney are obtained, to person under inspection's intravenous injection contrast medium, contrast medium used herein is iodine pair
Compare agent.
(enhanced ct scans)
As iodine preparation is filtered and the accumulation in renal tubules, collecting tubule after intravenouss injection iodine contrast medium, arterial phase, medulla renis
Matter is in gradual enhancing, and renal medulla CT values are shown as apparently higher than cortex, part kidney to excretion early stage phase (2-3 minutes after injection)
Also the iodine preparation of high concentration is filled in small cup and renal pelvis, at this moment the CT of kidney strengthens and mainly caused by the filtration of glomerule.
About 2-3 minutes after injection, the iodine preparation that Jing glomerule is leached is gathered in renal medulla and renal pelvis, is not yet substantially entered defeated
Urinary catheter, which provides a suitable CT scan time window for measuring the filtration quantity of iodine preparation, the quantity can thus when
CT reinforcement values (relative to the incrementss of unenhanced CT values) and kidney volume (contain renal pelvis) product describing.
For this purpose, from start injection of contrast agent through the stipulated time t1Moment, this enforcement is selected 2.5 after injection of contrast agent
Minute, Enhanced CT is carried out to kidney, obtain volume V of enhanced CT data and kidney.Here CT scan pattern can be appearance
Product or spiral scanning mode, scope causes the double kidneys of few covering.
With reference to Fig. 2, t is being carried out1During moment Enhanced CT, the contrast medium overwhelming majority of Jing glomerular filtrations is confined to kidney
In dirty and renal pelvis.
In the first embodiment, Enhanced CT adopts tube voltage identical with plain CT, and scope is both needed to covering two
Side kidney, can arrange fixing pipe current scanning, it would however also be possible to employ tube current modulates (tube current modulation).
(calculating the accumulation amount of contrast medium in kidney)
T is obtained according to the CT data of the unilateral kidney and volume V of unilateral kidney1The moment contrast medium of one side kidney gathers
Accumulated amount Q.
In the present embodiment, the t1Moment contrast medium accumulation amount Q of one side kidney is by the t1Moment one side kidney
CT data in all voxels mean CT-number HU1Deduct the mean CT-number of kidney all voxels in one side in the plain CT data
HU0, then it is multiplied by what volume V of unilateral kidney was obtained.
Or be also an option that by the t1CT value (the HU of each voxel in the enhanced CT data of moment kidneyi) and its
Volume (Vvoxel) the sum of products, deduct the CT value (HU of each voxel in the plain CT datak) and its volume (Vvoxel)
The sum of products is obtained.
(calculating the accumulation fraction of contrast medium in kidney)
By accumulation amount Q of contrast medium in the kidney divided by total contrast medium injection rate (I), the accumulation fraction of kidney is obtained
FRA。
Total contrast medium injection rate be by the volume of contrast medium be multiplied by contrast medium iodine concentration be multiplied by again contrast medium and
What the conversion coefficient k of CT data was obtained.
In the present embodiment, as shown in figure 4, the line between the iodine contrast medium concentrations obtained for phantom study experiment and CT scan value
Sexual intercourse figure.As can be seen from the figure：When it is determined that the bulb voltage of CT scan is 100KV, the conversion of iodine contrast medium and CT values
Coefficient k is 30HU*ml/mgI.
(calculating the glomerular filtration rate GFR per side kidney)
As shown in figure 5, setting up linear regression model (LRM) based on the accumulation fraction of kidney to calculate glomerular filtration rate (GFR)
Model is as shown in Equation 1：
GFR=1416.42 × FRA+2.58 ... (formula 1).
The CT data gathered in this enforcement suffer from the nephrydrotic patient of right kidney from one, and Fig. 3 is based on the present invention
Detection kidney filtering function CT measuring methods t1Double kidney three-dimensional reconstruction images of the moment Enhanced CT patient, specifically
Data and calculating process are as follows：
Iodine contrast medium injected slurry volume=64mL
Iodine contrast medium concentrations=370mgI/mL
Bulb voltage=100kV
Conversion coefficient=30.00 (HU*ml/mgI) of iodine contrast medium and CT values
Left kidney
Volume=179.7ml
HU0=40.3Hu
HU1=189.8Hu
Right kidney
Volume=144.3ml
HU0=32.5Hu
HU1=111.8Hu
The calculating process of FRA is as follows：
FRA (left side)=(volume × HU1- volume × HU0)/(iodine contrasts injected slurry volume × iodine contrast medium concentrations × iodine contrast medium
With the conversion coefficient of CT values)
=(179.7 × 189.8-179.7 × 40.3)/(64 × 370 × 30.00)
=0.0378
FRA (right side)=(volume × HU1- volume × HU0)/(iodine contrasts injected slurry volume × iodine contrast medium concentrations × iodine contrast medium
With the conversion coefficient of CT values)
=(144.3 × 111.8-144.3 × 32.5)/(64 × 370 × 30.00)
=0.0161
The computing formula of GFR is：CT-GFR=1416.42 × FRA+2.58
The calculating process of left kidney GFR is as follows：
CT-GFR (left side)=1416.42 × FRA (left side)+2.58
=1416.42 × 0.0378+2.58
=56.12 (ml/min)
The calculating process of right kidney GFR is as follows：CT-GFR (right side)=1416.42 × FRA (right side)+2.58
=1416.42 × 0.0161+2.58
=25.38 (ml/min)
【Contrast embodiment】
With the 99mTc-DTPA agent as a comparison of radiosiotope medicine, kidney dynamic scan is carried out with SPECT, used
Gates methods are calculating renal function (Gates-GFR).Specific operating process is referring to article Determination of
Glomerular Filtration Rate with CT Measurement of Renal Clearance of
Iodinated Contrast Material versus 99mTc-DTPA Dynamic Imaging"Gates"Method:A
Validation Study in Asymmetrical Renal Disease.Radiology 2016Aug 24:
160425.Yuan X, Zhang J, Tang K are waited.
【First embodiment and the interpretation for contrasting embodiment】
Compare the CT measuring methods (CT-GFR) and clinic of detection kidney (filtration) function of the invention referring to Fig. 6-7
Conventional use of ING Gates methods (Gates-GFR) are verifying the accuracy of CT-GFR.
Fig. 6 is based on the correlation analysiss between list kidney CT-GFR and single kidney Gates-GFR.The scatterplot confirms between the two
Dependency (r=0.924, p with height<0.001).
Can draw from the analysis result of Fig. 6：The result and Gates-GFR height correlations of CT-GFR.
Fig. 7 is CT-GFR and Gates-GFR result consistency analysis figures (Bland-Altman figure).Abscissa represents two kinds
The mean of measurement result, vertical coordinate represents the difference of two kinds of measurement results, and the figure shows that the difference of two kinds of measurement results is controlled
(it is distributed between two horizontal dotted lines) in rational scope.
Can draw from the analysis result of Fig. 7：CT-GFR and Gates-GFR survey reason deviations are little, and concordance is high.
The present invention passes through to propose a kind of CT measuring methods of new detection kidney filtering function, using conventional CT enhanced ct scans
Data information assessment renal function.The method can obtain accumulation fraction FRA of the iodine contrast medium in kidney, and the FRA can be straight
The pointer as assessment kidney filtering function GFR height is connect, or linear regression model (LRM) is set up on the basis of FRA calculating kidney
Glomerular filtration rate GFR.The method simple, intuitive, while the conventional CT of being expert at is diagnosed can quantitative evaluation renal function, and not
The x beta radiation amounts of patient can additionally be increased, clinically there is very important value.
Claims (10)
1. it is a kind of detection kidney filtering function CT measuring methods, it is characterised in that comprise the following steps：
To start to start contrast medium occur as time zero, Jing Guogui in injection of contrast agent or injection of contrast agent aorta posterior
The moment t for fixing time1, CT scan is carried out to kidney and obtains t1The enhanced CT data of moment kidney, the data include constituting kidney three
All voxels of dimension image, can therefrom obtain the CT values of each voxel, that is, represent contrast medium concentrations, and all voxel volumes
Sum, that is, represent the volume (V) of kidney；
T is obtained according to the enhanced CT data of the kidney1Contrast medium accumulation amount (Q) of moment kidney；
By the accumulation amount (Q) of contrast medium in the kidney divided by total contrast medium injection rate (I), the accumulation fraction of kidney is obtained
(FRA), i.e. FRA=Q/I；
The accumulation fraction (FRA) of the kidney can be directly as the index of the unilateral kidney filtering function height of assessment.
2. it is according to claim 1 detection kidney filtering function CT measuring methods, it is characterised in that
Methods described also includes setting up linear regression model (LRM) to calculate glomerular filtration rate based on the accumulation fraction (FRA) of kidney
(GFR)。
3. it is according to claim 2 detection kidney filtering function CT measuring methods, it is characterised in that
When injection of contrast agent will be started as time zero, set up linear regression model (LRM) to calculate based on the accumulation fraction of kidney
The computing formula of unilateral renal glomerulus filtration rate (GFR) is as shown in Equation 1：
GFR=1416.42 × FRA+2.58 ... (formula 1).
4. it is according to claim 1 detection kidney filtering function CT measuring methods, it is characterised in that
The moment t1Span be after injection of contrast agent, or to start 2-3 minutes after contrast medium occur in aorta；
Preferably, the moment t1Value be or to start to occur 2.5 minutes after contrast medium in aorta after injection of contrast agent.
5. it is according to claim 1 detection kidney filtering function CT measuring methods, it is characterised in that
The injection control of contrast medium was completed in 20 seconds；Preferably, the injection control of contrast medium was completed in 15 seconds；
Before injection of contrast agent, also include carrying out kidney the unenhanced steps of CT, to obtain the plain CT data of kidney, the number
According to including all voxels for constituting kidney 3-D view, it is zero that can therefrom obtain the CT values of each voxel, i.e. contrast medium concentrations
When CT values background；
The t1Contrast medium accumulation amount (Q) of moment kidney is by the t1All voxels in the enhanced CT data of moment kidney
Mean CT-number (HU1) deduct the mean CT-number (HU of all voxels in the plain CT data0), then the volume (V) for being multiplied by kidney
Obtain,
Or by the t1CT value (the HU of each voxel in the enhanced CT data of moment kidneyi) and its volume (Vvoxel)
The sum of products, deducts the CT value (HU of each voxel in the plain CT datak) and its volume (Vvoxel) the sum of products obtain.
6. it is according to claim 1 detection kidney filtering function CT measuring methods, it is characterised in that
When the CT scan is dual intensity CT scan or power spectrum CT scan,
The t1Contrast medium accumulation amount (Q) of moment kidney subtracts shadow and obtains by dual energy subtraction or power spectrum, including following mistake
Journey：
When dual intensity CT scan or power spectrum CT scan, it is not necessary to injection contrast medium before kidney is carried out it is unenhanced, directly by t1Moment
Dual energy subtraction or power spectrum subtract the incrementss (HU that shadow obtains each voxel CT values in kidneyy), it is multiplied by the volume (V of each voxelvoxel)
And summation obtains contrast medium accumulation amount (Q) of kidney,
Or by the incrementss (HU of each voxel CT valuesy) obtain the average increase (HU of all voxel CT values in kidneym), then
The volume (V) for being multiplied by kidney obtains contrast medium accumulation amount (Q) of kidney.
7. it is according to claim 1 detection kidney filtering function CT measuring methods, it is characterised in that
Total contrast medium injection rate (I) is to be multiplied by contrast medium concentrations by the volume of contrast medium to be multiplied by contrast medium and CT values again
Conversion coefficient (k) obtain.
8. it is according to claim 1 detection kidney filtering function CT measuring methods, it is characterised in that
Total contrast medium injection rate can also be represented by the volume of contrast medium and the product of contrast medium concentrations, calculate poly-
In order to unify dimension during fraction (FRA), contrast medium accumulation amount (Q) of the kidney should be divided by the conversion of contrast medium and CT values
Coefficient (k).
9. according to claim 7 or 8 detection kidney filtering function CT measuring methods, it is characterised in that
The contrast medium is iodine, and the height of the bulb magnitude of voltage of CT scan affects the conversion coefficient of contrast medium iodine concentration and CT values
K (), when bulb voltage is a certain fixed value, the conversion coefficient relative constancy can be determined by the scanning of CT body mould.
10. it is according to claim 9 detection kidney filtering function CT measuring methods, it is characterised in that
When the bulb voltage of CT scan is 120KV, the conversion coefficient (k) of contrast medium iodine concentration and CT values is 25, and unit is HU*
ml/mgI；
When the bulb voltage of CT scan is 100KV, the conversion coefficient (k) of contrast medium iodine concentration and CT values is 30, and unit is HU*
ml/mgI。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611159034.0A CN106667513A (en) | 2016-12-15 | 2016-12-15 | CT measurement method for detecting renal filtration function |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611159034.0A CN106667513A (en) | 2016-12-15 | 2016-12-15 | CT measurement method for detecting renal filtration function |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106667513A true CN106667513A (en) | 2017-05-17 |
Family
ID=58868997
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611159034.0A Pending CN106667513A (en) | 2016-12-15 | 2016-12-15 | CT measurement method for detecting renal filtration function |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106667513A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109350102A (en) * | 2017-11-28 | 2019-02-19 | 北京昆仑医云科技有限公司 | For determining the method, apparatus and system of Vascular capacitance and coronary blood flow |
| CN109545377A (en) * | 2018-10-22 | 2019-03-29 | 中南大学湘雅三医院 | Obtain method for building up and the application of the model of glomerular filtration rate |
| CN110702888A (en) * | 2019-09-19 | 2020-01-17 | 中国人民解放军总医院第八医学中心 | Liver and kidney clearance detection method for two isomers of Pumei |
| CN112700875A (en) * | 2021-01-29 | 2021-04-23 | 中山大学肿瘤防治中心（中山大学附属肿瘤医院、中山大学肿瘤研究所） | System for predicting kidney function after kidney tumor operation and computer readable storage medium |
| CN113069146A (en) * | 2021-03-22 | 2021-07-06 | 莆田学院附属医院 | Evaluation method for unilateral hydronephrosis of children patients |
| US11369278B2 (en) | 2017-11-28 | 2022-06-28 | Beijing Keya Medical Technology Co., Ltd. | Methods and systems for determining lumen volume and coronary blood flow |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1923140A (en) * | 2006-07-27 | 2007-03-07 | 沈阳东软医疗***有限公司 | Image measuring device, method and image instrumentation system of glomerular filtration rate |
| CN102065904A (en) * | 2008-04-18 | 2011-05-18 | 药物影像股份有限公司 | Renal function analysis method and apparatus |
-
2016
- 2016-12-15 CN CN201611159034.0A patent/CN106667513A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1923140A (en) * | 2006-07-27 | 2007-03-07 | 沈阳东软医疗***有限公司 | Image measuring device, method and image instrumentation system of glomerular filtration rate |
| CN102065904A (en) * | 2008-04-18 | 2011-05-18 | 药物影像股份有限公司 | Renal function analysis method and apparatus |
Non-Patent Citations (5)
| Title |
|---|
| XIAODONG YUAN ET AL.: "A Simplified Whole-Organ CT Perfusion Technique with Biphasic Acquisition: Preliminary Investigation of Accuracy and Protocol Feasibility in Kidneys", 《RADIOLOGY》 * |
| 周科峰等: "16 排螺旋CT 测定单侧肾小球滤过率的初步探讨", 《医学影像学杂志》 * |
| 汪依民: "《生理学》", 30 September 2005 * |
| 秦将均等: "心电门控模式下能谱扫描和传统扫描CT 值与碘对比剂浓度相关性的比较", 《放射学实践》 * |
| 钟建国等: "螺旋CT 多期扫描测定肾小球滤过率不同扫描技术的比较", 《医学影像学杂志》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109350102A (en) * | 2017-11-28 | 2019-02-19 | 北京昆仑医云科技有限公司 | For determining the method, apparatus and system of Vascular capacitance and coronary blood flow |
| US11369278B2 (en) | 2017-11-28 | 2022-06-28 | Beijing Keya Medical Technology Co., Ltd. | Methods and systems for determining lumen volume and coronary blood flow |
| CN109545377A (en) * | 2018-10-22 | 2019-03-29 | 中南大学湘雅三医院 | Obtain method for building up and the application of the model of glomerular filtration rate |
| CN109545377B (en) * | 2018-10-22 | 2021-11-23 | 中南大学湘雅三医院 | Method for establishing model for obtaining glomerular filtration rate |
| CN110702888A (en) * | 2019-09-19 | 2020-01-17 | 中国人民解放军总医院第八医学中心 | Liver and kidney clearance detection method for two isomers of Pumei |
| CN112700875A (en) * | 2021-01-29 | 2021-04-23 | 中山大学肿瘤防治中心（中山大学附属肿瘤医院、中山大学肿瘤研究所） | System for predicting kidney function after kidney tumor operation and computer readable storage medium |
| CN113069146A (en) * | 2021-03-22 | 2021-07-06 | 莆田学院附属医院 | Evaluation method for unilateral hydronephrosis of children patients |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106667513A (en) | CT measurement method for detecting renal filtration function | |
| CN100463655C (en) | Image measuring device, method and image instrumentation system of glomerular filtration rate | |
| Ferrier et al. | Validation of dynamic contrast-enhanced magnetic resonance imaging-derived vascular permeability measurements using quantitative autoradiography in the RG2 rat brain tumor model | |
| JP5775310B2 (en) | Computer-based method and system for imaging-based assessment of organ dynamic function | |
| Hwang et al. | Noninvasive and direct measures of kidney size in kidney donors | |
| Gomberg et al. | Method for Cortical Bone Structural Analysis From Magnetic Resonance Images1 | |
| WO2014167325A1 (en) | Methods and apparatus for quantifying inflammation | |
| Satkunasingham et al. | Liver fat quantification: comparison of dual-echo and triple-echo chemical shift MRI to MR spectroscopy | |
| Wang et al. | Reliability and validity of abdominal skeletal muscle area measurement using magnetic resonance imaging | |
| Zakhari et al. | Simplified method to measure renal volume: the best correction factor for the ellipsoid formula volume calculation in pretransplant computed tomographic live donor | |
| Nilsson et al. | Measuring split renal function in renal donors: can computed tomography replace renography? | |
| Smith et al. | Obesity and visceral fat in Indonesia: An unseen epidemic? A study using iDXA and surrogate anthropometric measures | |
| Zeng et al. | Measurement of single-kidney glomerular filtration function from magnetic resonance perfusion renography | |
| Jones et al. | MRU post-processing | |
| Pasha et al. | Quantitative imaging of the spine in adolescent idiopathic scoliosis: shifting the paradigm from diagnostic to comprehensive prognostic evaluation | |
| Hamed | New advances in assessment of the individual renal function in chronic unilateral renal obstruction using functional CT compared to 99mTc-DTPA renal scan | |
| US7146204B2 (en) | Method and apparatus for quantitatively evaluating a kidney | |
| Evrimler et al. | Quantitative assessment of disease severity of primary sclerosing cholangitis with T1 mapping and extracellular volume imaging | |
| Eckstein | Noninvasive study of human cartilage structure by MRI | |
| Tipirneni‐Sajja et al. | Measurement of glomerular filtration rate by dynamic contrast‐enhanced magnetic resonance imaging using a subject‐specific two‐compartment model | |
| Yap et al. | Bio-conductivity characteristics of chronic kidney disease stages examined by portable frequency-difference electrical impedance tomography | |
| Fonseca-Matheus et al. | Contrast-enhanced dynamic magnetic resonance nephrography in healthy dogs | |
| Michalski et al. | QCT Reconstruction Kernel has Important Quantitative Effects on Finite Element Estimated Bone Strength | |
| Hopper et al. | Comparison of high-resolution MRI, optical microscopy and SEM for quantitation of trabecular architecture in the rat femur | |
| Cai et al. | Semi-automatic segmentation of renal cortex and medulla based on dynamic magnetic resonance images |
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 |