for(i=0; i<num-1; ++i) { p_latent_minutia_1 = & latent_template.m_minutiae[get<1>(corr[i])]; p_rolled_minutia_1 = & rolled_template.m_minutiae[get<2>(corr[i])]; for(j=i+1; j<num;++j) { p_latent_minutia_2 = & latent_template.m_minutiae[get<1>(corr[j])]; p_rolled_minutia_2 = & rolled_template.m_minutiae[get<2>(corr[j])]; dx_1 = p_latent_minutia_1->x-p_latent_minutia_2->x; dx_2 = p_rolled_minutia_1->x-p_rolled_minutia_2->x; dy_1 = p_latent_minutia_1->y-p_latent_minutia_2->y; dy_2 = p_rolled_minutia_1->y-p_rolled_minutia_2->y; dist_1 = (dx_1*dx_1)+(dy_1*dy_1); dist_1 = sqrt(dist_1); dist_2 = (dx_2*dx_2)+(dy_2*dy_2); dist_2 = sqrt(dist_2); dist = fabs(dist_1-dist_2); H[i*num+j] = (30-dist)/(25.0); if(H[i*num+j]>1) H[i*num+j] = 1.0; else if(H[i*num+j]<0) H[i*num+j] = 0.0; H[j*num+i] = H[i*num+j]; } }
时间: 2024-04-13 13:27:15 浏览: 105
这段代码是一个嵌套的循环,用于计算两个模板中指纹特征点之间的距离,并将计算结果存储在矩阵`H`中。
外层循环`for(i=0; i<num-1; ++i)`遍历`corr`向量中的元素(即匹配的特征点对),从第一个元素开始到倒数第二个元素。
内层循环`for(j=i+1; j<num;++j)`从外层循环的当前元素的下一个元素开始,遍历到最后一个元素。
在内层循环中,首先通过`get<1>(corr[i])`和`get<2>(corr[i])`获取到当前匹配对的特征点在`latent_template`和`rolled_template`中的索引,然后通过索引获取到对应的特征点。
接下来,计算两个特征点在x轴和y轴上的位移,分别存储在`dx_1`、`dx_2`、`dy_1`和`dy_2`变量中。
然后,分别计算两个特征点之间的欧氏距离`dist_1`和`dist_2`,并通过平方根函数求得真实的距离值。
接下来,计算两个距离之间的差异,并存储在`dist`变量中。
然后,根据一定的公式将距离差异转换为相似度,并存储在矩阵`H`中的相应位置。
最后,将矩阵`H`中对称位置的元素进行设置,保持对称性。
整个循环的目的是计算两个模板中所有特征点对之间的距离,并将相似度结果存储在矩阵`H`中,用于后续的匹配过程。
相关问题
int num = corr.size(); float *H = new float [num*num](); vector<short> flag_latent(latent_template.m_nrof_minu),flag_rolled(rolled_template.m_nrof_minu); register int i,j,k; MinuPoint *p_latent_minutia_1, *p_latent_minutia_2, *p_rolled_minutia_1, *p_rolled_minutia_2; float dist_1, dist_2, dist; float dx_1, dy_1, dx_2, dy_2;
这段代码是在Matcher::LSS_R_Fast2_Dist_eigen函数内部定义了一些变量和数组,并进行了初始化。下面是对这段代码的解释:
- `int num = corr.size();`:获取变量`corr`的大小,即元素的个数,并将其赋值给变量`num`。
- `float *H = new float [num*num]();`:动态分配一个大小为`num*num`的一维浮点型数组,并初始化为0。这个数组被命名为`H`。
- `vector<short> flag_latent(latent_template.m_nrof_minu),flag_rolled(rolled_template.m_nrof_minu);`:创建了两个名为`flag_latent`和`flag_rolled`的短整型向量,并分别初始化为`latent_template.m_nrof_minu`和`rolled_template.m_nrof_minu`。
- `register int i,j,k;`:定义了三个用于循环迭代的整型变量`i`,`j`和`k`。
- `MinuPoint *p_latent_minutia_1, *p_latent_minutia_2, *p_rolled_minutia_1, *p_rolled_minutia_2;`:定义了四个指向MinuPoint类型的指针变量,分别是`p_latent_minutia_1`,`p_latent_minutia_2`,`p_rolled_minutia_1`和`p_rolled_minutia_2`。
- `float dist_1, dist_2, dist;`:定义了三个浮点型变量`dist_1`,`dist_2`和`dist`。
- `float dx_1, dy_1, dx_2, dy_2;`:定义了四个浮点型变量`dx_1`,`dy_1`,`dx_2`和`dy_2`。
这段代码主要是声明和初始化了一些变量和数组,后续的代码可能会使用这些变量和数组来进行计算和存储数据。
if (save_corr){ ofstream output; output.open(corr_file); for(int i = 0; i < corr3.size(); i++){ output<<latent_minu_template.m_minutiae[get<1>(corr3[i])].x<<","<<latent_minu_template.m_minutiae[get<1>(corr3[i])].y <<","<<rolled_minu_template.m_minutiae[get<2>(corr3[i])].x<<","<<rolled_minu_template.m_minutiae[get<2>(corr3[i])].y<<endl; } output.close(); }
这段代码中的代码块是一个条件语句,判断变量save_corr是否为真。如果是真的话,执行以下操作:
1. 创建一个ofstream对象output。
2. 使用open函数打开名为corr_file的文件。
3. 使用循环遍历corr3容器中的元素。
4. 在每次迭代中,将一对坐标值分别从latent_minu_template和rolled_minu_template中取出,并将它们以逗号分隔的形式写入到output文件中。
5. 在循环结束后,关闭output文件。
这段代码的作用是将corr3容器中的一些坐标信息写入到文件中。如果save_corr为真,则将相应的坐标信息按照一定格式写入到名为corr_file的文件中。
阅读全文