This article has Open Peer Review reports available.

# Masked smoothing using separable kernels for CT perfusion images

- David S Wack
^{1}Email author, - Kenneth V Snyder
^{2}, - Kevin F Seals
^{3}and - Adnan H Siddiqui
^{2}

**14**:28

https://doi.org/10.1186/1471-2342-14-28

© Wack et al.; licensee BioMed Central Ltd. 2014

**Received: **26 July 2013

**Accepted: **7 August 2014

**Published: **21 August 2014

## Abstract

### Background

CT perfusion images have a high contrast ratio between voxels representing different anatomy, such as tissue or vessels, which makes image segmentation of tissue and vascular regions relatively easy. However, grey and white matter tissue regions have relatively low values and can suffer from poor signal to noise ratios. While smoothing can improve the image quality of the tissue regions, the inclusion of much higher valued vascular voxels can skew the tissue values. It is thus desirable to smooth tissue voxels separately from other voxel types, as has been previously implemented using mean filter kernels. We created a novel Masked Smoothing method that performs Gaussian smoothing restricted to tissue voxels. Unlike previous methods, it is implemented as a combination of separable kernels and is therefore fast enough to consider for clinical work, even for large kernel sizes.

### Methods

We compare our Masked Smoothing method to alternatives using Gaussian smoothing on an unaltered image volume and Gaussian smoothing on an image volume with vascular voxels set to zero. Each method was tested on simulation data, collected phantom data, and CT perfusion data sets. We then examined tissue voxels for bias and noise reduction.

### Results

Simulation and phantom experiments demonstrate that Masked Smoothing does not bias the underlying tissue value, whereas the other smoothing methods create significant bias. Furthermore, using actual CT perfusion data, we demonstrate significant differences in the calculated CBF and CBV values dependent on the smoothing method used.

### Conclusion

The Masked Smoothing is fast enough to allow eventual clinical usage and can remove the bias of tissue voxel values that neighbor blood vessels. Conversely, the other Gaussian smoothing methods introduced significant bias to the tissue voxels.

## Background

CT perfusion imaging uses many high resolution scans in a dynamic series to determine parametric image maps of Cerebral Blood Flow (CBF), Cerebral Blood Volume (CBV), and Time to Peak (TTP), among other data types. A characteristic of CT image volumes is the high contrast ratio of voxel intensity values located in skull (or calcified regions) versus tissue regions, which can exceed 15:1. Furthermore, with the injection of a tracer, voxels representing vascular regions may have intensity values greater than four times higher than neighboring tissue regions. Kudo et al. demonstrated that the inclusion of vascular voxels could overestimate CBF [1]. The SNR within tissue regions is relatively low. Spatial smoothing is often applied to trade high spatial resolution for improved SNR characteristics. However, regular smoothing overestimates many tissue voxels due to nearby, high-valued vascular voxels.

While the importance of smoothing has been noted in the literature, it usually receives little discussion [2–4]. Klotz and König gave a brief but important description of their smoothing method as a “running mean smoothing procedure that operates separately on brain and vascular pixels” [5] (pg 173). As such, their approach operated in 2D and avoided blurring from smoothing high valued vascular pixels into tissue regions. Our method also operates separately on brain and vascular pixels, however we use a Gaussian kernel. Furthermore, our Masked Smoothing method can execute quickly, even when applied as 3D, by utilizing a combination of separable kernels. This offers an improvement in execution time of a few orders of magnitude relative to what could be achieved otherwise.

A “separable” 3D smoothing kernel can be expressed as the outer product of three vectors, and 3D smoothing can be applied as three successive 1D smoothings in the x, y, and z directions. While Gaussian kernels and mean kernels are separable, they do not remain separable if they must exclude vascular voxels. Our method overcomes this hurdle.

### Related methods

