Add files via upload

Author: liweiwu

README.MD

@@ -0,0 +1,30 @@
+Integrate Python Code into Simulink for Simulation
+
+Using Python Human Detection Algorithm in Simulink Example
+Copyright 2012-2021 The MathWorks, Inc.
+
+This example shows users how to integrate Python Human detection code into Simulink for simulation. This example reads a pre-recorded video, applies the human detection algorithm defined in Python, and generates the output video with human marked. This example includes the following files:
+
+runme.m: This file provides the prep steps when using this example. Since the Python algorithm is based on OpenCV, you will need to install the needed packages if they are not available on your computer. This file also helps set the "out of Process" execution mode to avoid possible MATLAB crashes.   
+
+
+detectHuman.py: this file contains the Python algorithm using OpenCV Histogram of Oriented Gradients (HOG) for human detection. 
+
+livedata.mp4: this file is a pre-recorded video showing several people walking.  
+
+python_HumanDetector.m: this file contains the MATLAB System Object integrating the Python human detection algorithm.
+
+videoReader.m: this file contains the MATLAB System Object reading the pre-recorded video of livedata.mp4. This can be replaced if you have DSP System Toolbox available (see below). 
+
+base_python_example_21a.slx: this file shows how Simulink can integrate Python code for simulation by using either MATLAB Function block or MATLAB System block. You could use the manual switch to select either method to bring in Python code. When you run this file, the output video will show blue boxes drawn around detected people.  
+
+If you have DSP System Toolbox available, you could use the "From Multimedia File" block to replace the videoReader System block for simpler video reading; If you have Computer Vision Toolbox available, you can use the "To video Display" block to replace the Video_with_human_detection MATLAB function block for simpler video play.  
+
+To use this example, open and run the runme.m first in MATLAB; Then start the base_python_example_21a.slx and run the simulation.   
+
+Products needed for using this example :
+
+MATLAB; 
+Simulink; 
+DSP System Toolbox (optional); 
+Computer Vision Toolbox (optional); 

SECURITY.md

@@ -0,0 +1,6 @@
+# Reporting Security Vulnerabilities 
+
+If you believe you have discovered a security vulnerability, please report it to 
+[security@mathworks.com](mailto:security@mathworks.com). Please see 
+[MathWorks Vulnerability Disclosure Policy for Security Researchers](https://www.mathworks.com/company/aboutus/policies_statements/vulnerability-disclosure-policy.html) 
+for additional information.  

base_python_example_21a.slx

@@ -0,0 +1,53 @@
+import time
+import imutils
+import numpy as np
+from cv2 import cv2
+
+
+class hogObject:
+    detector = cv2.HOGDescriptor()
+    def __init__(self):
+        self.detector.setSVMDetector(cv2.HOGDescriptor_getDefaultPeopleDetector())
+
+
+def getHogObject():
+    return hogObject()
+
+
+def detectHumanFromFrame(image, hog):
+    image = np.asarray(image)
+    image = imutils.resize(image, width=min(400, image.shape[1]))
+
+    # Detecting all the regions in the Image that has a pedestrians inside it
+    (regions, _) = hog.detector.detectMultiScale(
+        image, winStride=(4, 4), padding=(4, 4), scale=1.05)
+
+    # Drawing the regions in the image
+    for (x, y, w, h) in regions:
+        cv2.rectangle(image, (x, y), (x + w, y + h), (0, 0, 255), 2)
+        
+    return  image 
+
+
+if __name__ == "__main__":
+
+    hog = hogObject()
+    
+    cap = cv2.VideoCapture('livedata.mp4')
+    ret, image = cap.read()
+
+    t_end = time.time() + 10
+    while time.time() < t_end:
+        # Reading the video stream
+        ret, image = cap.read()
+        if ret:
+            image = detectHumanFromFrame(image, hog)
+            cv2.imshow("Image", image)
+            cv2.waitKey(5)
+            if cv2.waitKey(5) & 0xFF == ord('q'):
+                exit()
+        else:
+            break
+
+    cap.release()
+    cv2.destroyAllWindows()

detectHuman.py

@@ -0,0 +1,11 @@
+Copyright (c) 2020, The MathWorks, Inc.
+All rights reserved.
+Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
+1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
+2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
+3. In all cases, the software is, and all modifications and derivatives of the software shall be, licensed to you solely for use in conjunction with MathWorks products and service offerings. 
+THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+
+
+

license.txt

@@ -0,0 +1,36 @@
+classdef python_HumanDetector < matlab.System
+    % Pre-computed constants
+    properties(Access = private, Nontunable)
+        hog
+    end
+
+    methods(Access = protected)
+        function setupImpl(obj)
+            % Perform one-time calculations, such as computing constants
+            obj.hog = py.detectHuman.getHogObject();
+        end
+
+        function y = stepImpl(obj,u)
+            % Calculate y as a function of input u and discrete states.
+            out = py.detectHuman.detectHumanFromFrame(u, obj.hog);
+            y = uint8(out);
+        end
+
+        function out = getOutputSizeImpl(obj)
+            out = [300 400 3];
+        end
+        
+        function y1 = getOutputDataTypeImpl(obj)
+            y1 = 'uint8';
+        end
+        
+        function y1 = isOutputComplexImpl(~)
+            y1 = false;
+        end
+        
+        function out = isOutputFixedSizeImpl(obj)
+            out = true;
+        end
+
+    end
+end

livedata.mp4

@@ -0,0 +1,25 @@
+%%%%%%%%%%%%%%%%%
+% Setting up python installer and getting packages 
+% 1. install pip
+% 2. install opencv
+% 3. install imutils
+%
+% 1. intall pip (on Linux console)
+% sudo apt install python-pip
+%
+% 2. install opencv
+% python -m pip install opencv-python
+%
+% 3. install imutils
+% python -m pip install imutils
+
+% MATLAB may crash with python process - use "OutOfProcess" execution mode
+% could avoid that
+pyenv("ExecutionMode","OutOfProcess");
+    
+py.importlib.import_module('detectHuman');
+
+pathToPyfun = fileparts(which('detectHuman.py'));
+if count(py.sys.path,pathToPyfun) == 0
+insert(py.sys.path,int32(0),pathToPyfun);
+end

python_HumanDetector.m

@@ -0,0 +1,41 @@
+classdef videoReader < matlab.System
+    % Pre-computed constants
+    properties(Access = private, Nontunable)
+        v
+    end
+
+    methods(Access = protected)
+        function setupImpl(obj)
+            % Perform one-time calculations, such as computing constants
+            obj.v = VideoReader('livedata.mp4');
+        end
+
+        function y = stepImpl(obj)
+            % Calculate y as a function of input u and discrete states.
+            if hasFrame(obj.v)
+                out = readFrame(obj.v);
+            else
+                out = zeros(480, 640, 3);
+            end
+
+            y = uint8(out);
+        end
+
+        function out = getOutputSizeImpl(obj)
+            out = [480 640 3];
+        end
+        
+        function y1 = getOutputDataTypeImpl(obj)
+            y1 = 'uint8';
+        end
+        
+        function y1 = isOutputComplexImpl(~)
+            y1 = false;
+        end
+        
+        function out = isOutputFixedSizeImpl(obj)
+            out = true;
+        end
+
+    end
+end