/// Put the origin on a voxel to avoid small skips
binary_image->SetOrigin(floor((bounds[0]-samp_origin[0])/spacing[0]-2)*spacing[0]+samp_origin[0],
floor((bounds[2]-samp_origin[1])/spacing[1]-2)*spacing[1]+samp_origin[1],
/// Put the origin on a voxel to avoid small skips
binary_image->SetOrigin(floor((bounds[0]-samp_origin[0])/spacing[0]-2)*spacing[0]+samp_origin[0],
floor((bounds[2]-samp_origin[1])/spacing[1]-2)*spacing[1]+samp_origin[1],