Many smoothing methods are “adaptive” [6] and arrive at an optimal solution through the progressive refinement of an initial solution. Some methods preserve edges [7, 8], similar to our desire to separate vessel and tissue voxels. Other methods consider the first or second spatial derivatives [9–12] or use the Discreet Cosine Transform [13]. A strength is these methods do not need a priori knowledge, such as voxel classifications [8]. A 4D extension of bilateral filters varies the weight of neighboring voxels according to distance and intensity, or “similarity” differences [14, 15], and has been applied to CT perfusion scans [16]. The TIPS (Time Intensity Profile Similarity) bilateral filter method [17] calculates the similarity of neighboring pixels across all image frames. While this reduces processing time to some degree, the TIPS bilateral smoothing kernel is not strictly separable. While TIPS offers great flexibility in expressing the smoothing formulation, its execution time [17], even applied as a 2D filter, is much slower than what can be achieved using separable 3D kernels.

There are two advantages of CT perfusion imaging over most other image smoothing problems. First, there are multiple image volumes such that voxels in the same spatial location will have the same classification. The second is that there are extreme voxel intensity differences between voxels of different classifications for some image volumes. While thresholding the mean image and the difference of the maximum and minimum images is a simple but powerful way of identifying vascular voxels, more sophisticated methods have been presented for the identification of arteries and veins [18]. Hence Masked Smoothing makes use of the easy access to a mask image of the tissue regions that adaptive smoothing or bilateral filters fail to utilize.

### Masked smoothing algorithm

Smoothing methods typically use a weighted sum of voxels within the smoothing neighborhood of a given tissue voxel, V_{x}, to assign a new value to V_{x}. The weights are all nonnegative and sum to one. The smoothing neighborhood for a given tissue voxel will, in general, include voxels of different segmentation classes--such as a high valued vessel voxel. This inclusion will have a tendency to artificially increase the smoothed value found at V_{x} from the true underlying tissue value. Our goal is to apply smoothing by only using voxels of like classes. Excluding voxels of a different class could be achieved by setting their weight values to zero, while rescaling the weights of same class voxels so they sum to one.

We define sum of weights (SW) for V_{x} as the sum of all weights of voxels that are both within the area of the smoothing kernel and the mask of the same tissue region (without rescaling). SW will equal one if the all the voxels within the smoothing neighborhood of V_{x} are all of the same class as V_{x}. Otherwise SW will be less than one. The reciprocal of SW (1/SW) can be used to rescale the weights so that they sum to one.

Setting some weights to zero and rescaling the remaining weights associated with each voxel within the smoothing neighborhood of Vx is computationally cumbersome. We can simplify the computation by making two changes: 1) Rather than resetting the smoothing kernel weight values of voxels outside of our tissue mask to zero, we instead set voxel values outside of our tissue mask to zero. 2) Rather than rescaling the individual weights within our mask by 1/SW, we rescale the weighted sum of voxel values by (1/SW), employing the distributive property. That is, if SW is known for each voxel, then smoothing the image with non-tissue voxels set to zero and dividing voxel by voxel by SW will result in the desired with-in class masked smoothing. Post-smoothing, non-tissue voxels can be set to zero, replaced by their original values, or smoothed separately.

Fortunately, calculating SW for each tissue voxel is easy. *SW for each tissue voxel is the result of applying the smoothing kernel to the binary mask that designates voxels classified as tissue with a 1 and non-tissue voxels with 0.* This is true since the within-tissue class weights get multiplied by the mask image value of one, whereas the weights for voxels outside our mask are multiplied by zero. Smoothing the binary image is then simply the sum of the weights that are within class, i.e. SW. While we made the above argument for voxels classified as tissue, the same argument can be made in general for any classification.

