拍個自拍,讓Python告訴你,軍訓過后你黑了幾度?
大數據文摘出品
作者:曹培信、寧靜
一年一度的大學開學季,一年一度的軍訓季。
在六月中旬高考結束之后,萬千學子迎來了他們人生中最長也是最無憂無慮的假期,到了八月底九月初,他們將踏上人生的一段重要旅程——大學。
然而等待他們的第一課,便是軍訓。
每所高校對軍訓的要求不同,從時間來看,短的只有5天,長的多達一個月,根據網上一份軍訓時間排行榜,比如清華大學,以34天穩居第二。(警校排第一也是無可厚非了,不過吉林的院校確實熱衷軍訓,前十中兩所吉林的院校上榜)
除了一些開始又晚、時間又長的院校(比如重慶大學),大部分院校的軍訓應該已經結束了,然而經歷了軍姿、齊步、正步、閱兵甚至拉練的“摧殘”后,軍訓歲月在身上留下的最深的痕跡便是——曬!黑!了!
圖片來自網絡
軍訓前和軍訓后,就是“白古”和“黑古”的差別啊!想知道自己軍訓后究竟曬黑了多少么?下面文摘菌就帶你用Python看看,自己究竟軍訓后黑了幾個度。
基于RGB和YCbCr顏色空間的混合膚色檢測
想知道自己的皮膚顏色,首先要將皮膚檢測出來。
膚色檢測有很多方法,比如:
- 基于RGB的顏色空間模型;
- 基于橢圓皮膚模型的皮膚檢測;
- YCrCb顏色空間Cr分量+Otsu法閾值分割;
- 基于YCrCb顏色空間Cr,Cb范圍篩選法;
- HSV顏色空間H范圍篩選法;
- opencv自帶膚色檢測類AdaptiveSkinDetector;
相關鏈接:https://blog.csdn.net/qq_22527639/article/details/81501565
2004年,Georgy Kukharev和Adam Nowosielski為了提高模型的穩定性,將多個顏色空間結合,提出RGB顏色空間和YCbCr顏色空間的混合膚色檢測器。像素值滿足如下條件:
實現的代碼也很簡單,首先引入必要的包:
- import cv2
- import numpy as np
- from matplotlib import pyplot as plt
然而操縱圖像,將RGB顏色空間3通道的值和YCbCr顏色空間3通道的值結合起來,然后根據判別條件進行膚色檢測:
- def skin_color(imgFile):
- # load an original image
- img = cv2.imread(imgFile)
- rows,cols,channels = img.shape
- # convert color space from rgb to ycbcr
- imgYcc = cv2.cvtColor(img, cv2.COLOR_BGR2YCR_CB)
- # convert color space from bgr to rgb
- img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB)
- # prepare an empty image space
- imgSkin = np.zeros(img.shape, np.uint8)
- # copy original image
- imgimgSkin = img.copy()
- s=0
- sum_R=sum_G=sum_B=0
- for r in range(rows):
- for c in range(cols):
- # non-skin area if skin equals 0, skin area otherwise
- skin = 0
- # get values from rgb color space
- R = img.item(r,c,0)
- G = img.item(r,c,1)
- B = img.item(r,c,2)
- # get values from ycbcr color space
- Y = imgYcc.item(r,c,0)
- Cr = imgYcc.item(r,c,1)
- Cb = imgYcc.item(r,c,2)
- # skin color detection
- if R > G and R > B:
- if (G >= B and 5 * R - 12 * G + 7 * B >= 0) or (G < B and 5 * R + 7 * G - 12 * B >= 0):
- if Cr > 135 and Cr < 180 and Cb > 85 and Cb < 135 and Y > 80:
- # print 'Skin detected!'
參考鏈接:https://blog.csdn.net/shadow_guo/article/details/43635181
對比色卡,看看你黑了幾度
檢測出了皮膚的區域,我們就需要將皮膚區域的顏色RGB值計算出來,然后與色卡進行對比。
說到色卡,不得不提到Pantone(潘通)——一家專注于研究顏色的公司,以其Pantone顏色匹配系統(PMS)而聞名,該系統被廣泛用于各種行業,特別是平面設計,時裝設計,產品設計,印刷和制造,并支持從設計到生產的色彩管理。
2013年,潘通發布了一款彩通膚色指南(PANTONE SkinTone Guide),這個指南根據科學測量各種人類皮膚類型中數千種實際膚色而建立,為再現實體膚色而配制,被稱為人類膚色的完整視覺參考。
110種人類膚色被從1Y01 SP至4R15 SP編號,適用于各個人種。
文摘菌從中選取了比符合中國人膚色的1Y01-1Y13作為對比色卡。
從左至右依次為0-12度
然后將膚色顏色RGB與色卡的RGB數值進行對比,因為人眼對于R、G、B的敏感程度不同,在轉換的時候需要給予不同的權重。
人類視覺對綠色最敏感,因此它具有最大的系數值(0.7152),對藍色最不敏感,因此具有最小的系數(0.0722)。
- def Compare(list,color):
- min=100
- count=len(list)
- for i in range(count):
- value=list[i]
- error=abs(color[0]-value[0])*0.2126+abs(color[1]-value[1])*0.7152+abs(color[2]-value[2])*0.0722
- if(error<min):
- min=error
- iindex=i
- return index+1
最后我們輸入“白古”和“黑古”的圖片,與比色卡的RGB數值進行對比,輸出結果:
- if __name__ == '__main__':
- img_before= 'images/before1.jpg'
- img_after='images/after.jpg'
- color_before=skin_color(img_before)
- color_after=skin_color(img_after)
- print("before:the extracted RGB value of the skin color is {0}".format(color_before))
- print("after:the extracted RGB value of the color is {0}".format(color_after))
- list=[(200, 172, 153), (200, 169, 149), (197, 166, 145), (194, 163, 142), (190, 157, 134), (187, 152, 129), (182, 146, 121), (177, 136, 108), (168, 127, 100), (160, 118, 90), (148, 108, 81), (135, 98, 73), (119, 87, 65)]
- #color_dir="skin_color"
- #skincolor.load_color(color_dir,list)
- #print(list)
- #print(list)
- degree_before=Compare(list,color_before)
- degree_after=Compare(list,color_after)
- print("before:the degree of the skin color is {0}".format(degree_before))
- print("after:the degree of the skin color is {0}".format(degree_after))
輸出的結果如下圖所示:
也就是說,“黑古”比“白古”整整黑了5個度,雖然文摘菌認為結果可能還有偏差,因為12是比色表的最高值,“黑古”很可能已經爆表了。
軍訓結束,寫代碼的日子開始了
也許當你拍下自己軍訓后的照片,用文摘菌的代碼和自己一個月前的自拍對比了一下,然后留下了傷心的眼淚。
文摘菌想告訴你:你不是一個人在戰斗!
今年參加軍訓的大學新生有數百萬,短則5天,長則一個月的軍訓,風吹日曬,大家的皮膚變黑在所難免,然而,就算是1個月的軍訓,也并不能真正掌握多少軍事技能,其實重要的還是原本初次見面還陌生的同班同學,通過軍訓開始熟絡起來;原本還陌生的校園,也齊步走了幾個來回。
可以說,軍訓是大學的預科班,結束了軍訓,就正式開始了大學生活。
文摘菌本著行業內人士的身份,自然要對那些報讀計算機相關專業的新生們說一句:軍訓結束,你們寫代碼的日子就要開始了!
也許在某個深夜,當你停下敲擊鍵盤的雙手,撫摸著額前稀疏的頭發,你可能會意識到,原來軍訓并不是大學最痛苦的時候。
【本文是51CTO專欄機構大數據文摘的原創譯文,微信公眾號“大數據文摘( id: BigDataDigest)”】