用哈夫曼编码压缩图像

最近在上图像处理的课，要完成的一些作业就顺便分享到这上面来了。本来是编辑好了，就贴代码就好了，但由于新手操作，一放代码就出问题导致心态崩了。所以就不讲原理了，试试看代码怎么放吧。枯了。

原理就看下图吧。（感觉变糊了）

代码思路大概就是，我认为只有出现在最后两位最小的概率值需要加码字0或者1。所以就可以正向直接求每一个概率值的码字。比如说第一个概率a1，排序后判断是否是最后一个，是的话就加一个码字0；如果是倒数第二个，就加码字1；如果两种情况都不是，就保持码字不变。

值得注意的是代码中出现的（find(S(j,:)==loc)==L）、（find(S(j,:)==loc)==L-1）其实是想要找到当前求的概率的位置，其实就是想要实现（find(ai)==L-1）。但是由于每次最后两个概率相加之后都要重新排序，所以a1的位置不可能总是在第一个（即ai不可能总是在第i个），所以每次排序都要记下排序前ai的位置loc（第一次排序的位置就是i），然后根据这个loc来判断下一次排序后ai的位置是否在最后两个。

loc：第一次排序时是i，因为哈夫曼编码先将概率从大到小排序；当不是第一次时，就是上一次排序的位置

find(S(j,:)==loc)：S保存的是每次排序前的位置，所以这一命令可以找到排序后，即当前该概率值ai的位置

find(S(j,:)==loc)==L：判断是否在最后一个

下面就是写的代码（尝试了几次，不知道怎么办我就只能贴图了）：