Python实现mmseg分词算法和吐嘈

http://blog.csdn.net/acceptedxukai/article/details/7390300#comments

前两天一直在写爬虫，也是一直有心写个简单的搜索引擎吧，写出来肯定是没有Web界面的，本人不知道怎么地，对Web的代码一直有抵触心理。

搜索引擎嘛，我想要写出来必须要懂五个部分：

1、基础的多线程爬虫

2、分词算法，对用户输入的查找文本进行切割

3、MapReduce，用来统计词语出现的次数，这个直接关系到URL在队列中的等级

4、PageRank，用来给网页排名

5、Web界面

上述的五个部分，打算全部用自己的代码实现，不用任何第三方类库，用Python语言实现

爬虫大家都会写，这里就不说了，待会把自己花了1个点编码+调试的优先队列放出来吧，好久没写算法了，写个最大顶堆都花了老久。

下面继续说分词吧，咱是白手学习，啥也不会，就去Google了一下，先学了个正向最大匹配算法

啥正向最大匹配啊，说白了，就是拿着你的字典然后对文本进行暴力切割，这种算法切割出来的错误率肯定超级高，算法原理满大街都是，就不说了，时间复杂度倒挺低O(n)。

接着就换算法，学习了个也比较简单的mmseg算法，这是台湾同胞提出来的，其实就是对正向最大匹配算法进行了优化，这里的优化不是指优化时间复杂度，是优化降低错误率的，但是也会出现很奇葩的错误，让人看了都蛋疼。

算法的原文：http://technology.chtsai.org/mmseg/，中文有翻译的。

其实就是用了四个判断的法则：

这些我也就不去复制粘贴了，算法的说明也是满大街都是

简单给个博客：http://www.byywee.com/page/M0/S602/602088.html

斯坦福大学有个具体的实现包：http://nlp.stanford.edu/software/segmenter.shtml

这里举个正向最大匹配和mmseg算法对同样的文本切割后的效果吧：

1、南京市长江大桥欢迎您

正向最大匹配：南京/市长/江/大桥/欢迎您

mmseg：南京市/长江大桥/欢迎您

2、研究生命科学

正向最大匹配：研究生/命/科学

mmseg：研究/生命/科学

3、长春市长春药店

正向最大匹配：长春/市长/春药/店

mmseg：长春市/长春/药店

4、台湾国中学生

正向最大匹配：台湾国/中学生

mmseg：台湾/国中/学生

当然这些匹配的结果都取决于您的字典

大家在看了上面的博客，理解了什么是mmseg算法再往下看

具体Python代码如下：PyMmseg.py

此代码参照google code的一个代码，经过了我自己的修改和自己去其的理解

代码结构说明：

Word类就相当于C语言中的结构体，用来存单词和词频的

Chunk类是用来实现具体的切割判断方法的前期预处理计算的

ComplexCompare类是用来具体实现mmseg算法的四种评估方法的

有几个全局变量和全局函数是用来加载字典的，用全局是为了不让字典多次的加载

Analysis类是用来具体实现切割算法的

代码里的注释已经比较清楚了，在理解正向最大匹配算法的基础上再理解mmseg算法的原理就应该很容易明白下面的源码，具体不清楚的可以留言。

[python] view plaincopyprint?

