ОШИБКА: градиентная норма ниже точности пикселей - PullRequest
0 голосов
/ 03 апреля 2020

Я использую elastix для регистрации медицинских изображений с преобразованием (Transform "MultiBSplineTransformWithNormal"), чтобы исправить скользящую регистрацию движения. В базе данных параметров elastix пример param0016 использовал это преобразование, которое должно обеспечить изображение метки для определения двух различных скользящих объектов в фиксированном изображении. Но файлы меток автора не могут быть получены.

Я использовал ITK-snap для сегментирования изображения метки 0/1 скользящего объекта на фиксированном изображении, все файлы изображений имеют формат .nii с примерно 0,6 * 0,6 * 3 мм. И команда, которую я использую: p param-bspline-multi.txt

Но я получил ошибку при запуске:

itk::ERROR: IsoContourDistanceImageFilter(000002004EBFB380): Gradient norm is lower than pixel precision

Кто-то сказал, что это потому, что вы должны использовать тип пикселя с плавающей точкой, а не короткий. Но то, что я использовал, - это просто плавание. Так в чем проблема? Я также удаляю изображение fixMask, чтобы увидеть, так ли это, потому что маска исправления и метка исправления конфликтуют, но эта ошибка все еще возникает. Кто-нибудь может дать совет? Большое спасибо!

полезная информация: fixImage изображение метки

Файл параметров, который я использовал

(FixedInternalImagePixelType "float")
(MovingInternalImagePixelType "float")

(FixedImageDimension 3)
(MovingImageDimension 3)

(UseDirectionCosines "true")

// **************** Main Components **************************

(Registration "MultiResolutionRegistration")
(Interpolator "BSplineInterpolator")
(ResampleInterpolator "FinalBSplineInterpolator")
(Resampler "DefaultResampler")

(FixedImagePyramid "FixedRecursiveImagePyramid")
(MovingImagePyramid "MovingRecursiveImagePyramid")

(Optimizer "AdaptiveStochasticGradientDescent")
(Transform "MultiBSplineTransformWithNormal")
(Metric "AdvancedMattesMutualInformation")

// ***************** Transformation **************************

(FinalGridSpacingInPhysicalUnits 7)
(HowToCombineTransforms "Compose")

// ******************* Similarity measure *********************

(NumberOfHistogramBins 64)
(ErodeMask "false")

// ******************** Multiresolution **********************

(NumberOfResolutions 4)

// ******************* Optimizer ****************************

(MaximumNumberOfIterations 800)

// **************** Image sampling **********************

(NumberOfSpatialSamples 2048)
(NewSamplesEveryIteration "true")
(ImageSampler "Random")

// ************* Interpolation and Resampling ****************

(BSplineInterpolationOrder 3)
(FinalBSplineInterpolationOrder 3)
(DefaultPixelValue 0)
(WriteResultImage "true")
(ResultImagePixelType "short")
(ResultImageFormat "nii")
...