u010409517的博客

我们用ITK对图像进行处理的时候往往先会将16位的原始图像转换为0～255的8位无符号整数，但是要想输出为原始连续的序列，还要将这些掩码结果重新转换为16位的原始图像格式。这个过程用ITK实现起来比较复杂，我会一步步进行讲解。一、读取dicom图像将读取的signed short图像归一化到0～255的unsigned char图像ImageType::Pointer readdicom(s...

三维医学图像分别可从Axial、Coronal、Sagittal三个方向进行观察，常用的格式也有很多，最常见的就是dcm格式，代表一张切片。为了方便处理，通常在三维上直接对图像进行处理，常用的三维体数据格式有比如nii.gz、mhd、nrrd。但有时我们对三维体数据处理完后不可避免的要对每一张切片进行观察或者后处理，这样就要将一个三维体数据转换成二维的序列。下面我们用ITK将nii.gz体数据转...

均值滤波和和中值滤波都可以起到平滑图像，虑去噪声的功能。均值滤波采用线性的方法，平均整个窗口范围内的像素值，均值滤波本身存在着固有的缺陷，即它不能很好地保护图像细节，在图像去噪的同时也破坏了图像的细节部分，从而使图像变得模糊，不能很好地去除噪声点。均值滤波对高斯噪声表现较好，对椒盐噪声表现较差。中值滤波采用非线性的方法，它在平滑脉冲噪声方面非常有效,同时它可以保护图像尖锐的边缘，选择适当的点来...

由于ITK不支持可视化，处理完后必须先写入文件才能查看结果，我们可以结合VTK将处理结果可视化，快速查看图像处理结果。这里举一个简单的例子，读取dicom图像，并rescale到0～255，然后VTK显示结果。有两种方式显示结果： viewer.AddImage(image.GetPointer()); viewer.AddImage(rescaleFilter->GetOutput...

由于dicom通常都是一系列的图片，为方便后续处理，需要一次性将许多dicom图像一起读取写入。维数变成了三维，除了长宽，多了切片数。#include "itkGDCMImageIO.h"#include "itkGDCMSeriesFileNames.h"#include "itkImageSeriesReader.h"#include "itkImageSeriesWriter.h"...

ITK毕竟是医学图像处理专用，主要还是处理dicom医学格式的图像。dicom数据类型为signed short型。#include "itkImageFileReader.h"#include "itkImageFileWriter.h"#include "itkGDCMImageIO.h"#include <list>#include <fstream>...

RGB图像读取写入属于ITK基本操作，比较简单，直接套用下面代码即可，由于ITK不支持可视化，所以需要将处理完后的图片写入文件中打开查看，所以只需要改下面代码中的输入输出路径即可。#include "itkRGBPixel.h"#include "itkImage.h"#include "itkImageFileReader.h"#include "itkImageFileWriter.h...

原理介绍：Ostu大律法方法又名最大类间差方法，通过统计整个图像的直方图特性来实现全局阈值T的自动选取，其算法步骤为：先计算图像的直方图，将图像所有的像素点按照0~255共256个bin，统计落在每个bin的像素点数量；计算每个像素值的比例；i表示分类的最佳阈值，从0开始迭代；以i为阈值分类，计算第一类和第二类总的概率w0,w1，计算第一类的平均灰度和第二类的平均灰度以及整幅图像的平均...

ITK提供了完整的读取DICOM图像的方法，但不能实现可视化，VTK是一套进行数据可视化但通用开发工具包，所以在应用ITK进行图像读取和处理后，往往再利用VTK进行可视化。