标签：index 同态滤波 python I3 shape range 图像处理 np chanel

基于网友的方案，手动改的方案。传送门链接

'''
同态滤波
'''


def homofilter(img):
    I = np.float32(img)
    I3 = I
    m, n, chanels = I.shape
    rL = 0.5
    rH = 2
    c = 2
    d0 = 20
    for chanel in range(chanels):
        I1 = np.log(I[:, :, chanel] + 1)
        FI = np.fft.fft2(I1)
        n1 = np.floor(m / 2)
        n2 = np.floor(n / 2)
        D = np.zeros((m, n))
        H = np.zeros((m, n))
        for i in range(m):
            for j in range(n):
                D[i, j] = ((i - n1) ** 2 + (j - n2) ** 2)
                H[i, j] = (rH - rL) * (np.exp(c * (-D[i, j] ** 2 / (d0 ** 2)))) + rL
        I2 = np.fft.ifft2(H * FI)
        I3[:, :, chanel] = np.real(np.exp(I2))
    max_index = np.unravel_index(I3.argmax(), I3.shape)
    maxV = I3[max_index[0], max_index[1], max_index[2]]
    min_index = np.unravel_index(I3.argmin(), I3.shape)
    minV = I3[min_index[0], min_index[1], min_index[2]]
    print(I3.shape)
    for chanel in range(chanels):
        for i in range(m):
            for j in range(n):
                I3[i, j, chanel] = 255 * (I3[i, j, chanel] - minV) / (maxV - minV)
    return np.int32(I3)

##实现效果

 同态滤波之后，颜色有点变得离谱。大神路过有建议的话，麻烦评论一下，我会关注的5555

可能有些问题，大家积极评论，我积极改进。更新正文内容。

标签：index,同态滤波,python,I3,shape,range,图像处理,np,chanel
来源： https://blog.csdn.net/QuickservedNoodle/article/details/120420255

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。