# Local maximum
    local_maximas = []
    max_h = 2 * max_dimension
    max_w = theta_dimension

    for i in range(0, max_h-1):
        for j in range(0, max_w-1):
            local_maxima = H[i][j]
            local_maxima_pair = (i, j)
            if(j<max_w and H[i][j+1] > local_maxima): #right
                local_maxima = H[i][j+1]
                local_maxima_pair= (i,j+1)
            if(j>0 and H[i][j-1] > local_maxima): #left
                local_maxima = H[i][j-1]
                local_maxima_pair = (i,j-1)
            if(i<max_h and H[i+1][j] > local_maxima): #up
                local_maxima = H[i+1][j]
                local_maxima_pair = (i+1,j)
            if(i>0 and H[i-1][j] > local_maxima):   #down
                local_maxima = H[i-1][j]
                local_maxima_pair = (i-1,j)
            if(i<max_h and j<max_w and H[i+1][j+1] > local_maxima): #S-V
                local_maxima = H[i+1][j+1]
                local_maxima_pair = (i+1,j+1)
            if(i<max_h and j<0 and H[i+1][j-1] > local_maxima): #S-E
                local_maxima = H[i+1][j-1]
                local_maxima_pair = (i+1,j-1)
            if(i<0 and j<max_w and H[i-1][j+1] > local_maxima): #N-V
                local_maxima = H[i-1][j+1]
                local_maxima_pair = (i-1,j+1)
            if(i<0 and j<0 and H[i-1][j-1] > local_maxima): #N-E
                local_maxima = H[i-1][j-1]
                local_maxima_pair = (i-1,j-1)
            local_maximas.append(local_maxima_pair)