ITK/Examples/Smoothing/BilateralImageFilter: Difference between revisions
From KitwarePublic
Jump to navigationJump to search
(Show difference.) |
No edit summary |
||
Line 1: | Line 1: | ||
<div class=" | <div class="floatcenter">[[File:ITK_Examples_Baseline_Smoothing_TestBilateralImageFilter.png]]</div> | ||
Blur an image while preserving edges. | Blur an image while preserving edges. | ||
Revision as of 19:42, 19 December 2010
Blur an image while preserving edges.
BilateralImageFilter.cxx
<source lang="cpp">
- include "itkImage.h"
- include "itkImageFileReader.h"
- include "itkBilateralImageFilter.h"
- include "itkSubtractImageFilter.h"
- include "QuickView.h"
int main(int argc, char * argv[]) {
// Verify command line arguments if( argc < 2 ) { std::cerr << "Usage: " << std::endl; std::cerr << argv[0] << " inputImageFile [domainSigma] [rangeSigma]" << std::endl; return EXIT_FAILURE; } double rangeSigma = 2.0; double domainSigma = 2.0; if (argc > 2) { domainSigma = atof(argv[2]); } if (argc > 3) { rangeSigma = atof(argv[3]); }
// Parse command line arguments std::string inputFilename = argv[1]; // Setup types typedef itk::Image< float, 2 > ImageType; typedef itk::ImageFileReader< ImageType > ReaderType; typedef itk::BilateralImageFilter< ImageType, ImageType > FilterType; typedef itk::SubtractImageFilter< ImageType > SubtractImageFilterType; // Create and setup a reader ReaderType::Pointer reader = ReaderType::New(); reader->SetFileName( inputFilename.c_str() ); // Create and setup a derivative filter FilterType::Pointer bilateralFilter = FilterType::New(); bilateralFilter->SetInput( reader->GetOutput() ); bilateralFilter->SetDomainSigma(domainSigma); bilateralFilter->SetRangeSigma(rangeSigma);
SubtractImageFilterType::Pointer diff = SubtractImageFilterType::New(); diff->SetInput1(reader->GetOutput()); diff->SetInput2(bilateralFilter->GetOutput());
QuickView viewer; viewer.AddImage( reader->GetOutput(),true, itksys::SystemTools::GetFilenameName(argv[1]));
std::stringstream desc; desc << "Bilateral\ndomainSigma = " << domainSigma << " rangeSigma = " << rangeSigma; viewer.AddImage( bilateralFilter->GetOutput(), true, desc.str());
std::stringstream desc2; desc2 << "Original - Bilateral"; viewer.AddImage( diff->GetOutput(), true, desc2.str());
viewer.Visualize(); return EXIT_SUCCESS;
} </source>
CMakeLists.txt
<source lang="cmake"> cmake_minimum_required(VERSION 2.6)
PROJECT(BilateralImageFilter)
include_directories(/home/doriad/src/ITK/Wrapping/WrapITK/ExternalProjects/ItkVtkGlue/src/)
FIND_PACKAGE(VTK REQUIRED) INCLUDE(${VTK_USE_FILE})
FIND_PACKAGE(ITK REQUIRED) INCLUDE(${ITK_USE_FILE})
ADD_EXECUTABLE(BilateralImageFilter BilateralImageFilter.cxx) TARGET_LINK_LIBRARIES(BilateralImageFilter vtkHybrid ITKBasicFilters ITKIO ITKCommon)
</source>