?? minutia.c
字號:
/*############################################################################# * 文件名：minutia.c * 功能： 細節點的一些函數接口 * modified by PRTsinghua@hotmail.com#############################################################################*/#include <string.h>#include <stdlib.h>#include <math.h>#include "minutia.h"typedef struct iFvsMinutiaSet_t{ FvsInt_t nbminutia; FvsInt_t tablesize; FvsMinutia_t ptable[1]; } iFvsMinutiaSet_t;/****************************************************************************** * 功能：創建一個細節點集合 * 參數：size 集合的大小 * 返回：若失敗，返回空；否則返回新的對象句柄******************************************************************************/FvsMinutiaSet_t MinutiaSetCreate(const FvsInt_t size){ iFvsMinutiaSet_t* p = NULL; p = (iFvsMinutiaSet_t*)malloc(sizeof(iFvsMinutiaSet_t) +size*sizeof(FvsMinutia_t)); if (p!=NULL) { /* 表中無細節點 */ p->nbminutia = 0; p->tablesize = size; } return (FvsMinutiaSet_t)p;}/****************************************************************************** * 功能：銷毀細節點集合。 * 一旦銷毀，該對象不再可以為任何函數所用，直到重新申請。 * 參數：minutia 細節點集合 * 返回：無******************************************************************************/void MinutiaSetDestroy(FvsMinutiaSet_t minutia){ iFvsMinutiaSet_t* p = NULL; if (minutia==NULL) return; p = minutia; free(p);}/****************************************************************************** * 功能：獲得細節點集合的大小 * 參數：minutia 細節點集合 * 返回：細節點集合大小******************************************************************************/FvsInt_t MinutiaSetGetSize(const FvsMinutiaSet_t min){ const iFvsMinutiaSet_t* minutia = (iFvsMinutiaSet_t*)min; FvsInt_t nret = 0; if (minutia!=NULL) nret = minutia->tablesize; return nret;}/****************************************************************************** * 功能：細節點集合的實際元素個數 * 參數：minutia 細節點集合 * 返回：元素個數******************************************************************************/FvsInt_t MinutiaSetGetCount(const FvsMinutiaSet_t min){ const iFvsMinutiaSet_t* minutia = (iFvsMinutiaSet_t*)min; FvsInt_t nret = 0; if (minutia!=NULL) nret = minutia->nbminutia; return nret;}/****************************************************************************** * 功能：返回細節點集合的數據緩沖區指針 * 參數：minutia 細節點集合 * 返回：指針******************************************************************************/FvsMinutia_t* MinutiaSetGetBuffer(FvsMinutiaSet_t min){ iFvsMinutiaSet_t* minutia = (iFvsMinutiaSet_t*)min; FvsMinutia_t* pret = NULL; if (minutia!=NULL) pret = minutia->ptable; return pret;}/****************************************************************************** * 功能：清空細節點集合 * 參數：minutia 細節點集合 * 返回：錯誤編號******************************************************************************/FvsError_t MinutiaSetEmpty(FvsMinutiaSet_t min){ iFvsMinutiaSet_t* minutia = (iFvsMinutiaSet_t*)min; FvsError_t nRet = FvsOK; if (minutia!=NULL) minutia->nbminutia = 0; else nRet = FvsMemory; return nRet;}/****************************************************************************** * 功能：在集合中添加一個細節點，如果滿了，返回一個錯誤 * 參數：minutia 細節點集合 * x 細節點的X坐標 * y 細節點的Y坐標 * type 細節點類型 * angle 角度 * 返回：錯誤編號******************************************************************************/FvsError_t MinutiaSetAdd(FvsMinutiaSet_t min, const FvsFloat_t x, const FvsFloat_t y, const FvsMinutiaType_t type, const FvsFloat_t angle){ iFvsMinutiaSet_t* minutia = (iFvsMinutiaSet_t*)min; FvsError_t nRet = FvsOK; if (minutia->nbminutia < minutia->tablesize) { minutia->ptable[minutia->nbminutia].x = x; minutia->ptable[minutia->nbminutia].y = y; minutia->ptable[minutia->nbminutia].type = type; minutia->ptable[minutia->nbminutia].angle = angle; minutia->nbminutia++; } else /* 表中無空間 */ nRet = FvsMemory; return nRet;}static FvsError_t MinutiaSetCheckClean(FvsMinutiaSet_t min){ iFvsMinutiaSet_t* minutia = (iFvsMinutiaSet_t*)min; FvsError_t nRet = FvsOK; FvsFloat_t tx, ty, ta; FvsInt_t i, j; FvsMinutia_t* mi, *mj; tx = 4.0; ty = 4.0; ta = 0.5; if (minutia!=NULL && minutia->nbminutia > 1) { /* 檢驗表中是否已有該細節點 */ for (j = 0; j < minutia->nbminutia; j++) for (i = j+1; i < minutia->nbminutia; i++) { mi = minutia->ptable + i; mj = minutia->ptable + j; /* 比較細節點i，j */ /* 規則 1: 相似的細節點彼此靠近 -> 刪除一個 */ if ( (fabs(mi->x - mj->x ) < tx) && (fabs(mi->y - mj->y ) < ty) && (fabs(mi->angle - mj->angle) < ta) ) { minutia->nbminutia--; *mi = minutia->ptable[minutia->nbminutia]; } /* 規則 2: 方向相反，距離靠近 -> 同時刪除 */ } } else /* 表中無空間 */ nRet = FvsMemory; return nRet;}#define P(x,y) p[(x)+(y)*pitch]/****************************************************************************** * 功能：在圖像中畫出細節點，不改變背景 * 參數：minutia 細節點集合 * image 指紋圖像 * 返回：錯誤編號******************************************************************************/FvsError_t MinutiaSetDraw(const FvsMinutiaSet_t min, FvsImage_t image){ FvsInt_t w = ImageGetWidth(image); FvsInt_t pitch = ImageGetPitch(image); FvsByte_t* p = ImageGetBuffer(image); FvsInt_t n, x, y; FvsFloat_t fx, fy; FvsMinutia_t* minutia = MinutiaSetGetBuffer(min); FvsInt_t nbminutia = MinutiaSetGetCount(min); if (minutia==NULL || p==NULL) return FvsMemory; /* 畫出每個細節點 */ for (n = 0; n < nbminutia; n++) { x = (FvsInt_t)minutia[n].x; y = (FvsInt_t)minutia[n].y; if (x<w-5 && x>4) { if (y<h-5 && y>4) { switch (minutia[n].type) { case FvsMinutiaTypeEnding: P(x,y) = 0xFF; P(x-1, y) = 0xA0; P(x+1, y) = 0xA0; P(x, y-1) = 0xA0; P(x, y+1) = 0xA0; break; case FvsMinutiaTypeBranching: P(x,y) = 0xFF; P(x-1, y-1) = 0xA0; P(x+1, y-1) = 0xA0; P(x-1, y+1) = 0xA0; P(x+1, y+1) = 0xA0; break; default: continue; } fx = sin(minutia[n].angle); fy = -cos(minutia[n].angle); P(x+(int32_t)(fx) ,y+(int32_t)(fy) ) = 0xFF; P(x+(int32_t)(fx*2.0),y+(int32_t)(fy*2.0)) = 0xFF; P(x+(int32_t)(fx*3.0),y+(int32_t)(fy*3.0)) = 0xFF; P(x+(int32_t)(fx*4.0),y+(int32_t)(fy*4.0)) = 0xFF; P(x+(int32_t)(fx*5.0),y+(int32_t)(fy*5.0)) = 0xFF; } } } return FvsOK;}/* 宏定義 */#define P1 P(x ,y-1)#define P2 P(x+1,y-1)#define P3 P(x+1,y )#define P4 P(x+1,y+1)#define P5 P(x ,y+1)#define P6 P(x-1,y+1)#define P7 P(x-1,y )#define P8 P(x-1,y-1)/****************************************************************************** * 功能：從細化圖像中提取細節點，并儲存到集合中。 * 申請的細節點集合必須足夠大，如果太小了，滿了后會停止搜索細節點。 * 參數：minutia 細節點集合，用來保存細節點 * image 細化后的圖像 * direction 用來計算方向用 * mask 用來表示有效的指紋區域 * 返回：錯誤編號******************************************************************************/FvsError_t MinutiaSetExtract ( FvsMinutiaSet_t minutia, const FvsImage_t image, const FvsFloatField_t direction, const FvsImage_t mask ){ FvsInt_t w = ImageGetWidth(image); FvsInt_t h = ImageGetHeight(image); FvsInt_t pitch = ImageGetPitch(image); FvsInt_t pitchm = ImageGetPitch(mask); FvsByte_t* p = ImageGetBuffer(image); FvsByte_t* m = ImageGetBuffer(mask); FvsInt_t x, y; FvsFloat_t angle = 0.0; FvsInt_t whitecount; if (m==NULL || p==NULL) return FvsMemory; (void)MinutiaSetEmpty(minutia); /* 遍歷圖像，提取細節點 */ for (y = 1; y < h-1; y++) for (x = 1; x < w-1; x++) { if (m[x+y*pitchm]==0) continue; if (P(x,y)==0xFF) { whitecount = 0; if (P1!=0) whitecount++; if (P2!=0) whitecount++; if (P3!=0) whitecount++; if (P4!=0) whitecount++; if (P5!=0) whitecount++; if (P6!=0) whitecount++; if (P7!=0) whitecount++; if (P8!=0) whitecount++; switch(whitecount) { case 0: /* 孤立點，忽略 */ break; case 1: /* 檢測角度 */ angle = FloatFieldGetValue(direction, x, y); (void)MinutiaSetAdd(minutia, (FvsFloat_t)x, (FvsFloat_t)y, FvsMinutiaTypeEnding, (FvsFloat_t)angle); break; case 2: break; default: { angle = FloatFieldGetValue(direction, x, y); (void)MinutiaSetAdd(minutia, (FvsFloat_t)x, (FvsFloat_t)y, FvsMinutiaTypeBranching, (FvsFloat_t)angle); } break; } } } (void)MinutiaSetCheckClean(minutia); return FvsOK;}?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -