# # Unless required by applicable law or agreed to in writing, software # distributed under the License is distributed on an "AS IS" BASIS, # WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. # See the License for the specific language governing permissions and # limitations under the License. # # ==========================================================================*/ import InsightToolkit as itk import sys reader = itk.itkImageFileReaderF2_New() reader.SetFileName(sys.argv[1]) diffusion = itk.itkGradientAnisotropicDiffusionImageFilterF2F2_New() diffusion.SetInput(reader.GetOutput()) diffusion.SetTimeStep(0.0625) diffusion.SetConductanceParameter(9.0) diffusion.SetNumberOfIterations(5) gradient = itk.itkGradientMagnitudeImageFilterF2F2_New() gradient.SetInput(diffusion.GetOutput()) watershed = itk.itkWatershedImageFilterF2_New() watershed.SetInput(gradient.GetOutput()) watershed.SetThreshold(0.01) watershed.SetLevel(0.2) writer = itk.itkImageFileWriterUL2_New() writer.SetFileName(sys.argv[2])
import InsightToolkit as itk import sys reader = itk.itkImageFileReaderF2_New() reader.SetFileName( sys.argv[1] ) diffusion = itk.itkGradientAnisotropicDiffusionImageFilterF2F2_New() diffusion.SetInput(reader.GetOutput()) diffusion.SetTimeStep(0.0625) diffusion.SetConductanceParameter(9.0) diffusion.SetNumberOfIterations( 5 ); gradient = itk.itkGradientMagnitudeImageFilterF2F2_New() gradient.SetInput(diffusion.GetOutput()) watershed = itk.itkWatershedImageFilterF2_New() watershed.SetInput(gradient.GetOutput()) watershed.SetThreshold(0.01) watershed.SetLevel(0.2) # not UL writer available # will wait for LabelToRGBImageFIlter to be integrated in the toolkit # to add this test writer = itk.itkImageFileWriterUL2_New() writer.SetFileName( sys.argv[2] ) writer.SetInput( watershed.GetOutput() ) writer.Update()