1. class Word:  
2.     def __init__(self,text = '',freq = 0):  
3.         self.text = text  
4.         self.freq = freq  
5.         self.length = len(text)  
6.   
7. class Chunk:  
8.     def __init__(self,w1,w2 = None,w3 = None):  
9.         self.words = []  
10.         self.words.append(w1)  
11.         if w2:  
12.             self.words.append(w2)  
13.         if w3:  
14.             self.words.append(w3)  
15.       
16.     #计算chunk的总长度  
17.     def totalWordLength(self):  
18.         length = 0  
19.         for word in self.words:  
20.             length += len(word.text)  
21.         return length  
22.       
23.     #计算平均长度  
24.     def averageWordLength(self):  
25.         return float(self.totalWordLength()) / float(len(self.words))  
26.       
27.     #计算标准差  
28.     def standardDeviation(self):  
29.         average = self.averageWordLength()  
30.         sum = 0.0  
31.         for word in self.words:  
32.             tmp = (len(word.text) - average)  
33.             sum += float(tmp) * float(tmp)  
34.         return sum  
35.       
36.     #自由语素度  
37.     def wordFrequency(self):  
38.         sum = 0  
39.         for word in self.words:  
40.             sum += word.freq  
41.         return sum  
42.   
43. class ComplexCompare:  
44.       
45.     def takeHightest(self,chunks,comparator):  
46.         i = 1  
47.         for j in range(1, len(chunks)):  
48.             rlt = comparator(chunks[j], chunks[0])  
49.             if rlt > 0:  
50.                 i = 0  
51.             if rlt >= 0:  
52.                 chunks[i], chunks[j] = chunks[j], chunks[i]  
53.                 i += 1  
54.         return chunks[0:i]  
55.       
56.     #以下四个函数是mmseg算法的四种过滤原则，核心算法  
57.     def mmFilter(self,chunks):  
58.         def comparator(a,b):  
59.             return a.totalWordLength() - b.totalWordLength()  
60.         return self.takeHightest(chunks,comparator)  
61.       
62.     def lawlFilter(self,chunks):  
63.         def comparator(a,b):  
64.             return a.averageWordLength() - b.averageWordLength()  
65.         return self.takeHightest(chunks,comparator)  
66.       
67.     def svmlFilter(self,chunks):  
68.         def comparator(a,b):  
69.             return b.standardDeviation() - a.standardDeviation()  
70.         return self.takeHightest(chunks, comparator)  
71.       
72.     def logFreqFilter(self,chunks):  
73.         def comparator(a,b):  
74.             return a.wordFrequency() - b.wordFrequency()  
75.         return self.takeHightest(chunks, comparator)  
76.    
77.    
78. #加载词组字典和字符字典  
79. dictWord = {}  
80. maxWordLength = 0  
81.       
82. def loadDictChars(filepath):  
83.     global maxWordLength  
84.     fsock = file(filepath)  
85.     for line in fsock.readlines():  
86.         freq, word = line.split(' ')  
87.         word = unicode(word.strip(), 'utf-8')  
88.         dictWord[word] = (len(word), int(freq))  
89.         maxWordLength = maxWordLength < len(word) and len(word) or maxWordLength  
90.     fsock.close()  
91.       
92. def loadDictWords(filepath):  
93.     global maxWordLength  
94.     fsock = file(filepath)  
95.     for line in fsock.readlines():  
96.         word = unicode(line.strip(), 'utf-8')  
97.         dictWord[word] = (len(word), 0)  
98.         maxWordLength = maxWordLength < len(word) and len(word) or maxWordLength  
99.     fsock.close()  
100.   
101. #判断该词word是否在字典dictWord中      
102. def getDictWord(word):  
103.     result = dictWord.get(word)  
104.     if result:  
105.         return Word(word,result[1])  
106.     return None  
107.       
108. #开始加载字典  
109. def run():  
110.     from os.path import join, dirname  
111.     loadDictChars(join(dirname(__file__), 'data', 'chars.dic'))  
112.     loadDictWords(join(dirname(__file__), 'data', 'words.dic'))  
113.   
114. class Analysis:  
115.       
116.     def __init__(self,text):  
117.         if isinstance(text,unicode):  
118.             self.text = text  
119.         else:  
120.             self.text = text.encode('utf-8')  
121.         self.cacheSize = 3  
122.         self.pos = 0  
123.         self.textLength = len(self.text)  
124.         self.cache = []  
125.         self.cacheIndex = 0  
126.         self.complexCompare = ComplexCompare()  
127.           
128.         #简单小技巧，用到个缓存，不知道具体有没有用处  
129.         for i in range(self.cacheSize):  
130.             self.cache.append([-1,Word()])  
131.           
132.         #控制字典只加载一次  
133.         if not dictWord:  
134.             run()  
135.       
136.     def __iter__(self):  
137.         while True:  
138.             token = self.getNextToken()  
139.             if token == None:  
140.                 raise StopIteration  
141.             yield token  
142.               
143.     def getNextChar(self):  
144.         return self.text[self.pos]  
145.           
146.     #判断该字符是否是中文字符（不包括中文标点）    
147.     def isChineseChar(self,charater):  
148.         return 0x4e00 <= ord(charater) < 0x9fa6  
149.           
150.     #判断是否是ASCII码  
151.     def isASCIIChar(self, ch):  
152.         import string  
153.         if ch in string.whitespace:  
154.             return False  
155.         if ch in string.punctuation:  
156.             return False  
157.         return ch in string.printable  
158.       
159.     #得到下一个切割结果  
160.     def getNextToken(self):  
161.         while self.pos < self.textLength:  
162.             if self.isChineseChar(self.getNextChar()):  
163.                 token = self.getChineseWords()  
164.             else :  
165.                 token = self.getASCIIWords()+'/'  
166.             if len(token) > 0:  
167.                 return token  
168.         return None  
169.       
170.     #切割出非中文词  
171.     def getASCIIWords(self):  
172.         # Skip pre-word whitespaces and punctuations  
173.         #跳过中英文标点和空格  
174.         while self.pos < self.textLength:  
175.             ch = self.getNextChar()  
176.             if self.isASCIIChar(ch) or self.isChineseChar(ch):  
177.                 break  
178.             self.pos += 1  
179.         #得到英文单词的起始位置      
180.         start = self.pos  
181.           
182.         #找出英文单词的结束位置  
183.         while self.pos < self.textLength:  
184.             ch = self.getNextChar()  
185.             if not self.isASCIIChar(ch):  
186.                 break  
187.             self.pos += 1  
188.         end = self.pos  
189.           
190.         #Skip chinese word whitespaces and punctuations  
191.         #跳过中英文标点和空格  
192.         while self.pos < self.textLength:  
193.             ch = self.getNextChar()  
194.             if self.isASCIIChar(ch) or self.isChineseChar(ch):  
195.                 break  
196.             self.pos += 1  
197.               
198.         #返回英文单词  
199.         return self.text[start:end]  
200.       
201.     #切割出中文词，并且做处理，用上述4种方法  
202.     def getChineseWords(self):  
203.         chunks = self.createChunks()  
204.         if len(chunks) > 1:  
205.             chunks = self.complexCompare.mmFilter(chunks)  
206.         if len(chunks) > 1:  
207.             chunks = self.complexCompare.lawlFilter(chunks)  
208.         if len(chunks) > 1:  
209.             chunks = self.complexCompare.svmlFilter(chunks)  
210.         if len(chunks) > 1:  
211.             chunks = self.complexCompare.logFreqFilter(chunks)  
212.         if len(chunks) == 0 :  
213.             return ''  
214.           
215.         #最后只有一种切割方法  
216.         word = chunks[0].words  
217.         token = ""  
218.         length = 0  
219.         for x in word:  
220.             if x.length <> -1:  
221.                 token += x.text + "/"  
222.                 length += len(x.text)  
223.         self.pos += length  
224.         return token  
225.       
226.     #三重循环来枚举切割方法，这里也可以运用递归来实现  
227.     def createChunks(self):  
228.         chunks = []  
229.         originalPos = self.pos  
230.         words1 = self.getMatchChineseWords()  
231.           
232.         for word1 in words1:  
233.             self.pos += len(word1.text)  
234.             if self.pos < self.textLength:  
235.                 words2 = self.getMatchChineseWords()  
236.                 for word2 in words2:  
237.                     self.pos += len(word2.text)  
238.                     if self.pos < self.textLength:  
239.                         words3 = self.getMatchChineseWords()  
240.                         for word3 in words3:  
241.                             print word3.length,word3.text  
242.                             if word3.length == -1:  
243.                                 chunk = Chunk(word1,word2)  
244.                                 print "Ture"  
245.                             else :  
246.                                 chunk = Chunk(word1,word2,word3)  
247.                             chunks.append(chunk)  
248.                     elif self.pos == self.textLength:  
249.                         chunks.append(Chunk(word1,word2))  
250.                     self.pos -= len(word2.text)  
251.             elif self.pos == self.textLength:  
252.                 chunks.append(Chunk(word1))  
253.             self.pos -= len(word1.text)  
254.                                   
255.         self.pos = originalPos  
256.         return chunks  
257.       
258.     #运用正向最大匹配算法结合字典来切割中文文本    
259.     def getMatchChineseWords(self):  
260.         #use cache,check it   
261.         for i in range(self.cacheSize):  
262.             if self.cache[i][0] == self.pos:  
263.                 return self.cache[i][1]  
264.               
265.         originalPos = self.pos  
266.         words = []  
267.         index = 0  
268.         while self.pos < self.textLength:  
269.             if index >= maxWordLength :  
270.                 break  
271.             if not self.isChineseChar(self.getNextChar()):  
272.                 break  
273.             self.pos += 1  
274.             index += 1  
275.               
276.             text = self.text[originalPos:self.pos]  
277.             word = getDictWord(text)  
278.             if word:  
279.                 words.append(word)  
280.                   
281.         self.pos = originalPos  
282.         #没有词则放置个‘X’，将文本长度标记为-1  
283.         if not words:  
284.             word = Word()  
285.             word.length = -1  
286.             word.text = 'X'  
287.             words.append(word)  
288.           
289.         self.cache[self.cacheIndex] = (self.pos,words)  
290.         self.cacheIndex += 1  
291.         if self.cacheIndex >= self.cacheSize:  
292.             self.cacheIndex = 0  
293.         return words  

