The problem is not with shape_index any more, x, y and image have transposed dimensions. If you replace ax2.plot_surface(x, y, image, linewidth=0, alpha=0.5) with ax2.plot_surface(x.T, y.T, image, linewidth=0, alpha=0.5) it works. Cheers Emma On Wed, Jan 02, 2019 at 03:54:29PM +0100, Stefanie Lück wrote:
As Emma pointed out, an image and the code would be helpful to solve the problem. I can't sent the image because of confidentiality reasons but I have the problem with other images too...
Cheers Stefanie
Am Mi., 2. Jan. 2019 um 12:11 Uhr schrieb Stefanie Lück <luecks@gmail.com>:
Hi Emma,
thanks for you reply, I tried already to load the image as gray scale and after your suggestion to convert the image (should be the same) but I got the following error message:
RuntimeWarning: divide by zero encountered in true_divide return (2.0 / np.pi) * np.arctan((l2 + l1) / (l2 - l1)) D:\Anaconda\lib\site-packages\skimage\feature\corner.py:421: RuntimeWarning: invalid value encountered in true_divide return (2.0 / np.pi) * np.arctan((l2 + l1) / (l2 - l1)) C:/Users/stefanie/PycharmProjects/BluVision/Testing/shape_index.py:27: RuntimeWarning: invalid value encountered in less point_y, point_x = np.where(np.abs(s - target) < delta) C:/Users/stefanie/PycharmProjects/BluVision/Testing/shape_index.py:36: RuntimeWarning: invalid value encountered in less point_y_s, point_x_s = np.where(np.abs(s_smooth - target) < delta) Traceback (most recent call last): File "C:/Users/stefanie/PycharmProjects/BluVision/Testing/ shape_index.py", line 59, in <module> ax2.plot_surface(x, y, image, linewidth=0, alpha=0.5) File "D:\Anaconda\lib\site-packages\mpl_toolkits\mplot3d\axes3d.py", line 1612, in plot_surface X, Y, Z = np.broadcast_arrays(X, Y, Z) File "D:\Anaconda\lib\site-packages\numpy\lib\stride_tricks.py", line 252, in broadcast_arrays shape = _broadcast_shape(*args) File "D:\Anaconda\lib\site-packages\numpy\lib\stride_tricks.py", line 187, in _broadcast_shape b = np.broadcast(*args[:32]) ValueError: shape mismatch: objects cannot be broadcast to a single shape
I would love to open an issue or even to try a pull request. Since this would be the first time, I would need some help and might come back to you if it is OK for you. I notice quite often an uncertainty in the documentation about loading images, so I am glad to help to improve it.
Cheers Stefanie
Am Di., 1. Jan. 2019 um 22:50 Uhr schrieb Emmanuelle Gouillart < emmanuelle.gouillart@normalesup.org>:
Hi Stefanie,
happy new year to you too :-). I had a quick look at the code of feature.shape_index, and the function requires a 2D (hence grayscale image). Can you check that the function works on color.rgb2gray(image)?
This also means that the docstring is incorrect, we should say
image: 2D array
instead of
image : ndarray
Would you mind opening an issue for this? (or even submitting a pull request?)
All the best Emma
On Tue, Jan 01, 2019 at 08:48:52PM +0100, Stefanie Lück wrote: > Hi everybody,
> first of all I wish you a happy new year :)
> I have a problem with the Shape Index example. When I try with my image (shape > (778, 1802, 3)) I am getting the following error message:
> Traceback (most recent call last): > File "C:/Users/stefanie/PycharmProjects/BluVision/Testing/ shape_index.py", > line 16, in <module> > s = shape_index(image) > File "D:\Anaconda\lib\site-packages\skimage\feature\corner.py", line 419, in > shape_index > l1, l2 = hessian_matrix_eigvals(H) > ValueError: too many values to unpack (expected 2)
> System: > scikit-image 0.14.1 > Windows 10 > Python 3.6
> Any suggestions? > Thanks > Stefanie
> _______________________________________________ > scikit-image mailing list -- scikit-image@python.org > To unsubscribe send an email to scikit-image-leave@python.org _______________________________________________ scikit-image mailing list -- scikit-image@python.org To unsubscribe send an email to scikit-image-leave@python.org
_______________________________________________ scikit-image mailing list -- scikit-image@python.org To unsubscribe send an email to scikit-image-leave@python.org