c++ - How to use region-growing algorithms to define a region of interest? -


i working on dicom images (ct scans) & isolate structures of interest in picture such human organs (like aorta, cf image enclosed). coding in c++ of itk & vtk.

the bright tube describes sick aorta , corresponds structure isolate.

let's assume these organs have particular brightness intensity, therefore can automatically identify them using region-growing algorithm (code below). in order so, computed threshold values based on mean & standard deviation values of voxels belonging organ.

how can keep aorta in image of itk/vtk features? guess i'm looking filter exact opposite of itk mask image filter.

please find (pseudo) code corresponding organ isolation below. computed 5 voxels dilation on result of region-growing sure include voxels of organ , have sufficient margin around organ after cropping.

typedef short inputpixeltype; typedef unsigned char outputpixeltype;  const int dimension = 3;  typedef itk::image< inputpixeltype, dimension > inputimagetype; typedef itk::image< outputpixeltype, dimension > outputimagetype;  // region growing typedef itk::connectedthresholdimagefilter< inputimagetype,  outputimagetype > connectedfiltertype;  connectedfiltertype::pointer connectedthreshold = connectedfiltertype::new();  connectedthreshold->setinput(input); connectedthreshold->setupper(upperthreshold); connectedthreshold->setlower(lowerthreshold);  //initializing seed internalimagetype::indextype index; index[0] = seed_x;  index[1] = seed_y; connectedthreshold->setseed(index);  // dilate resulting region-growing of 5 voxels safety typedef itk::binaryballstructuringelement< outputimagetype,  dimension > structuringelementtype; typedef itk::binarydilateimagefilter< outputimagetype,  outputimagetype, struturingelementtype > dilatefiltertype;  structuringelementtype structuringelement; structuringelement.setradius(5); structuringelement.createstructuringelement();  dilatefiltertype::pointer dilatefilter = dilatefiltertype::new(); dilatefilter->setinput(connectedthreshold->getoutput()); dilatefilter->setkernel(structuringelement);  // saving results of rg+dilation typedef itk::imagefilewriter< outputimagetype > writertype; writertype::pointer writer = writertype::new(); writer->setinput(dilatefilter->getoutput()); writer->setfilename("organ-segmented-with-dilation.mhd"); try {     writer->update(); } catch(itk::exceptionobject& err) {     std::cerr << "exception caught! " << err.what() << std::endl;     return exit_failure; }  // next crop input image region-growing?

any or remark welcomed.

mask filter can opposite of does. default, masking value 0, , outside value. means parts of image correspond non-zero part of mask kept, , rest zeroed out. if not want, can invert logic setting different masking , outside values.


Comments

Post a Comment

Popular posts from this blog

Is there a better way to structure post methods in Class Based Views -

i'm using class based views in django , never specify url in ajax calls , omit action parameter in form because know these requests go through post method of class based view. notice view - profile page make many requests depending on user updating on page - post method become convoluted , if/else mess example: def post(self, request, *args, **kwargs): # handle user changing available date if request.post.get('availabledate') != none: self.updateavailabledate(request) return jsonresponse({'result':'success'}) if request.post.get('newprojectname') != none: creative_user = creativeuserprofile.objects.get(id = request.user.id) project = project.create_project(creative_user, request.post['newprojectname']) project.save() return jsonresponse({'projectid': project.id }) if request.post.get('projectimage') != none , request.files['file']: proj...
Read more

performance - Why is XCHG reg, reg a 3 micro-op instruction on modern Intel architectures? -

i'm doing micro-optimization on performance critical part of code , came across sequence of instructions (in at&t syntax): add %rax, %rbx mov %rdx, %rax mov %rbx, %rdx i thought had use case xchg allow me shave instruction , write: add %rbx, %rax xchg %rax, %rdx however, dimay found agner fog's instruction tables , xchg 3 micro-op instruction 2 cycle latency on sandy bridge, ivy bridge, broadwell, haswell , skylake. 3 whole micro-ops , 2 cycles of latency! 3 micro-ops throws off 4-1-1-1 cadence , 2 cycle latency makes worse original in best case since last 2 instructions in original might execute in parallel. now... cpu might breaking instruction micro-ops equivalent to: mov %rax, %tmp mov %rdx, %rax mov %tmp, %rdx where tmp anonymous internal register , suppose last 2 micro-ops run in parallel latency 2 cycles. given register renaming occurs on these micro-architectures, though, doesn't make sense me done way. why wouldn't register renamer ...
Read more

c# - Asp.net web api : redirect unauthorized requst to forbidden page -

im trying redirect unauthorized request forbidden page instead i'm getting forbidden page in response body , how can fix ? here's startup class : app.createperowincontext(storecontext.create); app.createperowincontext<applicationusermanager>(applicationusermanager.create); app.createperowincontext<applicationsigninmanager>(applicationsigninmanager.create); app.usecookieauthentication(new cookieauthenticationoptions { authenticationtype = defaultauthenticationtypes.applicationcookie, expiretimespan = timespan.fromdays(30), }); app.useexternalsignincookie(defaultauthenticationtypes.externalcookie); app.usetwofactorsignincookie(defaultauthenticationtypes.twofactorcookie, timespan.fromminutes(5)); app.usetwofactorrememberbrowsercookie(defaultauthenticationtypes.twofactorrememberbrowsercookie); the method im trying reach : [httpget] [authorize(roles = "admin")] public string getcurrentusername() { return userman...
Read more