CN103126704A - Equipment and method of image data correction - Google Patents

Equipment and method of image data correction Download PDF

Info

Publication number
CN103126704A
CN103126704A CN2011104054219A CN201110405421A CN103126704A CN 103126704 A CN103126704 A CN 103126704A CN 2011104054219 A CN2011104054219 A CN 2011104054219A CN 201110405421 A CN201110405421 A CN 201110405421A CN 103126704 A CN103126704 A CN 103126704A
Authority
CN
China
Prior art keywords
image detecting
view data
detecting element
corrected
detecting unit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2011104054219A
Other languages
Chinese (zh)
Other versions
CN103126704B (en
Inventor
曹蹊渺
王学礼
H·江
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to CN201110405421.9A priority Critical patent/CN103126704B/en
Priority claimed from CN201110405421.9A external-priority patent/CN103126704B/en
Publication of CN103126704A publication Critical patent/CN103126704A/en
Application granted granted Critical
Publication of CN103126704B publication Critical patent/CN103126704B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

The invention relates to a system and a method of imaging. A method of image data correction comprises obtaining image data of an adjacent image detecting unit of a to-be-corrected image detecting unit in an image detecting unit array; carrying out interpolation on the to-be-corrected image detecting unit according to the image data of the adjacent image detecting unit of the to-be-corrected image detecting unit, and obtaining estimated image data of the to-be-corrected image detecting unit; judging whether a step change point of the image date exists in the image data of the adjacent image detecting unit of the to-be-corrected image detecting unit; and selecting various correcting modes on the basis of a judgment result so as to carry out correction on the estimated image data.

Description

Image data correction equipment and method
Technical field
Embodiments of the invention relate generally to imaging system and method, are specifically related to equipment and method on the picture quality impact for reducing fault picture detecting unit in image detecting element array in imaging system.
Background technology
Medical diagnostic imaging system comprises various imaging forms, such as plane x ray, ultrasonic, magnetic resonance (MR), electron beam x line laminagraphy (EBT), positron emission tomography take a picture (PET), single photon emission computed tomography is taken a picture (SPECT), pico computer tomography (micro CT), grand computed tomography (macro CT) imaging system etc.Medical diagnostic imaging system carrys out the image of formation object by being exposed to such as patient's object energy source (such as the x ray by the patient).The image generated can be for many purposes.For example, internal flaw that can detected object.In addition, can determine the variation in aligning or internal structure.Fluid stream in all right expressive object.And, image can display object in parts existence or do not exist.The information obtained from diagnosis imaging is applied to many fields, comprises medical treatment and manufacture field.
Typical imaging system adopts the image detecting element array to carry out detected object, the image that then reconstruct and demonstration detect.Above-mentioned array comprises a plurality of image detecting elements capable (row).The capable a plurality of image detecting elements (cell) that comprise of each image detecting element, each image detecting element is connected to different data collecting system (DAS) passage (channel).That is to say, the DAS passage can be mapped to image detecting element.Each image detecting element produces signal.Jumbo array comprises image detecting element and the DAS passage of large quantity.Due to the quantity increase of image detecting element and DAS passage, the probability of image detecting element, DAS channel failure or DAS special IC (ASIC) fault increases.In addition, due to the quantity increase of image detecting element and DAS passage, wish that image detecting element array and imaging system assembly become more integrated.
The artifact (artifact) of problem (problem) in image detecting element array or DAS passage in any one image detecting element in may causing reconstructed image.In image detecting element array and/or DAS passage, problematic image detecting element is called fault picture detecting unit (malfunctioning cell).The fault picture detecting unit may be with several different modes faults, such as fault off and on, cause the signal lower than a certain percentage ratio of other image detecting element signal accuracy, produce the signal that is weaker than a certain percentage ratio of adjacent image detecting unit signal, and fully inoperative.Any mistake or " artifact " in the image that image detecting element produces in imaging system or imaging system may cause doctor, medical science practitioner or other observer to take action based on error message.
The capacity that increases the image detecting element array allows the user to carry out quickly imaging to larger object, because once scan (sweep) of larger array carries out imaging to the greater part of object.The capacity that increases the image detecting element array also carries out imaging to object more accurately, because disappear the less time during imaging when adopting less number of times to carry out imaging to object.Along with the pursuit to larger capacity, image detecting element array and DAS number of channels increase sharply.As a result, increase the probability of fault picture detecting unit.Replacing the fault picture detecting unit has increased the cost of system significantly.All faulty components that replacement has on the system of great amount of images array of detection units and DAS passage may be uneconomical.In addition, faulty components has been interrupted the workflow in the hospital.Therefore, wish very much the deterioration of image quality or the minimized system of image artifacts that cause due to the fault picture detecting unit.
Summary of the invention
An aspect, provide a kind of method for the treatment of view data according to an embodiment of the invention, comprising: the image data acquisition step, obtain the view data of the adjacent image detecting unit of image detecting element to be corrected in the image detecting element array; The view data estimating step, the identical line adjacency graph of going together mutually with described image detecting element to be corrected in described view data based on described adjacent image detecting unit carries out interpolation as the view data of detecting unit to described image detecting element to be corrected, obtains the estimated image data of described image detecting element to be corrected; And view data correction step, the characteristic of the view data based on described adjacent image detecting unit is revised described estimated image data.
Another aspect, provide a kind of equipment for the treatment of view data according to an embodiment of the invention, comprising: the image data acquisition device, obtain the view data of the adjacent image detecting unit of image detecting element to be corrected in the image detecting element array; The view data estimation unit, the identical line adjacency graph of going together mutually with described image detecting element to be corrected in described view data based on described adjacent image detecting unit carries out interpolation as the view data of detecting unit to described image detecting element to be corrected, obtains the estimated image data of described image detecting element to be corrected; And image data correction apparatus, the characteristic of the view data based on described adjacent image detecting unit is revised described estimated image data.
Further, the method and apparatus according to the above-mentioned aspect of embodiments of the invention also comprises: judge the step point that whether has view data in the view data of adjacent image detecting unit of described image detecting element to be corrected; And based on described judged result, described estimated image data are proofreaied and correct.
Further, while having the step point of view data in the view data of described adjacent image detecting unit, according to the first correction mode, described estimated image data are proofreaied and correct.
Further, described the first correction mode comprises: along the adjacent image detecting unit of the set direction described to be corrected image detecting element different from the row of described image detecting element array, described direction has reflected the distribution arrangement of step point in the view data of described adjacent image detecting unit; Calculate the deviation of directivity Δ of the view data of selected adjacent image detecting element, and proofread and correct described estimated image data based on described deviation of directivity Δ, obtain the view data after directivity is proofreaied and correct.
Further, the deviation of directivity Δ of view data that calculates the adjacent image detecting unit of described image detecting element to be corrected comprises: the difference of calculating the estimated image data of the view data of selected adjacent image detecting element and described image detecting element to be corrected and calculate described difference
Figure BSA00000632350600032
meansigma methods as described deviation of directivity Δ.
Further, described the first correction mode comprises: described image detecting element to be corrected is carried out to monolateral extrapolation, obtain the view data after the monolateral extrapolation of described image detecting element to be corrected.
Further, described image detecting element to be corrected being carried out to monolateral extrapolation comprises: calculate the left side weight w of described image detecting element to be corrected in described image detecting element array leftwith right side weight w right, and according to described left side weight w leftwith right side weight w rightcalculate the view data X:x=w after monolateral extrapolation left* x left+ w right* x right, wherein, x leftand x rightit is respectively the view data of described image detecting element left side adjacent image detecting unit to be corrected and right side adjacent image detecting unit.
Further, calculate in such a way described left side weight w leftwith described right side weight w right:
For described image detecting element to be corrected, the extrapolated value x of lastrow left side extrapolation estimation lastrow up_left, the extrapolated value x of lastrow right side extrapolation estimation lastrow up_left, the extrapolated value x of next line left side extrapolation estimation lastrow up_left, the extrapolated value x of next line right side extrapolation estimation lastrow up_left, wherein, described extrapolated value is to adopt the view data of adjacent image detecting unit of the respective side of corresponding line to carry out Lagrange's interpolation, and
By described extrapolated value x up_left, x up_left, x up_left, x up_leftwith the upper adjacent image detecting unit of described image detecting element to be corrected and the view data x of lower adjacent image detecting unit up, x downsolve following equation and obtain right outside slotting weight w up_left, w down_leftwith right outside slotting weight w up_right, w down right:
x up=w up_left*x up_left+w up_right*x up_right
x down=w down_left*x down_left+w down_right*x down_right
W up_left+ w up_right=1, and
w down_left+w down_right=1
Pass through obtained left outside slotting weight w up_left, w down_left, and right outside slotting weight w up_right, w down rightcalculate described left side weight w leftwith described right side weight w right:
w left=a*w up_left+(1-a)*w down_left
w right=bw up_right+(1-b)*w down?right
Wherein, a and b are the constant of value between 0 to 1.
Further, the view data of view data after the correction of described image detecting element to be corrected after described directivity is proofreaied and correct.
Further, the view data after adopting view data after described monolateral extrapolation as the correction of described image detecting element to be corrected; Perhaps calculate the deviation of directivity Δ of view data of the selected adjacent image detecting element of described image detecting element to be corrected, be greater than at 0 o'clock at described deviation of directivity Δ, view data after selecting view data after described directivity is proofreaied and correct and the maximum in the view data after described monolateral extrapolation as the correction of described image detecting element to be corrected, otherwise, the view data after selecting view data after described directivity is proofreaied and correct and the reckling in the view data after described monolateral extrapolation as the correction of described image detecting element to be corrected.
Alternatively, while not having the step point of view data in the view data of described adjacent image detecting unit, according to the second correction mode, described estimated image data are proofreaied and correct.
Further, described the second correction mode comprises: the meansigma methods of calculating the deviation of the view data of adjacent image detecting unit of described image detecting element to be corrected and described estimated image data, meansigma methods based on described deviation is proofreaied and correct described estimated image data, obtains the view data after the correction of described image detecting element to be corrected.
Further, the meansigma methods of the view data of the adjacent image detecting unit of the described image detecting element to be corrected of described calculating and the deviation of described estimated image data comprises: calculate the meansigma methods of the arithmetic difference of the view data of the upper and lower, left and right of described image detecting element to be corrected, upper left, upper right, bottom left, bottom right adjacent image detecting unit and described estimated image data, as the meansigma methods of described deviation.
Further, the meansigma methods of the view data of the adjacent image detecting unit of the described image detecting element to be corrected of described calculating and the deviation of described estimated image data comprises: calculate the meansigma methods of the arithmetic difference of the view data of upper and lower adjacent image detecting unit of described image detecting element to be corrected and described estimated image data, as the meansigma methods of described deviation.
Further, the meansigma methods of the deviation of the view data of the adjacent image detecting unit of the described image detecting element to be corrected of described calculating and described estimated image data comprises: judge whether the size of image detecting element column direction in described image detecting element array and the ratio R of the size of line direction are more than or equal to predetermined threshold V thresholdif R is more than or equal to described predetermined threshold V thresholdcalculate the meansigma methods of the arithmetic difference of the view data of the upper and lower, left and right of described image detecting element to be corrected, upper left, upper right, bottom left, bottom right adjacent image detecting unit and described estimated image data, meansigma methods as described deviation, otherwise, calculate the meansigma methods of the arithmetic difference of the view data of upper and lower adjacent image detecting unit of described image detecting element to be corrected and described estimated image data, as the meansigma methods of described deviation.
Further, wherein said described predetermined threshold V thresholdthe numerical value between 0.5 to 2.
Further, the estimated image data that obtain described image detecting element to be corrected comprise: the view data of the adjacent image detecting unit of being expert at based on described image detecting element to be corrected is carried out Lagrange's interpolation to described image detecting element to be corrected, obtains described estimated image data.
Below in conjunction with accompanying drawing, embodiments of the invention are described.
The accompanying drawing explanation
In the accompanying drawings:
Fig. 1 is the diagram that can apply the imaging system of embodiments of the invention.
Fig. 2 is the schematic block diagram of the system shown in Fig. 1.
Fig. 3 is the diagram that the image detecting element array that comprises fault picture detecting unit and fault-free image detecting element is described.
Fig. 4 is the flow chart for the method for proofreading and correct fault picture detecting unit of explanation according to one embodiment of the present of invention.
Fig. 5 a and 5b are the diagrams of the step point of explanation view data.
Fig. 6 is the flow chart of determination methods that whether has the step point of view data in the view data of explanation according to the adjacent image detecting unit of the fault picture detecting unit of one embodiment of the present of invention.
Fig. 7 is the diagram of directivity of the step point of explanation view data.
Fig. 8 is for the diagram of monolateral extrapolation is described.
Fig. 9 is the flow chart of explanation according to the picture data correction method of the fault picture detecting unit of a preferred embodiment of the present invention.
Figure 10 is the block diagram for the equipment of proofreading and correct fault picture detecting unit of explanation according to one embodiment of the present of invention.
Figure 11 shows according to embodiments of the invention and carries out image ratio before and after image data correction.
The specific embodiment
Some embodiments of the present invention are provided for the system and method that the mistake in image detecting element, DAS passage and/or DAS ASIC is carried out to error correction.
For the purpose of illustrating, below describe some embodiment that relate to the CT imaging system in detail.One of skill in the art will appreciate that embodiments of the invention be applicable to equally other imaging system (such as plane x ray, ultrasonic, magnetic resonance (MR), electron beam x line laminagraphy (EBT), positron emission tomography take a picture (PET), single photon emission computed tomography is taken a picture (SPECT), pico computer tomography (micro CT), grand computed tomography (macro CT) imaging system etc.)
Fig. 1 is the diagram that can apply the imaging system 10 of embodiments of the invention.Those skilled in the art will recognize that the present invention be applicable to such as plane x ray, ultrasonic, magnetic resonance (MR), electron beam tomography photography (EBT), positron emission tomography take a picture (PET), single photon emission computed tomography is taken a picture (SPECT), pico computer tomography (micro CT), grand computed tomography (macro CT) imaging system etc.The following discussion of x ray system 10 is only the example of such realization and is unexpectedly restrictive from the form aspect.
With reference to Fig. 1, computer tomography (CT) imaging system 10 comprises frame 12.Frame 12 has x ray tube assembly or x radiographic source assembly 14, and it is projected in the cone beam of x ray on the opposite side of this frame 12 towards image detecting element assembly or collimator 16.Referring now to Fig. 2, image detecting element assembly 16 is formed by image detecting element array 18 and data collecting system (DAS) 20.These image detecting element array 18 sensings are by the x ray 22 of the projection of medical patient 24, and DAS 20 becomes digital signal for subsequent treatment data transaction.Each image detecting element 18 generation analog electrical signal, the intensity of its representative collision x beam and the intensity that therefore represents the attenuated beam when it passes through this patient 24.During scanning gathers the x ray projection data, frame 12 and parts mounted thereto rotate around center of rotation 26.
The operation of the rotation of frame 12 and x radiographic source assembly 14 is by controlling organization 28 management of CT system 10.Controlling organization 28 comprises: x ray controller 30, and it provides electric power and timing signal to x radiographic source assembly 14; With frame motor controller 32, it controls rotating speed and the position of frame 12.Image reconstructor 34 from DAS 20 receive sampling with digitized x ray data and carry out high speed reconstruction.The image of rebuilding offers computer 36 (it is stored in this image in mass storage device 38).Computer 36 also has the software corresponding to positioning of beam and magnetic field control stored thereon, as described in detail hereinafter.
Computer 36 also receives order and sweep parameter via control station 40 from the operator, and this control station 40 has such as the controller of keyboard, mouse, voice-activated or the operator interface of any other suitable certain form such as input equipment.Associated display 42 allows this operator to observe image and other data from the reconstruction of computer 36.The order that this operator provides and parameter are used for providing control signal and information to DAS 20, x ray controller 30 and frame motor controller 32 by computer 36.In addition, computer 36 operation stand motor controllers 44, it is controlled electronic stand 46 and settles patient 24 and frame 12.Especially, stand 46 makes the frame openings 48 of patient 24 all or part of Fig. 1 of moving through.
Fig. 3 is the diagram that the image detecting element array 300 that comprises fault picture detecting unit (or image detecting element to be corrected) and fault-free image detecting element is described.Image detecting element array 300 has image detecting element capable (in figure, horizontal direction is extended) and image detecting element row (in figure, vertical direction extends).Image detecting element array 300 comprises fault picture detecting unit 330, is positioned at the cross point of row 310 and row 320.Fault picture detecting unit 330 may not can be reported view data accurately.The view data of the image detecting element around fault picture detecting unit 330 can be used for the view data of suspected fault image detecting element 330.
What Fig. 4 was explanation according to one embodiment of the present of invention is positioned at for correction the flow chart of method 400 of view data that the image detecting element array is positioned at the fault picture detecting unit 330 of row 310 and row 320.
In step 410, input is corresponding to the view data of each image detecting element in image detecting element array 300.
In step 420, the adjacent image detecting unit based on fault picture detecting unit 330 carries out interpolation, obtains the estimated image data (Xbase) of fault picture detecting unit 330.According to one embodiment of the present of invention, based on fault picture detecting unit 330, with the view data of the adjacent image detecting unit of the fault picture detecting unit 330 of a line 310, carrying out interpolation.In one embodiment of the invention, can based on the view data of fault picture detecting unit 330 with the adjacent image detecting unit of the fault picture detecting unit 330 of a line 310, adopt Lagrange's interpolation, obtain the estimated image data of fault picture detecting unit 330.
In an embodiment of the present invention, the adjacent image detecting unit of fault picture detecting unit comprises the adjacent image detecting unit be close to most of fault picture detecting unit, the adjacent image detecting unit that perhaps also comprises further time next-door neighbour of fault picture detecting unit, the adjacent image detecting unit that perhaps also comprises further time time next-door neighbour of fault picture detecting unit, etc.
Can adopt several different methods to locate the fault picture detecting unit.For example, fault picture detecting unit 330 can adopt such as at least one method of following several method and locate: (1) is exposed to all image detecting elements in the image detecting element array under the x ray, then measuring with the adjacent image detecting unit of image detecting element of each image detecting element compared to find signal intensity (variations), and (2) find the change (changes) (capable all emissions are observed for each image detecting element) in the average signal (data) between the adjacent image detecting unit that can identify the fault picture detecting unit.
First kind of way is called air calibration scanning (air-calibration scan).In one embodiment, each imaging system experiences one group of quickly calibrated processing (being called " air calibration ").Air calibration is processed image detecting element directly is exposed to x ray photons beam, and without any object in beam.By checking the different reading obtained from adjacent image detecting unit 340, fault picture detecting unit 330 can identify soon.Along with the past of time, fault picture detecting unit 330 may some the time suitably move, and other the time inappropriate operation.Fault picture detecting unit 330 can identify into the image detecting element as one group of percentage ratio below the reading of adjacent image detecting unit 340,350.Above-mentioned percentage ratio can also change according to the position of image detecting element.
For the second way that identifies the fault picture detecting unit, be to use from start to finish the patient scan data.The advantage of the second way is: the operator who is not all hospitals or imaging system carries out air calibration every day.Therefore, if mistake is introduced image detecting element between air calibration scanning, image detecting element may suitably not identified.For example, by using average reading (capable all emissions are observed for each image detecting element), relatively the average signal (data) between the adjacent image detecting unit can identify fault picture detecting unit 330.The emission of image detecting element is observed and on average can be kept the reading of each signal to approach another, unless the image detecting element fault.
Air calibration mode and patient scan mode can be implemented by the hardware separated or the software that are positioned at each position, such as at image detecting element array 18, DAS 32, image reconstructor 34, computer 36 etc.The image detecting element array worsens in time and can be checked through.Carry out each the inspection with sign fault picture detecting unit 330.Change in the image detecting element reading can be stored.Change in can check image detecting unit reading significantly changes to search in image detecting element array 18.Can identify change, thus the image detecting element array 18 that identifies fault picture detecting unit 330 and will replace potentially.Can keep the daily record over all variations that detected of one group of threshold value, thus sign fault picture detecting unit 330.This daily record can keep by the hardware separated or the software that is positioned at each position, for example, at arrays 18, DAS 32, image reconstructor 34, computer 36, or for example remains in the hardware for imaging system 10 preexists.
According to embodiments of the invention, the characteristic of the view data of the adjacent image detecting unit of failure judgement image detecting element 330, the characteristic based on this view data is proofreaied and correct the estimated image data (Xbase) of fault picture detecting unit 330.For example, the step response of the view data of the adjacent image detecting unit based on fault picture detecting unit 330 is proofreaied and correct the estimated image data (Xbase) of fault picture detecting unit 330.Embodiments of the invention are not limited to the step response of the view data of the adjacent image detecting unit based on fault picture detecting unit 330 the estimated image data (Xbase) of fault picture detecting unit 330 are proofreaied and correct, and other characteristics of the view data of adjacent image detecting unit that also can be based on fault picture detecting unit 330 are proofreaied and correct the estimated image data (Xbase) of fault picture detecting unit 330.
In step 430, whether there is the step point of view data in the view data of the adjacent image detecting unit of failure judgement image detecting element.The step point feature of view data is: in the image detecting element array, the relevant view data of the left-side images detecting unit of a certain image detecting element of a certain image detecting element in capable and image right detecting unit is relative steady, but the view data of left and right sides has obvious numerical value difference.In the CT projecting image data, the step point of typical view data appears at the edge of high-density matter.If the shape of high-density matter slowly changes, projecting image data also can present the variation of " slowly continuously ", referring to the situation shown in Fig. 5 a.If the shape of high-density matter has corner angle, the step point of view data appears possibly in projecting image data, referring to the situation shown in Fig. 5 b.When there is the step point of view data in the view data of the adjacent image detecting unit of fault picture detecting unit, mean the step point that the actual image data of this fault picture detecting unit may be in view data.The present inventor has insight into, and for the step point of view data, traditional interpolation method is difficult to recover its real view data.The low order interpolation is to running into the problem that accuracy is inadequate, and high-order interpolation can be introduced the interpolation overshoot, all can affect the accuracy of interpolation.This inaccurate artifact that is presented as on image, the most common with streak artifact.The introducing meeting of artifact is to clinical generation mistaken diagnosis, and this is on medical image, to do one's utmost to be avoided.The method of embodiments of the invention is considered the step point factor of view data, to the estimated image data of fault picture detecting unit 330, takes different modes to be proofreaied and correct, with the deficiency of the method that overcomes prior art.
Fig. 6 is whether explanation exists the diagram of the step point determination methods of view data according to the view data of the adjacent image detecting unit of the fault picture detecting unit of one embodiment of the present of invention.In Fig. 6, X1, X2, the X3 representative is positioned at the view data of image detecting element C with the adjacent image detecting unit in the left side of a line, X4, X5, the X6 representative is positioned at the view data of image detecting element C with the adjacent image detecting unit on the right side of a line, wherein X3 is the view data of the image detecting element C left side image detecting element the most adjacent with image detecting element C, the view data of the X2 image detecting element that to be image detecting element C left side adjacent with image detecting element C time, X1 is the view data at C time adjacent image detecting element of image detecting element C left side and image detecting element, X4 is the view data of the image detecting element C right side image detecting element the most adjacent with image detecting element C, X5 is the view data of C adjacent image detecting element of fault image detecting element C right side and image detecting element, X6 is the view data of C time adjacent image detecting element of image detecting element C right side and image detecting element.δ-value between the view data of the relevant adjacent image detecting unit that computed image detecting unit C is expert at: δ 1=|X2-X3|, δ 2=|X3-X4|, δ 3=|X4-X5|, if δ is 2>3 * max, (δ 1, δ 3) and δ 2>0.8 * (max (Xi)-min (Xi)), i ∈ [2,4] wherein, think the step point of image detecting element C in view data.The method that judges the step point of view data more than using is carried out the whether step point in view data of the adjacent image detecting unit of failure judgement image detecting element, the step point methods of judgement view data of the present invention is not limited to said method, and those skilled in the art can modify or the view data that adopts other method the to judge image detecting element step point in view data whether to above method.
During step point when the view data of the adjacent image detecting unit of determining fault picture detecting unit 330 in step 420 in view data, enter step 450 and carry out the first bearing calibration, otherwise, enter step 470 and carry out the second bearing calibration.
In step 490, provide from step 450 or from the correction estimated image data of step 470.
In one embodiment of the invention, the first bearing calibration comprises, whether the step point of view data of adjacent image detecting unit that judges fault picture detecting unit in the image detecting element array deposits certain trend in array, for example, shows certain directivity.
Fig. 7 is the diagram of directivity of the step point of explanation view data.The distribution of the step point of the view data of adjacent image detecting unit that can be by fault picture detecting unit in the check image array of detection units determines whether it exists directivity.Fig. 7 a means that lower-left and the upper right adjacent pattern detecting unit of fault graphical detecting unit (central point in figure, the i.e. the darkest piece of gray scale) are step points (piece that gray scale is the most shallow), the directivity of the dotted line reflection step point in Fig. 7 a.Fig. 7 b means that upper right and the below adjacent pattern detecting unit of fault graphical detecting unit (central point in figure, the i.e. the darkest piece of gray scale) are step points (piece that gray scale is the most shallow), the directivity of the dotted line reflection step point in Fig. 7 b.Fig. 7 c means that the upper and lower adjacent pattern detecting unit of fault graphical detecting unit (central point in figure, the i.e. the darkest piece of gray scale) is step point (piece that gray scale is the most shallow), the directivity of the dotted line reflection step point of the connection step point in Fig. 7 c.Embodiments of the invention are not limited to illustrated situation in Fig. 7, the also applicable the present invention of the situation of other reflection step point directivity.
If there is directivity in the step point of the view data of the adjacent image detecting unit of definite fault picture detecting unit in the image detecting element array, centered by the fault picture detecting unit, for example, according to the view data of the adjacent image detecting unit of the set direction fault picture detecting unit that represents described directivity (the adjacent image detecting unit that in Fig. 7, dotted line connects), the view data based on selected adjacent image detecting unit is proofreaied and correct the estimated image data travel direction of the fault picture detecting unit that obtains in step 420.According to one embodiment of the present of invention, calculate the deviation (Δ of the estimated image data of the view data of each image detecting element (i) in selected adjacent image detecting element and fault picture detecting unit i), calculate the deviation (Δ of the estimated image data of the view data of all image detecting elements in selected adjacent image detecting element and fault picture detecting unit i) the average deviation
Figure BSA00000632350600131
by above-mentioned average deviation
Figure BSA00000632350600132
with estimated image data (Xbase) addition of fault picture detecting unit, the correction estimated image data (X) after the directivity that obtains the fault picture detecting unit is proofreaied and correct.
Alternatively, if there is directivity in the step point of the view data of the adjacent image detecting unit of definite fault picture detecting unit in the image detecting element array, for the fault picture detecting unit, carry out monolateral extrapolation correction.Aligning step according to one embodiment of the present of invention is as follows.
(1) adjacent row twocouese extrapolation estimation
Referring to Fig. 8 a, adjacent row twocouese extrapolation estimation comprises the extrapolated value x of the lastrow left side extrapolation estimation lastrow of fault picture detecting unit up_left, the extrapolated value x of lastrow right side extrapolation estimation lastrow up_left, the extrapolated value x of next line left side extrapolation estimation lastrow up_left, the extrapolated value x of next line right side extrapolation estimation lastrow up_left.
x up_left=extra_lagrange(x up_left_1,x up_left_2,x up_left_3);
x up_right=extra_lagrange(x up_right_1,x up_right_2,x up_right_3);
x down_left=extra_lagrange(x down_left_1,x down_left_2,x down_left_3);
x down_right=extra_lagrange(x down_right_1,x down_right_2,x down_right_3);
In formula, extra_lagrange () means that carrying out Lagrange for the parameter in bracket calculates.
(2) weight analysis of adjacent row twocouese
Referring to Fig. 8 b, the view data by extrapolated value and fault picture detecting unit vertical direction point, can calculate the left and right extrapolation and put to obtain weight.
That is, by by extrapolated value x up_left, x up_left, x up_left, x up_left, and the view data x of the adjacent image detecting unit of fault picture detecting unit vertical direction up, x downthe substitution following formula carrys out solving equation,
x up=w up_left*x up_left+w up_right*x up_right
x down=w down_left*x down_left+w down_right*x down_right
w up_left+w up_right=1
w down_left+w down_right=1
Thereby obtain weight w up_left, w down_left, w up_right, and w down right.
(3) one's own profession twocouese weight calculation
Calculate one's own profession twocouese weight w according to following formula leftand w right.
w left=a*w up_left+(1-a)*w down_left
w right=bw up_right+(1-b)*w down?right
Wherein, constant a and b are the numerical value between 0 to 1.
Alternatively, calculate one's own profession twocouese weight w according to following formula leftand w right.
w left=0.5*w up_left+0.5*w down_left
w right=0.5*w up_right+0.5*w down_right
(4) monolateral extrapolation result
The result of monolateral extrapolation considers the extrapolated value of both direction, by following formula, calculates X.
x=w left*x left+w right*x right
Then, use the correction estimated image data (X) after tried to achieve X proofreaies and correct as the monolateral extrapolation of fault picture detecting unit.
According to one embodiment of the present of invention, if there is directivity in the step point of the view data of the adjacent image detecting unit of definite fault picture detecting unit in the image detecting element array, for the travel direction correction simultaneously of fault picture detecting unit and monolateral extrapolation, proofread and correct.Then, the judgement directivity in proofreading and correct, calculate the average deviation
Figure BSA00000632350600151
whether be greater than zero, if the average deviation be greater than zero, get in the correction estimated image data after correction estimated image data after directivity is proofreaied and correct and monolateral extrapolation are proofreaied and correct the maximum as the correction estimated image data (X) of fault picture detecting unit, otherwise, get in the correction estimated image data after correction estimated image data after directivity is proofreaied and correct and monolateral extrapolation are proofreaied and correct reckling as the correction estimated image data (X) of fault picture detecting unit.
In one embodiment of the invention, the second bearing calibration comprises, by the view data of fault picture detecting unit up-downgoing adjacent image detecting unit, carries out the delta correction.Alternatively, calculate fault picture detecting unit delta of the view data of each adjacent image detecting unit in the neighbouring image detecting element on vertical direction (with the difference of the estimated image data Xbase of fault picture detecting unit) in the image detecting element array, to all these calculated delta calculating mean values, by estimated image data (Xbase) addition of this meansigma methods and fault picture detecting unit, obtain the correction estimated image data (X) of fault picture detecting unit.Alternatively, the delta of the view data of each adjacent image detecting unit in the adjacent image detecting unit of calculating fault picture detecting unit upper and lower, left and right, upper left, upper right, lower-left, bottom right, to all these calculated delta calculating mean values, by estimated image data (Xbase) addition of this meansigma methods and fault picture detecting unit, obtain the correction estimated image data (X) of fault picture detecting unit.
According to one embodiment of the present of invention, the second bearing calibration also comprises whether the size that judges image detecting element vertical direction in the image detecting element array is more than or equal to predetermined threshold V with the ratio R of the size of line direction threshold.If R is more than or equal to predetermined threshold V thresholdcalculate the delta of the view data of each adjacent image detecting unit in the adjacent image detecting unit of fault picture detecting unit upper and lower, left and right, upper left, upper right, lower-left, bottom right, to all these calculated delta calculating mean values, by estimated image data (Xbase) addition of this meansigma methods and fault picture detecting unit, obtain the correction estimated image data (X) of fault picture detecting unit.If R is less than predetermined threshold V thresholdcalculate fault picture detecting unit delta of the view data of each adjacent image detecting unit in the neighbouring image detecting element on vertical direction in the image detecting element array, to all these calculated delta calculating mean values, by estimated image data (Xbase) addition of this meansigma methods and fault picture detecting unit, obtain the correction estimated image data (X) of fault picture detecting unit.Predetermined threshold V thresholdcan be in 0.5 to 2 numerical range value.
Fig. 9 is the diagram of explanation according to the picture data correction method of the fault picture detecting unit of a preferred embodiment of the present invention.In step 910, receive the view data of image detecting element array, wherein in the image detecting element array, include fault picture detecting unit or image detecting element to be corrected.In step 920, the view data of the adjacent image detecting unit of being expert at the image detecting element array based on the fault picture detecting unit is carried out interpolation calculation to the fault picture detecting unit, obtains the estimated image data (Xbase) of fault picture detecting unit.In step 930, whether the view data of the adjacent image detecting unit of failure judgement image detecting element exists the step point of view data, if there is the step point of view data, enters step 940, otherwise, enter step 970.In step 940, whether the step point of the view data of the adjacent image detecting unit of failure judgement image detecting element has directivity in the image detecting element array, if having directivity, enters step 970, otherwise, enter step 952 and step 953.In step 952, calculate the difference (Δ of the estimated image data of the view data of the adjacent image detecting unit in the fault picture detecting unit on the represented direction of above-mentioned directivity and fault picture detecting unit i).In step 954, all differences of calculating are asked to (Δ in step 952 i) average, obtain mean difference
Figure BSA00000632350600161
and by mean difference
Figure BSA00000632350600162
estimated image data (Xbase) addition with the fault picture detecting unit
Figure BSA00000632350600163
correction estimated image data (Xin) after the directivity that obtains fault picture is proofreaied and correct.In step 953, for the view data of the adjacent image detecting unit of fault picture detecting unit, solve extrapolation weights W l and Wr.In step 957, extrapolation weights W l and Wr based on obtaining in step 953, calculate the monolateral extrapolated value of fault picture detecting unit, obtains the correction estimated image data (Xextra) after the monolateral extrapolation of fault picture detecting unit.In step 960, judge above-mentioned mean difference
Figure BSA00000632350600171
whether be greater than 0, if be greater than 0, enter step 962, otherwise, step 963 entered.In step 962, get middle the maximum of correction estimated image data (Xextra) after correction estimated image data (Xin) after directivity is proofreaied and correct and monolateral extrapolation are proofreaied and correct correction estimated image data (X) as the fault picture detecting unit.In step 963, get the middle reckling of correction estimated image data (Xextra) after correction estimated image data (Xin) after directivity is proofreaied and correct and monolateral extrapolation the are proofreaied and correct correction estimated image data (X) as the fault picture detecting unit.In step 990, provide and proofread and correct the estimated image data.
Figure 10 is the block diagram of explanation according to the equipment of the detecting unit of the fault picture for the correcting image detection arrays 330 of one embodiment of the present of invention.This equipment comprises: view data input equipment 1010, the view data for input corresponding to image detecting element array 300 each image detecting elements; View data estimation unit 1020, for the adjacent image detecting unit based on this fault picture detecting unit 330 carries out interpolation for the fault picture detecting unit, obtain the estimated image data (Xbase) of fault picture detecting unit 330; View data step point judgment means 1030; Whether the view data that is used for the adjacent image detecting unit of failure judgement image detecting element exists the step point of view data; The first correcting unit 1050, be used in view data step point judgment means, when there is the step point of view data in the view data of the adjacent image detecting unit of failure judgement image detecting element, adopting the first correction mode to be proofreaied and correct the estimated image data of fault picture detecting unit; The second correcting unit 1070, be used in view data step point judgment means, when there is the step point of view data in the view data of the adjacent image detecting unit of failure judgement image detecting element, adopting the second correction mode to be proofreaied and correct the estimated image data of fault picture detecting unit; Image-data output device 1090, for exporting the correction estimated image data of fault picture detecting unit.
View data step point judgment means 1030 can comprise further: whether the first judgment means exists the step point of view data for the view data of the adjacent image detecting unit of failure judgement image detecting element; With the second judgment means, while for the view data of the adjacent image detecting unit at the failure judgement image detecting element, having the step point of view data, judge whether the step point of view data exists directivity in the distribution of image detecting element array.
The first correcting unit 1050 may further include: the directivity correcting unit, for the step point of the view data of the adjacent image detecting unit determining the fault picture detecting unit when there is directivity in the image detecting element array, centered by the fault picture detecting unit, view data according to the adjacent image detecting unit of the set direction fault picture detecting unit that represents described directivity, the estimated image data travel direction correction of view data based on selected adjacent image detecting unit to the fault picture detecting unit, correction estimated image data after the directivity that obtains the fault picture detecting unit is proofreaied and correct, monolateral extrapolation correcting unit, for for the fault picture detecting unit, carrying out monolateral extrapolation, the correction estimated image data after proofreading and correct with the monolateral extrapolation that obtains the fault picture detecting unit, directivity-monolateral extrapolation is proofreaied and correct selecting arrangement, for selection, is the correction estimated image data that adopt after the correction estimated image data of the directivity of fault picture detecting unit correction or the monolateral extrapolation of selection fault picture detecting unit are proofreaied and correct.
The second correcting unit may further include judgment means, for judging image detecting element, in the size of image detecting element array vertical direction, with the ratio R of the size of line direction, whether is more than or equal to predetermined threshold V threshold.When the judgment means judgement, R is more than or equal to predetermined threshold V threshold, the second correcting unit is calculated on the fault picture detecting unit, under, left, right, upper left, upper right, lower-left, the delta of the view data of each adjacent image detecting unit in the adjacent image detecting unit of bottom right, to all these calculated delta calculating mean values, estimated image data (Xbase) addition by this meansigma methods and fault picture detecting unit, obtain the correction estimated image data (X) of fault picture detecting unit, otherwise, calculate fault picture detecting unit delta of the view data of each adjacent image detecting unit in the neighbouring image detecting element on vertical direction in the image detecting element array, to all these calculated delta calculating mean values, estimated image data (Xbase) addition by this meansigma methods and fault picture detecting unit, obtain the correction estimated image data (X) of fault picture detecting unit.
The ingredient that each ingredient for the equipment 1000 of proofreading and correct fault picture detecting unit 330 of embodiments of the invention can integrally be arranged in system shown in Figure 2 for example (for example, be positioned at DAS 20 or be positioned at image reconstructor 34 etc.), also can be distributed to the different ingredients (for example, being distributed in DAS 20 and image reconstructor 34 etc.) of system shown in Figure 2.Each ingredient for the equipment 1000 of proofreading and correct fault picture detecting unit 330 of embodiments of the invention can be not only hardware, can also be the combination of firmware or hardware and firmware, or the combination of hardware, firmware and software.
Figure 11 shows according to embodiments of the invention and carries out image ratio before and after image data correction.The image on the right is to carry out image data correction image (image result that adopts existing alignment technique to be proofreaied and correct) before according to embodiments of the invention, has " artifact " in figure, striped as indicated as arrow in figure; The image on the left side is to carry out the image after image data correction according to embodiments of the invention, " artifact " wherein not resembling in right image.
Abovely by specific embodiment, describe the present invention, but the present invention is not limited to these specific embodiments.Those skilled in the art will understand, and can also carry out various modifications to the present invention, be equal to replacement, variation etc.For example, by a step in above-described embodiment or module is divided into two or more steps or module realizes, or contrary, the function of two or more steps in above-described embodiment or module or device is placed in a step or module and realizes.But these conversion, all should be within protection scope of the present invention as long as do not deviate from spirit of the present invention.In addition, some terms that present specification and claims are used not are restriction, and are only for convenience of description.

Claims (35)

1. the method for the treatment of view data, described method comprises:
The image data acquisition step, obtain the view data of the adjacent image detecting unit of image detecting element to be corrected in the image detecting element array;
The view data estimating step, the identical line adjacency graph of going together mutually with described image detecting element to be corrected in described view data based on described adjacent image detecting unit carries out interpolation as the view data of detecting unit to described image detecting element to be corrected, obtains the estimated image data of described image detecting element to be corrected; And
View data correction step, the characteristic of the view data based on described adjacent image detecting unit is revised described estimated image data.
2. the method for claim 1, wherein described view data correction step comprises:
Step point determining step, judge the step point that whether has view data in the view data of adjacent image detecting unit of described image detecting element to be corrected; And
The image data correction step, select different correction modes to be proofreaied and correct described estimated image data based on described judged result.
3. method as claimed in claim 2, wherein, described image data correction step also comprises:
While having the step point of view data in the view data of described adjacent image detecting unit, according to the first correction mode, described estimated image data are proofreaied and correct, otherwise, according to the second correction mode, described estimated image data are proofreaied and correct.
4. method as claimed in claim 3, wherein, described the first correction mode comprises:
Along the adjacent image detecting unit of the set direction described to be corrected image detecting element different from the row of described image detecting element array, described direction has reflected the distribution arrangement of step point in the view data of described adjacent image detecting unit;
Calculate the deviation of directivity Δ of the view data of selected adjacent image detecting element, and
Proofread and correct described estimated image data based on described deviation of directivity Δ, obtain the view data after directivity is proofreaied and correct.
5. method as claimed in claim 4, wherein, the deviation of directivity Δ of view data that calculates the adjacent image detecting unit of described image detecting element to be corrected comprises:
Calculate the difference of the estimated image data of the view data of selected adjacent image detecting element and described image detecting element to be corrected and
Calculate described difference
Figure FSA00000632350500022
meansigma methods
Figure FSA00000632350500023
as described deviation of directivity Δ.
6. method as described as claim 3-5, wherein, described the first correction mode comprises:
Described image detecting element to be corrected is carried out to monolateral extrapolation, obtain the view data after the monolateral extrapolation of described image detecting element to be corrected.
7. method as claimed in claim 6, wherein, described image detecting element to be corrected is carried out to monolateral extrapolation and comprise:
Calculate the left side weight w of described image detecting element to be corrected in described image detecting element array leftwith right side weight w right, and
According to described left side weight w leftwith right side weight w rightcalculate the view data X after monolateral extrapolation:
x=w left*x left+w right*x right
Wherein, x leftand x rightit is respectively the view data of described image detecting element left side adjacent image detecting unit to be corrected and right side adjacent image detecting unit.
8. method as claimed in claim 7, wherein, calculate described left side weight w in such a way leftwith described right side weight w right:
For described image detecting element to be corrected, the extrapolated value x of lastrow left side extrapolation estimation lastrow up_left, the extrapolated value x of lastrow right side extrapolation estimation lastrow up_left, the extrapolated value x of next line left side extrapolation estimation lastrow up_left, the extrapolated value x of next line right side extrapolation estimation lastrow up_left, wherein, described extrapolated value is to adopt the view data of adjacent image detecting unit of the respective side of corresponding line to carry out Lagrange's interpolation, and
By described extrapolated value x up_left, x up_left, x up_left, x up_leftwith the upper adjacent image detecting unit of described image detecting element to be corrected and the view data x of lower adjacent image detecting unit up, x downsolve following equation and obtain right outside slotting weight w up_left, w down_leftwith right outside slotting weight w up_right, w down right:
x up=w up_left*x up_left+w up_right*x up_right
x down=w down_left*x down_left+w down_right*x down_right
W up_left+ w up_right=1, and
w down_left+w down_right=1
Pass through obtained left outside slotting weight w up_left, w down_left, and right outside slotting weight w up_right, w down rightcalculate described left side weight w leftwith described right side weight w right:
w left=a*w up_left+(1-a)*w down_left
w right=bw up_right+(1-b)*w down?right
Wherein, a and b are the constant of value between 0 to 1.
9. method as described as claim 4 or 5, wherein, the view data of view data after the correction of described image detecting element to be corrected after described directivity is proofreaied and correct.
10. method as described as any one in claim 6-8, wherein,
View data after adopting view data after described monolateral extrapolation as the correction of described image detecting element to be corrected; Perhaps
Calculate the deviation of directivity Δ of view data of the selected adjacent image detecting element of described image detecting element to be corrected, be greater than at 0 o'clock at described deviation of directivity Δ, view data after selecting view data after described directivity is proofreaied and correct and the maximum in the view data after described monolateral extrapolation as the correction of described image detecting element to be corrected, otherwise, the view data after selecting view data after described directivity is proofreaied and correct and the reckling in the view data after described monolateral extrapolation as the correction of described image detecting element to be corrected.
11. method as described as any one in claim 3-10, wherein, described the second correction mode comprises:
Calculate the meansigma methods of the deviation of the view data of adjacent image detecting unit of described image detecting element to be corrected and described estimated image data, meansigma methods based on described deviation is proofreaied and correct described estimated image data, obtains the view data after the correction of described image detecting element to be corrected.
12. method as claimed in claim 11, wherein, the meansigma methods of the view data of the adjacent image detecting unit of the described image detecting element to be corrected of described calculating and the deviation of described estimated image data comprises:
Calculate the meansigma methods of the arithmetic difference of the view data of the upper and lower, left and right of described image detecting element to be corrected, upper left, upper right, bottom left, bottom right adjacent image detecting unit and described estimated image data, as the meansigma methods of described deviation.
13. method as claimed in claim 11, wherein, the meansigma methods of the view data of the adjacent image detecting unit of the described image detecting element to be corrected of described calculating and the deviation of described estimated image data comprises:
Calculate the meansigma methods of the arithmetic difference of the view data of upper and lower adjacent image detecting unit of described image detecting element to be corrected and described estimated image data, as the meansigma methods of described deviation.
14. method as claimed in claim 11, wherein, the meansigma methods of the view data of the adjacent image detecting unit of the described image detecting element to be corrected of described calculating and the deviation of described estimated image data comprises:
Whether the size that judges image detecting element column direction in described image detecting element array is more than or equal to predetermined threshold V with the ratio R of the size of line direction thresholdif R is more than or equal to described predetermined threshold V thresholdcalculate the meansigma methods of the arithmetic difference of the view data of the upper and lower, left and right of described image detecting element to be corrected, upper left, upper right, bottom left, bottom right adjacent image detecting unit and described estimated image data, meansigma methods as described deviation, otherwise, calculate the meansigma methods of the arithmetic difference of the view data of upper and lower adjacent image detecting unit of described image detecting element to be corrected and described estimated image data, as the meansigma methods of described deviation.
15. method as claimed in claim 14, wherein said described predetermined threshold V thresholdthe numerical value between 0.5 to 2.
16. method as described as any one in claim 1-15, wherein, described view data estimating step also comprises:
The view data of the adjacent image detecting unit of being expert at based on described image detecting element to be corrected is carried out Lagrange's interpolation to described image detecting element to be corrected, obtains described estimated image data.
17. the equipment for the treatment of view data, described equipment comprises:
The image data acquisition device, obtain the view data of the adjacent image detecting unit of image detecting element to be corrected in the image detecting element array;
The view data estimation unit, the identical line adjacency graph of going together mutually with described image detecting element to be corrected in described view data based on described adjacent image detecting unit carries out interpolation as the view data of detecting unit to described image detecting element to be corrected, obtains the estimated image data of described image detecting element to be corrected; And
Image data correction apparatus, the characteristic of the view data based on described adjacent image detecting unit is revised described estimated image data.
18. equipment as claimed in claim 17, wherein, described image data correction apparatus comprises:
Step point judgment means, judge the step point that whether has view data in the view data of adjacent image detecting unit of described image detecting element to be corrected; And
The image data correction device, select different correction modes to be proofreaied and correct described estimated image data based on described judged result.
19. equipment as claimed in claim 18, wherein, described image data correction device also comprises:
The first mode correcting unit, while for the view data at described adjacent image detecting unit, having the step point of view data, proofreaied and correct described estimated image data according to the first correction mode, and
The second pattern correcting unit, while for the view data at described adjacent image detecting unit, not having the step point of view data, proofreaied and correct described estimated image data according to the second correction mode.
20. equipment as claimed in claim 19, wherein, described first mode correcting unit also comprises:
Direction selection device, for the adjacent image detecting unit along the set direction described to be corrected image detecting element different from the row of described image detecting element array, described direction has reflected the distribution arrangement of step point in the view data of described adjacent image detecting unit;
Deviation of directivity accountant, for the deviation of directivity Δ of the view data of calculating selected adjacent image detecting element, and
The directivity correcting unit, for based on described deviation of directivity Δ, proofreading and correct described estimated image data, obtain the view data after directivity is proofreaied and correct.
21. equipment as claimed in claim 20, wherein, described deviation of directivity accountant also comprises:
Difference for the estimated image data of calculating the view data of selected adjacent image detecting element and described image detecting element to be corrected
Figure FSA00000632350500061
device, and
For calculating described difference
Figure FSA00000632350500062
meansigma methods
Figure FSA00000632350500063
and described this meansigma methods device as described deviation of directivity Δ.
22. equipment as described as claim 19-21, wherein, described first mode correcting unit also comprises:
For described image detecting element to be corrected being carried out to the device of the view data of monolateral extrapolation after with the monolateral extrapolation that obtains described image detecting element to be corrected.
23. equipment as claimed in claim 22 wherein, describedly also comprises for the device that described image detecting element to be corrected is carried out to the view data of monolateral extrapolation after with the monolateral extrapolation that obtains described image detecting element to be corrected:
For calculating described image detecting element to be corrected at the left side of described image detecting element array weight w leftwith right side weight w rightdevice, and
Be used for according to described left side weight w leftwith right side weight w rightcalculate the device of the view data X after monolateral extrapolation, wherein,
x=w left*x left+w right*x right
Wherein, x leftand x rightit is respectively the view data of described image detecting element left side adjacent image detecting unit to be corrected and right side adjacent image detecting unit.
24. equipment as claimed in claim 23 is wherein, described for calculating described image detecting element to be corrected at the left side of described image detecting element array weight w leftwith right side weight w rightdevice:
For described image detecting element to be corrected, the extrapolated value x of lastrow left side extrapolation estimation lastrow up_left, the extrapolated value x of lastrow right side extrapolation estimation lastrow up_left, the extrapolated value x of next line left side extrapolation estimation lastrow up_left, the extrapolated value x of next line right side extrapolation estimation lastrow up_left, wherein, described extrapolated value is to adopt the view data of adjacent image detecting unit of the respective side of corresponding line to carry out Lagrange's interpolation, and
By described extrapolated value x up_left, x up_left, x up_left, x up_leftwith the upper adjacent image detecting unit of described image detecting element to be corrected and the view data x of lower adjacent image detecting unit up, x downsolve following equation and obtain right outside slotting weight w up_left, w down_leftwith right outside slotting weight w up_right, w down right:
x up=w up_left*x up_left+w up_right*x up_right
x down=w down_left*x down_left+w down_right*x down_right
W up_left+ w up_right=1, and
w down_left+w down_right=1
Pass through obtained left outside slotting weight w up_left, w down_left, and right outside slotting weight w up_right, w down_rightcalculate described left side weight w leftwith described right side weight w right:
w left=a*w up_left+(1-a)*w down_left
w right=bw up_right+(1-b)*w down?right
Wherein, a and b are the constant of value between 0 to 1.
25. equipment as described as claim 20 or 21, wherein, the view data view data as the correction of described to be corrected image detecting element after of described image data correction device after described directivity is proofreaied and correct.
26. equipment as described as any one in claim 22-24, wherein, described image data correction device:
View data after adopting view data after described monolateral extrapolation as the correction of described image detecting element to be corrected; Perhaps
Calculate the deviation of directivity Δ of view data of the selected adjacent image detecting element of described image detecting element to be corrected, be greater than at 0 o'clock at described deviation of directivity Δ, view data after selecting view data after described directivity is proofreaied and correct and the maximum in the view data after described monolateral extrapolation as the correction of described image detecting element to be corrected, otherwise, the view data after selecting view data after described directivity is proofreaied and correct and the reckling in the view data after described monolateral extrapolation as the correction of described image detecting element to be corrected.
27. equipment as described as any one in claim 19-26, wherein, described the second pattern correcting unit also comprises:
For the device of the meansigma methods of the deviation of the view data of adjacent image detecting unit of calculating described image detecting element to be corrected and described estimated image data, and
Proofread and correct the device of the view data of described estimated image data after with the correction that obtains described image detecting element to be corrected for meansigma methods based on described deviation.
28. equipment as claimed in claim 27, wherein, the device of the meansigma methods of the deviation of the described view data for the adjacent image detecting unit that calculates described image detecting element to be corrected and described estimated image data comprises:
For the meansigma methods of the arithmetic difference of the view data of the upper and lower, left and right of calculating described image detecting element to be corrected, upper left, upper right, bottom left, bottom right adjacent image detecting unit and described estimated image data the device of the meansigma methods using the meansigma methods of described arithmetic difference as described deviation.
29. equipment as claimed in claim 27, wherein, the device of the meansigma methods of the deviation of the described view data for the adjacent image detecting unit that calculates described image detecting element to be corrected and described estimated image data comprises:
For the meansigma methods of the arithmetic difference of the view data of upper and lower adjacent image detecting unit of calculating described image detecting element to be corrected and described estimated image data the device of the meansigma methods using the meansigma methods of described arithmetic difference as described deviation.
30. equipment as claimed in claim 27, wherein, the device of the meansigma methods of the deviation of the described view data for the adjacent image detecting unit that calculates described image detecting element to be corrected and described estimated image data also comprises:
Whether the threshold decision device, be more than or equal to predetermined threshold V in the size of described image detecting element array column direction with the ratio R of the size of line direction for judging image detecting element threshold,
For at threshold decision device judgement R, being more than or equal to described predetermined threshold V thresholdsituation under calculate the meansigma methods of arithmetic difference of the view data of the upper and lower, left and right of described image detecting element to be corrected, upper left, upper right, bottom left, bottom right adjacent image detecting unit and described estimated image data the device of the meansigma methods using the meansigma methods of described arithmetic difference as described deviation, and
For at threshold decision device judgement R, being less than described predetermined threshold V thresholdsituation under calculate the meansigma methods of arithmetic difference of the view data of upper and lower adjacent image detecting unit of described image detecting element to be corrected and described estimated image data the device of the meansigma methods using the meansigma methods of described arithmetic difference as described deviation.
31. equipment as claimed in claim 30, wherein, described predetermined threshold V thresholdthe numerical value between 0.5 to 2.
32. equipment as described as any one in claim 17-31, wherein, described view data estimation unit also comprises:
View data for the adjacent image detecting unit of being expert at based on described image detecting element to be corrected is carried out Lagrange's interpolation to obtain the device of described estimated image data to described image detecting element to be corrected.
33. the product for the treatment of view data, comprise for executing claims the unit of the step of 1-16 in any one.
34. an imaging device, is characterized in that comprising the equipment for the treatment of view data as described as any one in claim 17-32 or the described product for the treatment of view data of claim 33.
35. imaging device as claimed in claim 34, wherein, described imaging device is selected from one of CT imaging device, plane x radiation imaging apparatus, supersonic imaging apparatus, magnetic resonance (MR) imaging device, electron beam x line laminagraphy (EBT) imaging device, positron emission tomography photograph (PET) imaging device, single photon emission computed tomography photograph (SPECT) imaging device, pico computer tomography (micro CT) imaging device, grand computed tomography (macro CT) imaging device.
CN201110405421.9A 2011-11-30 Image data correction equipment and method Active CN103126704B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110405421.9A CN103126704B (en) 2011-11-30 Image data correction equipment and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110405421.9A CN103126704B (en) 2011-11-30 Image data correction equipment and method

Publications (2)

Publication Number Publication Date
CN103126704A true CN103126704A (en) 2013-06-05
CN103126704B CN103126704B (en) 2016-12-14

Family

ID=

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4787098A (en) * 1987-04-10 1988-11-22 Kabushiki Kaisha Toshiba Method for obtaining calibrated tomographic image data to correct for collimator width differences
JPH09275490A (en) * 1996-04-06 1997-10-21 Ricoh Co Ltd Dot area separator
JPH10243300A (en) * 1997-02-27 1998-09-11 Matsushita Electric Ind Co Ltd Pixel defect detector and pixel defect corrector
US20010038706A1 (en) * 2000-04-20 2001-11-08 Kai Eck X-ray examination apparatus and method for forming an X-ray image
CN1551744A (en) * 2001-09-03 2004-12-01 西门子公司 Method for recovering fault channel in multi-channel optical beam detector
US20050117708A1 (en) * 2003-08-08 2005-06-02 Young-Bin Cho Method and system for calibrating a source and detector instrument
US7142636B2 (en) * 2003-09-23 2006-11-28 General Electric Company System and method for defective detector cell and DAS channel correction
CN101305620A (en) * 2005-11-09 2008-11-12 辉达公司 Graphics processing unit use and device
WO2010032497A1 (en) * 2008-09-19 2010-03-25 コニカミノルタエムジー株式会社 Defective pixel determining method, defective pixel determining program, radiological image detector, and defective pixel determining system
US20110075909A1 (en) * 2009-09-30 2011-03-31 Fujifilm Corporation X-ray imaging system, imaging method and computer readable media including imaging program

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4787098A (en) * 1987-04-10 1988-11-22 Kabushiki Kaisha Toshiba Method for obtaining calibrated tomographic image data to correct for collimator width differences
JPH09275490A (en) * 1996-04-06 1997-10-21 Ricoh Co Ltd Dot area separator
JPH10243300A (en) * 1997-02-27 1998-09-11 Matsushita Electric Ind Co Ltd Pixel defect detector and pixel defect corrector
US20010038706A1 (en) * 2000-04-20 2001-11-08 Kai Eck X-ray examination apparatus and method for forming an X-ray image
CN1551744A (en) * 2001-09-03 2004-12-01 西门子公司 Method for recovering fault channel in multi-channel optical beam detector
US20050117708A1 (en) * 2003-08-08 2005-06-02 Young-Bin Cho Method and system for calibrating a source and detector instrument
US7142636B2 (en) * 2003-09-23 2006-11-28 General Electric Company System and method for defective detector cell and DAS channel correction
CN101305620A (en) * 2005-11-09 2008-11-12 辉达公司 Graphics processing unit use and device
WO2010032497A1 (en) * 2008-09-19 2010-03-25 コニカミノルタエムジー株式会社 Defective pixel determining method, defective pixel determining program, radiological image detector, and defective pixel determining system
US20110075909A1 (en) * 2009-09-30 2011-03-31 Fujifilm Corporation X-ray imaging system, imaging method and computer readable media including imaging program

Similar Documents

Publication Publication Date Title
JP6253426B2 (en) Method and apparatus for geometric calibration of CT scanner
CN101454801B (en) Local motion compensation based on list mode data
CN106232007B (en) X ray CT device and processing unit
CN107330949A (en) A kind of artifact correction method and system
US7142636B2 (en) System and method for defective detector cell and DAS channel correction
US20180061045A1 (en) Systems and methods for adaptive imaging systems
CN105377141B (en) X ray CT device
US20120243759A1 (en) Image processing apparatus, x-ray ct apparatus, and image processing method
US9538974B2 (en) Methods and systems for correcting table deflection
CN101911119B (en) Image registration alignment metric
CA2768296A1 (en) Extended low contrast detectability for radiographic imaging systems
CN110246580B (en) Cranial image analysis method and system based on neural network and random forest
JP4035320B2 (en) Method and apparatus for pixel replacement in panel detector
CN103026379A (en) Method for estimating image noise level
US20140198892A1 (en) X-ray ct device, and method
JP2012115659A (en) Diagnostic imaging apparatus and method
CN106539590A (en) Computer tomography bearing calibration and computed tomography imaging system
Pretorius et al. Diminishing the impact of the partial volume effect in cardiac SPECT perfusion imaging
Johnston et al. Temporal and spectral imaging with micro‐CT
JP2012088300A (en) Positron radiation tomography (pet) imaging system, method and nuclear medicine imaging system
JP2000083945A (en) Device and method for forming tomographic image
US7187750B1 (en) Method and apparatus for compensating non-uniform detector collimator plates
JP5588697B2 (en) X-ray CT system
US9858688B2 (en) Methods and systems for computed tomography motion compensation
CN104000618B (en) The true number of photons gate control method of one ring carries out the bearing calibration of respiratory movement gate

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant