Iterative filter framework for medical images

[0002] The field of the invention is medical imaging, and particularly, the filtering and enhancement of medical images to improve their visual quality.

[0003] The quality of medical images is a function of the imaging modality used and the particular method used to acquire the image data. As a general matter, regardless of the imaging modality used, the quality of medical images can be improved by acquiring redundant data which can be averaged to reduce the effects of random noise. Unfortunately, this solution results in an increased scan time that inconveniences the patient and is expensive. Another approach is to increase the power of the imaging system. In MRI this may take the form of a stronger polarizing field (e.g. 1.5 Tesla instead of 0.5 Tesla), in x-ray systems this can take the form of a higher power x-ray beam, and in ultrasound this may take the form of a more powerful rf amplifier and transducer. These measures drive up the cost of the equipment, and in the case of x-ray, increase the dose of ionizing radiation to the patient.

[0004] Attempts have been made to improve the imagability of the subject by injecting contrast agents into the patient. See for example, U.S. Pat. No. 4,834,964 of Rosen. However, injected contrast agents only improve a limited range of image characteristics, and because it is as an invasive technique, it is sometimes inappropriate for medical reasons.

[0005] Acquired medical images can also be processed to enhance their clinical value by modifying with the histogram or distribution of signal values on a global or local basis as described, for example, in U.S. Pat. No. 5,063,607. In other enhancement techniques, the gray scale range of each subimage region is stretched such that it covers the entire display range as, for example, in U.S. Pat. No. 4,991,092. However, histogram modifications which expand the dynamic range of the data, also increase the noise in the image. Local histogram modifications cause a blocking effect that results in a lack of uniformity over the entire image.

[0006] Images may be enhanced using convolution or filtering techniques. Such techniques include the amplification of selected frequency bands as illustrated in U.S. Pat. No. 5,072,314 of Chang. Others have used a combination of high and low pass filtering to enhance images as illustrated for example in U.S. Pat. No. 5,081,692 to Kwon or U.S. Pat. No. 4,972,256 to Hirosawa. However, global filtering techniques tend to blur the images and eliminate the lower frequency regions. This makes clinical evaluation of the images difficult.

[0007] To eliminate some of the drawbacks of global filtering or convolution, locally adjusted filtering may be used. See for example, U.S. Pat. No. 4,761,819 of Denison, et al.; U.S. Pat. No. 4,991,092 of Swon; U.S. Pat. No. 5,050,227 of Furusawa and “Adaptive Smoothing: A General Tool For Early Vision”, IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. 13, No. 6, June 1991, Saint-Marc, et al. However, the local filtering techniques have difficulty distinguishing between sudden image variations attributable to edges and sudden image variations attributable to noise. These techniques also fail to account for differences in edge direction and regional variance, producing an image which is overly smooth without consideration of the edges.

[0008] Restoration approaches may also be used in which the acquisition process is modeled and the degradation of the imaging process described mathematically. These techniques then attempt to invert the degradations using restoration techniques such as least squares, Bayesian, or Kalman filtering. However, the restoration methods require an accurate model for the acquisition process. Complicated acquisition processes, such as MRI imaging are too difficult to model accurately, and the parameters of a complicated model for a given image can require lengthy iterative computations.

[0009] U.S. Pat. No. 4,691,366 of Fenster, et al. uses filters which are adjusted to enhance long edges and attenuate noise and points. However, this technique requires an analysis of the imaging system in order to produce the appropriate filters. Such analysis is computationally intensive, time consuming, and prone to errors.

[0010] The present invention is a multi-pass filter framework which is employed to improve the quality of medical images. More particularly, the multi-pass filter framework includes storage means for storing the acquired image data, a plurality of separate filter functions for receiving the acquired image data in a pass therethrough to produce corresponding filtered image data, scaling means for multiplying the filtered image data by a constant (α), summing means for combining filtered and original image data by appropriate scaling to produce an output image; and iteration means for passing the image data through selected filter functions a plurality of times. The filter functions can take many different forms to tailor the multi-pass filter to many different imaging modalities and clinical applications.

[0011] A general object of the invention is to enhance the quality of medical images without the addition of special hardware or significantly increasing the image processing time. By using the multi-pass filtering of the image data, the separate filtering functions may be considerably simplified. For example, rather than filtering with 15×15 element filter kernels, the same result can be accomplished in three passes through the framework which uses a 3×1 element filter kernel. This translates into simplified and faster processing. Also, by using the multi-pass technique the filter can be changed for each pass to provide a much greater range of image effects.

[0012] A more specific object of the invention is to improve the quality of images acquired with an MRI system. Noise may be filtered from the image data using a directional smoothing filter that is sensitive to edges in the image. The edges are smoothed along them, but their sharpness is not reduced. Segmentation may be used to differentially filter regions with strong and weak edges. These different filter functions may be employed in different combinations during the multiple passes through the filter framework.

[0013] The foregoing and other objects and advantages of the invention will appear from the following description. In the description, reference is made to the accompanying drawings which form a part hereof, and in which there is shown by way of illustration a preferred embodiment of the invention. Such embodiment does not necessarily represent the full scope of the invention, however, and reference is made therefore to the claims herein for interpreting the scope of the invention.

[0014] FIG. 1 is a block diagram of an MRI system which employs the present invention; and

[0015] FIG. 2 is a flow chart of the functions performed by the MRI system of FIG. 1 to filter acquired images according to the preferred embodiment of the invention.

[0016] FIG. 3 is a block diagram of a one-pass filter according to the present invention;

[0017] FIG. 4 is a block diagram of a multi-pass first order filter using the one-pass filter of FIG. 3; and

[0018] FIG. 5 is a block diagram of a multi-pass second order filter using the filter of FIG. 4.

[0019] The present invention is a filter framework which can be configured by the user to implement a wide range of filtering operations on medical images. Referring particularly to FIG. 3, the framework is formed around a basic filter block F₁₁which is comprised of a filter F₁that operates on input image data input at 10. The output of filter F₁is scaled by a factor C₁and is combined at summer 12 with scaled input image data. The scaling factor for the input image data is (1-C). The output 14 of the filter block F₁₁is intermediate image data that has undergone one pass of first degree filtering. The filter F₁may be relatively simple in construction such as a smoothing or sharpening filter using a 3×1 filter kernel, or it may be a combination of different filters that can be selectively employed during each pass.

[0020] The filter framework uses the filter block in FIG. 3 to form a multi-pass filter F₂₁such as that shown in FIG. 4. The intermediate image data is fed back into the filter block for another pass, and the particular filter function and the scaling factor C used may be different during each pass. Any number (n) of passes may be performed, as illustrated in FIG. 4, and after the last pass, the filtered output image data is scaled by a scaling factor (d) and combined with scaled original input image data at summer point 16. The filtered image data is output at 18. As will be described in detail below, the preferred embodiment of the invention is a three-pass (i.e. n=3) filter for use with image data acquired with an MRI system.

[0021] With particularly noisy images it may be necessary to provide a second degree filter. Referring particularly to FIG. 5, this can be achieved with the filter framework of the present invention by concatenating multi-pass filter blocks such as filter F₂₁in FIG. 4. The output 18 of each multi-pass filter block F₂₁, F₂₂. . . F_2mis fed back through the framework to provide the input to the next filter block in the chain. The output 18 of the last filter block F_2mis scaled by a scaling factor (e) and combined with scaled original input image data at summing point 20. The filtered image data is output at 22. It should be apparent that the second order filter F₃₁can itself be concatenated to form still higher order filters as needed to provide the desired image effects.

[0022] While the summing points 12, 16 and 20 may be as simple as the arithmetic sum of the scaled pixel intensities, it is also possible to segment the image into different regions and combine the regions with the original image. Regions with strong edges, for example, may not be combined, whereas regions with weak edges may be combined with corresponding data in the original image.

[0023] Another aspect of the present invention is the ability to configure the filter framework to provide the desired filter functions. A set of parameters control the filter framework to provide the filter blocks and concatenate them as described above. These parameters select the particular filter function to be performed during each pass, and select the value of the scaling factors c₁-c_n, d₁-d_n, etc. Stored sets of parameters may be used to configure the filter framework for particular medical applications.

[0024] Referring first to FIG. 1, there is shown the major components of a preferred MRI system which incorporates the present invention. The operation of the system is controlled from an operator console 100 which includes a keyboard and control panel 102 and a display 104. The console 100 communicates through a link 116 with a separate computer system 107 that enables an operator to control the production and display of images on the screen 104. The computer system 107 includes a number of modules which communicate with each other through a backplane. These include an image processor module 106, a CPU module 108 and a memory module 113, known in the art as a frame buffer for storing image data arrays. The computer system 107 is linked to a disk storage 111 and a tape drive 112 for storage of image data and programs, and it communicates with a separate system control 122 through a high speed serial link 115.

