c++ - OpenCV Active Contour doesn't work well for Biomed data compared it's Matlab version -


i working on translating matlab version active contour opencv. seems opencv has own version of active contour deprecated then. curious how different between version matlab version. did comparison experiment. given mri brain image , initial boundary , applied both matlab , opencv version active contour segment white matter in brain.

for both matlab , opencv, choose gradient(edge) method of active contour. labeled results on image. red curve initial boundary while green curve final boundary.

the matlab code is: 

clear clc = imread('brain.png'); imasks = imread('shrinkmask.png');  figure(1) imshow(i) hold on visboundaries( imasks, 'color', 'r');  bwsnake = activecontour(i,  imasks, 250, 'edge','contractionbias',-0.45,'smoothfactor',0.7); finalboundary = bwboundaries(bwsnake); visboundaries(finalboundary, 'color', 'g');

enter image description here

it seems matlab job when choosing appropriate parameters. "contractionbias" parameter (refer here) useful controlling curve shrink or expand.

however, opencv version active contour seems doesn't work @ all. use slider bar change parameters alpha, beta , gamma cannot descend result. seems curve doesn't expand no matter how set parameters. use opencv 2.4.13.3 , active contour can found in legacy folder snakes.cpp. chek opencv github or more snakes.cpp. opencv test code is:

#include <iostream> #include "opencv2/core.hpp" #include "opencv2/highgui.hpp" #include "opencv2/imgproc.hpp" #include "opencv2/legacy/legacy.hpp"   using namespace std; using namespace cv;  iplimage *image = 0;  iplimage *image2 = 0;  iplimage *imagemask = 0; using namespace std;  int ialpha = 10; int ibeta = 40; int igamma = 50;   void onchange(int pos) {      if (image2) cvreleaseimage(&image2);     if (image) cvreleaseimage(&image);     if (imagemask) cvreleaseimage(&imagemask);      image2 = cvloadimage("brain.png", 1);      image = cvloadimage("brain.png", 0);     imagemask = cvloadimage("shrinkmask.png", 0);        cvmemstorage* storage = cvcreatememstorage(0);     cvseq* contours = 0;      cvfindcontours(imagemask, storage, &contours, sizeof(cvcontour),         cv_retr_external, cv_chain_approx_simple);       if (!contours) return;     int length = contours->total;     if (length<10) return;     cvpoint* point = new cvpoint[length];       cvseqreader reader;     cvpoint pt = cvpoint(0, 0);;     cvseq *contour2 = contours;      cvstartreadseq(contour2, &reader);     (int = 0; < length; i++)     {         cv_read_seq_elem(pt, reader);         point[i] = pt;     }     cvreleasememstorage(&storage);       (int = 0; i<length; i++)     {         int j = (i + 1) % length;         cvline(image2, point[i], point[j], cv_rgb(255, 0, 0), 1, 8, 0);     }      float alpha = ialpha / 100.0f;     float beta = ibeta / 100.0f;     float gamma = igamma / 100.0f;      cvsize size;     size.width = 3;     size.height = 3;     cvtermcriteria criteria;     criteria.type = cv_termcrit_iter;     criteria.max_iter = 500;     criteria.epsilon = 0.1;     cvsnakeimage(image, point, length, &alpha, &beta, &gamma, cv_value, size, criteria, 1);       (int = 0; i<length; i++)     {         int j = (i + 1) % length;         cvline(image2, point[i], point[j], cv_rgb(0, 255, 0), 1, 8, 0);     }     delete[]point;  }  int main(int argc, char* argv[]) {       cvnamedwindow("win1", 0);     cvcreatetrackbar("alpha", "win1", &ialpha, 100, onchange);     cvcreatetrackbar("beta", "win1", &ibeta, 100, onchange);     cvcreatetrackbar("gamma", "win1", &igamma, 100, onchange);     //cvresizewindow("win1", 696, 520);     cvresizewindow("win1", 256, 256);     onchange(0);      (;;)     {         if (cvwaitkey(40) == 27) break;         cvshowimage("win1", image2);     }      return 0; }

enter image description here

is there did wrong or version of active contour in opencv useless? checked snake.cpp, source code opencv version active contour, seems implemented based on theory in paper. because of file protection, cannot access matlab version active contour. therefore, cannot compare them directly. since eventually, use opencv implement stuffs, suggestion improving or changing opencv code appreciated.

i attached original image , mask in case wants repeat test.

thanks!

brain image mask


Comments

Post a Comment

Popular posts from this blog

android - InAppBilling registering BroadcastReceiver in AndroidManifest -

i'm implementing inapp billing in android app. works well, however, i'm trying decouple broadcast receiver activity manifest. suggestion in android's trivialdrive sample: // important: dynamically register broadcast messages updated purchases. // register receiver here instead of <receiver> in manifest // because call getpurchases() @ startup, therefore can ignore // broadcasts sent while app isn't running. // note: registering listener in activity bad idea, done here // because sample. regardless, receiver must registered after // iabhelper setup, before first call getpurchases(). currently there's class extends broadcastreceiver : public class iabbroadcastreceiver extends broadcastreceiver { /** * intent action receiver should filter for. */ public static final string action = "com.android.vending.billing.purchases_updated"; private final iabbroadcastlistener mlistener; public iabbroadcastreceiver(iabbroadcastlistener listener) { m...
Read more

python Tkinter Capturing keyboard events save as one single string -

is possible store keyboard events 1 string? code below stores 1 char , prints it. card reader or bar code reader, contains collection of character/ string not 1 character @ time. goal save char pressed text variable. from tkinter import * root = tk() def key(event): text= event.char text+= event.char print ("pressed", text) def callback(event): frame.focus_set() print ("clicked at", event.x, event.y) frame = frame(root, width=100, height=100) frame.bind("<key>", key) frame.bind("<button-1>", callback) frame.pack() root.mainloop() currently, creating text variable , everytime key() function gets called, , text stores last character typed. you can define text module level variable , , use module level text inside key function - text = '' def key(event): global text text+= event.char print("pressed", text)
Read more

sql server - Why does Linq-to-SQL add unnecessary COUNT()? -

here tables: create table customer ( cid int primary key, name varchar(32) ) create table ord ( oid int primary key, cid int foreign key references customer, address varchar(20) ) here's linq-to-sql statement: var aaa = c in db.customer select new { c, o = c.ord.tolist() }; and here's inexplicable query generated linq-to-sql: select [t0].[cid] [cid], [t0].[name] [name], [t1].[oid] [oid], [t1].[cid] [cid2], [t1].[address] [address], (select count(*) [dbo].[ord] [t2] [t2].[cid] = [t0].[cid]) [value] [dbo].[customer] [t0] left outer join [dbo].[ord] [t1] on [t1].[cid] = [t0].[cid] order [t0].[cid], [t1].[oid] i want understand how rid of count(*) part. it's entirely uncalled for! linq needs count(*) in order determine number of entries returned o = c.ord.tolist() .
Read more