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

分段线性变换

# -*- 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=30elif pixel<=150:ant=(200-30)/(150-30)*(pixel-30)+30else:ant=(255-200)/(255-150)*(pixel-150)+200return 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 npdef 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 cv2def 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_imgimg = 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)

原图：

对比图：

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！