Merge pull request #215 from thewtex/resample-a-scalar-image-python

Author: thewtex

.gitignore

@@ -27,3 +27,6 @@ CMakeLists.txt.user*
 *.orig
 
 *.DS_Store
+
+# Jupyter notebooks
+.ipynb_checkpoints/

src/Filtering/ImageGrid/CMakeLists.txt

@@ -44,9 +44,9 @@ compare_to_baseline(
   BASELINE_PREFIX BrainProtonDensitySlice
 )
 
-add_example( UpsampleOrDownsampleScalarImage )
+add_example( ResampleAScalarImage )
 compare_to_baseline(
-  EXAMPLE_NAME UpsampleOrDownsampleScalarImage
+  EXAMPLE_NAME ResampleAScalarImage
   BASELINE_PREFIX OutputBaseline
 )
 
@@ -105,4 +105,4 @@ compare_to_baseline(EXAMPLE_NAME RunImageFilterOnRegionOfImage
 add_example(CropImageBySpecifyingRegion2)
 compare_to_baseline(EXAMPLE_NAME CropImageBySpecifyingRegion2
         BASELINE_PREFIX CropImageBySpecifyingRegion2
-        )
+        )

src/Filtering/ImageGrid/ResampleAScalarImage/CMakeLists.txt

@@ -0,0 +1,38 @@
+cmake_minimum_required(VERSION 3.10.2)
+
+project(ResampleAScalarImage )
+
+find_package(ITK REQUIRED)
+include(${ITK_USE_FILE})
+
+add_executable(ResampleAScalarImage Code.cxx)
+target_link_libraries(ResampleAScalarImage ${ITK_LIBRARIES})
+
+install(TARGETS ResampleAScalarImage
+  DESTINATION bin/ITKExamples/Filtering/ImageGrid
+  COMPONENT Runtime
+  )
+
+install(FILES Code.cxx CMakeLists.txt Code.py
+  DESTINATION share/ITKExamples/Code/Filtering/ImageGrid/ResampleAScalarImage
+  COMPONENT Code
+  )
+
+enable_testing()
+add_test(NAME ResampleAScalarImageTest
+  COMMAND ${CMAKE_RUNTIME_OUTPUT_DIRECTORY}/ResampleAScalarImage
+    ${CMAKE_CURRENT_BINARY_DIR}/Gourds.png
+    Output.png
+    160
+    120
+  )
+
+if(ITK_WRAP_PYTHON)
+  add_test(NAME ResampleAScalarImageTestPython
+    COMMAND ${PYTHON_EXECUTABLE} ${CMAKE_CURRENT_SOURCE_DIR}/Code.py
+      ${CMAKE_CURRENT_BINARY_DIR}/Gourds.png
+      OutputPython.png
+      160
+      120
+    )
+endif()

src/Filtering/ImageGrid/UpsampleOrDownsampleScalarImage/Code.cxx renamed to src/Filtering/ImageGrid/ResampleAScalarImage/Code.cxx

@@ -0,0 +1,45 @@
+#!/usr/bin/env python
+
+# Copyright NumFOCUS
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#        http://www.apache.org/licenses/LICENSE-2.0.txt
+#
+# 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 sys
+import itk
+
+if len(sys.argv) != 5:
+    print("Usage: " + sys.argv[0] + " <input_filename> <output_filename> <size_x> <size_y>")
+    sys.exit(1)
+
+input_filename = sys.argv[1]
+output_filename = sys.argv[2]
+output_size = [int(sys.argv[3]), int(sys.argv[4])]
+
+input_image = itk.imread(input_filename)
+
+input_spacing = itk.spacing(input_image)
+input_size = itk.size(input_image)
+dimension = input_image.GetImageDimension()
+output_spacing = [input_spacing[dim] * input_size[dim] / output_size[dim] for
+        dim in range(dimension)]
+
+transform = itk.IdentityTransform[itk.D, dimension].New()
+
+output_image = itk.resample_image_filter(input_image,
+        size=output_size,
+        output_spacing=output_spacing,
+        output_origin=itk.origin(input_image),
+        transform=transform
+        )
+
+itk.imwrite(output_image, output_filename)

src/Filtering/ImageGrid/ResampleAScalarImage/Code.py

@@ -1,13 +1,13 @@
-Resample Scalar Image
-=====================
+Resample a Scalar Image
+=======================
 
 .. index::
    single: ResampleImageFilter
 
 Synopsis
 --------
 
-Resample one scalar image
+Resample a scalar image.
 
 
 Results
@@ -35,6 +35,12 @@ C++
 .. literalinclude:: Code.cxx
    :lines: 18-
 
+Python
+......
+
+.. literalinclude:: Code.py
+   :language: python
+   :lines: 1, 16-
 
 Classes demonstrated
 --------------------

src/Filtering/ImageGrid/UpsampleOrDownsampleScalarImage/Documentation.rst renamed to src/Filtering/ImageGrid/ResampleAScalarImage/Documentation.rst

@@ -24,5 +24,5 @@ ImageGrid
   Stack2DImagesInto3DImage/Documentation.rst
   TileImagesSideBySide/Documentation.rst
   UpsampleAnImage/Documentation.rst
-  UpsampleOrDownsampleScalarImage/Documentation.rst
+  ResampleAScalarImage/Documentation.rst
   WarpAnImageUsingADeformationField/Documentation.rst