function Rn=refinement(N,LSP_Old)
% 函數 REFINEMENT()為精細編碼程序，對上一級編碼產生的重要系數列表LSP_Old，讀取每個
% 表項相應小波系數絕對值的二進制表示，輸出其中第N個重要的位，即相應于 2^N 處的碼數
% 輸入參數：N —— 本級編碼閾值的指數
%          LSP_Old —— 上一級編碼產生的重要系數列表
% 輸出參數：Rn —— 精細掃描輸出位流

global Mat
% Mat是輸入的小波分解系數矩陣，作為全局變量，在編碼的相關程序中使用

Rn=[];
% 每級精細掃描開始時，Rn 均為空表
% LSP_Old 非空時才執行精細掃描程序
if ~isempty(LSP_Old)
	rlsp=size(LSP_Old,1);
% 獲取 LSP_Old 的表項個數，對每個表項進行掃描
	for r=1:rlsp
        tMat=Mat(LSP_Old(r,1),LSP_Old(r,2));
% 讀取該表項對應的小波系數值
        [biLSP,Np]=fracnum2bin(abs(tMat),N);
% 函數 FRACNUM2BIN() 根據精細掃描對應的權位 N ，將任意的十進制正數轉換為二進制數，
% 輸出參數為二進制表示列表 biLSP 和 權位N與最高權位的距離 Np 。
        Rn=[Rn,biLSP(Np)];
% biLSP(Np)即為小波系數絕對值的二進制表示中第N個重要的位
	end
end