Segmentation fault in 64 bit systems

[hgrecco]

I have a working script in my 32 bit python that fails when running in a 64 bit python with segmentation fault. The function that fails is:

def _shift(image, template):
# Form the affine registration instance.
shift = register.Register(
model.Shift,
metric.Residual,
sampler.CubicConvolution
)

# Coerce the image data into RegisterData.
image = register.RegisterData(image)
template = register.RegisterData(template)

# Register.
params, warp, img, error = shift.register(
    image,
    template,
    verbose=False,
    )

return params, warp, img, error

In addition to this: should I use https://github.com/pyimreg/imreg or https://github.com/pyimreg/python-register

Thanks

[hgrecco]

I also get segfault when the sampler is Nearest, Bilinear or Spline

[hgrecco]

To my surprise, the problem was in np.gradient. Everything got fixed when I changed to float32. Maybe the images are too big (1024, 1344)?

diff --git a/imreg/metrics/metric.py b/imreg/metrics/metric.py
index b1fd29b..f141853 100644
--- a/imreg/metrics/metric.py
+++ b/imreg/metrics/metric.py
@@ -104,7 +104,7 @@ class Residual(Metric):
| p = number of parameters.
"""

•    grad = np.gradient(warpedImage)
  

•    grad = np.gradient(np.asarray(warpedImage, dtype=np.float32))
  
   dIx = grad[1].flatten()
   dIy = grad[0].flatten()