HEVC学习（十四） —— SAO函数解析之二

这篇博客写得不错，觉得对我对大家刚开始学习时会有帮助，于是转载之。原文地址：http://blog.youkuaiyun.com/hevc_cjl/article/details/8288572

1. /** rate distortion optimization of all SAO units
2. * \param saoParam SAO parameters
3. * \param lambda
4. * \param lambdaChroma
5. */
6. #if SAO_ENCODING_CHOICE
7. Void TEncSampleAdaptiveOffset::rdoSaoUnitAll(SAOParam *saoParam, Double lambda, Double lambdaChroma, Int depth)
8. #else
9. Void TEncSampleAdaptiveOffset::rdoSaoUnitAll(SAOParam *saoParam, Double lambda, Double lambdaChroma)
10. #endif
11. {
13. Int idxY;
14. Int idxX;
15. Int frameHeightInCU = saoParam->numCuInHeight;
16. Int frameWidthInCU = saoParam->numCuInWidth;
17. Int j, k;
18. Int addr = 0;
19. Int addrUp = -1;
20. Int addrLeft = -1;
21. Int compIdx = 0;
22. SaoLcuParam mergeSaoParam[3][2];
23. Double compDistortion[3];
25. saoParam->bSaoFlag[0] = true;
26. saoParam->bSaoFlag[1] = true;
27. saoParam->oneUnitFlag[0] = false;
28. saoParam->oneUnitFlag[1] = false;
29. saoParam->oneUnitFlag[2] = false;
31. #if SAO_ENCODING_CHOICE
32. #if SAO_ENCODING_CHOICE_CHROMA
33. Int numNoSao[2];
34. numNoSao[0] = 0;// Luma
35. numNoSao[1] = 0;// Chroma
36. if( depth > 0 && m_depthSaoRate[0][depth-1] > SAO_ENCODING_RATE )
37. {
38. saoParam->bSaoFlag[0] = false;
39. }
40. if( depth > 0 && m_depthSaoRate[1][depth-1] > SAO_ENCODING_RATE_CHROMA )
41. {
42. saoParam->bSaoFlag[1] = false;
43. }
44. #else
45. Int numNoSao = 0;
47. if( depth > 0 && m_depth0SaoRate > SAO_ENCODING_RATE )
48. {
49. saoParam->bSaoFlag[0] = false;
50. saoParam->bSaoFlag[1] = false;
51. }
52. #endif
53. #endif
54. //!< 以LCU为单位对图像中的每个LCU进行遍历
55. for (idxY = 0; idxY< frameHeightInCU; idxY++)
56. {
57. for (idxX = 0; idxX< frameWidthInCU; idxX++)
58. {
59. addr = idxX + frameWidthInCU*idxY; //!< 当前LCU地址
60. addrUp = addr < frameWidthInCU ? -1:idxX + frameWidthInCU*(idxY-1); //!< 当前LCU上邻块地址
61. addrLeft = idxX == 0 ? -1:idxX-1 + frameWidthInCU*idxY; //!< 当前LCU左邻块地址
62. Int allowMergeLeft = 1;
63. Int allowMergeUp = 1;
64. UInt rate;
65. Double bestCost, mergeCost;
66. if (idxX!=0) //!< 非第1列
67. {
68. // check tile id and slice id //! 检查当前LCU与其左邻块是否属于同一个tile以及是否属于同一个slice，不同的话该邻块不可用
69. if ( (m_pcPic->getPicSym()->getTileIdxMap(addr-1) != m_pcPic->getPicSym()->getTileIdxMap(addr)) || (m_pcPic->getCU(addr-1)->getSlice()->getSliceIdx() != m_pcPic->getCU(addr)->getSlice()->getSliceIdx()))
70. {
71. allowMergeLeft = 0; //!< 左邻块不可用
72. }
73. }
74. else
75. {
76. allowMergeLeft = 0; //!< 第1列的左邻块均不可用
77. }
78. if (idxY!=0) //!< 非第1行
79. {//! 检查当前LCU与其上邻块是否属于同一个tile以及是否属于同一个slice，不同的话该邻块不可用
80. if ( (m_pcPic->getPicSym()->getTileIdxMap(addr-m_iNumCuInWidth) != m_pcPic->getPicSym()->getTileIdxMap(addr)) || (m_pcPic->getCU(addr-m_iNumCuInWidth)->getSlice()->getSliceIdx() != m_pcPic->getCU(addr)->getSlice()->getSliceIdx()))
81. {
82. allowMergeUp = 0; //!< 上邻块不可用
83. }
84. }
85. else
86. {
87. allowMergeUp = 0; //!< 第1行的上邻块均不可用
88. }
90. compDistortion[0] = 0; //!< Y distortion
91. compDistortion[1] = 0; //!< Cb distortion
92. compDistortion[2] = 0; //!< Cr distortion
93. m_pcRDGoOnSbacCoder->load(m_pppcRDSbacCoder[0][CI_CURR_BEST]);
94. if (allowMergeLeft)
95. {
96. m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(0); //!< 编码句法元素sao_merge_left_flag
97. }
98. if (allowMergeUp)
99. {
100. m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(0); //!< 编码句法元素sao_merge_up_flag
101. }
102. m_pcRDGoOnSbacCoder->store( m_pppcRDSbacCoder[0][CI_TEMP_BEST] );
103. // reset stats Y, Cb, Cr
104. for ( compIdx=0;compIdx<3;compIdx++)
105. {
106. for ( j=0;j<MAX_NUM_SAO_TYPE;j++)
107. {
108. for ( k=0;k< MAX_NUM_SAO_CLASS;k++)
109. {
110. m_iOffset [compIdx][j][k] = 0;
111. if( m_saoLcuBasedOptimization && m_saoLcuBoundary ){ //!< true && false
112. m_iCount [compIdx][j][k] = m_count_PreDblk [addr][compIdx][j][k];
113. m_iOffsetOrg[compIdx][j][k] = m_offsetOrg_PreDblk[addr][compIdx][j][k];
114. }
115. else
116. {
117. m_iCount [compIdx][j][k] = 0;
118. m_iOffsetOrg[compIdx][j][k] = 0;
119. }
120. }
121. }
122. saoParam->saoLcuParam[compIdx][addr].typeIdx = -1;
123. saoParam->saoLcuParam[compIdx][addr].mergeUpFlag = 0;
124. saoParam->saoLcuParam[compIdx][addr].mergeLeftFlag = 0;
125. saoParam->saoLcuParam[compIdx][addr].subTypeIdx = 0;
126. #if SAO_ENCODING_CHOICE
127. if( (compIdx ==0 && saoParam->bSaoFlag[0])|| (compIdx >0 && saoParam->bSaoFlag[1]) )
128. #endif
129. {//! 统计BO和EO各个模式下，对应classIdx下滤波前的重建像素值与原始像素值的差值的总和，以及对classIdx的计数
130. calcSaoStatsCu(addr, compIdx, compIdx);
131. }
132. }
133. //!< Y分量最佳滤波模式的选择
134. saoComponentParamDist(allowMergeLeft, allowMergeUp, saoParam, addr, addrUp, addrLeft, 0, lambda, &mergeSaoParam[0][0], &compDistortion[0]);
135. //!< CbCr分量最佳滤波模式的选择
136. sao2ChromaParamDist(allowMergeLeft, allowMergeUp, saoParam, addr, addrUp, addrLeft, lambdaChroma, &mergeSaoParam[1][0], &mergeSaoParam[2][0], &compDistortion[0]);
137. if( saoParam->bSaoFlag[0] || saoParam->bSaoFlag[1] )
138. {
139. // Cost of new SAO_params
140. m_pcRDGoOnSbacCoder->load(m_pppcRDSbacCoder[0][CI_CURR_BEST]);
141. m_pcRDGoOnSbacCoder->resetBits();
142. if (allowMergeLeft)
143. {
144. m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(0);
145. }
146. if (allowMergeUp)
147. {
148. m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(0);
149. }
150. for ( compIdx=0;compIdx<3;compIdx++)
151. {
152. if( (compIdx ==0 && saoParam->bSaoFlag[0]) || (compIdx >0 && saoParam->bSaoFlag[1]))
153. {
154. m_pcEntropyCoder->encodeSaoOffset(&saoParam->saoLcuParam[compIdx][addr], compIdx);
155. }
156. }
158. rate = m_pcEntropyCoder->getNumberOfWrittenBits();
159. bestCost = compDistortion[0] + (Double)rate;
160. m_pcRDGoOnSbacCoder->store(m_pppcRDSbacCoder[0][CI_TEMP_BEST]);
162. // Cost of Merge
163. for(Int mergeUp=0; mergeUp<2; ++mergeUp)
164. {
165. if ( (allowMergeLeft && (mergeUp==0)) || (allowMergeUp && (mergeUp==1)) )
166. {
167. m_pcRDGoOnSbacCoder->load(m_pppcRDSbacCoder[0][CI_CURR_BEST]);
168. m_pcRDGoOnSbacCoder->resetBits();
169. if (allowMergeLeft)
170. {
171. m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(1-mergeUp);
172. }
173. if ( allowMergeUp && (mergeUp==1) )
174. {
175. m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(1);
176. }
178. rate = m_pcEntropyCoder->getNumberOfWrittenBits();
179. mergeCost = compDistortion[mergeUp+1] + (Double)rate;
180. if (mergeCost < bestCost)
181. {
182. bestCost = mergeCost;
183. m_pcRDGoOnSbacCoder->store(m_pppcRDSbacCoder[0][CI_TEMP_BEST]);
184. for ( compIdx=0;compIdx<3;compIdx++)
185. {
186. mergeSaoParam[compIdx][mergeUp].mergeLeftFlag = 1-mergeUp;
187. mergeSaoParam[compIdx][mergeUp].mergeUpFlag = mergeUp;
188. if( (compIdx==0 && saoParam->bSaoFlag[0]) || (compIdx>0 && saoParam->bSaoFlag[1]))
189. {
190. copySaoUnit(&saoParam->saoLcuParam[compIdx][addr], &mergeSaoParam[compIdx][mergeUp] );
191. }
192. }
193. }
194. }
195. }
196. #if SAO_ENCODING_CHOICE
197. #if SAO_ENCODING_CHOICE_CHROMA
198. if( saoParam->saoLcuParam[0][addr].typeIdx == -1) //!< Y分量不存在SAO参数
199. {
200. numNoSao[0]++;
201. }
202. if( saoParam->saoLcuParam[1][addr].typeIdx == -1) //!< CbCr分量不存在SAO参数
203. {
204. numNoSao[1]+=2;
205. }
206. #else
207. for ( compIdx=0;compIdx<3;compIdx++)
208. {
209. if( depth == 0 && saoParam->saoLcuParam[compIdx][addr].typeIdx == -1)
210. {
211. numNoSao++;
212. }
213. }
214. #endif
215. #endif
216. m_pcRDGoOnSbacCoder->load(m_pppcRDSbacCoder[0][CI_TEMP_BEST]);
217. m_pcRDGoOnSbacCoder->store(m_pppcRDSbacCoder[0][CI_CURR_BEST]);
218. } //!< if( saoParam->bSaoFlag[0] || saoParam->bSaoFlag[1] )
219. } //!< for (idxX = 0; idxX< frameWidthInCU; idxX++)
220. } //!< for (idxY = 0; idxY< frameHeightInCU; idxY++)
221. #if SAO_ENCODING_CHOICE
222. #if SAO_ENCODING_CHOICE_CHROMA
223. #if SAO_ENCODING_CHOICE_CHROMA_BF
224. if( !saoParam->bSaoFlag[0])
225. {
226. m_depthSaoRate[0][depth] = 1.0;
227. }
228. else
229. {
230. m_depthSaoRate[0][depth] = numNoSao[0]/((Double) frameHeightInCU*frameWidthInCU);
231. }
232. if( !saoParam->bSaoFlag[1])
233. {
234. m_depthSaoRate[1][depth] = 1.0;
235. }
236. else
237. {
238. m_depthSaoRate[1][depth] = numNoSao[1]/((Double) frameHeightInCU*frameWidthInCU*2);
239. }
240. #else
241. m_depthSaoRate[0][depth] = numNoSao[0]/((Double) frameHeightInCU*frameWidthInCU);
242. m_depthSaoRate[1][depth] = numNoSao[1]/((Double) frameHeightInCU*frameWidthInCU*2);
243. #endif
244. #else
245. if( depth == 0)
246. {
247. // update SAO Rate
248. m_depth0SaoRate = numNoSao/((Double) frameHeightInCU*frameWidthInCU*3);
249. }
250. #endif
251. #endif
253. }