[0025] The system control 122 includes a set of modules connected together by a backplane. These include a CPU module 119 and a pulse generator module 121 which connects to the operator console 100 through a serial link 125. It is through this link 125 that the system control 122 receives commands from the operator which indicate the scan sequence that is to be performed. The pulse generator module 121 operates the system components to carry out the desired scan sequence. It produces data which indicates the timing, strength and shape of the RF pulses which are to be produced, and the timing of and length of the data acquisition window. The pulse generator module 121 connects to a set of gradient amplifiers 127, to indicate the timing and shape of the gradient pulses to be produced during the scan. The pulse generator module 121 also receives patient data from a physiological acquisition controller 129 that receives signals from a number of different sensors connected to the patient, such as ECG signals from electrodes or respiratory signals from a bellows. And finally, the pulse generator module 121 connects to a scan room interface circuit 133 which receives signals from various sensors associated with the condition of the patient and the magnet system. It is also through the scan room interface circuit 133 that a patient positioning system 134 receives commands to move the patient to the desired position for the scan.

[0026] The gradient waveforms produced by the pulse generator module 121 are applied to a gradient amplifier system 127 comprised of G_x, G_yand G_z, amplifiers. Each gradient amplifier excites a corresponding gradient coil in an assembly generally designated 139 to produce the magnetic field gradients used for position encoding acquired signals. The gradient coil assembly 139 forms part of a magnet assembly 141 which includes a polarizing magnet 140 and a whole-body RF coil 152. A transceiver module 150 in the system control 122 produces pulses which are amplified by an RF amplifier 151 and coupled to the RF coil 152 by a transmit/receive switch 154. The resulting signals radiated by the excited nuclei in the patient may be sensed by the same RF coil 152 and coupled through the transmit/receive switch 154 to a preamplifier 153. The amplified NMR signals are demodulated, filtered, and digitized in the receiver section of the transceiver 150. The transmit/receive switch 154 is controlled by a signal from the pulse generator module 121 to electrically connect the RF amplifier 151 to the coil 152 during the transmit mode and to connect the preamplifier 153 during the receive mode. The transmit/receive switch 154 also enables a separate RF coil (for example, a head coil or surface coil) to be used in either the transmit or receive mode.

[0027] The NMR signals picked up by the RF coil 152 are digitized by the transceiver module 150 and transferred to a memory module 160 in the system control 122. When the scan is completed and an entire array of data has been acquired in the memory module 160, an array processor 161 operates to Fourier transform the data into an array of image data. This image data is conveyed through the serial link 115 to the computer system 107 where it is stored in the disk memory 111. In response to commands received from the operator console 100, this image data may be archived on the tape drive 112, or it may be further processed by the image processor 106 in accordance with the present invention and conveyed to the operator console 100 and presented on the display 104.

[0028] For a more detailed description of the transceiver 150, reference is made to U.S. Pat. Nos. 4,952,877 and 4,922,736 which are incorporated herein by reference.

[0029] The present invention is implemented in the image processor 106 shown in FIG. 1 by carrying out the sequence of steps indicated by the flow chart in FIG. 2 on an image stored in digital form. The filtering process provides a multi-pass (iterative) framework with options for changing the parameters used during each pass. As shown by process block 200, eighteen parameters that control the filter process are input by the operator for the particular application. As will be described in more detail below, these parameters control the filtering effect which can range from complete smoothing of the image to complete sharpening of image features, as well as segmentation-based adaptive smoothing and sharpening. In research settings, the separate parameters may be specified by the operator, whereas in clinical settings, a group of stored parameters may be selected by the operator depending on the particular application.

[0030] As indicated at process block 202, the next step is to initialize the system for the multi-pass filter operation to be performed, and to scale the input image data. The main reason for scaling is to make the filtering functions independent of the dynamic range of the image data or its DC offset value. More specifically, the following steps are performed:

[0031] 1) Internal parameters which will be discussed in the next section are initialized and the iteration counter is set to 1;

[0032] 2) Maximum (MAX_ORIGINAL)and minimum (MIN_ORIGINAL) image intensity values are determined;

[0033] 3) If the iteration number is 1, the original data is save as I₂;

[0034] 4) scale=4095.0/MAX_ORIGINAL;

[0035] 5) Scaled image I is obtained using the relation: I=(I-MIN_ORIGINAL)* scale;

[0036] 6) Save the prefilteration image as I₁=I;

[0037] 7) Compute the average image intensity value (MEAN_BEFORE) of I before filtering.

[0038] Because the MR image data has 12-bit precision, 4095/MAX_ORIGINAL is used for scaling. The image I₁is used in an intermediate blending operation described below, while the stored image I₂is used for a final blending operation also to be described below.

[0039] As indicated by decision block 204, either of two directional smoothing functions is performed on the image data. The choice is one of the 18 control parameters described above. One choice (“1”) is a directional smoothing function indicated by process block 206 in which a simultaneous direction smoothing operation is applied to the image I. This simultaneous smoothing operation is described in L. S. Davis, “A Survey of Edge-Detection Techniques,” Computer Graphics and Image Processing, volume 4, pp. 248-270, 1975. It employs the following steps:

[0040] 1) Scan the image in the raster format one time.

[0041] 2) At each image pixel, consider four 1-dimensional 3×1 filters in all the four directions (45°, 135°, 90°, 0°).

[0042] 3) At each direction, compute the mean value and the variance of the 3 corresponding pixels.

[0043] 4) Determine the direction which gives the minimum variance and choose the corresponding mean value for the subsequent calculations.

[0044] 5) This new pixel value at the given location (x,y) is thus modified according to the equation:

[0045] I(x,y)=mean.

[0046] By design, this process smooths along the line structures and edges but not across them. The border pixels in a given image have to be treated differently and are set to zero in the processed image.

[0047] The other choice (“2”) at decision block 204 is a sequential smoothing function indicated by process block 208. The smoothing function is applied to the image I and its function is described in P. Chan and J. S. Lim, “One-Dimensional Processing for Adaptive Image Restoration,” IEEE Transactions In Acoustics, Speech and Signal Processing, volume 33, pp. 117-126, Feb. 1985. This sequential smoothing function employs the following steps:

[0048] 1) Scan the image in the raster format four times.

[0049] 2) At each pixel, consider four 1-dimensional 3×1 filters in all the four directions (45°, 135°, 90°, 0°) sequentially.

[0050] 3) For one direction, compute the mean value and the variance of the 3 corresponding pixels and use these for the subsequent calculations.

[0051] 4) If the variance (var) is greater than a preselected noise variance (nvar), then set var=var−nvar, otherwise, set var=0.

[0052] 5) The new pixel value at the given location (x,y) is thus modified according to the equation:

[0053] 0(x,y)=mean+(I(x,y)−mean)*(var/var+nvar),

[0054] where 0 (x,y) is a temporary image. After all the points in the image are processed, set I(x,y)=0(x,y) at every point (x,y).

[0055] 6) Repeat steps 3-5 for all the remaining direction angles.

[0056] With this filter function the amount of smoothing can be controlled using a single parameter, nvar. Since this is a sequential direction filtering, it is not cumulative. However, this approach has been found to produce minimal smoothing across edges while producing an overall stronger smoothing effect while preserving edges. The border pixels in a given image have to be treated differently and are set equal to zero in the processed image.

[0057] Regardless of which smoothing function is performed, the next step in the filter framework is to determine if the image edges are to be sharpened as indicated at process block 210. If they are, a directional sharpening and segmentation function is performed at process block 212. The functions performed by this filter function include the following:

[0058] 1) The 1-D Laplacian of the image at every pixel in each of the 4 directions (0°, 45°, 90°, 135°) is obtained using the equation:

[0059] E(k)=2.0*I(k)=I(k−1)−I(k+1);

[0060] where the index x refers to the current location along a given direction, E(k) is the edge strength and I(k) is the intensity value at the pixel.

[0061] 2) After computing all the four edge strengths at a given pixel, determine the maximum directional edge strength and use it in the subsequent steps as the edge strength E(x,y) at that location.

[0062] 3) Continue this process for all the pixels in the image. Border pixels in a given image have to be treated differently and are set equal to zero for the subsequent steps.

[0063] 4) Smooth the edges of image E(x,y) to produce a smoothed edge strength image ES(x,y) by performing the following steps:

[0064] a) Scan the image E(x,y) in the raster format one time.

[0065] b) At each pixel, consider four 1-dimensional 3×1 filters in all the four directions (45°; 135°, 90°, 0°).

[0066] c) At each direction, compute the mean value and the variance of the 3 corresponding pixels.

[0067] d) Determine the direction which gives the minimum variance and choose the corresponding mean value for the subsequent calculations.

[0068] e) The new pixel value at the given location (x,y) is thus modified according to the equation:

[0069] ES(x,y)=mean.

[0070] 5) Determine the maximum and the minimum values of ES(x,y). Select the greater of the absolute value of the maximum and the minimum as parameter MAX.

[0071] 6) Compute the threshold value using the equation:

[0072] threshold=MAX/threshold_parameter,

[0073] where threshold_parameter is one of the eighteen selectable parameters.

[0074] 7) In this step computations are performed based on decisions made using another one of the eighteen selectable parameters called segment₁₃parameter. The decision logic also creates a mask called mask(x,y). The parameters HPF_attenuated_parameter and HPF_unattenuated_parameter are also user selectable. These are used to essentially control the amount of edge strength image ES(x,y) added back to the smoothed image I(x,y) based on the segmentation criterion as determined in the following:

[0075] a) If segment₁₃parameter=1, perform the following:

[0076] if (ES(x,y)≧=threshold)

[0077] set mask(x,y)=1; and

[0078] set I(x,y)=I(x,y)+ES(x,y)*HPF_attenuated_parameter; or

[0079] if ES(x,y)>threshold)

[0080] set mask(x,y)=0; and

[0081] set I(x,y)=I(x,y)+ES(x,y)*HPF_unattenuated_parameter;

[0082] b) If segment_parameter=2, perform the following:

[0083] if (ES(x,y))>=-threshold)

[0084] set mask(x,y)=1; and

[0085] set I(x,y)=I(x,y)+ES(x,y)*HPF_attenuated_parameter; or

[0086] if (ES (x,y)<-threshold)

[0087] set mask(x,y)=0; and

[0088] set I(x,y)=I(x,y)+ES(x,y)*HPF_unattenuated_parameter;

[0089] c) If segment_parameter=3, perform the following:

[0090] if(abs(ES(x,y))<=threshold)

[0091] set mask(x,y)=1; and

[0092] set I(x,y)=I(x,y)+ES(x,y)*HPF_attenuated_parameter; or

[0093] if(abs ES(x,y))=threshold)

[0094] set mask(x,y)=0; and

[0095] set I(x,y)=I(x,y)+ES(x,y)*HPF_unattenuated_parameter.

[0096] At the completion of the sharpening and segmentation process 212 we have differentially enhanced strong edges and weak edges and the image is separated into strong edge regions and weak edge regions.

[0097] As indicated by process block 214, the next function to be performed is a blending and rescaling operation in which the filtered imaged I(x,y) is combined with the stored prefiltered image I₁(x,y) and the combined image is rescaled. This operation includes the following steps.

[0098] 1) The filtered image I(x,y) and the prefilteration image I₁(x,y) are blended using the equation:

[0099] I(x,y)=alpha*(I(x,y)−I₁,(x,y))+I₁(x,y);

[0100] where alpha is one of the selectable parameters such that 0<alpha<1.

[0101] 2) Compute the average pixel intensity (MEAN_AFTER) in the blended image I(x,y).

[0102] 3) Compute the normalization factor NORM_FACTOR=MEAN_BEFORE/MEAN_AFTER.

[0103] 4) Compute the normalized image using: I(x,y)=(I(x,y)*NORM_FACTOR)+MIN_ORIGINAL, where MIN_ORIGINAL is the minimum intensity of the original image.

[0104] 5) Set the value of image I₁(x,y) to be equal to image I(x,y).

[0105] At the completion of process 214, a test is made at decision block 216 to determine if all the passes through the filter framework are complete. In the preferred embodiment three passes are performed by branching back to process block 202 to repeat all the selected filter functions on the image I(x,y). The filter functions performed need not be the same during each pass. For example, it has been found that the best results are achieved when the simultaneous directional smoothing function 202 is performed during the first pass and the directional sharpening function 212 is performed during the final pass. Other than these, the operator can choose any functions during each iteration to achieve the desired effects.

[0106] The final function performed by the filter framework is a segmentation-based blending as indicated at process block 218. The purpose of this function is to maintain the texture of the original image in regions where there are no edges (i.e. low edge strength). The final image I_f(x,y) is obtained by blending the filtered image I(x,y) and the prefilteration image I₂(x,y) using the equation:

[0107] if mask (x,y)=1,

[0108] set I_f(x,y)=alpha_f*I(x,y)+beta_f*I₂(x,y); or

[0109] if mask(x,y)=0,

[0110] set I_f(x,y)=I(x,y);

[0111] where alpha_f is the scaling factor, and alpha_f, and beta_f are two of the eighteen selectable parameters and mask (x,y) is produced by process block 212 described above. Depending upon the selection of the parameters alpha_f and beta_f used in this process, one can obtain adaptive grey level compression and expansion. Various visual effects can thus be obtained as follows:

[0112] if (alpha_f+beta_f)=1, no grey scale compression or expansion is performed, (i.e. beta_f=1−alpha_f);

[0113] if (alpha_f+beta_f)>1, grey scale compression occurs to provide an overall smoothing effect.

[0114] If (alpha_f+beta_f)<1, contrast enhancement occurs to provide improved contrast.

[0115] There are 18 user selectable parameters for use in the preferred embodiment of the filter framework. Various visual effects are produced by different combinations of these parameters. In the following description the parameters are indicated as P[i], where i indicates the parameter number.

[0116] P[1]:

[0117] Description: Sets the upper limit for the sharpening process. This can be set equal to 4095 for 12-bit images.

[0118] Range: 4000-6000;

[0119] P[2]:

[0120] Description: Dampening factor to prevent weak edge amplification. This is set to 0.1 in the preferred embodiment.

[0121] Range: 0<P[2]<1.

[0122] P[3] and P[4]:

[0123] Description: Fractional texture blended from original image to weak edge regions of filtered image in the final stage of processing as determined by the segmentation mask. Set P[3 ] and P[4] appropriately to avoid “smushiness” while providing overall smoothness.

[0124] Range: 0<=P[3]<=10, 0<=P[4]<=10, where “=” stands for “less than or equal to.”

[0125] P[5]:

[0126] Description: Multiplication factor (HPF_unattenuated_factor) to the edge image for adding edges back to the filtered image in the final pass. Select this parameter high (>0.2) to increase overall sharpness of strong edges. Never set this parameter equal to zero to avoid missing out on the segmentation step.

[0127] Range: 0.00001 <=P[5]<=5.

[0128] P[6]:

[0129] Description: Multiplication factor (HPF_unattenuated_factor) to the edge image for adding edges back to the filtered image in the penultimate iteration (pass number 2 in the preferred embodiment). Use this parameter if additional edge enhancement boost is needed for tasks such as making the dark structures darker and white edges whiter. This parameter can also be used to put some texture back in images while providing higher overall smoothness to the image.

[0130] Range: 0.0 <=P[6]<=5.

[0131] P[7], P[8], P[9]:

[0132] Description: Variance threshold for the first, the second and the third passes to control the amount of smoothing. These values are useful if and only if directional smoothing 2 is used. Setting these values too high tends to blur the image.

[0133] Range: 0.1<=P[7]<=10000000, 0.1<=P[8]<=10000000, 0.l<=P[9]<=10000000.

[0134] P[10], P[11], P[12]:

[0135] Description: These parameters represent the fractional post-pass filtered image blended in with pre-pass filtered image for each of the three passes.

[0136] Range: Range: 0<=P[10]<=1, 0<=P[11]<=1, 0<=P[12]<=1.

[0137] P[13], P[14], P[15]:

[0138] Description: Dampening factors (HPF_attenuated_factors) to prevent weak edges amplification in each of three iterations;

[0139] Range: Range: 0<=P[10]<−1, 0<=P[11]<=1, 0<=P[12]<=1.

[0140] P[16]:

[0141] Description: Control parameter to decide during which iteration directional smoothing 2 (block E) should be done instead of directional smoothing 1 (block D). Always use directional smoothing 1 for the first pass and use either of the smoothing schemes for the remaining passes.

[0142] Range: 1<=P[16]<=4.

[0143] P[17]:

[0144] Description: Reciprocal of this parameter (threshold_parameter) is multiplied to the maximum edge strength to determine the threshold for segmentation of the image into strong and weak edge regions.

[0145] Range 1<=P[17]<=200.

[0146] P[18]:

