{"id":4058,"date":"2014-08-19T18:00:21","date_gmt":"2014-08-19T23:00:21","guid":{"rendered":"https:\/\/blogs.mathworks.com\/seth\/?p=4058"},"modified":"2014-08-15T18:39:32","modified_gmt":"2014-08-15T23:39:32","slug":"a-matlab-vs-simulink-battle-im-in","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/simulink\/2014\/08\/19\/a-matlab-vs-simulink-battle-im-in\/","title":{"rendered":"A MATLAB vs. Simulink battle? I’m in!"},"content":{"rendered":"

When I saw the latest post by Matt Tearle<\/a> on Steve's blog<\/a>, I had to do something.<\/p>\r\n\r\n

Since the video at the center of this post is about a rap battle between MATLAB and Simulink, I thought it would be interesting to bring that to a higher level and see how the video processing could have been implemented in Simulink instead of MATLAB.<\/p>\r\n\r\n

Seriously... who wants to write code when you can simply connect blocks?<\/p>\r\n\r\n

Overview<\/strong><\/p>\r\n\r\n

Since I have absolutely zero knowledge in image processing, I simply re-implemented the algorithm suggested by Matt in Simulink. The resulting model looks like:<\/p>\r\n\r\n

\"Green<\/p>\r\n\r\n

Let's look at the details.<\/p>\r\n\r\n

Importing the video<\/strong><\/p>\r\n\r\n

To begin, I used the From Multimedia File<\/a> from the DSP System Toolbox<\/a> to directly read the AVI-file with our rappers in front of a green background. Then I use the resize<\/a> block from the Computer Vision System Toolbox<\/a>. This allows me to work on images of a relatively small size. After, I normalize the data, as suggested by Matt.<\/p> \r\n\r\n

\"Reading<\/p>\r\n\r\n

The Greenness factor<\/strong><\/p>\r\n\r\n

I really like the simple equation the Matt figured out to calculate the \"greenness\" of each pixel in the image. In MATLAB, it looks like:<\/p>\r\n\r\n

% Greenness = G*(G-R)*(G-B)<\/span>\r\ngreenness = yd(:,:,2).*(yd(:,:,2)-yd(:,:,1)).*(yd(:,:,2)-yd(:,:,3));<\/code><\/p>\r\n\r\nIn Simulink, I can use Selector<\/a> blocks to separate the three colors (third dimension of the image signal) and implement Matt's equation:\r\n\r\n

\"greenness\"<\/p>\r\n\r\n

Thresholding<\/strong><\/p>\r\n\r\n

The next piece of code is something that many users struggle to figure out:<\/p>\r\n\r\n

thresh = 0.3*mean(greenness(greenness>0));<\/code><\/p>\r\n\r\n

To implement logical indexing (e.g. A(A>0)<\/tt>) in Simulink, we need to combine the Find<\/a> block with a Selector<\/a> block. This results in a variable-size signal<\/a><\/p>\r\n\r\n

\"thresholding\"<\/p>\r\n\r\n\r\n

Combining the images<\/strong><\/p>\r\n\r\n

The last line we need to implement is to replace the pixels identified as green in the foreground image but their values in the background image. In MATLAB, this is accomplished using this line:<\/p>\r\n\r\n

rgb2(isgreen) = rgb1(isgreen);<\/code><\/p>\r\n\r\n

In Simulink, we need to extend the logic seen at the previous step. We first select the desired pixels from the background image, and then assign them to the foreground image:<\/p>\r\n\r\n

\"combining\"<\/p>\r\n\r\n

Note that the Computer Vision System Toolbox<\/a> has a Compositing<\/a> block implementing this functionality.<\/p>\r\n\r\n

Here is what the final result looks like:<\/p>\r\n\r\n

\"battle\"<\/p>\r\n\r\n\r\n

Now it's your turn<\/strong><\/p>\r\n\r\n

What do you think? Is it more convenient to implement this algorithm in MATLAB or Simulink? How do you typically decide between a code-based<\/em> or a block-based<\/em> implementation? Let us know by leaving a comment here<\/a>.<\/p>","protected":false},"excerpt":{"rendered":"

<\/div>

When I saw the latest post by Matt Tearle on Steve's blog, I had to do something.\r\n\r\nSince the video at the center of this post is about a rap battle between MATLAB and Simulink, I thought it would... read more >><\/a><\/p>","protected":false},"author":41,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":[],"categories":[50,71],"tags":[402,194],"_links":{"self":[{"href":"https:\/\/blogs.mathworks.com\/simulink\/wp-json\/wp\/v2\/posts\/4058"}],"collection":[{"href":"https:\/\/blogs.mathworks.com\/simulink\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blogs.mathworks.com\/simulink\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.mathworks.com\/simulink\/wp-json\/wp\/v2\/users\/41"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.mathworks.com\/simulink\/wp-json\/wp\/v2\/comments?post=4058"}],"version-history":[{"count":15,"href":"https:\/\/blogs.mathworks.com\/simulink\/wp-json\/wp\/v2\/posts\/4058\/revisions"}],"predecessor-version":[{"id":4083,"href":"https:\/\/blogs.mathworks.com\/simulink\/wp-json\/wp\/v2\/posts\/4058\/revisions\/4083"}],"wp:attachment":[{"href":"https:\/\/blogs.mathworks.com\/simulink\/wp-json\/wp\/v2\/media?parent=4058"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blogs.mathworks.com\/simulink\/wp-json\/wp\/v2\/categories?post=4058"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blogs.mathworks.com\/simulink\/wp-json\/wp\/v2\/tags?post=4058"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}