int main(int, char *[])
{
ImageType::Pointer image = ImageType::New();
const ImageType::SizeType size = {{ 200, 200, 200}};
const ImageType::IndexType start = {{ 0, 0, 0 }};
ImageType::RegionType region;
region.SetSize( size );
region.SetIndex( start );
image->SetRegions( region );
image->Allocate(true);
const ImageType::IndexType pixelIndex = {{27,29,37}};
ImageType::PixelType pixelValue = image->GetPixel( pixelIndex );
image->SetPixel( pixelIndex, pixelValue+1 );
return 0;
}