首先讀取圖片資料，轉化成矩陣格式：

from PIL import Image
import numpy as np

def getimage(path = path):
    #path為需要讀取圖片的路徑
    img = Image.open(path)
    M,N = img.size
    r,g,b = img.split()
    rd = np.asarray(r)
    gd = np.asarray(g)
    bd = np.asarray(b)
    return rd,gd,bd

分別計算一二三階顏色矩：

#r通道的一階顏色矩
rd_1 = rd.mean()

#r通道的二階顏色矩
rd_2 = rd.std()

#r通道的三階顏色矩
#定義一個求三階顏色矩的函式
def var(x=None):
    mid = np.mean(((x - x.mean()) ** 3))
    return np.sign(mid) * abs(mid) ** (1/3)
 

批量求一個資料夾下所有圖片的各階顏色矩：

import os

def getimagedata(path=path):
    filename = os.listdir(path)
    n = len(filename)
    data = np.zeros([n,9])
    for i in range(n):
            img = Image.open(path+'\\'+imagenames[i])
            M,N = img.size
            r,g,b = img.split()
            rd = np.asarray(r)
            gd = np.asarray(g)
            bd = np.asarray(b)
            data[i,0] = rd.mean();data[i,1] = gd.mean();data[i,2] = bd.mean()
            data[i,3] = rd.std();data[i,4] = gd.std();data[i,5] = bd.std()
            data[i,6] = self.var(rd);data[i,7] = self.var(gd);data[i,8] = self.var(bd)
        return data