Skalierungsproblem contour

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

benjamin.spunkt

Forum-Newbie


	Beiträge: 9

	Anmeldedatum: 17.06.14

	Wohnort: Berlin

	Version: R2013b

Verfasst am: 18.11.2014, 15:18 Titel: Skalierungsproblem contour

Hallo liebes Forum,
ich möchte einen contour plot einer matrix erstellen.
Leider kommt es manchmal vor, dass in den Daten ein (relativ zu den restlichen Werten) sehr hoher Wert enthalten ist. Ist dies der Fall so stimmt die Skalierung nicht mehr, bzw. werden einem Wert falsche Farben zugeordnet:

Code:

clear;

%zufällige matrix erstellen
x =rand(6,7)*400;

%nun schreiben wir an irgend eine Position der Matrix einen (vergleichsweiße) hohen Wert:
x(1,1)=10000;

n = [1 2 3 4 5 6 7];
dz = [1 2 3 4 5 6];

%Skalierung Z Achse
ZAchseuntereGrenze = 0;
ZAchseobereGrenze = 300;

figure4 = figure('PaperType','a4letter',...
'PaperSize',[20.98404194812 29.67743169791],...
'Name','test',...
'Colormap',[1 0 0;1 0.0266666673123837 0;1 0.0533333346247673 0;1 0.0799999982118607 0;1 0.106666669249535 0;1 0.133333340287209 0;1 0.159999996423721 0;1 0.186666667461395 0;1 0.213333338499069 0;1 0.240000009536743 0;1 0.266666680574417 0;1 0.293333351612091 0;1 0.319999992847443 0;1 0.346666663885117 0;1 0.373333334922791 0;1 0.400000005960464 0;1 0.426666676998138 0;1 0.453333348035812 0;1 0.480000019073486 0;1 0.506666660308838 0;1 0.533333361148834 0;1 0.560000002384186 0;1 0.586666703224182 0;1 0.613333344459534 0;1 0.639999985694885 0;1 0.666666686534882 0;1 0.693333327770233 0;1 0.720000028610229 0;1 0.746666669845581 0;1 0.773333370685577 0;1 0.800000011920929 0;0.969696998596191 0.790849685668945 0;0.939393937587738 0.781699359416962 0;0.909090936183929 0.772549033164978 0;0.878787875175476 0.763398706912994 0;0.848484873771667 0.754248380661011 0;0.818181812763214 0.745098054409027 0;0.787878811359406 0.735947728157043 0;0.757575750350952 0.72679740190506 0;0.727272748947144 0.717647075653076 0;0.69696968793869 0.708496749401093 0;0.666666686534882 0.699346423149109 0;0.636363625526428 0.690196096897125 0;0.60606062412262 0.681045770645142 0;0.575757563114166 0.671895444393158 0;0.545454561710358 0.662745118141174 0;0.515151500701904 0.653594791889191 0;0.484848499298096 0.644444465637207 0;0.454545468091965 0.635294139385223 0;0.424242436885834 0.62614381313324 0;0.393939405679703 0.616993486881256 0;0.363636374473572 0.607843160629272 0;0.333333343267441 0.598692834377289 0;0.30303031206131 0.589542508125305 0;0.272727280855179 0.580392181873322 0;0.242424249649048 0.571241855621338 0;0.212121218442917 0.562091529369354 0;0.181818187236786 0.552941203117371 0;0.151515156030655 0.543790876865387 0;0.121212124824524 0.534640550613403 0;0.0909090936183929 0.52549022436142 0;0.060606062412262 0.516339898109436 0;0.030303031206131 0.507189571857452 0;0 0.498039215803146 0],...
'Color',[1 1 1]);
% Create axes
axes1 = axes('Parent',figure4,'YGrid','on','XGrid','on',...
'Position',[0.13 0.11 0.749699248120301 0.815],...
'Layer','top',...
'CLim',[ZAchseuntereGrenze ZAchseobereGrenze], 'CLimMode','manual');

box(axes1,'on');
hold(axes1,'all');

