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

Grayscale DICOM Bild in NxMx3-Matrix wandeln
 

tteichmann
Forum-Newbie

Forum-Newbie

Beiträge: 1
Anmeldedatum: 18.11.15
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 18.11.2015, 14:54     Titel: Grayscale DICOM Bild in NxMx3-Matrix wandeln
  Antworten mit Zitat      
Hi Leute Smile

Ich lese in meinem Code ein DICOM-Bild ein, um dieses weiter zu verarbeiten (Winkeldektektion, danach Kantenerkennung). Zur Verarbeitung nutze ich unter Anderem den rgb2gray-Befehl aus einem Skript aus Foren (notwendig, damit die Winkelkorrektur funktioniert), welcher mit Grayscale-DICOM-Bildern leider nicht funktioniert, weil diese nur eine NxM-Matrix beim Einlesen erzeugen (NxMx3-Matrix wird benötigt).

Hat jemand von euch eine Idee, wie ich die Konversion innerhalb von Matlab schaffen könnte, ohne das DICOM-Format zu verlassen?

Ideen meinerseits:
1) Grayscale DICOM in Truecolor DICOM wandeln
2) DICOM in anderes Bildformat abändern (bereits in Teilen probiert, ohne Erfolg)

Code:

Image1 = dicomread('SPI01130.dcm');

Image1invert = imcomplement(Image1);

Image1norm = double(Image1);
Image1norm = (Image1norm - min(Image1norm(:))) / (max(Image1norm(:)) - min(Image1norm(:))); %// Normalisieren 0-1
Image1norm = uint16(255*Image1norm); %// Skalieren bis 255

d = rgb2gray(image);

IndXVec=[];
IndYVec=[];
 
[n,m]=size(d);
 
% Minimum jeder Spalte bestimmen
for i=1:m
    IndY=max(find(d(:,i)==0));
    if(~isempty(IndY))
        IndXVec(end+1)=i;
        IndYVec(end+1)=n-IndY+1;        
    end
end
 
% Wiederanstieg filtern
IndXVec=IndXVec(diff(IndYVec)<1);
IndYVec=IndYVec(diff(IndYVec)<1);
 
% Lineare Regression der Kante
p=polyfit(IndXVec,IndYVec,1);
 
% Winkel berechnen
alpha=atan(p(1))*180/pi;

% Rotationskorrektur
image_rotiert = imrotate(d,-alpha_rund,'bilinear','crop');


 


Grüße
Private Nachricht senden Benutzer-Profile anzeigen


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.