机器视觉实验灰度图像分段、指数、幂律变换

分段线性变换

# -*- coding: utf-8 -*-
"""
Created on Wed Oct 14 20:02:40 2020

@author: chenpeng
"""

import cv2
import math
import numpy as np
def XX(pixel):
    if pixel<=30:
        ant=30
    elif pixel<=150:
        ant=(200-30)/(150-30)*(pixel-30)+30
    else:
        ant=(255-200)/(255-150)*(pixel-150)+200
    return ant
def logTransform(img):

    #3通道RGB
    '''h,w,d = img.shape[0],img.shape[1],img.shape[2]
    new_img = np.zeros((h,w,d))
    for i in range(h):
        for j in range(w):
            for k in range(d):
                new_img[i,j,k] = c*(math.log(1.0+img[i,j,k]))'''

    #灰度图专属
    h,w = img.shape[0], img.shape[1]
    new_img = np.zeros((h, w))
    for i in range(h):
        for j in range(w):
            new_img[i, j] = XX(img[i, j])
            
    new_img = cv2.normalize(new_img,new_img,0,255,cv2.NORM_MINMAX)
    return new_img
    
#替换为你的图片路径
img = cv2.imread('cc.jpg',0)
cv2.imshow("img",img)
log_img = logTransform(img)
log_img=log_img.astype(np.uint8)#类型转换
cv2.imshow('xian_xing_img',log_img)
#保存图片
#cv2.imwrite('xian_xing_img',log_img)
#cv2.waitKey(0)

对数变换

"""
Created on Wed Oct 14 20:19:31 2020

@author: chenpeng
"""
import cv2
import math
import numpy as np

def logTransform(c,img):

    #3通道RGB
    '''h,w,d = img.shape[0],img.shape[1],img.shape[2]
    new_img = np.zeros((h,w,d))
    for i in range(h):
        for j in range(w):
            for k in range(d):
                new_img[i,j,k] = c*(math.log(1.0+img[i,j,k]))'''

    #灰度图专属
    h,w = img.shape[0], img.shape[1]
    new_img = np.zeros((h, w))
    for i in range(h):
        for j in range(w):
            new_img[i, j] = c * (math.log(1.0 + img[i, j]))
    new_img = cv2.normalize(new_img,new_img,0,255,cv2.NORM_MINMAX)
    return new_img
    
#替换为你的图片路径
img = cv2.imread('cc.jpg',0)
cv2.imshow("img",img)
log_img = logTransform(41.0,img)
log_img=log_img.astype(np.uint8)#类型转换
cv2.imshow('log_img',log_img)
#保存图片
#cv2.imwrite('log_img.jpg',log_img)
#cv2.waitKey(0)

幂律变换

# -*- coding: utf-8 -*-
"""
Created on Wed Oct 14 20:19:31 2020

@author: chenpeng
"""

import math
import numpy as np
import cv2

def gammaTranform(c,gamma,img):
    h,w = img.shape[0],img.shape[1]
    new_img = np.zeros((h,w))
    for i in range(h):
        for j in range(w):
            new_img[i,j] = c*math.pow(img[i, j], gamma)
    new_img = cv2.normalize(new_img,new_img,0,255,cv2.NORM_MINMAX)
    return new_img

img = cv2.imread('cc.jpg',0)
cv2.imshow('img',img)
log_img = gammaTranform(1.5,2.5,img)
log_img=log_img.astype(np.uint8)#类型转换
cv2.imshow('milv_img',log_img)
#cv2.imwrite(r'C:\Users\xxx\Desktop\gray_2.5.jpg',new_img)
#cv2.waitKey(0)

原图：

对比图：