|
|
|
@ -67,10 +67,6 @@ min4( float a, float b, float c, float d ) |
|
|
|
|
#define INSIDE 2 //unknown
|
|
|
|
|
#define CHANGE 3 //servise
|
|
|
|
|
|
|
|
|
|
// Processing data types
|
|
|
|
|
typedef double data_type; |
|
|
|
|
int data_type_cv = CV_64F; |
|
|
|
|
|
|
|
|
|
typedef struct CvHeapElem |
|
|
|
|
{ |
|
|
|
|
float T; |
|
|
|
@ -281,7 +277,7 @@ icvCalcFMM(const CvMat *f, CvMat *t, CvPriorityQueueFloat *Heap, bool negate) { |
|
|
|
|
} |
|
|
|
|
} |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
template <typename data_type> |
|
|
|
|
static void |
|
|
|
|
icvTeleaInpaintFMM(const CvMat *f, CvMat *t, CvMat *out, int range, CvPriorityQueueFloat *Heap ) { |
|
|
|
|
int i = 0, j = 0, ii = 0, jj = 0, k, l, q, color = 0; |
|
|
|
@ -513,7 +509,7 @@ icvTeleaInpaintFMM(const CvMat *f, CvMat *t, CvMat *out, int range, CvPriorityQu |
|
|
|
|
} |
|
|
|
|
} |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
template <typename data_type> |
|
|
|
|
static void |
|
|
|
|
icvNSInpaintFMM(const CvMat *f, CvMat *t, CvMat *out, int range, CvPriorityQueueFloat *Heap) { |
|
|
|
|
int i = 0, j = 0, ii = 0, jj = 0, k, l, q, color = 0; |
|
|
|
@ -644,28 +640,28 @@ icvNSInpaintFMM(const CvMat *f, CvMat *t, CvMat *out, int range, CvPriorityQueue |
|
|
|
|
|
|
|
|
|
if (CV_MAT_ELEM(*f,uchar,k+1,l)!=INSIDE) { |
|
|
|
|
if (CV_MAT_ELEM(*f,uchar,k-1,l)!=INSIDE) { |
|
|
|
|
gradI.x=(float)(abs(CV_MAT_ELEM(*out,data_type,kp+1,lm)-CV_MAT_ELEM(*out,data_type,kp,lm))+ |
|
|
|
|
abs(CV_MAT_ELEM(*out,data_type,kp,lm)-CV_MAT_ELEM(*out,data_type,km-1,lm))); |
|
|
|
|
gradI.x=(float)(std::abs(CV_MAT_ELEM(*out,data_type,kp+1,lm)-CV_MAT_ELEM(*out,data_type,kp,lm))+ |
|
|
|
|
std::abs(CV_MAT_ELEM(*out,data_type,kp,lm)-CV_MAT_ELEM(*out,data_type,km-1,lm))); |
|
|
|
|
} else { |
|
|
|
|
gradI.x=(float)(abs(CV_MAT_ELEM(*out,data_type,kp+1,lm)-CV_MAT_ELEM(*out,data_type,kp,lm)))*2.0f; |
|
|
|
|
gradI.x=(float)(std::abs(CV_MAT_ELEM(*out,data_type,kp+1,lm)-CV_MAT_ELEM(*out,data_type,kp,lm)))*2.0f; |
|
|
|
|
} |
|
|
|
|
} else { |
|
|
|
|
if (CV_MAT_ELEM(*f,uchar,k-1,l)!=INSIDE) { |
|
|
|
|
gradI.x=(float)(abs(CV_MAT_ELEM(*out,data_type,kp,lm)-CV_MAT_ELEM(*out,data_type,km-1,lm)))*2.0f; |
|
|
|
|
gradI.x=(float)(std::abs(CV_MAT_ELEM(*out,data_type,kp,lm)-CV_MAT_ELEM(*out,data_type,km-1,lm)))*2.0f; |
|
|
|
|
} else { |
|
|
|
|
gradI.x=0; |
|
|
|
|
} |
|
|
|
|
} |
|
|
|
|
if (CV_MAT_ELEM(*f,uchar,k,l+1)!=INSIDE) { |
|
|
|
|
if (CV_MAT_ELEM(*f,uchar,k,l-1)!=INSIDE) { |
|
|
|
|
gradI.y=(float)(abs(CV_MAT_ELEM(*out,data_type,km,lp+1)-CV_MAT_ELEM(*out,data_type,km,lm))+ |
|
|
|
|
abs(CV_MAT_ELEM(*out,data_type,km,lm)-CV_MAT_ELEM(*out,data_type,km,lm-1))); |
|
|
|
|
gradI.y=(float)(std::abs(CV_MAT_ELEM(*out,data_type,km,lp+1)-CV_MAT_ELEM(*out,data_type,km,lm))+ |
|
|
|
|
std::abs(CV_MAT_ELEM(*out,data_type,km,lm)-CV_MAT_ELEM(*out,data_type,km,lm-1))); |
|
|
|
|
} else { |
|
|
|
|
gradI.y=(float)(abs(CV_MAT_ELEM(*out,data_type,km,lp+1)-CV_MAT_ELEM(*out,data_type,km,lm)))*2.0f; |
|
|
|
|
gradI.y=(float)(std::abs(CV_MAT_ELEM(*out,data_type,km,lp+1)-CV_MAT_ELEM(*out,data_type,km,lm)))*2.0f; |
|
|
|
|
} |
|
|
|
|
} else { |
|
|
|
|
if (CV_MAT_ELEM(*f,uchar,k,l-1)!=INSIDE) { |
|
|
|
|
gradI.y=(float)(abs(CV_MAT_ELEM(*out,data_type,km,lm)-CV_MAT_ELEM(*out,data_type,km,lm-1)))*2.0f; |
|
|
|
|
gradI.y=(float)(std::abs(CV_MAT_ELEM(*out,data_type,km,lm)-CV_MAT_ELEM(*out,data_type,km,lm-1)))*2.0f; |
|
|
|
|
} else { |
|
|
|
|
gradI.y=0; |
|
|
|
|
} |
|
|
|
@ -748,11 +744,13 @@ cvInpaint( const CvArr* _input_img, const CvArr* _inpaint_mask, CvArr* _output_i |
|
|
|
|
if( !CV_ARE_SIZES_EQ(input_img,output_img) || !CV_ARE_SIZES_EQ(input_img,inpaint_mask)) |
|
|
|
|
CV_Error( CV_StsUnmatchedSizes, "All the input and output images must have the same size" ); |
|
|
|
|
|
|
|
|
|
if( (CV_MAT_CN(input_img->type) != 1 && |
|
|
|
|
if( (CV_MAT_TYPE(input_img->type) != CV_8U && |
|
|
|
|
CV_MAT_TYPE(input_img->type) != CV_16U && |
|
|
|
|
CV_MAT_TYPE(input_img->type) != CV_32F && |
|
|
|
|
CV_MAT_TYPE(input_img->type) != CV_8UC3) || |
|
|
|
|
!CV_ARE_TYPES_EQ(input_img,output_img) ) |
|
|
|
|
CV_Error( CV_StsUnsupportedFormat, |
|
|
|
|
"Any 1-channel and 8-bit 3-channel input/output images are supported" ); |
|
|
|
|
"8-bit, 16-bit unsigned or 32-bit float 1-channel and 8-bit 3-channel input/output images are supported" ); |
|
|
|
|
|
|
|
|
|
if( CV_MAT_TYPE(inpaint_mask->type) != CV_8UC1 ) |
|
|
|
|
CV_Error( CV_StsUnsupportedFormat, "The mask must be 8-bit 1-channel image" ); |
|
|
|
@ -802,10 +800,36 @@ cvInpaint( const CvArr* _input_img, const CvArr* _inpaint_mask, CvArr* _output_i |
|
|
|
|
cvSub(out,band,out,NULL); |
|
|
|
|
SET_BORDER1_C1(out,uchar,0); |
|
|
|
|
icvCalcFMM(out,t,Out,true); |
|
|
|
|
icvTeleaInpaintFMM(mask,t,output_img,range,Heap); |
|
|
|
|
switch(CV_MAT_DEPTH(output_img->type)) |
|
|
|
|
{ |
|
|
|
|
case CV_8U: |
|
|
|
|
icvTeleaInpaintFMM<uchar>(mask,t,output_img,range,Heap); |
|
|
|
|
break; |
|
|
|
|
case CV_16U: |
|
|
|
|
icvTeleaInpaintFMM<ushort>(mask,t,output_img,range,Heap); |
|
|
|
|
break; |
|
|
|
|
case CV_32F: |
|
|
|
|
icvTeleaInpaintFMM<float>(mask,t,output_img,range,Heap); |
|
|
|
|
break; |
|
|
|
|
default: |
|
|
|
|
CV_Error( cv::Error::StsBadArg, "Unsupportedformat of the input image" ); |
|
|
|
|
} |
|
|
|
|
} |
|
|
|
|
else if (flags == cv::INPAINT_NS) { |
|
|
|
|
icvNSInpaintFMM(mask,t,output_img,range,Heap); |
|
|
|
|
switch(CV_MAT_DEPTH(output_img->type)) |
|
|
|
|
{ |
|
|
|
|
case CV_8U: |
|
|
|
|
icvNSInpaintFMM<uchar>(mask,t,output_img,range,Heap); |
|
|
|
|
break; |
|
|
|
|
case CV_16U: |
|
|
|
|
icvNSInpaintFMM<ushort>(mask,t,output_img,range,Heap); |
|
|
|
|
break; |
|
|
|
|
case CV_32F: |
|
|
|
|
icvNSInpaintFMM<float>(mask,t,output_img,range,Heap); |
|
|
|
|
break; |
|
|
|
|
default: |
|
|
|
|
CV_Error( cv::Error::StsBadArg, "Unsupported format of the input image" ); |
|
|
|
|
} |
|
|
|
|
} else { |
|
|
|
|
CV_Error( cv::Error::StsBadArg, "The flags argument must be one of CV_INPAINT_TELEA or CV_INPAINT_NS" ); |
|
|
|
|
} |
|
|
|
@ -819,17 +843,6 @@ void cv::inpaint( InputArray _src, InputArray _mask, OutputArray _dst, |
|
|
|
|
Mat src = _src.getMat(), mask = _mask.getMat(); |
|
|
|
|
_dst.create( src.size(), src.type() ); |
|
|
|
|
Mat dst = _dst.getMat(); |
|
|
|
|
|
|
|
|
|
Mat src_temp; |
|
|
|
|
if (src.channels() == 1) { |
|
|
|
|
src.convertTo(src_temp, data_type_cv); |
|
|
|
|
} else |
|
|
|
|
src_temp = src.clone();
|
|
|
|
|
|
|
|
|
|
Mat dst_temp = Mat::zeros(src_temp.size(), src_temp.type()); |
|
|
|
|
CvMat c_src = src_temp, c_mask = mask, c_dst = dst_temp; |
|
|
|
|
CvMat c_src = src, c_mask = mask, c_dst = dst; |
|
|
|
|
cvInpaint( &c_src, &c_mask, &c_dst, inpaintRange, flags ); |
|
|
|
|
|
|
|
|
|
if (src.channels() == 1) dst_temp.convertTo(dst_temp, src.type()); |
|
|
|
|
dst_temp.copyTo(dst); |
|
|
|
|
} |
|
|
|
|