CN102842114A - Filling length determination device and method and image processing method and device - Google Patents
Filling length determination device and method and image processing method and device Download PDFInfo
- Publication number
- CN102842114A CN102842114A CN2011101901209A CN201110190120A CN102842114A CN 102842114 A CN102842114 A CN 102842114A CN 2011101901209 A CN2011101901209 A CN 2011101901209A CN 201110190120 A CN201110190120 A CN 201110190120A CN 102842114 A CN102842114 A CN 102842114A
- Authority
- CN
- China
- Prior art keywords
- image
- aliasing effect
- filtering
- filling length
- length
- 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
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000003672 processing method Methods 0.000 title claims abstract description 16
- 230000000694 effects Effects 0.000 claims abstract description 239
- 238000001914 filtration Methods 0.000 claims abstract description 167
- 238000012545 processing Methods 0.000 claims abstract description 90
- 238000012360 testing method Methods 0.000 claims description 53
- 230000007480 spreading Effects 0.000 claims description 26
- 238000003745 diagnosis Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 105
- 230000008569 process Effects 0.000 description 23
- 238000010586 diagram Methods 0.000 description 10
- 230000009466 transformation Effects 0.000 description 9
- 238000004590 computer program Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000002224 dissection Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000005055 memory storage Effects 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration by non-spatial domain filtering
-
- G06T5/70—
Abstract
The invention discloses a filling length determination device and method and an image processing method and device. The device for determining a filling length of an image in an image filtering direction comprises a aliasing effect function acquisition unit and a filling length determination unit, wherein the aliasing effect function acquisition unit is configured to acquire an aliasing effect function of a frequency domain filter according to the image, and the aliasing effect function is a function representing the change of an aliasing effect in a spatial domain caused by frequency domain filter processing; and the filling length determination unit is configured to determine the filling length of the image in the image filtering direction based on the aliasing effect function. Thus, the aliasing effect can be substantially eliminated in the filter processing, and the calculation amount is reduced.
Description
Technical field
Present invention relates in general to computer realm, relate in particular to a kind of apparatus and method that are used for when carrying out image filtering, confirming the filling length of image on the image filtering direction, and the image processing method and the device that utilize these apparatus and method.
Background technology
In existing image processing techniques,, need carry out Filtering Processing to image for the edge effect in the removal of images.
In a kind of common Filtering Processing, utilize discrete Fourier transformation (DFT) that image is transformed from a spatial domain to frequency domain earlier.Then, in frequency domain, utilize filter function to come the image after the conversion is carried out filtering.Specifically, Filtering Processing can comprise filter function is acted on the image, for example, in frequency domain, with the mode (be equivalent to and in spatial domain, carry out circular convolution) that filter function and image function are multiplied each other filter function acted on the image.Afterwards, spatial domain is returned in the image transformation after utilizing inversefouriertransform with Filtering Processing.Thus, eliminated the edge effect in the image.
Wherein, When utilizing wave filter that image is carried out Filtering Processing; For fear of in Filtering Processing because the wave filter aliasing effect that (being equivalent to the circular convolution in the spatial domain) produces when acting on image (that this aliasing effect comes down to is adjacent with image header owing to image afterbody when the circular convolution, mutual interference mutually between image afterbody and the image header cause); Can expand and reasonable value (for example, be set to 0 image value) is filled in the image expanding position eliminate or avoid aliasing effect image.
For example; In the length of image on the filtering direction is A and filter length is under the situation of B; In order to avoid aliasing effect fully; Then at least need be on the filtering direction begin image length is extended to the length of A+B-1, and expanding location is filled rational image value from the position of image end.For example, be that the image spreading part of B-1 all is set under the situation of 0 image value to filling length, when in frequency domain, carrying out Filtering Processing, these are that 0 picture position also is 0 to the interference that image header forms for the image afterbody, have avoided aliasing effect thus.In addition, in a concrete example, be under the situation of A in length on the filtering direction and filter length, then need on the filtering direction, begin the length of image spreading the position from image end to (2*A-1) at image.In this case, the filling length of image is A-1.
Because in common Filtering Processing, filter length is consistent with image length often.Therefore, in the Filtering Processing of above description, owing to need to carry out Fourier transform and Fourier inversion to the image (for example, the length of above-mentioned 2*A-1) that length almost is expanded to twice, the calculating quantitative change is big.
Summary of the invention
Provided hereinafter about brief overview of the present invention, so that the basic comprehension about some aspect of the present invention is provided.Should be appreciated that this general introduction is not about exhaustive general introduction of the present invention.It is not that intention is confirmed key of the present invention or pith, neither be intended to limit scope of the present invention.Its purpose only is to provide some notion with the form of simplifying, with this as the preorder in greater detail of argumentation after a while.
In technical scheme of the present invention; Provide to be used for when image filtering confirming to eliminate basically the aliasing effect that Filtering Processing produces and can not cause method and apparatus simultaneously, and utilized this method and apparatus to carry out the method and apparatus that image filtering is handled than the filling length of intensive.
Handle owing to adopt long filling length (for example, the above-mentioned filling length that reaches A-1 or B-1) to carry out image filtering, can under the little situation of calculated amount, eliminate the influence of the aliasing item that produces in the Filtering Processing basically picture quality.
For this reason, on the one hand, embodiments of the invention provide the device of the filling length of a kind of definite image on the image filtering direction, comprising: aliasing effect function acquiring unit is configured to obtain the aliasing effect function of frequency domain filter to said image; Fill length and confirm the unit, be configured to confirm the filling length of said image based on said aliasing effect function.
On the other hand, embodiments of the invention also provide a kind of image processing apparatus, comprising: the device of the filling length of above-mentioned definite image on the image filtering direction is configured to confirm the filling length of image; The image spreading unit is configured in spatial domain according to said filling length expanded images on the filtering direction; Filter unit is configured to the image after the expansion is carried out filtering.
Correspondingly, embodiments of the invention also provide the method for the filling length of a kind of definite image on the image filtering direction, comprising: obtain the aliasing effect function of frequency domain filter to this image; Confirm the filling length of this image based on said aliasing effect function.
In addition, embodiments of the invention also provide a kind of image processing method, comprising: the filling length of coming to confirm image according to the method described above; In spatial domain according to said filling length expanded images on the filtering direction; And to the expansion after image carry out filtering.
Another object of the present invention provides when being carried out by computing equipment and makes computing equipment can carry out the computer program code of the processing procedure of said method, with and on store the computer-readable recording medium and the computer program of this computer program code.
Through below in conjunction with the detailed description of accompanying drawing to most preferred embodiment of the present invention, above-mentioned and other advantage of the present invention will be more obvious.
Description of drawings
The present invention can wherein use same or analogous Reference numeral to represent identical or similar parts in institute's drawings attached through with reference to hereinafter combining the given description of accompanying drawing to be better understood.Said accompanying drawing comprises in this manual and forms the part of this instructions together with following detailed description, and is used for further illustrating the preferred embodiments of the present invention and explains principle and advantage of the present invention.In the accompanying drawings:
Fig. 1 shows the method that is used for when carrying out image filtering, confirming the filling length of image on the image filtering direction according to an embodiment of the invention;
Fig. 2 is the process flow diagram that the processing of the aliasing effect function that obtains frequency domain filter according to an embodiment of the invention is shown;
Fig. 3 shows the aliasing effect function that obtains to Butterworth (Butterworth) wave filter according to an embodiment of the invention;
Fig. 4 shows the process flow diagram of confirming the processing of filling length according to an embodiment of the invention according to the aliasing effect function;
Fig. 5 shows the process flow diagram of image processing method according to an embodiment of the invention;
Fig. 6 shows the particular flow sheet of Filtering Processing according to an embodiment of the invention;
Fig. 7 is the process flow diagram that the image processing method that comprises the processing that need to judge whether expanded images according to an embodiment of the invention is shown;
Fig. 8 shows the process flow diagram that judges whether the processing of needs expanded images according to an embodiment of the invention according to filling length;
Fig. 9 illustrates the device that is used for confirming the filling length of image on the image filtering direction according to an embodiment of the invention;
Figure 10 illustrates the structural representation of aliasing effect function acquiring unit according to an embodiment of the invention;
Figure 11 shows and fills the structural representation that length is confirmed unit 1100 according to an embodiment of the invention;
Figure 12 shows the structural representation according to the image processing apparatus 1200 of this embodiment;
Figure 13 shows the structural representation according to the image processing apparatus 1300 of this embodiment;
Figure 14 shows the structural drawing of judging part according to an embodiment of the invention;
Figure 15 shows the block scheme that wherein can realize according to the exemplary configurations of the general purpose personal computer of the method for the embodiment of the invention and/or device.
Embodiment
To combine accompanying drawing that example embodiment of the present invention is described hereinafter.In order to know and for simplicity, in instructions, not describe all characteristics of actual embodiment.Yet; Should understand; In the process of any this practical embodiments of exploitation, must make a lot of decisions, so that realize developer's objectives, for example specific to embodiment; Meet and system and professional those relevant restrictive conditions, and these restrictive conditions may change along with the difference of embodiment to some extent.In addition, might be very complicated and time-consuming though will also be appreciated that development, concerning the those skilled in the art that have benefited from present disclosure, this development only is customary task.
At this; What also need explain a bit is; For fear of having blured the present invention, only show in the accompanying drawings and closely-related apparatus structure of scheme according to the present invention and/or treatment step, and omitted other details little with relation of the present invention because of unnecessary details.
Fig. 1 shows the method that is used for when carrying out image filtering, confirming the filling length of image on the image filtering direction according to an embodiment of the invention.
As shown in Figure 1, can obtain the aliasing effect function at step S102 place.
Specifically, at step S102 place, can be to treating that the image of filtering obtains the aliasing effect function of frequency domain filter on the image filtering direction.
Wherein, the aliasing effect function can be to characterize because the function that the aliasing effect that wave filter acts on the image to be caused changes on spatial domain.
For example; When the filtering direction of image be one dimension (for example; On directions X in 3-D view, Y direction or the Z direction during filtering) time, can be the one dimension function on this image filtering direction (for example, the directions X in the 3-D view, Y direction or Z direction) to the aliasing effect function of this image.Perhaps; When the image filtering direction be the two dimension (for example; On directions X in 3-D view and the Y direction during filtering) or more during multidimensional; To the aliasing effect function of this image can be the two-dimensional function or the function of multidimensional more on this image filtering direction (for example, directions X in 3-D view and Y direction).
As a preferred exemplary, under the application scenarios of medical image, the image filtering direction for example can be the Z-direction in the 3-D view.In this case, the aliasing effect function can be the function on the Z of this medical image axle.Wherein, said medical image can be the image that the data according to the detected person who utilizes the medical diagnosis imaging device to obtain form.Described medical diagnosis device includes but not limited to: x-ray imaging diagnostic device, ultrasonic diagnostic imaging device, computed tomography (CT) device, magnetic resonance imaging (MRI) diagnostic imaging device, positron emission computerized tomography (Positron Emission Tomography, PET) device etc.
In addition, the method for acquisition aliasing effect function can have a variety of.
For example, can obtain the aliasing effect function through the variation before and after the analysis image filtering and through analyzing this variation.
Perhaps, also can directly obtain the aliasing effect function of wave filter from outside (for example, from the memory storage of the aliasing effect function that stores each wave filter).
Then, can confirm the filling length of image at step S104 place.
Specifically, can confirm the filling length of image based on the aliasing effect function.
Because the aliasing effect that the aliasing effect function causes in the time of can characterizing the frequency domain filtering processing is in the variation of spatial domain, the filling length in the time of therefore can confirming can aliasing effect be restricted to the degree that is enough to ignore according to the aliasing effect function on spatial domain.
Through confirming to fill length according to the aliasing effect function, and needn't be on spatial domain with image spreading to long length.Therefore, can reduce on the one hand in the image filtering processing because the long bigger operand that image length brought; On the other hand, owing to considered the influence of aliasing effect,, guarantee picture quality so also can eliminate the influence of aliasing effect basically.
In one embodiment of the invention, can obtain the aliasing effect function in the variation before and after the filtering through analysis image.
Specifically; Can be through wave filter (be for example acted on image; On the filtering direction, have the image of identical length with treating filtering image), obtain to characterize the big or small aliasing effect function of aliasing effect that wave filter causes according to this image in the variation before and after the Filtering Processing.
Fig. 2 is the process flow diagram that the processing of the aliasing effect function that obtains frequency domain filter according to an embodiment of the invention is shown.
As shown in Figure 2, at step S202 place, make up a test pattern.
Specifically, can make up an image that on spatial domain, has and treat the filtering image equal length as test pattern.
Said test pattern can be the image of arbitrary form, as long as it has on spatial domain and treats the filtering image equal length.For example, in some cases, can be directly with treating that filtering image is as test pattern.
In a preferred embodiment of the invention, said test pattern can be that the image that on the filtering direction, is rendered as the form of impulse function calculates to simplify.
At one more in the preferred embodiment; Said test pattern can be following this image: wherein; On the filtering direction (for example; Z axle under the one-dimensional filtering situation), the image value of first position of said image can be for unit strength (for example, 1) other position image value be 0.
Then, as shown in Figure 2, at step S204 place, carry out Filtering Processing to the test pattern that constructs.
Specifically, can through Fourier transform this test pattern be transformed to frequency domain earlier.
Then, in frequency domain, test pattern is carried out Filtering Processing.
For example, can multiply each other at frequency domain through function filter function and test pattern (be equivalent in spatial domain and carry out circular convolution) realize Filtering Processing.
Then, can be through Fourier inversion again with the test pattern inverse transformation spatial domain after the Filtering Processing.
Realized thus the Filtering Processing of test pattern in frequency domain.
Then, at step S206 place, can obtain the aliasing effect function based on the variation of said test pattern before and after the image filtering.
Obviously, if in frequency domain, do not carry out Filtering Processing, the test pattern after the inverse transformation should not change (can still be the form of impulse function for example) then.
But if test pattern carries out Filtering Processing in frequency domain, on the one hand, Filtering Processing can be eliminated the edge effect of test pattern.And on the other hand, the aliasing effect that produces in the Filtering Processing can be so that the image value at the diverse location place of the test pattern after the inverse transformation in spatial domain changes when inverse transformation is returned spatial domain.
Thereby, can confirm to be used to characterize frequency domain filtering according to the conversion of the image value at the diverse location place of the test pattern after the Filtering Processing in spatial domain and handle the aliasing effect function that the aliasing effect that causes changes on spatial domain.
As stated; At test pattern be under the situation of form of impulse function (wherein; The intensity of first pixel of this test pattern be unit strength 1 and the value of other pixel is 0), the variation of the diverse location that the test pattern after the inverse transformation can reflect image comparatively intuitively before and after image filtering.Thereby, can obtain the aliasing effect function according to this variation simple, intuitive ground.
More than utilize the test pattern of impulse function form to confirm that the description of aliasing effect function is merely example, the invention is not restricted to this.In fact, said test pattern can be arbitrary form.
For example, under the situation of the test pattern that adopts other form, can confirm the aliasing effect function in the variation before and after the Filtering Processing according to each images of positions value in the test pattern.
Specifically, can calculate the variation of each images of positions value before and after Filtering Processing in the test pattern earlier; Then, confirm the aliasing effect function through mathematics manipulation (for example, dissection process etc.).
Fig. 3 shows the aliasing effect function that obtains to Butterworth (Butterworth) wave filter according to an embodiment of the invention.
It is to be noted; Aliasing effect function shown in Figure 3 is an example only, the invention is not restricted to this, in fact; Aliasing effect function to Butterworth filter also can adopt other suitable form, as long as it can reflect the distribution of the size of aliasing effect in spatial domain.
Aliasing effect function shown in Figure 3 can illustrate the distribution of size in spatial domain of aliasing effect; Wherein transverse axis representative on spatial domain apart from by the image section at the different distance place of the image section of filtering; The longitudinal axis is represented the aliasing effect corresponding with each image section size, and the area that curve and transverse axis are surrounded is relevant with the size of aliasing effect.
Therefore; Can determine the size of the aliasing effect corresponding based on Fig. 3 with different distance; Come correspondingly to eliminate because frequency domain filtering is handled the aliasing effect that causes through image spreading is filled suitable value (for example, be 0 value) to respective and in corresponding position.
In the prior art, through filtering image being expanded to the length identical and expanding location being filled to 0 image value, come to eliminate fully the aliasing item that possibly occur and eliminate aliasing effect with filtering image.
Although this mode can be eliminated aliasing effect,, brought bigger calculated amount to the Fourier transform and the Fourier inversion of the image after the expansion.
Thereby, in an embodiment according to the present invention, can determine the filling length that can eliminate aliasing effect basically according to the aliasing effect function that obtains.Through utilizing this short filling length to come expanded images and carry out Filtering Processing, can eliminate the aliasing effect that Filtering Processing is brought basically on the one hand, can reduce calculated amount effectively on the other hand.
About confirm to fill the processing of length according to the aliasing effect function, can realize with various suitable manner.
Fig. 4 shows the process flow diagram of confirming the processing of filling length according to an embodiment of the invention according to the aliasing effect function.
As shown in Figure 4, at step S402 place, can calculate the size of total aliasing effect according to the aliasing effect function.
For example, be under the situation of form shown in Figure 3 at the aliasing effect function, can be through the aliasing effect function be obtained the size of total aliasing effect at entire image length range (for example, length range) corresponding with the length of filtering image integration.
Then, at step S404 place, can calculate the size of the aliasing effect of diverse location place in image (promptly distance is by the image section different distance place of filtering) generation according to the aliasing effect function.
For example, be under the situation of form shown in Figure 3 at the aliasing effect function, can calculate the aliasing effect that diverse location place in image (promptly distance is by the image section different distance place of filtering) produces.
Specifically, can obtain size through the aliasing effect function is carried out integration about diverse location in the image in the aliasing effect of corresponding position.
Then, at step S406 place, the ratio that can the basis aliasing effect corresponding with each picture position accounts for total aliasing effect is confirmed said filling length.
Specifically, can calculate the ratio of aliasing effect corresponding and total aliasing effect earlier, and the distance when ratio surpassed predetermined threshold is as said filling length with each picture position.Wherein, predetermined threshold can be according to circumstances and predetermined numerical value, for example can be 85%, 90%, 95%, 96%, 97%, 98%, 99% etc.
For example, in the embodiment shown in fig. 3, can find out that the aliasing effect that 6 positions produce accounts for 95% of total aliasing effect.Therefore; Be set in said threshold value under 95% the situation; Can the filling length of image be confirmed as 6, afterwards 6 of image spreading and each extension bits are set to 0 value, the aliasing effect that caused these 6 positions when the image transformation after expanding was to frequency domain like this is 0.Therefore, 95% in the aliasing effect that produces in the time of can eliminating Filtering Processing eliminated the aliasing effect that when frequency domain filtering is handled, causes thus basically.
Through above-mentioned mode, can determine the filling length that to eliminate aliasing effect (for example, 95% in the foregoing description) basically, and need not can reduce the calculated amount in the follow-up Filtering Processing thus with image spreading to whole complete image length.
Through utilizing above-mentioned definite image to fill the method for length, embodiments of the invention further provide a kind of image processing method that can effectively reduce calculated amount and can eliminate aliasing effect basically.
Fig. 5 shows the process flow diagram of image processing method according to an embodiment of the invention.
As shown in Figure 5, at step S502 place, can confirm the filling length of image.
For example, can be through confirm the filling length of image according to the method for the filling length of definite image of the foregoing description.
Then, at step S504 place, can come expanded images according to said filling length.
Specifically, in spatial domain, can on the image filtering direction, begin expanded images and fill according to said filling length from the position of image end.
Wherein, the image value of the expanding location of image can be filled to be 0.In this case, the size of image expanding position aliasing effect of generation in image filtering is handled is 0.
Although the image value of the expanding location of image is set to 0 in above description, above description is merely example, the invention is not restricted to this.In fact, also can be for expanding location is provided with other suitable image value, as long as it can reduce the influence of aliasing effect effectively.
For example, in another embodiment of the present invention, the image value of the expanding location of image also can be set to the image value of image spreading boundary.In this case, the aliasing effect of image spreading position generation is not 0.But because the image spreading position has been filled the boundary value at extended boundary place, so the exercising result of aliasing effect can be to make the image boundary place be able to strengthen and image subject is not caused interference.Thereby, in this case, also can reduce the influence of aliasing effect effectively.
In addition, above-mentioned image value can be the image parameter value of any appropriate.For example, can be gray-scale value etc.
Then, at step S506 place, can carry out filtering to the image after the expansion.
Owing to targetedly image is expanded,, on the other hand, can also be reduced calculated amount so when image expanding is carried out filtering, can reduce aliasing effect basically on the one hand.
Fig. 6 shows the particular flow sheet of Filtering Processing according to an embodiment of the invention.
As shown in Figure 6, at step S602 place, can will arrive frequency domain by the image transitions after the said filling extended length through Fourier transform.
Wherein, owing to confirmed in advance than short filling length and by said filling length image to be expanded according to embodiments of the invention, thereby in Fourier transform is handled, can reduce calculated amount.
Then, at step S604 place, can in frequency domain, utilize frequency domain filter that image is carried out Filtering Processing.
Because confirmed in advance than short filling length and by said filling length image to be expanded according to embodiments of the invention, therefore, in frequency domain filtering is handled, most aliasing effect will be eliminated.That is to say, in Filtering Processing, do not produce aliasing effect basically.
Then, at step S606 place, can return spatial domain through the image transitions of Fourier inversion after with Filtering Processing.
Obviously, owing to confirmed in advance image to be expanded, therefore, in Fourier inversion is handled, can reduce calculated amount than short filling length and by said filling length according to embodiments of the invention.
Thereby, can not produce basically under the lower situation of aliasing effect and calculated amount according to the image processing method of the foregoing description, realize Filtering Processing to image.
More than describe the image processing method that is used for image filtering and be merely example, the invention is not restricted to this.For example, the image processing method that is used for image filtering according to another embodiment of the invention can further include the processing that need to judge whether expanded images.
Specifically, filling length that can basis is determined after the processing of the filling length of confirming image judges whether the processing of needs expanded images.
Fig. 7 is the process flow diagram that the image processing method that comprises the processing that need to judge whether expanded images according to an embodiment of the invention is shown.
As shown in Figure 7, the filling length of definite image at step S702 place.
For example, can confirm the filling length on the image filtering direction according to the method for any the foregoing description.
Then, step S704 place judges whether that according to the filling length of determining the needs expanded images avoids aliasing effect.
Specifically, aliasing effect comes down to that the image afterbody caused the phase mutual interference of image header because image (is equivalent to and carries out circular convolution in the spatial domain) when frequency domain carries out Filtering Processing.Therefore; Can and whether be surpassed threshold value according to each image section that causes aliasing effect and judge whether that needs avoid aliasing effect through expanded images by the similarity between each image section that disturbs; Wherein, can rule of thumb be worth with actual and require to be provided with said threshold value.
As shown in Figure 7, if judgement need not expanded at step S704 place, then can directly carry out Filtering Processing at step S708 place.
If need expand, then can come the blank map picture according to said extension length at step S706 place in step S704 place's judgement.
Then, at step S708 place, can carry out Filtering Processing to image.
Wherein, step S502, S504, S506 among the embodiment that describes with reference to Fig. 5 before the processing of step S702, S706, S708 can be similar to no longer carry out repetition so that instructions keeps succinct at this.
Can find out,,, can after determining filling length, further not judge and need carry out image spreading, can reduce operand further thus for some images that need not fill through above-mentioned method.
Fig. 8 shows the process flow diagram that judges whether the processing of needs expanded images according to an embodiment of the invention according to filling length.
As shown in Figure 8, at step S802 place, the similarity between can two parts of computed image.Specifically, can calculate the image section that causes aliasing effect and the image section that disturbed by aliasing effect between similarity.
More specifically, in this embodiment, can with determine what cause most of aliasing effect is the image section that causes aliasing effect with filling that the corresponding image section of length regards as.In addition; Frequency domain is handled the median filter function and is caused with multiply each other (be equivalent in spatial domain and carry out circular convolution) of image function because aliasing effect is, so the image section that had the greatest impact by aliasing effect is first image section of when filtering of Z axle (for example Z axle) on spatial domain on the filtering direction.Thereby in this embodiment, can first image section of when filtering of Z axle (for example Z axle) on the filtering direction be regarded as is the image section that influenced by aliasing effect.
Obviously, if confirm to cause the image section and the similarity degree between the quilt image section that disturbs very high (surpassing threshold value) of aliasing effect, show that then (being aliasing effect) interfering with each other between them is very little, thereby can not judge and need carry out extension process; On the other hand; If confirm to cause the image section of aliasing effect and the similarity degree not high (less than threshold value) between the quilt image section that disturbs; Then show the stronger interference (being that aliasing effect is higher) of existence between them, thereby can judge need be through eliminating aliasing effect in the spatial domain extension process.
For example; In the embodiment shown in fig. 3; The image section of determining 6 positions in the image causes aliasing effect to account for 95% of total aliasing effect to image, that is to say, can regard the image section of these 6 positions as be to cause aliasing effect image section.Therefore, can judge whether the needs expanded images according to the similarity between the image section of the image section of these 6 positions and first position.
Referring to Fig. 8,, judge whether the needs expanded images according to similarity at step S804 place.Specifically, whether judge whether the needs expanded images according to similarity greater than threshold value.
For example, before carrying out image spreading, can whether surpass threshold value according to the similarity between the image section of each position among the image section of first position in the image and 6 positions of determining and judge whether the needs expanded images.
Owing to handle to be equivalent to the circular convolution on the spatial domain at frequency domain filtering, thereby on spatial domain, these 6 positions of causing 95% aliasing effect of determining among Fig. 3 come down to last 6 positions of image in spatial domain.
Thereby; If the similarity of the image section of each position in the image section of first position in the image and the image among last 6 positions all surpasses threshold value; The image section of first position is very similar and aliasing effect that they cause is very little in the image section that then shows last 6 positions in the image and the image, so can not judge at step S806 place and need carry out extension process.
On the other hand; If the similarity of a position in the image section of first position in the image and the image among last 6 positions or the image section of a plurality of positions is less than threshold value; Then show the part that has the aliasing effect that image is caused in last 6 positions in the image, therefore can judge and to carry out extension process at step S808 place.
In addition; In above-mentioned description; Compare according in the image section of first position in the image and the image section that causes most aliasing effect (with filling the corresponding image section of length) each, above description is merely example, the invention is not restricted to this.
For example; In the application of medical image, consider what medical image slowly changed often continuously on the filtering direction, therefore; Can be (in example shown in Figure 3 with the image section that causes most aliasing effect; Image section for last 6 positions) approximately regard as identically, thereby any part in the image section that only 6 picture positions of the image section of first position and this are located in the computed image (preferably; Can be the image section of last position) between similarity, judge whether the needs expanded images according to comparative result.Calculate similarity in this way, can further reduce calculated amount, raise the efficiency.
In above-mentioned processing, can adopt variety of way to come the similarity of two parts of computed image.
For example, be two dimensional image and the filtering direction is under the situation of one dimension (for example, in the situation of directions X filtering) at image, two image sections to be calculated can be two image strips of two dimensional image intercepting on directions X.
Again for example, be that 3-D view and filtering direction are that two image sections to be calculated can be two image sheets of 3-D view intercepting on the Z direction under the situation of one dimension (for example, in the situation of Z trend pass filtering) (for example, the situation of medical image) at image.
The method of calculating similarity to two image sections can have a variety of.In one embodiment of the invention, can calculate the similarity of two image sheets through following formula:
In above-mentioned formula, S representes the similarity of two image sheets, and N representes the number of the pixel in each image sheet, and gi representes the parameter value of i pixel in first image sheet, and si representes the parameter value of second i pixel in the image sheet.As stated, said parameter value for example can be the gray-scale value of pixel.
In above description, be that example is described to calculate similarity, but the invention is not restricted to this that above-mentioned formula can wait to image strip similarly and calculate similarity, repeats no more so that instructions keeps succinct at this to image sheet.
Similar with said method, embodiments of the invention also provide corresponding device thereof.
Fig. 9 illustrates the device that is used for confirming the filling length of image on the image filtering direction according to an embodiment of the invention.
As shown in Figure 9, can comprise aliasing effect function acquiring unit 902 and fill length and confirm unit 904 according to the device that is used for confirming the filling length of image on the image filtering direction of this embodiment.
Specifically, aliasing effect function acquiring unit 902 can obtain the aliasing effect function of frequency domain filter to said image; Confirm that unit 904 confirms the filling length of image on the image filtering direction based on said aliasing effect function and fill length.Wherein, as stated, the aliasing effect function can be to characterize because the function that the aliasing effect that the frequency domain filtering processing causes image changes on spatial domain.
The method that aliasing effect function acquiring unit 902 obtains the aliasing effect function can have a variety of.
For example, aliasing effect function acquiring unit 902 can obtain the aliasing effect function through the variation before and after the analysis image filtering and through analyzing this variation.
Perhaps, aliasing effect function acquiring unit 902 also can directly obtain the aliasing effect function of wave filter from outside (for example, from the memory storage of the aliasing effect function that stores each wave filter).
Fill length and confirm that unit 904 can confirm the filling length of image based on the aliasing effect function.More particularly, fill length filling length on spatial domain when confirming that unit 904 can be confirmed can aliasing effect be restricted to the degree that is enough to ignore according to the aliasing effect function.
Can find out, can confirm to fill length according to the aliasing effect function according to the device that is used for the definite filling length of image on the image filtering direction of the foregoing description, and needn't on spatial domain, image spreading be arrived long length.Therefore, can reduce on the one hand in the image filtering processing because the long bigger operand that image length brought; On the other hand, owing to considered the influence of aliasing effect,, guarantee picture quality so also can eliminate the influence of aliasing effect basically.
In one embodiment of the invention, aliasing effect function acquiring unit 902 can obtain the aliasing effect function in the variation before and after the filtering through analysis image.
Figure 10 illustrates the structural representation of aliasing effect function acquiring unit according to an embodiment of the invention.
Shown in figure 10, can comprise test pattern construction unit 1002, aliasing effect function generation unit 1004 according to the aliasing effect function acquiring unit 1000 of this embodiment.
Specifically, test pattern construction unit 1002 can make up the length test pattern identical with the length of this image, and aliasing effect function generation unit 1004 can obtain the aliasing effect function based on the variation of said test pattern before and after the Filtering Processing.
Wherein, said test pattern can be the image of arbitrary form, as long as it has on spatial domain and treats the filtering image equal length.For example, in some cases, can be directly with treating that filtering image is as test pattern.
In a preferred embodiment of the invention, said test pattern can be that the image that on the filtering direction, is rendered as the form of impulse function calculates to simplify.
At one more in the preferred embodiment; Said test pattern can be following this image: wherein; On the filtering direction (for example; Z axle under the one-dimensional filtering situation), the image value of first position of said image can be for unit strength (for example, 1) other position image value be 0.
To test pattern, aliasing effect function generation unit 1004 can obtain the aliasing effect function based on the variation of said test pattern before and after the image filtering.More particularly, aliasing effect function generation unit 1004 can be confirmed to be used to characterize frequency domain filtering according to the variation of the image value at the diverse location place of the test pattern after the Filtering Processing in spatial domain handle the aliasing effect function that the aliasing effect that causes changes on spatial domain.
As stated; At test pattern be under the situation of form of impulse function (wherein; The intensity of first pixel of this test pattern be unit strength 1 and the value of other pixel is 0), the variation of the diverse location that the test pattern after the inverse transformation can reflect image comparatively intuitively before and after image filtering.Thereby aliasing effect function generation unit 1004 can obtain the aliasing effect function according to this variation simple, intuitive ground.
More than utilize the test pattern of impulse function form to confirm that the description of aliasing effect function is merely example, the invention is not restricted to this.In fact, said test pattern can be arbitrary form.
For example, under the situation of the test pattern that adopts other form, aliasing effect function generation unit 1004 can be confirmed the aliasing effect function in the variation before and after the Filtering Processing according to each images of positions value in the test pattern.
Specifically, can calculate the variation of each images of positions value before and after Filtering Processing in the test pattern earlier; Then, confirm the aliasing effect function through mathematics manipulation (for example, dissection process etc.).
Filling length among Fig. 9 confirms that the implementation of unit also can have a variety of.
Figure 11 shows and fills the structural representation that length is confirmed unit 1100 according to an embodiment of the invention.
Shown in figure 11, fill length and confirm that unit 1100 can comprise aliasing effect size computing unit 1102, part aliasing effect computing unit 1104, fill length and confirm unit 1106.Specifically, aliasing effect size computing unit 1102 can calculate the size of total aliasing effect according to the aliasing effect function; Part aliasing effect computing unit 1104 can according to the aliasing effect function calculate with image in the size of the corresponding aliasing effect of diverse location; The ratio of aliasing effect that length confirms that unit 1106 is can basis corresponding with diverse location and total aliasing effect of filling is confirmed said filling length.
For example; At the aliasing effect function is under the situation of form shown in Figure 3; Aliasing effect size computing unit 1102 can obtain the size of total aliasing effect at entire image length range (for example, length range) corresponding with the length of filtering image integration with the aliasing effect function.
Part aliasing effect computing unit 1104 can obtain the size in the aliasing effect of corresponding position through the aliasing effect function is carried out integration about diverse location in the image.
Fill length and confirm that unit 1106 can calculate the ratio of aliasing effect corresponding with each picture position and total aliasing effect earlier, and with ratio the distance during above predetermined threshold as said filling length.Wherein, predetermined threshold can be according to circumstances and predetermined numerical value, for example can be 85%, 90%, 95%, 96%, 97%, 98%, 99% etc.
For example, in the embodiment shown in fig. 3, can find out that distance is that the aliasing effect that produced accounts for 95% of total aliasing effect by the distance of the image section of filtering at 6 o'clock.Therefore, being set in said threshold value under 95% the situation, filling length and confirm that unit 1100 can confirm as 6 with the filling length of image, is 0 so that eliminate the aliasing effect that causes these 6 positions.Therefore, can eliminate aliasing effect, eliminate the aliasing effect that when frequency domain filtering is handled, causes thus basically up to 95%.
In the above embodiments, can determine the filling length that to eliminate aliasing effect (for example, 95% in the foregoing description) basically, and need not can reduce the calculated amount in the follow-up Filtering Processing thus with image spreading to whole complete image length.
In addition; Based on the above-mentioned device that is used for confirming the filling length of image on the image filtering direction; Embodiments of the invention further provide a kind of image processing apparatus that comprises this device, and it can effectively reduce calculated amount and can eliminate aliasing effect basically.
Figure 12 shows the structural representation according to the image processing apparatus 1200 of this embodiment.
Shown in figure 12, image processing apparatus 1200 can comprise that filling length confirms device 1202, image spreading unit 1204 and wave filter 1206.
Specifically, fill length and confirm that device 1202 can be the device according to the filling length that is used for definite image of any the foregoing description.
Wherein, the image value of the expanding location of image can be filled to be 0.In this case, the size of image expanding position aliasing effect of generation in image filtering is handled is 0.
Although the image value of the expanding location of image is set to 0 in above description, above description is merely example, the invention is not restricted to this.In fact, also can be for expanding location is provided with other suitable image value, as long as it can reduce the influence of aliasing effect effectively.
For example, in another embodiment of the present invention, the image value of the expanding location of image also can be set to the image value of image spreading boundary.In this case, the aliasing effect of image spreading position generation is not 0.But because the image spreading position has been filled the boundary value at extended boundary place, so the exercising result of aliasing effect can be to make frontier point be able to strengthen and image subject is not caused interference.Thereby, in this case, also can reduce the influence of aliasing effect effectively.
In addition, above-mentioned image value can be the image parameter value of any appropriate.For example, can be gray-scale value etc.
1206 in wave filter can carry out filtering at frequency domain to the image after the above-mentioned expansion.
Owing to targetedly image is expanded,, on the other hand, can also be reduced calculated amount so when image expanding is carried out filtering, can reduce aliasing effect basically on the one hand.
More than describe the image processing apparatus that is used for image filtering and be merely example, the invention is not restricted to this.For example, at the image processing apparatus that is used for image filtering according to another embodiment of the invention, can further include the judging part that judges whether the needs expanded images according to said filling length.
Figure 13 shows the structural representation according to the image processing apparatus 1300 of this embodiment.
Shown in figure 13, image processing apparatus 1300 can comprise that filling length confirms device 1302, image spreading unit 1304, wave filter 1306 and judging part 1308.
Wherein, Fill length confirm device 1302, image spreading unit 1304, wave filter 1306 concrete ins and outs can referring to before combine Fig. 5 to describe filling length confirm device 1202, image spreading unit 1204, wave filter 1206, do not giving unnecessary details at this so that instructions keeps succinctly.
Judging part 1308 can confirm that the filling length that device 1302 provides judges whether the needs expanded images according to filling length.
If judging part 1308 judgements need expanded images, then notify image spreading unit 1304 to expand and the blank map picture according to this filling length; Do not need expanded images if judging part 1308 is judged, then notify image spreading unit 1304 image not to be expanded but image is offered wave filter 1306 and carry out Filtering Processing.
Can find out,, can judge further not after filling length that needs carry out image spreading determining, can reduce operand further thus according to the image processing apparatus 1300 of this embodiment for some images that need not fill
Judging part in the foregoing description can have multiple suitable embodiment.
Figure 14 shows the structural drawing of judging part according to an embodiment of the invention.
Shown in figure 14, comprise similarity calculating part 1402 and judgment processing portion 1404 according to the judging part 1400 of this embodiment.
Wherein, similarity calculating part 1402 can calculate the image section that causes aliasing effect and the image section that disturbed by aliasing effect between similarity.
In one embodiment of the invention; Can with determine what cause most of aliasing effect is the image section that causes aliasing effect with filling that the corresponding image section of length regards as, be the image section that influenced by aliasing effect and first image section of when filtering of Z axle (for example Z axle) on the filtering direction is regarded as.Thereby, similarity calculating part 1402 can the calculation of filtered direction on (and for example Z axle the during filtering of Z axle) first image section and corresponding to the similarity between each image section in the image section of filling length.
In above-mentioned description; Similarity calculating part 1402 compares according in the image section of first position in the image and the image section that causes most aliasing effect (with filling the corresponding image section of length) each; Above description is merely example, the invention is not restricted to this.
For example; In the application of medical image; Consider what medical image slowly changed often continuously on the filtering direction, therefore, can be (in example shown in Figure 3 with the image section that causes most aliasing effect; Image section for last 6 positions) approximate regard as identical; Thereby similarity calculating part 1402 is the similarity between any part (preferably, can be the image section of last position) in the image section at the image section of first position and this place, 6 picture positions in the computed image only.Calculate similarity in this way, can further reduce calculated amount, raise the efficiency.
After similarity was calculated end, judgment processing portion 1404 can carry out judgement according to the similarity that calculates.
Specifically; For example; To embodiment shown in Figure 3; If the similarity of the image section of each position in the image section of first position in the image and the image among last 6 positions all surpasses threshold value, the image section of first position is very similar and aliasing effect that they cause is very little in the image section that then shows last 6 positions in the image and the image, so can not judge at step S806 place and need carry out extension process.
On the other hand; If the similarity of a position in the image section of first position in the image and the image among last 6 positions or the image section of a plurality of positions is less than threshold value; Then show the part that has the aliasing effect that image is caused in last 6 positions in the image, therefore can judge and to carry out extension process at step S808 place.
Can adopt variety of way to come the similarity of two parts of computed image.
For example, be two dimensional image and the filtering direction is under the situation of one dimension (for example, in the situation of directions X filtering) at image, two image sections to be calculated can be two image strips of two dimensional image intercepting on directions X.
Again for example, be that 3-D view and filtering direction are that two image sections to be calculated can be two image sheets of 3-D view intercepting on the Z direction under the situation of one dimension (for example, in the situation of Z trend pass filtering) (for example, the situation of medical image) at image.
The method of calculating similarity to two image sections can have a variety of.In one embodiment of the invention, can calculate the similarity of two image sheets through following formula:
In above-mentioned formula, S representes the similarity of two image sheets, and N representes the number of the pixel in each image sheet, and gi representes the parameter value of i pixel in first image sheet, and si representes the parameter value of second i pixel in the image sheet.As stated, said parameter value for example can be the gray-scale value of pixel.
In above description, be that example is described to calculate similarity, but the invention is not restricted to this that above-mentioned formula can wait to image strip similarly and calculate similarity, repeats no more so that instructions keeps succinct at this to image sheet.
In addition, in the specific descriptions of device, omitted some and in method is described, carried out the content of describing in detail so that instructions keeps succinct.Thereby, about the ins and outs more specifically of device, can referring to before the description carried out of associated methods.
In addition, it will be appreciated that various examples as herein described and embodiment all are exemplary, the invention is not restricted to this.In this manual, statements such as " first ", " second " only are for described characteristic is distinguished on literal, clearly to describe the present invention.Therefore, should it be regarded as having any determinate implication.
Each forms module in the said apparatus, the unit can be configured through the mode of software, firmware, hardware or its combination.Dispose spendable concrete means or mode and be well known to those skilled in the art, repeat no more at this.Under situation about realizing through software or firmware; From storage medium or network the program that constitutes this software is installed to the computing machine with specialized hardware structure (multi-purpose computer 1500 for example shown in Figure 15); This computing machine can be carried out various functions etc. when various program is installed.
In Figure 15, CPU (CPU) 1501 carries out various processing according to program stored among ROM (read-only memory) (ROM) 1502 or from the program that storage area 1508 is loaded into random-access memory (ram) 1503.In RAM 1503, also store data required when CPU 1501 carries out various processing or the like as required.CPU 1501, ROM 1502 and RAM 1503 are connected to each other via bus 1504.Input/output interface 1505 also is connected to bus 1504.
Following parts are connected to input/output interface 1505: importation 1506 (comprising keyboard, mouse or the like), output 1507 (comprise display; Such as cathode ray tube (CRT), LCD (LCD) etc. and loudspeaker etc.), storage area 1508 (comprising hard disk etc.), communications portion 1509 (comprising that NIC is such as LAN card, modulator-demodular unit etc.).Communications portion 1509 is handled such as the Internet executive communication via network.As required, driver 1510 also can be connected to input/output interface 1505.Detachable media 1511 is installed on the driver 1510 such as disk, CD, magneto-optic disk, semiconductor memory or the like as required, makes the computer program of therefrom reading be installed to as required in the storage area 1508.
Realizing through software under the situation of above-mentioned series of processes, such as detachable media 1511 program that constitutes software is being installed such as the Internet or storage medium from network.
It will be understood by those of skill in the art that this storage medium is not limited to shown in Figure 15 wherein having program stored therein, distribute so that the detachable media 1511 of program to be provided to the user with equipment with being separated.The example of detachable media 1511 comprises disk (comprising floppy disk (registered trademark)), CD (comprising compact disc read-only memory (CD-ROM) and digital universal disc (DVD)), magneto-optic disk (comprising mini-disk (MD) (registered trademark)) and semiconductor memory.Perhaps, storage medium can be hard disk that comprises in ROM 1502, the storage area 1508 or the like, computer program stored wherein, and be distributed to the user with the equipment that comprises them.
The present invention also proposes a kind of program product that stores the instruction code of machine-readable.When said instruction code is read and carried out by machine, can carry out above-mentioned method according to the embodiment of the invention.
Correspondingly, the storage medium that is used for carrying the program product of the above-mentioned instruction code that stores machine-readable is also included within of the present invention open.Said storage medium includes but not limited to floppy disk, CD, magneto-optic disk, storage card, memory stick or the like.
At last; Also need to prove; Term " comprises ", " comprising " or its any other variant are intended to contain comprising of nonexcludability; Thereby make to comprise that process, method, article or the equipment of a series of key elements not only comprise those key elements, but also comprise other key element of clearly not listing, or also be included as this process, method, article or equipment intrinsic key element.In addition, under the situation that do not having much more more restrictions, the key element that limits by statement " comprising ... ", and be not precluded within process, method, article or the equipment that comprises said key element and also have other identical element.
Though more than combine accompanying drawing to describe embodiments of the invention in detail, should be understood that top described embodiment just is used to explain the present invention, and be not construed as limiting the invention.For a person skilled in the art, can make various modifications and change to above-mentioned embodiment and do not deviate from essence of the present invention and scope.Therefore, scope of the present invention is only limited appended claim and equivalents thereof.
Claims (22)
1. device that is used for confirming the filling length of image on the image filtering direction comprises:
Aliasing effect function acquiring unit is configured to obtain the aliasing effect function of frequency domain filter to said image, and wherein said aliasing effect function is to characterize frequency domain filtering to handle the function that the aliasing effect that causes changes on spatial domain;
Fill length and confirm the unit, be configured to confirm the filling length of said image on the image filtering direction based on said aliasing effect function.
2. device according to claim 1, wherein, said aliasing effect function acquiring unit comprises:
The test pattern construction unit is configured to make up the length test pattern identical with the length of this image;
Aliasing effect function generation unit is configured to obtain the aliasing effect function based on the variation of said test pattern before and after the Filtering Processing.
3. device according to claim 2, wherein, said test pattern is the image of following form: the image value in first position of image filtering direction is a unit value, and the image value of other position is 0.
4. device according to claim 2, wherein, said filling length confirms that the unit comprises:
Aliasing effect size computing unit is configured to calculate according to the aliasing effect function size of total aliasing effect;
Part aliasing effect computing unit, be configured to according to the aliasing effect function calculate with image in the size of the corresponding aliasing effect of diverse location;
Fill length and confirm the unit, the ratio that is configured to basis aliasing effect corresponding with diverse location and total aliasing effect is confirmed said filling length.
5. according to any one of claims 1 to 4 image processing apparatus, wherein, the medical image of said image for forming according to the data that obtain through medical diagnosis device.
6. image processing apparatus comprises:
Filling length that is used for the definite filling length of image on the image filtering direction according to any one among the claim 1-4 is confirmed device;
The image spreading unit is configured in spatial domain according to said filling length expanded images on the filtering direction;
Wave filter is configured to the image after the expansion is carried out filtering.
7. device according to claim 6 also comprises: judging part is configured to judge whether the needs expanded images according to said filling length.
8. device according to claim 7, wherein, judging part comprises:
The similarity calculating part, be configured to calculate receive each image section that aliasing effect influences and cause aliasing effect corresponding to the similarity between each image section in the image section of said filling length;
Judgment processing portion is configured to carry out judgement according to the similarity that calculates.
9. device according to claim 8, wherein, calculate the similarity of two parts according to following formula or its mathematic(al) manipulation:
Wherein, S representes similarity, and N representes the number of the pixel in each part, and gi representes the parameter value of an i pixel in the image section, and si representes the parameter value of i pixel in another part.
10. device according to claim 9, wherein said parameter value comprises gray-scale value.
11. like each described image processing apparatus in the claim 5 to 10, wherein, the medical image of said image for forming according to the data that obtain through medical diagnosis device.
12. the method for the filling length of a definite image on the image filtering direction comprises:
Obtain the aliasing effect function of frequency domain filter to said image, wherein said aliasing effect function is to characterize frequency domain filtering to handle the function that the aliasing effect that causes changes on spatial domain;
Confirm the filling length of said image on the image filtering direction based on said aliasing effect function.
13. method according to claim 12 wherein, is obtained frequency domain filter the processing of the aliasing effect function of this image is comprised:
Make up the length test pattern identical with the length of this image;
Said test pattern is carried out Filtering Processing;
Variation based on said test pattern before and after the Filtering Processing obtains the aliasing effect function.
14. method according to claim 13, wherein, said test pattern is the image of following form: the image value in first position of image filtering direction is a unit value, and the image value of other position is 0.
15. method according to claim 12 wherein, confirms that based on the aliasing effect function processing of the filling length of image comprises:
Calculate the size of total aliasing effect according to the aliasing effect function;
According to the aliasing effect function calculate with image in the size of the corresponding aliasing effect of diverse location;
Ratio according to aliasing effect corresponding with diverse location and total aliasing effect is confirmed said filling length.
16. like each described image processing method in the claim 12 to 15, wherein, the medical image of said image for forming according to the data that obtain through medical diagnosis device.
17. an image processing method comprises:
Confirm the filling length of image according to method any among the claim 12-15;
In spatial domain according to said filling length expanded images on the filtering direction;
Image to after the expansion carries out filtering.
18. method according to claim 17 also comprises: judge whether the needs expanded images according to said filling length.
19. method according to claim 18 wherein, judges whether that according to said filling length the processing of needs expanded images comprises:
Calculating receive each image section that aliasing effect influences and cause aliasing effect corresponding to the similarity between each image section in the image section of said filling length;
If all greater than threshold value, then judging, each similarity do not need expanded images;
Otherwise judgement needs expanded images.
20. method according to claim 19 wherein, is calculated the similarity of two parts according to following formula or its mathematic(al) manipulation:
Wherein, S representes similarity, and N representes the number of the pixel in each part, and gi representes the image value of an i pixel in the image section, and si representes the image value of i pixel in another part.
21. method according to claim 20, wherein said image value comprises gray-scale value.
22. like each described image processing method in the claim 15 to 21, wherein, the medical image of said image for forming according to the data that obtain through medical diagnosis device.
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2011101901209A CN102842114A (en) | 2011-06-21 | 2011-06-21 | Filling length determination device and method and image processing method and device |
| US13/529,173 US20120328172A1 (en) | 2011-06-21 | 2012-06-21 | Apparatus and method for determining padding length |
| JP2012140164A JP2013004105A (en) | 2011-06-21 | 2012-06-21 | Image processing system, method and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2011101901209A CN102842114A (en) | 2011-06-21 | 2011-06-21 | Filling length determination device and method and image processing method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN102842114A true CN102842114A (en) | 2012-12-26 |
Family
ID=47361898
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2011101901209A Pending CN102842114A (en) | 2011-06-21 | 2011-06-21 | Filling length determination device and method and image processing method and device |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20120328172A1 (en) |
| JP (1) | JP2013004105A (en) |
| CN (1) | CN102842114A (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2529157B (en) * | 2014-08-08 | 2017-06-28 | Marshall Wace Llp | Noise reduction method and system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6741643B1 (en) * | 1998-04-15 | 2004-05-25 | Telecommunications Research Laboratories | Asymmetric equalization system for data transmission |
| US20060110065A1 (en) * | 2004-11-23 | 2006-05-25 | Stmicroelectronics Asia Pacific Pte. Ltd. | Block filtering system for reducing artifacts and method |
| CN1806259A (en) * | 2003-06-17 | 2006-07-19 | 皇家飞利浦电子股份有限公司 | Primitive edge pre-filtering |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6990368B2 (en) * | 2002-04-04 | 2006-01-24 | Surgical Navigation Technologies, Inc. | Method and apparatus for virtual digital subtraction angiography |
| US8213687B2 (en) * | 2006-04-28 | 2012-07-03 | Hewlett-Packard Development Company, L.P. | Image processing methods, image processing systems, and articles of manufacture |
| US8014965B2 (en) * | 2007-02-21 | 2011-09-06 | ABB Bomen | Spectral estimation of spectrometers time-sampled signals using fast implementation of the reverse non-uniform discrete fourier transform |
| TW200933538A (en) * | 2008-01-31 | 2009-08-01 | Univ Nat Chiao Tung | Nursing system |
-
2011
- 2011-06-21 CN CN2011101901209A patent/CN102842114A/en active Pending
-
2012
- 2012-06-21 JP JP2012140164A patent/JP2013004105A/en active Pending
- 2012-06-21 US US13/529,173 patent/US20120328172A1/en not_active Abandoned
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6741643B1 (en) * | 1998-04-15 | 2004-05-25 | Telecommunications Research Laboratories | Asymmetric equalization system for data transmission |
| CN1806259A (en) * | 2003-06-17 | 2006-07-19 | 皇家飞利浦电子股份有限公司 | Primitive edge pre-filtering |
| US20060110065A1 (en) * | 2004-11-23 | 2006-05-25 | Stmicroelectronics Asia Pacific Pte. Ltd. | Block filtering system for reducing artifacts and method |
Non-Patent Citations (1)
| Title |
|---|
| R. VITULLI ET AL.: "Aliasing effects mitigation by optimised sampling grids and impact on image acquisition chains", 《IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM》, vol. 2, 31 December 2002 (2002-12-31), pages 979 - 981, XP010597897 * |
Also Published As
| Publication number | Publication date |
|---|---|
| US20120328172A1 (en) | 2012-12-27 |
| JP2013004105A (en) | 2013-01-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dosne et al. | Improving the estimation of parameter uncertainty distributions in nonlinear mixed effects models using sampling importance resampling | |
| LoVerde et al. | Effects of scale-dependent non-Gaussianity on cosmological structures | |
| Verhoosel et al. | Image-based goal-oriented adaptive isogeometric analysis with application to the micro-mechanical modeling of trabecular bone | |
| Torres-Roldan et al. | CSpace: An integrated workplace for the graphical and algebraic analysis of phase assemblages on 32-bit Wintel platforms | |
| Dietrich et al. | Comprehensive comparison of numerical relativity and effective-one-body results to inform improvements in waveform models for binary neutron star systems | |
| EP2972303B1 (en) | Systems and methods for improving direct numerical simulation of material properties from rock samples and determining uncertainty in the material properties | |
| CN107101714A (en) | A kind of transformer health Evaluation method based on multi-measuring point vibration signal characteristics | |
| CN105631090A (en) | Finite element model optimization device and method | |
| CN102755174B (en) | Beam-former, medical image system and the method for showing diagnostic image | |
| CN104990854B (en) | Determine the method and device of irreducible water saturation | |
| US20140270393A1 (en) | Systems and methods for improving direct numerical simulation of material properties from rock samples and determining uncertainty in the material properties | |
| Jung et al. | Enhancement of perivascular spaces using densely connected deep convolutional neural network | |
| CN102656581B (en) | Messaging device, information processing method and storage medium | |
| Ma et al. | Shearlet-based compressed sensing for fast 3D cardiac MR imaging using iterative reweighting | |
| Ismail et al. | Evaluation of quantitative MRCP (MRCP+) for risk stratification of primary sclerosing cholangitis: comparison with morphological MRCP, MR elastography, and biochemical risk scores | |
| Liravi et al. | A 2.5 D coupled FEM-BEM-MFS methodology for longitudinally invariant soil-structure interaction problems | |
| CN115201902A (en) | Fault intelligent identification method and system based on deep learning | |
| CN102842114A (en) | Filling length determination device and method and image processing method and device | |
| CN104463825A (en) | Apparatus and method for detecting objects in three-dimensional volumetric image | |
| US11656199B2 (en) | System and method for obtaining and analyzing flux leakage data in the inspection of oil and gas wells | |
| CN108135575A (en) | Diagnostic ultrasound equipment and method | |
| Bernal et al. | A multigrid-like algorithm for probabilistic domain decomposition | |
| Fata et al. | Elastic scatterer reconstruction via the adjoint sampling method | |
| Heinzl et al. | Fast estimation of optimal specimen placements in 3D X-ray computed tomography | |
| Toumi | A decomposition of orthogonally additive polynomials on Archimedean vector lattices |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20121226 |