Conv2 python

import numpy as np
from scipy.ndimage.filters import convolve

def conv2(x,y,mode='same'):
    """
    Emulate the function conv2 from Mathworks.

    Usage:

    z = conv2(x,y,mode='same')

    TODO: 
     - Support other modes than 'same' (see conv2.m)
    """

    if not(mode == 'same'):
        raise Exception("Mode not supported")

    # Add singleton dimensions
    if (len(x.shape) < len(y.shape)):
        dim = x.shape
        for i in range(len(x.shape),len(y.shape)):
            dim = (1,) + dim
        x = x.reshape(dim)
    elif (len(y.shape) < len(x.shape)):
        dim = y.shape
        for i in range(len(y.shape),len(x.shape)):
            dim = (1,) + dim
        y = y.reshape(dim)

    origin = ()

    # Apparently, the origin must be set in a special way to reproduce
    # the results of scipy.signal.convolve and Matlab
    for i in range(len(x.shape)):
        if ( (x.shape[i] - y.shape[i]) % 2 == 0 and
             x.shape[i] > 1 and
             y.shape[i] > 1):
            origin = origin + (-1,)
        else:
            origin = origin + (0,)

    z = convolve(x,y, mode='constant', origin=origin)

    return z

Source: stackoverflow.com

Similar codes

4 votes327 views1 answer
Plt add axis name
Python
0 vote349 views2 answers
From django.core.management import execute_from_command_line ImportError: No module named django.core.management
Python
0 vote319 views1 answer
How to make a label in python
Python
41 votes341 views5 answers
.replace python
Python
22 votes315 views5 answers
If else python
Python
1 vote314 views1 answer
Division euclidienne python
Python
0 vote353 views1 answer
Error 302 heroku django
Python Django
0 vote312 views1 answer
Beautifulsoup find get value
Python
0 vote319 views1 answer
Split and only grab first part of string
Python
0 vote315 views1 answer
How to open excel with more than one sheetpython
Python