WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Meerschweinchen Tracking

 

paul85
Forum-Newbie

Forum-Newbie


Beiträge: 4
Anmeldedatum: 27.06.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 27.06.2009, 16:13     Titel: Meerschweinchen Tracking
  Antworten mit Zitat      
Hallo,
ich soll ein Programm schreiben, dass Meerschweinchen trackt. leider habe ich keinen ansatz wie ich das machen soll.
meine matlab kenntnisse haben sich bis jetzt nur auf kantendetektion beschränkt.
kann mir jmd. erklären wie ich das machen soll?
das wäre echt super.


Ich denke, dass man die Tiere vom Hintergrund trennen sollte. Und dann mit einer for schleife jede Veränderung des Hintergrundes erfassen soll, um dann die Veränderung zu tracken. aber ich weiss nicht wie ich den Code niederschreiben soll.

video upload hat leider zu lange gedauert. ich hoffe, dass man sich das auch so vorstellen kann.

dankbar für jede hilfe.
Private Nachricht senden Benutzer-Profile anzeigen


paul85
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 4
Anmeldedatum: 27.06.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 27.06.2009, 16:25     Titel:
  Antworten mit Zitat      
Code:
meerschw = mmreader('Hausmeerschweinchen_1.mpg');
get(meerschw)
implay('Hausmeerschweinchen_1.mpg');

schwarzerMeerschweinchenWert = 50;
darkCar = rgb2gray(read(meerschw,25)); %zahl ist ich glaube die framerate
keinSchwarzesMeerschweinchen = imextendedmax(darkCar, schwarzerMeerschweinchenWert);
imshow(darkCar)
figure, imshow(keinSchwarzesMeerschweinchen)

%alles weiße wegmachen was nicht zu den meerschweinchen gehört
sedisk = strel('disk',8);
noSmallStructures = imopen(keinSchwarzesMeerschweinchen, sedisk);
imshow(noSmallStructures)


nframes = get(meerschw, 'NumberOfFrames');
I = read(meerschw, 71);
taggedCars = zeros([size(I,1) size(I,2)], class(I));

for k = 1 : nframes
    singleFrame = read(meerschw, k);

    % Convert to grayscale to do morphological processing.
    I = rgb2gray(singleFrame);

    % Remove dark cars.
    keinSchwarzesMeerschweinchen = imextendedmax(I, schwarzerMeerschweinchenWert);

    % Remove lane markings and other non-disk shaped structures.
    noSmallStructures = imopen(keinSchwarzesMeerschweinchen, sedisk);

    % Remove small structures.
    noSmallStructures = bwareaopen(noSmallStructures, 150);

 
    L = bwlabel(noSmallStructures);
    taggedmeerschweinchen(:,:,:,k) = singleFrame;
    if any(L(:))
       stats = regionprops(L, {'centroid','area'});
       areaArray = [stats.Area];
       [junk,idx] = max(areaArray);
       c = stats(idx).Centroid;
       c = floor(fliplr(c));
       width = 2;
       row = c(1)-width:c(1)+width;
       col = c(2)-width:c(2)+width;
       taggedmeerschweinchen(row,col,k) = 255;
       taggedmeerschweinchen(row,col,k) = 0;
       taggedmeerschweinchen(row,col,k) = 0;
    end
end
 


das ist mein code. der ist aus einem beispiel übernommen. das passt nicht ganz. habe es versucht zu ändern, aber es kommt immer der eine fehler:
??? Assignment has more non-singleton rhs dimensions than non-singleton
subscripts

Error in ==> meerschweinchen_track at 41
taggedCars(:,:,:,k) = singleFrame;

--> den verstehe ich nicht so ganz...
Private Nachricht senden Benutzer-Profile anzeigen
 
paul85
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 4
Anmeldedatum: 27.06.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 29.06.2009, 10:38     Titel:
  Antworten mit Zitat      
keiner ne idee?
Private Nachricht senden Benutzer-Profile anzeigen
 