[0147] Description: This parameter (segment_parameter) determines which of the three segmentation modes is to be applied for filtering. The parameter is set equal to 1 when the user needs the white edges but not the black edges to be emphasized. The parameter is set to 2 when the black edges but not the white edges have to be emphasized. The parameter is set to 3 when all edges have to be emphasized.

[0148] Range: 1<=P[18]<=3.

[0149] While there are millions of possible parameter combinations, all of them are not useful. Therefore, it is important to tune parameters to get the desired filtering effects. In order to achieve these effects, the following guidelines can be used:

[0150] Strong edge regions: In general, for noisy images, the strong edges must be emphasized or kept the same as in the original. Smoothing can be performed along a strong edge, but should never be performed across it. This is important to avoid blurring of details in the image which is undesirable in medical images.

[0151] Weak edge region: As a rule of thumb, weak edge regions must be attenuated. Sharpening the weak edge regions in general increases the chance of noise amplification and the appearance of “worms” in MR images. The texture in these regions must be maintained to have natural appearance in order to avoid “smudging” and “smushing” of pixel intensities.

[0152] The following are exemplary parameter sets for use in the preferred embodiment of the invention to achieve specific visual effects.

[0153] Short description: Sharpen with very little smoothing.

[0154] Detailed description: This is an overall filter with very little smoothing to enhance structural edges.

[0155] Suggested usage: High signal, very low noise 256×256 pixel images such as head, abdomens and extremities.

[0156] Parameter Combination

[0157] P[1]=4095;

[0158] P[2]=0.1;

[0159] P[3]=0.0;

[0160] P[4]=1.0;

[0161] P[5]=0.75;

[0162] P[6]=0.0;

[0163] P[7]=800000.;

[0164] P[8]=8000.;

[0165] P[9]=8000.;

[0166] P[10]=0.45;

[0167] P[11]=0.45;

[0168] P[12]=0.45;

[0169] P[13]=0;

[0170] P[14]=0;

[0171] P[15]=0.04;

[0172] P[16]=4;

[0173] P[17]=15;

[0174] P[18]=3;

[0175] Short description: Sharpen with smoothing.

[0176] Detailed Description: This is an overall sharpening filter with smoothing to enhance structural edges and suppress noise.

[0177] Suggested usage: Extremities.

[0178] Parameter Combination

[0179] P[1]=4095;

[0180] P[2]=0.1;

[0181] P[3]=0.7;

[0182] P[4]=0.3;

[0183] P[5]=0.6;

[0184] P[6]=0.45;

[0185] P[7]=800000.;

[0186] P[8]=8000.;

[0187] P[9]=8000.;

[0188] P[10]=0.45;

[0189] P[11]=0.45;

[0190] P[12]=0.45;

[0191] P[13]=0;

[0192] P[14]=0;

[0193] P[15]=0.04;

[0194] P[16]=3;

[0195] P[17]=15.;

[0196] P[18]=2;

[0197] Short description: High smooth and very little sharpen.

[0198] Detailed description: Overall this filter produces smoothing and very little sharpening to suppress noise and provides smoother structural edges.

[0199] Suggested usage: High resolution, low signal images.

[0200] Parameter Combination

[0201] P[1]=4095;

[0202] P[2]=0.1;

[0203] P[3]=0.8;

[0204] P[4]=0.45;

[0205] P[5]=0.01;

[0206] P[6]=0.1;

[0207] P[7]=800000.;

[0208] P[8]=8000.;

[0209] P[9]=8000.;

[0210] P[10]=1.0;

[0211] P[11]=1.0;

[0212] P[12]=1.0;

[0213] P[13]=0;

[0214] P[14]=0;

[0215] P[15]=0.001;

[0216] P[16]=2;

[0217] P[17]=1.0;

[0218] P[18]=1;

[0219] The invented filter provides a framework of selectable functions which when used in a multi-pass manner provide a variety of filtering functions that are performed in an efficient manner using 3=1 one-dimensional computations. This framework is capable of producing smoothing and sharpening depending upon edge strength of the region and it therefore produces varying filtration effects which are considered useful by radiologists. Since the weak edge regions are attenuated but never completely removed and strong edges are always smoothed along but never across, visually pleasing results can be obtained. Furthermore, the framework provides multilevel control to obtain sharpness and smoothness. Examples of the parameter combinations illustrate various visual effects that can be achieved using the present invention with high field MR images. Results indicate the proposed framework will also work well for all MRI systems with different strengths, but the parameter values should be separately tuned to obtain optimal results.