% Create contour
contour(dz,n,x','LineColor',[0 0 0],'Fill','on','Parent',axes1);

% Create xlabel
xlabel('Druck [bar(rel)]');

% Create ylabel
ylabel('Drehzahl [1/min]');

% Create title
title('test');

% Create colorbar
colorbar('peer',axes1,...
[0.885964912280702 0.250602157829869 0.0639097744360902 0.362043795620438]);

% Create textbox
annotation(figure4,'textbox',...
[0.896196800964674 0.641692826133855 0.04625 0.0733322830786442],...
'String',{'Verstell-','geschw.','[°/s]'},...
'HorizontalAlignment','center',...
'FitBoxToText','off',...
'LineStyle','none');

Funktion ohne Link?

Das Problem lässt sich teilweise lösen indem mehr Klassen (z.B. 100) einführt werden, also:

Code:

contour(dz,n,x',100,'LineColor',[0 0 0],'Fill','on','Parent',axes1);

Funktion ohne Link?

Dies ist aber unelegant. Hat jemand eine Idee wie die Skalierung auch so immer den richtigen Farbwert anzeigt? Die Klassen sollen sich lediglich innerhalb der Limits befinden.
Würde mich über Anregungen sehr freuen.

Viele Grüße,
Benjamin

Seban

Forum-Meister


	Beiträge: 600

	Anmeldedatum: 18.01.12

	Wohnort: ---

	Version: ab R2014b

Verfasst am: 19.11.2014, 01:00 Titel:

Hallo Benjamin,

Ist es eine Option, die zu großen Werte zu verändern? zB könnte man sie auf den Durchschnittswert des x-Arrays setzen oder auf den Durchschnitt der umliegenden Werte.

Als Schwelle habe ich beliebig 1000 gewählt.

Code:

%zufällige matrix erstellen
rng('default'); % so werden stets die selben random Werte erzeugt
x =rand(6,7)*400;

%nun schreiben wir an irgend eine Position der Matrix einen (vergleichsweiße) hohen Wert:
x(1,1)=10000;

x(x>1000)= mean(mean(x)); % zu große Werte ersetzen

n = [1 2 3 4 5 6 7];
dz = [1 2 3 4 5 6];

%Skalierung Z Achse
ZAchseuntereGrenze = 0;
ZAchseobereGrenze = 300;

figure4 = figure('PaperType','a4letter',...
'PaperSize',[20.98404194812 29.67743169791],...
'Name','test',...
'Colormap',[1 0 0;1 0.0266666673123837 0;1 0.0533333346247673 0;1 0.0799999982118607 0;1 0.106666669249535 0;1 0.133333340287209 0;1 0.159999996423721 0;1 0.186666667461395 0;1 0.213333338499069 0;1 0.240000009536743 0;1 0.266666680574417 0;1 0.293333351612091 0;1 0.319999992847443 0;1 0.346666663885117 0;1 0.373333334922791 0;1 0.400000005960464 0;1 0.426666676998138 0;1 0.453333348035812 0;1 0.480000019073486 0;1 0.506666660308838 0;1 0.533333361148834 0;1 0.560000002384186 0;1 0.586666703224182 0;1 0.613333344459534 0;1 0.639999985694885 0;1 0.666666686534882 0;1 0.693333327770233 0;1 0.720000028610229 0;1 0.746666669845581 0;1 0.773333370685577 0;1 0.800000011920929 0;0.969696998596191 0.790849685668945 0;0.939393937587738 0.781699359416962 0;0.909090936183929 0.772549033164978 0;0.878787875175476 0.763398706912994 0;0.848484873771667 0.754248380661011 0;0.818181812763214 0.745098054409027 0;0.787878811359406 0.735947728157043 0;0.757575750350952 0.72679740190506 0;0.727272748947144 0.717647075653076 0;0.69696968793869 0.708496749401093 0;0.666666686534882 0.699346423149109 0;0.636363625526428 0.690196096897125 0;0.60606062412262 0.681045770645142 0;0.575757563114166 0.671895444393158 0;0.545454561710358 0.662745118141174 0;0.515151500701904 0.653594791889191 0;0.484848499298096 0.644444465637207 0;0.454545468091965 0.635294139385223 0;0.424242436885834 0.62614381313324 0;0.393939405679703 0.616993486881256 0;0.363636374473572 0.607843160629272 0;0.333333343267441 0.598692834377289 0;0.30303031206131 0.589542508125305 0;0.272727280855179 0.580392181873322 0;0.242424249649048 0.571241855621338 0;0.212121218442917 0.562091529369354 0;0.181818187236786 0.552941203117371 0;0.151515156030655 0.543790876865387 0;0.121212124824524 0.534640550613403 0;0.0909090936183929 0.52549022436142 0;0.060606062412262 0.516339898109436 0;0.030303031206131 0.507189571857452 0;0 0.498039215803146 0],...
'Color',[1 1 1]);
% Create axes
axes1 = axes('Parent',figure4,'YGrid','on','XGrid','on',...
'Position',[0.13 0.11 0.749699248120301 0.815],...
'Layer','top',...
'CLim',[ZAchseuntereGrenze ZAchseobereGrenze], 'CLimMode','manual');

box(axes1,'on');
hold(axes1,'all');

% Create contour
contour(dz,n,x','LineColor',[0 0 0],'Fill','on','Parent',axes1);

% ...

Funktion ohne Link?

Grüße,
Seban
_________________

Richtig fragen
Debugging

benjamin.spunkt

Themenstarter

Forum-Newbie


	Beiträge: 9

	Anmeldedatum: 17.06.14

	Wohnort: Berlin

	Version: R2013b

Verfasst am: 19.11.2014, 14:37 Titel:

Hallo Seban,

vielen Dank für Deine Idee. Ich habe es nun wie folgt umgangen:

Code:

x(x>1800)= NaN; %zu hohe (unrealistische) Werte mit NaN überschreiben
x_range = 1:1:length(n); %X Vektor
y_range = 1:1:length(dz); %Y Vektor
valid = ~isnan(x); %Suche NaN
[X,Y] = meshgrid(x_range, y_range);
x = griddata(X(valid),Y(valid),x(valid),X,Y); %grid interpolieren

Funktion ohne Link?

D.h. der falsche Wert wird geköscht und interpoliert. Sollte der falsche Wert jedoch an den Ecken der Matrix auftreten so funktioniert dies nicht, das ist aber in diesem Fall aus anderen Gründen nicht Möglich.
Ich denke ich kann mit diesem Workaround gut leben...eleganter wäre natürlich eine wahre Skalierung gewesen.

Beste Grüße,
Benjamin


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 |

RSS

Hosted by:

Copyright © 2007 - 2025 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.