接下来给出测试代码：

[python] view plaincopyprint?

1. #-*- coding: utf-8 -*-  
2.   
3. from PyMmseg import Analysis  
4.   
5. def cuttest(text):  
6.     #cut =  Analysis(text)  
7.     wlist = [word for word in Analysis(text)]  
8.     tmp = ""  
9.     for w in wlist:  
10.         tmp += w  
11.     print tmp  
12.     print "================================"  
13.           
14. if __name__=="__main__":  
15. #    cuttest(u"研究生命来源")  
16. #    cuttest(u"南京市长江大桥欢迎您")  
17. #    cuttest(u"请把手抬高一点儿")  
18. #    cuttest(u"长春市长春节致词。")  
19. #    cuttest(u"长春市长春药店。")  
20. #    cuttest(u"我的和服务必在明天做好。")  
21. #    cuttest(u"我发现有很多人喜欢他。")  
22. #    cuttest(u"我喜欢看电视剧大长今。")  
23. #    cuttest(u"半夜给拎起来陪看欧洲杯糊着两眼半晌没搞明白谁和谁踢。")  
24. #    cuttest(u"李智伟高高兴兴以及王晓薇出去玩，后来智伟和晓薇又单独去玩了。")  
25. #    cuttest(u"一次性交出去很多钱。 ")  
26. #    cuttest(u"这是一个伸手不见五指的黑夜。我叫孙悟空，我爱北京，我爱Python和C++。")  
27. #    cuttest(u"我不喜欢日本和服。")  
28. #    cuttest(u"雷猴回归人间。")  
29. #    cuttest(u"工信处女干事每月经过下属科室都要亲口交代24口交换机等技术性器件的安装工作")  
30. #    cuttest(u"我需要廉租房")  
31. #    cuttest(u"永和服装饰品有限公司")  
32. #    cuttest(u"我爱北京天安门")  
33. #    cuttest(u"abc")  
34. #    cuttest(u"隐马尔可夫")  
35. #    cuttest(u"雷猴是个好网站")  
36. #    cuttest(u"“Microsoft”一词由“MICROcomputer（微型计算机）”和“SOFTware（软件）”两部分组成")  
37. #    cuttest(u"草泥马和欺实马是今年的流行词汇")  
38. #    cuttest(u"伊藤洋华堂总府店")  
39. #    cuttest(u"中国科学院计算技术研究所")  
40. #    cuttest(u"罗密欧与朱丽叶")  
41. #    cuttest(u"我购买了道具和服装")  
42. #    cuttest(u"PS: 我觉得开源有一个好处，就是能够敦促自己不断改进，避免敞帚自珍")  
43. #    cuttest(u"湖北省石首市")  
44. #    cuttest(u"总经理完成了这件事情")  
45. #    cuttest(u"电脑修好了")  
46. #    cuttest(u"做好了这件事情就一了百了了")  
47. #    cuttest(u"人们审美的观点是不同的")  
48. #    cuttest(u"我们买了一个美的空调")  
49. #    cuttest(u"线程初始化时我们要注意")  
50. #    cuttest(u"一个分子是由好多原子组织成的")  
51. #    cuttest(u"祝你马到功成")  
52. #    cuttest(u"他掉进了无底洞里")  
53. #    cuttest(u"中国的首都是北京")  
54. #    cuttest(u"孙君意")  
55. #    cuttest(u"外交部发言人马朝旭")  
56. #    cuttest(u"领导人会议和第四届东亚峰会")  
57. #    cuttest(u"在过去的这五年")  
58. #    cuttest(u"还需要很长的路要走")  
59. #    cuttest(u"60周年首都阅兵")  
60. #    cuttest(u"你好人们审美的观点是不同的")  
61. #    cuttest(u"买水果然后来世博园")  
62. #    cuttest(u"买水果然后去世博园")  
63. #    cuttest(u"但是后来我才知道你是对的")  
64. #    cuttest(u"存在即合理")  
65. #    cuttest(u"的的的的的在的的的的就以和和和")  
66. #    cuttest(u"I love你，不以为耻，反以为rong")  
67. #    cuttest(u" ")  
68. #    cuttest(u"")  
69. #    cuttest(u"hello你好人们审美的观点是不同的")  
70. #    cuttest(u"很好但主要是基于网页形式")  
71. #    cuttest(u"hello你好人们审美的观点是不同的")  
72. #    cuttest(u"为什么我不能拥有想要的生活")  
73. #    cuttest(u"后来我才")  
74. #    cuttest(u"此次来中国是为了")  
75. #    cuttest(u"使用了它就可以解决一些问题")  
76. #    cuttest(u",使用了它就可以解决一些问题")  
77. #    cuttest(u"其实使用了它就可以解决一些问题")  
78. #    cuttest(u"好人使用了它就可以解决一些问题")  
79. #    cuttest(u"是因为和国家")  
80. #    cuttest(u"老年搜索还支持")  
81.     cuttest(u"干脆就把那部蒙人的闲法给废了拉倒！RT @laoshipukong : 27日，全国人大常委会第三次审议侵权责任法草案，删除了有关医疗损害责任“举证倒置”的规定。在医患纠纷中本已处于弱势地位的消费者由此将陷入万劫不复的境地。 ")  

mmseg算法虽然比正向最大匹配算法在错误率上要好很多，但是也不是很完美，如果想让其比较好的运行需要一个很庞大的字典库外加词频库，其中词频库很重要，而且个人觉得判断依据3（标准差）的判断比较不靠谱，这里应该优化一下，能加上词频就更完美了。

关于mmseg的Python源码还是比较少的，除了google code 有一个外，其他的我没见过。倒是C++的代码比较多，但是代码写的比较乱，没有任何注释，STL乱用，也不说明用途，反正我是懒的看啦。

另外字典没有给出，我也不知道怎么上传压缩包，这类简单的字典网上到还是挺多的。

分词算法就研究到这里吧，这也只是用来自己娱乐，专业的分词包很多，更高级的分词算法也很多，当然最后直接拿来用就可以了，没有必要去纠结人家是怎么实现的，最基本的是要知道分词的原理和分词的算法，有算法还怕实现不了？