SimpleITK  
Python/ImageCreateAndSet.py
1#!/usr/bin/env python
2# =========================================================================
3#
4# Copyright NumFOCUS
5#
6# Licensed under the Apache License, Version 2.0 (the "License");
7# you may not use this file except in compliance with the License.
8# You may obtain a copy of the License at
9#
10# http://www.apache.org/licenses/LICENSE-2.0.txt
11#
12# Unless required by applicable law or agreed to in writing, software
13# distributed under the License is distributed on an "AS IS" BASIS,
14# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
15# See the License for the specific language governing permissions and
16# limitations under the License.
17#
18# =========================================================================
19
20""" A SimpleITK example that sets pixel values and computes a Gaussian blob. """
21
22import os
23
24import SimpleITK as sitk
25
26xImg = sitk.Image(256, 256, sitk.sitkFloat32)
27yImg = sitk.Image(256, 256, sitk.sitkFloat32)
28
29for y in range(0, xImg.GetSize()[1]):
30 for x in range(0, xImg.GetSize()[0]):
31 xImg.SetPixel(x, y, x)
32 yImg[x, y] = y
33
34sigma = 50
35
36xImg = sitk.Subtract(xImg, xImg.GetSize()[0] / 2)
37yImg = yImg - yImg.GetSize()[1] / 2
38
39gaussianImg = sitk.Exp(-1 * (xImg**2 + yImg**2) / (2.0 * sigma**2))
40
41if "SITK_NOSHOW" not in os.environ:
42 sitk.Show(gaussianImg, "Gaussian Blob")
The Image class for SimpleITK.
Definition sitkImage.h:77
void SITKIO_EXPORT Show(const Image &image, const std::string &title="", const bool debugOn=ProcessObject::GetGlobalDefaultDebug())
Image Exp(Image &&image1)
Computes the exponential function of each pixel.
Image Subtract(Image &&image1, const Image &image2)
Pixel-wise subtraction of two images.