paul85
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 4
Anmeldedatum: 27.06.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 29.06.2009, 13:52     Titel:
  Antworten mit Zitat      
so habe jetzt eine neue frage. kann mir jemand sagen, wie ich etwas geplottetes in ein Bild direkt einbinden und speichern kann bzw. in ein Videoframe?

hoffe jmd. antwortet...


Code:

function drawBox( rprops )
%DRAWBOX Summary of this function goes here
%   Detailed explanation goes here

centroids = cat(1, rprops.Centroid);
csize = size(centroids);

if (csize(1) == 4)
plot(centroids(1,1), centroids(1,2), 'y*');
plot(rectangle('Position', rprops(1).BoundingBox, 'EdgeColor', 'y'));
plot(centroids(2,1), centroids(2,2), 'g*');
plot(rectangle('Position', rprops(2).BoundingBox, 'EdgeColor', 'g'));
plot(centroids(3,1), centroids(3,2), 'r*');
plot(rectangle('Position', rprops(3).BoundingBox, 'EdgeColor', 'r'));
plot(centroids(4,1), centroids(4,2), 'b*');
plot(rectangle('Position', rprops(4).BoundingBox, 'EdgeColor', 'b'));
end;

if (csize(1) == 3)
plot(centroids(1,1), centroids(1,2), 'y*');
plot(rectangle('Position', rprops(1).BoundingBox, 'EdgeColor', 'y'));
plot(centroids(2,1), centroids(2,2), 'g*');
plot(rectangle('Position', rprops(2).BoundingBox, 'EdgeColor', 'g'));
plot(centroids(3,1), centroids(3,2), 'r*');
plot(rectangle('Position', rprops(3).BoundingBox, 'EdgeColor', 'r'));
end;

if (csize(1) == 2)
plot(centroids(1,1), centroids(1,2), 'y*');
plot(rectangle('Position', rprops(1).BoundingBox, 'EdgeColor', 'y'));
plot(centroids(2,1), centroids(2,2), 'g*');
plot(rectangle('Position', rprops(2).BoundingBox, 'EdgeColor', 'g'));
end;

if (csize(1) == 1)
plot(centroids(1,1), centroids(1,2), 'y*');
plot(rectangle('Position', rprops(1).BoundingBox, 'EdgeColor', 'y'));
end;

end


function [ binImage ] = processImage( rgbImage, grayBG )
% PROCESSIMAGE Summary of this function goes here
%   Detailed explanation goes here


grayImage = rgb2gray(rgbImage);
grayImage = threshDiff(grayImage, grayBG, 40);
grayImage = imfill(grayImage, 'holes');

level = graythresh(grayImage);
binImage = im2bw(grayImage,level);
binImage = bwareaopen(binImage, 15);

binImage = imfill(binImage, 'holes');
se = strel('disk',15);
binImage = imclose(binImage, se);

end


function [ msDiff ] = threshDiff( grayImage1, grayImage2, threshold )
%THRESHDIFF Summary of this function goes here
%   Detailed explanation goes here

msDiff1 = abs(grayImage1 - grayImage2);
msDiff2 = abs(grayImage2 - grayImage1);
msDiff = (msDiff1 + msDiff2 - threshold) * 3;

end


msvid = mmreader('msvid.mpg');
nframes = get(msvid, 'NumberOfFrames');
msBG = rgb2gray(imread('meersauBG.png'));
ms0 = read(msvid, 1);

for k = 1 : 100

   singleFrame = read(msvid, k);
   msBin = processImage(singleFrame, msBG);
   L = bwlabel(msBin, 4);
   s  = regionprops(L, 'centroid', 'BoundingBox', 'ConvexImage');
   imshow(singleFrame);

   hold on
       drawBox(s);
   hold off

end;
Private Nachricht senden Benutzer-Profile anzeigen
 
gast

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 01.07.2009, 22:59     Titel:
  Antworten mit Zitat      
hörst du CV an der WWU? Wink
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
.





 Impressum  | Nutzungsbedingungen  | Datenschutz | FAQ | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.