/** rate distortion optimization of all SAO units
 * \param saoParam SAO parameters
 * \param lambda 
 * \param lambdaChroma
 */
#if SAO_ENCODING_CHOICE
Void TEncSampleAdaptiveOffset::rdoSaoUnitAll(SAOParam *saoParam, Double lambda, Double lambdaChroma, Int depth)
#else
Void TEncSampleAdaptiveOffset::rdoSaoUnitAll(SAOParam *saoParam, Double lambda, Double lambdaChroma)
#endif
{

  Int idxY;
  Int idxX;
  Int frameHeightInCU = saoParam->numCuInHeight;
  Int frameWidthInCU  = saoParam->numCuInWidth;
  Int j, k;
  Int addr = 0;
  Int addrUp = -1;
  Int addrLeft = -1;
  Int compIdx = 0;
  SaoLcuParam mergeSaoParam[3][2];
  Double compDistortion[3];

  saoParam->bSaoFlag[0] = true;
  saoParam->bSaoFlag[1] = true;
  saoParam->oneUnitFlag[0] = false;
  saoParam->oneUnitFlag[1] = false;
  saoParam->oneUnitFlag[2] = false;

#if SAO_ENCODING_CHOICE
#if SAO_ENCODING_CHOICE_CHROMA
  Int numNoSao[2];
  numNoSao[0] = 0;// Luma 
  numNoSao[1] = 0;// Chroma 
  if( depth > 0 && m_depthSaoRate[0][depth-1] > SAO_ENCODING_RATE )
  {
    saoParam->bSaoFlag[0] = false;
  }
  if( depth > 0 && m_depthSaoRate[1][depth-1] > SAO_ENCODING_RATE_CHROMA )
  {
    saoParam->bSaoFlag[1] = false;
  }
#else
  Int numNoSao = 0;

  if( depth > 0 && m_depth0SaoRate > SAO_ENCODING_RATE )
  {
    saoParam->bSaoFlag[0] = false;
    saoParam->bSaoFlag[1] = false;
  }
#endif
#endif
  //!< 以LCU为单位对图像中的每个LCU进行遍历
  for (idxY = 0; idxY< frameHeightInCU; idxY++)
  {
    for (idxX = 0; idxX< frameWidthInCU; idxX++)
    {
      addr     = idxX  + frameWidthInCU*idxY; //!< 当前LCU地址
      addrUp   = addr < frameWidthInCU ? -1:idxX   + frameWidthInCU*(idxY-1); //!< 当前LCU上邻块地址
      addrLeft = idxX == 0               ? -1:idxX-1 + frameWidthInCU*idxY; //!< 当前LCU左邻块地址
      Int allowMergeLeft = 1;
      Int allowMergeUp   = 1;
      UInt rate;
      Double bestCost, mergeCost;
      if (idxX!=0) //!< 非第1列
      { 
        // check tile id and slice id //! 检查当前LCU与其左邻块是否属于同一个tile以及是否属于同一个slice，不同的话该邻块不可用
        if ( (m_pcPic->getPicSym()->getTileIdxMap(addr-1) != m_pcPic->getPicSym()->getTileIdxMap(addr)) || (m_pcPic->getCU(addr-1)->getSlice()->getSliceIdx() != m_pcPic->getCU(addr)->getSlice()->getSliceIdx()))
        {
          allowMergeLeft = 0; //!< 左邻块不可用
        }
      }
      else
      {
        allowMergeLeft = 0; //!< 第1列的左邻块均不可用
      }
      if (idxY!=0) //!< 非第1行
      {//! 检查当前LCU与其上邻块是否属于同一个tile以及是否属于同一个slice，不同的话该邻块不可用
        if ( (m_pcPic->getPicSym()->getTileIdxMap(addr-m_iNumCuInWidth) != m_pcPic->getPicSym()->getTileIdxMap(addr)) || (m_pcPic->getCU(addr-m_iNumCuInWidth)->getSlice()->getSliceIdx() != m_pcPic->getCU(addr)->getSlice()->getSliceIdx()))
        {
          allowMergeUp = 0; //!< 上邻块不可用
        }
      }
      else
      {
        allowMergeUp = 0; //!< 第1行的上邻块均不可用
      }

      compDistortion[0] = 0; //!< Y distortion
      compDistortion[1] = 0; //!< Cb distortion
      compDistortion[2] = 0; //!< Cr distortion
      m_pcRDGoOnSbacCoder->load(m_pppcRDSbacCoder[0][CI_CURR_BEST]);
      if (allowMergeLeft)
      {
        m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(0); //!< 编码句法元素sao_merge_left_flag
      }
      if (allowMergeUp)
      {
        m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(0); //!< 编码句法元素sao_merge_up_flag
      }
      m_pcRDGoOnSbacCoder->store( m_pppcRDSbacCoder[0][CI_TEMP_BEST] );
      // reset stats Y, Cb, Cr
      for ( compIdx=0;compIdx<3;compIdx++)
      {
        for ( j=0;j<MAX_NUM_SAO_TYPE;j++)
        {
          for ( k=0;k< MAX_NUM_SAO_CLASS;k++)
          {
            m_iOffset   [compIdx][j][k] = 0;
            if( m_saoLcuBasedOptimization && m_saoLcuBoundary ){ //!< true && false
              m_iCount    [compIdx][j][k] = m_count_PreDblk    [addr][compIdx][j][k];
              m_iOffsetOrg[compIdx][j][k] = m_offsetOrg_PreDblk[addr][compIdx][j][k];
            }
            else
            {
              m_iCount    [compIdx][j][k] = 0;
              m_iOffsetOrg[compIdx][j][k] = 0;
            }
          }  
        }
        saoParam->saoLcuParam[compIdx][addr].typeIdx       =  -1;
        saoParam->saoLcuParam[compIdx][addr].mergeUpFlag   = 0;
        saoParam->saoLcuParam[compIdx][addr].mergeLeftFlag = 0;
        saoParam->saoLcuParam[compIdx][addr].subTypeIdx    = 0;
#if SAO_ENCODING_CHOICE
  if( (compIdx ==0 && saoParam->bSaoFlag[0])|| (compIdx >0 && saoParam->bSaoFlag[1]) )
#endif
        {//! 统计BO和EO各个模式下，对应classIdx下滤波前的重建像素值与原始像素值的差值的总和，以及对classIdx的计数
          calcSaoStatsCu(addr, compIdx,  compIdx);
		}
      }
	  //!< Y分量最佳滤波模式的选择
      saoComponentParamDist(allowMergeLeft, allowMergeUp, saoParam, addr, addrUp, addrLeft, 0,  lambda, &mergeSaoParam[0][0], &compDistortion[0]);
	  //!< CbCr分量最佳滤波模式的选择
      sao2ChromaParamDist(allowMergeLeft, allowMergeUp, saoParam, addr, addrUp, addrLeft, lambdaChroma, &mergeSaoParam[1][0], &mergeSaoParam[2][0], &compDistortion[0]);
     if( saoParam->bSaoFlag[0] || saoParam->bSaoFlag[1] )
      {
        // Cost of new SAO_params
        m_pcRDGoOnSbacCoder->load(m_pppcRDSbacCoder[0][CI_CURR_BEST]);
        m_pcRDGoOnSbacCoder->resetBits();
        if (allowMergeLeft)
        {
          m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(0); 
        }
        if (allowMergeUp)
        {
          m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(0);
        }
        for ( compIdx=0;compIdx<3;compIdx++)
        {
        if( (compIdx ==0 && saoParam->bSaoFlag[0]) || (compIdx >0 && saoParam->bSaoFlag[1]))
          {
           m_pcEntropyCoder->encodeSaoOffset(&saoParam->saoLcuParam[compIdx][addr], compIdx);
          }
        }

        rate = m_pcEntropyCoder->getNumberOfWrittenBits();
        bestCost = compDistortion[0] + (Double)rate;
        m_pcRDGoOnSbacCoder->store(m_pppcRDSbacCoder[0][CI_TEMP_BEST]);

        // Cost of Merge
        for(Int mergeUp=0; mergeUp<2; ++mergeUp)
        {
          if ( (allowMergeLeft && (mergeUp==0)) || (allowMergeUp && (mergeUp==1)) )
          {
            m_pcRDGoOnSbacCoder->load(m_pppcRDSbacCoder[0][CI_CURR_BEST]);
            m_pcRDGoOnSbacCoder->resetBits();
            if (allowMergeLeft)
            {
              m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(1-mergeUp); 
            }
            if ( allowMergeUp && (mergeUp==1) )
            {
              m_pcEntropyCoder->m_pcEntropyCoderIf->codeSaoMerge(1); 
            }

            rate = m_pcEntropyCoder->getNumberOfWrittenBits();
            mergeCost = compDistortion[mergeUp+1] + (Double)rate;
            if (mergeCost < bestCost)
            {
              bestCost = mergeCost;
              m_pcRDGoOnSbacCoder->store(m_pppcRDSbacCoder[0][CI_TEMP_BEST]);              
              for ( compIdx=0;compIdx<3;compIdx++)
              {
                mergeSaoParam[compIdx][mergeUp].mergeLeftFlag = 1-mergeUp;
                mergeSaoParam[compIdx][mergeUp].mergeUpFlag = mergeUp;
                if( (compIdx==0 && saoParam->bSaoFlag[0]) || (compIdx>0 && saoParam->bSaoFlag[1]))
                {
                  copySaoUnit(&saoParam->saoLcuParam[compIdx][addr], &mergeSaoParam[compIdx][mergeUp] );             
                }
              }
            }
          }
        }
#if SAO_ENCODING_CHOICE
#if SAO_ENCODING_CHOICE_CHROMA
if( saoParam->saoLcuParam[0][addr].typeIdx == -1) //!< Y分量不存在SAO参数
{
  numNoSao[0]++;
}
if( saoParam->saoLcuParam[1][addr].typeIdx == -1) //!< CbCr分量不存在SAO参数
{
  numNoSao[1]+=2;
}
#else
        for ( compIdx=0;compIdx<3;compIdx++)
        {
          if( depth == 0 && saoParam->saoLcuParam[compIdx][addr].typeIdx == -1)
          {
            numNoSao++;
          }
        }
#endif
#endif
        m_pcRDGoOnSbacCoder->load(m_pppcRDSbacCoder[0][CI_TEMP_BEST]);
        m_pcRDGoOnSbacCoder->store(m_pppcRDSbacCoder[0][CI_CURR_BEST]);
      } //!< if( saoParam->bSaoFlag[0] || saoParam->bSaoFlag[1] )       
    } //!< for (idxX = 0; idxX< frameWidthInCU; idxX++)     
  } //!< for (idxY = 0; idxY< frameHeightInCU; idxY++)   
#if SAO_ENCODING_CHOICE
#if SAO_ENCODING_CHOICE_CHROMA
#if SAO_ENCODING_CHOICE_CHROMA_BF
  if( !saoParam->bSaoFlag[0]) 
  {
    m_depthSaoRate[0][depth] = 1.0;
  }
  else
  {
    m_depthSaoRate[0][depth] = numNoSao[0]/((Double) frameHeightInCU*frameWidthInCU);
  }
  if( !saoParam->bSaoFlag[1]) 
  {
    m_depthSaoRate[1][depth] = 1.0;
  }
  else 
  {
    m_depthSaoRate[1][depth] = numNoSao[1]/((Double) frameHeightInCU*frameWidthInCU*2);
  }
#else
m_depthSaoRate[0][depth] = numNoSao[0]/((Double) frameHeightInCU*frameWidthInCU);
m_depthSaoRate[1][depth] = numNoSao[1]/((Double) frameHeightInCU*frameWidthInCU*2);
#endif
#else
  if( depth == 0)
  {
    // update SAO Rate
    m_depth0SaoRate = numNoSao/((Double) frameHeightInCU*frameWidthInCU*3);
  }
#endif
#endif

}