Я проектирую морфологические операции в OpenCV. Я пытаюсь имитировать функции remove
а также bridge
в матлабе Для этого я сослался на определение функции bwmorph.m, там я получил массивы таблицы поиска для remove и bridge.
После этого шага процедура одинакова для Matlab и OpenCV.
lut(img,lutarray,img)
Проблема заключается в том, что Matlab использует схему просмотра таблицы из 512 элементов (9 бит), в то время как OpenCV использует схему поиска из 256 элементов (8 бит). Как использовать lutarray Matlab в OpenCV?
После некоторых исследований я наткнулся эта почта.
Что подразумевается под человеком, когда он говорит, что «разделяет» изображение с 0-512, а затем на две части?
Является ли вышеуказанный метод даже правильным? Есть ли альтернативы этому?
bwlookup(bw,lut)
http://se.mathworks.com/help/images/ref/bwlookup.html
или внутренне, applylut
оба выполняют операцию соседства 2 на 2 или 3 на 3 на двоичном & белый) изображение, тогда как OpenCV cv::LUT
выполняет преобразование уровня серого на пиксель (тесно связано с intlut
в MATLAB). Примером последнего является выполнение гамма-коррекции на изображении уровня серого.
//! transforms array of numbers using a lookup table: dst(i)=lut(src(i))
CV_EXPORTS_W void LUT(InputArray src, InputArray lut, OutputArray dst,
int interpolation=0);
Насколько мне известно, нет соседства bwlookup
реализация в OpenCV. Однако, следуя описанию MATLAB’s bwlookup
Вы можете написать это сами.
// performs 3-by-3 lookup on binary image
void bwlookup(
const cv::Mat & in,
cv::Mat & out,
const cv::Mat & lut,
int bordertype=cv::BORDER_CONSTANT,
cv::Scalar px = cv::Scalar(0) )
{
if ( in.type() != CV_8UC1 )
CV_Error(CV_StsError, "er");
if ( lut.type() != CV_8UC1 || lut.rows*lut.cols!=512 || !lut.isContinuous() )
CV_Error(CV_StsError, "lut size != 512" );
if ( out.type() != in.type() || out.size() != in.size() )
out = cv::Mat( in.size(), in.type() );
const unsigned char * _lut = lut.data;
cv::Mat t;
cv::copyMakeBorder( in,t,1,1,1,1,bordertype,px);
const int rows=in.rows+1;
const int cols=in.cols+1;
for ( int y=1;y<rows;++y)
{
for ( int x=1;x<cols;++x)
{
int L = 0;
const int jmax=y+1;
#if 0 // row-major order
for ( int j=y-1, k=1; j<=jmax; ++j, k<<=3 )
{
const unsigned char * p = t.ptr<unsigned char>(j) + x-1;
for ( unsigned int u=0;u<3;++u )
{
if ( p[u] )
L += (k<<u);
#else // column-major order (MATLAB)
for ( int j=y-1, k=1; j<=jmax; ++j, k<<=1 )
{
const unsigned char * p = t.ptr<unsigned char>(j) + x-1;
for ( unsigned int u=0;u<3;++u )
{
if ( p[u] )
L += (k<<3*u);
#endif
}
}
out.at<unsigned char>(y-1,x-1)=_lut[ L ];
}
}
}
Я проверил это против remove
а также bridge
так должно работать. Надеюсь, это поможет.
Изменить: После проверки по случайной таблице поиска,
lut = uint8( rand(512,1)>0.5 ); % @MATLAB
B = bwlookup( A, lut );
Я изменил порядок отображения индексов в таблице поиска (не имеет значения, является ли операция симметричной).
Других решений пока нет …