To summarize, the smoothing process for a given image, Im_{orig}, is: 1) create an image mask, Msk, with 1 values at voxel locations representing tissue, and 0 otherwise; 2) Create Im_{masked} by setting all non-tissue voxels of Im_{orig} to zero; 3) Apply the desired smoothing to Im_{masked} and Msk, creating Sm(Im_{masked}) and Sm (Msk); 4) Create the “Masked Smoothing” image by setting tissue voxels to the voxel-wise quotient Sm(Im_{masked})/Sm (Msk), and non-tissue voxels to their original values. By using a separable smoothing kernel in Step 3) the Masked Smoothing method will be orders of magnitude faster than directly using the 3D kernel for the calculation.

### Masked smoothing assertions

We developed Masked Smoothing as an alternative to basic Gaussian smoothing, which we term “Simple Smoothing”, and a method where the vessel voxels are set to zero prior to smoothing in an attempt to minimize the impact on tissue voxels, which we term “Removed Smoothing”. We believe that if Simple Smoothing, Removed Smoothing, and Masked Smoothing are used to process the same image, then tissue voxels that neighbor vessel voxels will best maintain their true value with Masked Smoothing. Furthermore, we believe that the differences between smoothing methods can lead to meaningful consequences in the determination of critical CT perfusion parameters. That is, if each smoothing method is applied to the individual time frames of a CT perfusion scan, then tissue voxels that are located near vessel voxels will have significant and meaningful differences in the resulting values of CBF, CBV and TTP depending on the smoothing method used.

## Methods

The Masked Smoothing method was tested against two smoothing methods (Simple and Removed Smoothing) that are similar, but which do not limit the smoothing to tissue voxels. The smoothing methods were tested using simulated data, phantom data, and anonymized CT perfusion data from patients. The simulations provide a framework for determining the noise reduction and bias for each method. The phantom data allows us to test for bias using real scanner data. The CT perfusion data from 23 patients allows an assessment of the impact of bias caused by smoothing the CT Perfusion time series images on the calculation of CBF, CBV, and TTP.

### Smoothing methods

- 1)
**Simple Smoothing**: The unmodified image volume is smoothed using a Gaussian kernel. - 2)
**Removed Smoothing**: The vascular voxels are set to zero, and the image is smoothed as in (1) above. - 3)
**Masked Smoothing**: First, Simple Smoothing is performed on the binary tissue mask. Second, the voxel-by-voxel ratio of the Removed Smoothing image and the smoothed tissue mask image (i.e. the result of the first step of this method) is returned as the Masked Smoothing image. Voxel values outside of the tissue mask are assigned the original image value, except for the phantom experiment. We used this case to also demonstrate that Masked Smoothing can be used to separately smooth both the object and object background.

### Simulation experiment

#### Parameters

- 1)
**Intensity ratio:**The intensity ratio is the ratio between the value assigned to vascular voxels and the tissue region. The tissue value was set to 50. The simulations used intensity ratios of: 2 to 1, 3 to 1, and 4 to 1, which correspond to vascular voxel values of 100, 150, and 200. - 2)
**SNR:**Gaussian noise was added to all simulation iterations. The standard deviation for the noise generator was set to 50, 25, and ~16.6, which corresponded to SNR values 1, 2, and 3 (lowest noise). - 3)
**Vessel width:**Is the cross-sectional width of the vascular region. Values used were: .8, 1.6, 2.4 mm. - 4)
**Smoothing kernel --FWHM**: The isotropic Gaussian smoothing kernel size was set to Full Width Half Max (FWHM) values of 1, 2, and 4 mm.

#### Iteration

- 1)
The simulated volume was formed with intensity values of tissue voxels set to 50. Vascular voxels were selected according to “Vessel Width”, and assigned an intensity value according to the variable “Intensity Ratio”.2) Gaussian noise was added at a level determined by the Signal to Noise Ratio (SNR), Figure 1, image a.3) All three smoothing methods were applied. All methods used the same Gaussian kernel with the kernel size determined from the variable “Kernel FWHM”, Figure 1, images b-d.4) The value at a tissue location located midway along, and directly next to, the vascular voxels was selected and the value with noise and smoothing applied was recorded. Figure 2 shows the intensity profile for the different smoothing methods for a vertical line passing through the images of Figure 1.

#### Assessment

500 iterations were used to determine the mean and standard deviation for a selected tissue voxel that neighbored the vessel voxels for different settings of parameters. Since in-tissue values in all cases were set to 50, the calculated mean value from 500 iterations even if a high level of noise is added is expected to be very close to 50 unless a bias is present. The four parameters (Intensity Ratio, SNR, Vessel Width, and Smoothing Kernel FWHM) were individually varied using the values given above for each. When one parameter was varied the other values were set to default values (Intensity Ratio = 2 to 1, SNR = 2, Vessel Width = 1.6 mm, and FWHM = 2 mm). One simulation run (500 iterations) was performed for each parameter set allowing for examination of deviation from the expected tissue value of 50 for each smoothing method.

#### CT phantom experiment

#### Influence of smoothing method on CBF, CBV, and TTP values

Twenty-three CT perfusion studies were selected from the Neurosurgery department’s stroke research database, at the University at Buffalo. Each dataset consisted of nineteen CT perfusion volumes from a Toshiba Aquilion ONE, 320 slice scanner (with voxel sizes of .42×.42×.5 mm) which were collected from patients presenting with symptoms of a stroke. Images were converted from Dicom to NifTI format, and corrected for motion using SPM8 (http://www.fil.ion.ucl.ac.uk/spm). Image volumes were “skull striped”, and vascular voxels were identified using in-house software written in Matlab. The middle cerebral artery was automatically identified and a center portion was segmented and used for the arterial input function. A similar procedure was used to select the sagittal sinus, and these values were used to ensure the proper scaling of the arterial input function. A parametric image of CBF values was calculated using the maximum slope method, while a CBV image was calculated using the integral of tracer activity divided by the integral of arterial activity.

Tissue voxels immediately adjacent to a selected artery were selected by performing a voxel-wise dilation of the voxels representing a selected artery followed by an intersection with the tissue masks resulting in the elimination of the vascular voxels. For each smoothing method we calculated CBF, CBV, and TTP values for the selected neighboring tissue voxels, using a Gaussian kernel size of 4 × 4 × 4 mm, FWHM. Mean and voxel-wise statistics were calculated to examine differences in CBF, CBV, and TTP due to the smoothing method used.

### Ethics

This project is approved by the University at Buffalo Health Sciences Institutional Review Board.

## Results

### Simulation experiment

Figure 2 displays line profiles, each corresponding to a vertical line across each image of Figure 1, to demonstrate the effects of the Simple, Removed, and Masked Smoothing methods.

#### Simulation experiment intensity ratio

**Simulation—intensity ratio: mean values and (standard deviation)**

Method\Intensity ratio | 1.5 to 1 | 2 to 1 | 4 to 1 |
---|---|---|---|

Noise – No smoothing | 49.7 (25.1) | 50.3 (25.0) | 49.0 (25.7) |

Simple smoothing | 58.2 (1.6) | 66.3 (1.6) | 99.3 (1.7) |

Removed smoothing | 33.6 (1.3) | 33.5 (1.3) | 33.5 (1.4) |

Masked smoothing | 50.0 (1.9) | 50.0 (1.9) | 50.0 (2.1) |

#### Simulation experiment -- SNR

**Simulation—SNR: mean values and (standard deviation)**

Method\SNR | 1 | 2 | 3 (least noise) |
---|---|---|---|

Noise – No smoothing | 53.6 (49.0) | 49.9 (27.1) | 50.8 (16.9) |

Simple smoothing | 66.4 (3.3) | 66.4 (1.6) | 60.4 (1.1) |

Removed smoothing | 33.6 (2.6) | 33.6 (1.3) | 33.6 (0.8) |

Masked smoothing | 50.0 (3.9) | 49.9 (1.9) | 50.0 (1.3) |

#### Simulation experiment: change of vessel diameter

**Simulation—vessel width: mean values and (standard deviation)**

Method\Vessel diameter | 1 | 2 | 3 mm |
---|---|---|---|

Noise – No smoothing | 51.2 (26.0) | 50.2 (24.5) | 48.5 (25.5) |

Simple smoothing | 61.3 (1.6) | 66.5 (1.7) | 68.6 (1.6) |

Removed smoothing | 38.8 (1.4) | 33.6 (1.4) | 31.6 (1.3) |

Masked smoothing | 50.0 (1.8) | 50.1 (2.0) | 50.1 (2.1) |

#### Simulation experiment: change of smoothing kernel FWHM

**Simulation—smoothing kernel FWHM: mean values and (standard deviation)**

Method\kernel FWHM | 1 | 2 | 4 mm |
---|---|---|---|

Noise – No smoothing | 51.7 (25.1) | 49.9 (25.9) | 50.0 (24.0) |

Simple smoothing | 64.6 (4.1) | 66.4 (1.6) | 60.4 (0.6) |

Removed smoothing | 35.7 (3.6) | 33.6 (1.3) | 39.6 (0.5) |

Masked smoothing | 50.2 (5.1) | 50.0 (2.0) | 50.0 (0.7) |

#### All reported biases

For all simulations the number of iterations was 500, and the standard deviation was relatively small compared to the size of the bias. Hence, all biases reported above were strongly significant (p < 0.0001).

#### Phantom data experiment

The line profiles passing through the small and large object show that each smoothing method is essentially identical for voxels away from the object (Figures 4 and 5). However, where the profiles cross through the object, the Simple Smoothing method has significantly lower valued voxels than the reference, i.e. mean across all image volumes. In contrast, for several voxels on either side of the object, the Simple Smoothing method has higher intensity values than the reference. The Masked Smoothing method has values close to the reference both for voxels located within the object, and outside of the object. We notice that there are a few voxels at either side of the object where the reference values lay between central values for the object and background. This is an indication of the limitation in the scanner resolution and reflects partial volume and an inherent smoothness of the raw data. Similar effects are seen for the line profile passing through the smaller object.

The mean and standard deviation for the 41 voxel line parallel and adjacent to the large object, for all 30 collected image volumes was 6.90 (13.08) HU. With Simple Smoothing, using 5 × 5 × 5 mm Gaussian kernel, the mean increased to 29.78 HU, but the standard deviation decreased to 1.76. With the corresponding Masked Smoothing applied the mean equaled 8.26 HU, i.e. much closer to the original. Further, the standard deviation equaled 1.96 HU, close to same value seen with Simple Smoothing.

#### Patient CT perfusion data – calculation of CBF, CBV, and TTP values

**Mean Parametric values for subjects’ CT perfusion scan data**

Method\Parameter | CBF (ml/(cc x min)) | CBV (ml/cc) | TTP (min) |
---|---|---|---|

Simple smooth | 1.48 | .112 | .52 |

Removed smooth | .35 | .029 | .54 |

Masked smooth | .67 | .048 | .54 |

Execution time of the Masked Smoothing method was 66 seconds for the 512×512×320×19 voxel CT perfusion image volume using a Gaussian Filter with size 2 × 2 × 2 mm FWHM, which required a 25 × 25 × 21 voxel kernel. Execution time using an 8 × 8 × 8 mm FWHM Gaussian filter, which required a kernel of 99×99×93 voxels, was 81 seconds. Execution time was measured using “tic” and “toc” Matlab functions, on a multi-user Dell PowerEdge R710 server with Dual 2.4 GHz processors, and 48 GB RAM.

## Discussion

Our simulation and phantom data show that the Simple (i.e. ordinary Gaussian smoothing) and Removed Smoothing introduce a significant bias to tissue voxels that neighbor vessels, whereas our Masked Smoothing method did not introduce a bias. Our experiment using patient data revealed that the bias of the Simple and Removed Smoothing methods had a large impact on the calculation of CBF and CBV. The Removed Smoothing method had the lowest values for CBF and CBV and were influenced by factoring in zero values in the place of neighboring vascular values. The Simple Smoothing method had increased CBF and CBV values for tissue voxels that neighbor vessels that were not physiologically reasonable. The Masked Smoothing method had physiologically reasonable values for CBF and CBV, between the extremes returned by the Removed and Simple smoothing methods (Table 5). Given that the Masked Smoothing method does not introduce bias (as opposed to the Removed and Simple Smoothing methods), is easy to implement, and executes fast enough to allow clinical use, we advocate its use over the other methods for the smoothing of CT perfusion images.

### Study design

We used simulated data to test the performance characteristics of the three smoothing methods in situations where the smoothing neighborhood for a tissue voxel included the much higher valued vessel voxels. In all simulations the tissue and vessel intensity values remained constant, allowing us to measure the effect of varying vessel characteristics (both vessel size and tracer concentration), the effect of varying SNR, and influence of the smoothing kernel FWHM. Using phantom imaging we were able to further show potential biases caused by the different smoothing methods. Using real world data from 23 patients, we also compared Simple, Masked, and Removed Smoothing to examine whether the theoretical improvement seen on simulations can have a real life impact in the calculation of CBF, CBV, and TTP. Using this approach we not only showed that Masked Smoothing did not have the bias of the other methods, but we also demonstrated the large practical impact this has on determining physiological parametric images for CBF and CBV.

### Change of intensity ratio/vessel width

Our simulation experiments indicate that the Simple Smoothing method has a large upward bias for tissue voxels surrounding a vessel that increases as the intensity of the vessel voxel increases. By setting the vessel voxels to zero for the purpose of smoothing, the Removed Smoothing method has a downward bias for tissue voxels neighboring a vessel voxel that is both fixed and independent of the vessel voxel’s intensity level. The Masked Smoothing method avoided bias by compensating for voxels set to zero. Increasing the vessel width increased the bias for the Simple Smoothing method, which reflects that a greater number of high intensity vessel voxels are within the smoothing neighborhood of the tissue voxel. The Removed Smoothing method also increased its bias (downward) with an increase in vessel size. This is reasonable, since for a given smoothing neighborhood the Removed Smoothing method would have a greater number of vascular voxels set to zero as the vessel width increases. Again, by compensating for voxels that were set to zero the Masked Smoothing method did not exhibit a bias.

### Change of SNR/smoothing kernel FWHM

All smoothing methods provided a large decrease in the noise level. Increasing the noise level caused an increase in the standard deviation measured for all methods, but had no effect on the calculated mean value. Increasing the filter kernel for all methods reduced the measured standard deviation. Increasing the filter size from 1 mm to 2 mm increased the bias for both Simple and Removed Smoothing. However, increasing the smoothing kernel further to 4 mm resulted in the smallest bias. The change in bias reflects the weighted proportion of voxels that are within the smoothing neighborhood. With the 4 mm smoothing kernel, the smoothing is incorporating a significant number of voxels from the “other-side” of the vessel, hence lessening the influence of the vessel itself. As in all cases, the Masked Smoothing exhibited no significant bias.

### Influence of smoothing method on the calculation of CBF, CBV, and TTP

Smoothing is a critical noise reduction pre-processing step prior to the calculation of physiologic parameters as we have demonstrated previously using simulations [19–21]. The CBF and CBV derived from the Simple and Removed Smoothing methods differed, approximately, by a factor of four for tissue voxels close to vessels, thus demonstrating the critical importance of the smoothing method. The CBF and CBV values, calculated using the Masked Smoothing method, were in-between and significantly different from the other smoothing methods, and closest to physiologically expected values. Since the Masked Smoothing approach showed no bias on the simulated data, we believe these CBF and CBV values are the most accurate. The TTP values for the Removed and Masked Smoothing were identical because the time activity curves for a given voxel will only differ by a scaling multiple, and were close to the Simple Smoothing method.

### Filter selection

We used a Gaussian smoothing kernel for our implementation because it is commonly used for medical images, and allows for fast implementations because it is separable. Our 3D execution times for an entire volume was significantly faster than a 2D TIPs bilateral filter on a single slice. Klotz and König [5] also applied smoothing separately on brain and vascular voxels. Their approach used multiple applications of a mean filter, whereas we utilized a Gaussian kernel. Our approach would also work with mean filters, since they are also separable. There are very fast methods for implementing mean filters; and furthermore, multiple passes of a mean filter can be used to approximate a Gaussian filter. However, internal timings during development favored our approach.

### Segmentation and segmentation

The Masked Smoothing method assumes that satisfactory segmentation is available. However, if the thickness between planes is high then partial volume effects may hinder segmentation. If a vessel voxel were to be classified as a tissue voxel, then neighboring tissue voxels will be biased upward, especially as the tracer concentration peaks in the vessel. However, this bias cannot exceed the bias from using Simple Smoothing. Because of the quantitation, some voxels partially represent both underlying tissue and vessel. This is not a problem in practice. If this voxel is excluded, the estimate for a nearby tissue voxel proceeds without using the value. If the voxel is included, then a neighboring voxel may be biased upward, but the effect will be minimal since the voxel partially represents tissue and thus will not reach especially high intensity levels. This is similar to the situation seen with the phantom data, where the mean of the raw data shows a gradual increase to the higher intensity object. In this case the Masked Smoothing best approximated the best estimate of the true value found by calculating the mean across 30 image volumes. Finally, our method allows both the arterial and tissue regions to be smoothed separately.

## Conclusion

We demonstrated that the Masked Smoothing method executes rapidly and can readily integrate into existing smoothing kernels. The Masked Smoothing method does not introduce a bias in situations where nearby voxels have a different classification and a large difference in intensity values. This accuracy, coupled with speed, gives the Masked Smoothing method the potential to significantly improve the clinical processing of perfusion imaging.

## Declarations

### Acknowledgements

The authors wish to thank Carmen Mieles, RTCT at WNY PET/CT, for her technical expertise and assistance for the phantom data collection.

## Authors’ Affiliations

## References

- Kudo K, Terae S, Katoh C, Oka M, Shiga T, Tamaki N, Miyasaka K: Quantitative cerebral blood flow measurement with dynamic perfusion CT using the vascular-pixel elimination method: comparison with H215O positron emission tomography. AJNR Am J Neuroradiol. 2003, 24 (3): 419-426.PubMedGoogle Scholar
- König M, Bültmann E, Bode-Schnurbus L, Koenen D, Mielke E, Heuser L: Image quality in CT perfusion imaging of the brain. Eur Radiol. 2007, 17 (1): 39-47. 10.1007/s00330-006-0277-3.View ArticlePubMedGoogle Scholar
- Kudo K, Sasaki M, Yamada K, Momoshima S, Utsunomiya H, Shirato H, Ogasawara K: Differences in CT Perfusion Maps Generated by Different Commercial Software: Quantitative Analysis by Using Identical Source Data of Acute Stroke Patients1. Radiology. 2010, 254 (1): 200-10.1148/radiol.254082000.View ArticlePubMedGoogle Scholar
- Sasaki M, Kudo K, Oikawa H: CT perfusion for acute stroke: current concepts on technical aspects and clinical applications. 2006, Elsevier, 30-36.Google Scholar
- Klotz E, Konig M: Perfusion measurements of the brain: using dynamic CT for the quantitative assessment of cerebral ischemia in acute stroke. Eur J Radiol. 1999, 30 (3): 170-184. 10.1016/S0720-048X(99)00009-1.View ArticlePubMedGoogle Scholar
- Saint-Marc P, Chen J, Medioni G: Adaptive smoothing: a general tool for early vision. IEEE Trans Pattern Anal Mach Intell. 1991, 13 (6): 618-624.View ArticleGoogle Scholar
- Nagao M, Matsuyama T: Edge preserving smoothing. Computer graphics and image processing. 1979, 9 (4): 394-407. 10.1016/0146-664X(79)90102-3.View ArticleGoogle Scholar
- Fang M, Qian J: Adaptive edge-preserving smoothing filter. Google Patents. 1998, US Patent: 08/672,194, Publication date: June 23Google Scholar
- Alvarez L, Guichard F, Lions PL, Morel JM: Axioms and fundamental equations of image processing. Archive for rational mechanics and analysis. 1993, 123 (3): 199-257. 10.1007/BF00375127.View ArticleGoogle Scholar
- Alvarez L, Lions PL, Morel JM: Image selective smoothing and edge detection by nonlinear diffusion. II SIAM Journal on numerical analysis. 1992, 29 (3): 845-866. 10.1137/0729052.View ArticleGoogle Scholar
- Angenent S, Pichon E, Tannenbaum A: Mathematical methods in medical image processing. Bulletin of the American Mathematical Society. 2006, 43 (3): 365-396. 10.1090/S0273-0979-06-01104-9.View ArticlePubMedPubMed CentralGoogle Scholar
- Chan TF, Shen J, Vese L: Variational PDE models in image processing. Notices AMS. 2003, 50: 14-26.Google Scholar
- Garcia D: Robust smoothing of gridded data in one and higher dimensions with missing values. Computational Statistics & Data Analysis. 2010, 54 (4): 1167-1178. 10.1016/j.csda.2009.09.020.View ArticleGoogle Scholar
- Paris S, Durand F: A fast approximation of the bilateral filter using a signal processing approach. International Journal of Computer Vision. 2009, 81 (1): 24-52. 10.1007/s11263-007-0110-8.View ArticleGoogle Scholar
- Paris S, Kornprobst P, Tumblin J: Bilateral filtering. Theory and applications, Foundations and Trends in Computer Graphics and Vision. 2009, 4 (1): 1-73.Google Scholar
- Mendrik A, Vonken E, Dankbaar JW, Prokop M, Van Ginneken B: Noise filtering in thin-slice 4D cerebral CT perfusion scans. 2010, SPIE Proceedings, 7623View ArticleGoogle Scholar
- Mendrik AM, Vonken E, van Ginneken B, de Jong HW, Riordan A, van Seeters T, Smit EJ, Viergever MA, Prokop M: TIPS bilateral noise reduction in 4D CT perfusion scans produces high-quality cerebral blood flow maps. Phys Med Biol. 2011, 56: 3857-10.1088/0031-9155/56/13/008.View ArticlePubMedGoogle Scholar
- Mendrik A, Vonken E, van Ginneken B, Smit E, Waaijer A, Bertolini G, Viergever MA, Prokop M: Automatic segmentation of intracranial arteries and veins in four-dimensional cerebral CT perfusion scans. Med Phys. 2010, 37: 2956-10.1118/1.3397813.View ArticlePubMedGoogle Scholar
- Fisher J: Improvements in Computed Tomography Perfusion Output using Complex Singular Value Decomposition and the Maximum Slope Algorithm. 2014, Master’s Thesis, Boston University, School of MedicineGoogle Scholar
- Wack DS, Badgaiyan RD: Complex singular value decomposition based noise reduction of dynamic PET images. Current Medical Imaging Reviews. 2011, 7 (2): 113-117. 10.2174/157340511795445685.View ArticleGoogle Scholar
- Snyder K, Seals K, Wack D: Using simulations to explore the characteristics of CT perfusion calculations in the assessment of stroke. Current Medical Imaging Reviews. 2014, 10 (3): In pressGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/14/28/prepub

### Pre-publication history

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.