Optimierung: was mach ich falsch?

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

MatLabNooB

Forum-Guru


	Beiträge: 262

	Anmeldedatum: 27.03.09

	Wohnort: ---

	Version: ---

Verfasst am: 21.11.2010, 05:26 Titel: Optimierung: was mach ich falsch?

Hallo,

hab jetzt schon ganzen Abend rumgetestet komme aber leider nicht zum Ziel, und zweifle langsam ob ich was falsch mach und die Optimierung so wie ich es versuch garnicht möglich ist?

Ich versuch quasi 2 bilder aus denen ich punktkoordinaten besitzte zu kalibrieren. Daher würd ich gern die Verzerrung rekonstruieren, indem ich die koordinatendifferenzen minimiere. Ich habe dazu zunächst einmal eine Simulation geschrieben:

Code:

% nonreflective spatial transformation

% input coordinate
samples = 0:10;
[X Y] = meshgrid(samples);
input = [X(:) Y(:) ones(numel(samples)^2,1)];

% rotation (rightsided)
angle_ = 27;

% scaling
scale = 1.3;

% translation
dx = -2;
dy = +1;

% transformationmatrix
tform = [scale*cos(angle_), -scale*sin(angle_);
scale*sin(angle_), scale*cos(angle_);
dx, dy];

% transformed coordinate
output = input * tform;

plot(X(:), Y(:), 'rx')
hold on
plot(output(:,1), output(:,2), 'bo')
grid on

Funktion ohne Link?

leider scheitere ich deutlich dabei die benutze Transformationsmatrix zu rekonstruieren -.- sogar eine globale suche war vergebens

Code:

x0=[1 20 -2.4 1.5] % [scale, angle_ dx dy]
[x,fval] =...
solveSimilarityTform(input(:,1), input(:,2), input(:,3), output(:,1)', output(:,2)', ones(1,size(output,1)), x0)

function [x,fval] =...
solveSimilarityTform(xIn, yIn, nIn, xOut, yOut, nOut, x0)

[x,fval] = fminunc(@SimilarityTform,x0);
% [x,fval] = ga(@SimilarityTform,4);

function SSE = SimilarityTform(x)
t = [xIn, yIn, nIn] * [x(1)*cos(x(2)), -x(1)*sin(x(2));
x(1)*sin(x(2)), x(1)*cos(x(2));
x(3), x(4)];

% minimize sum squared error
SSE =sum(sum(...
(nIn*xOut-t(:,1)*nOut).^2+...
(nIn*yOut-t(:,2)*nOut).^2));
end
end

Funktion ohne Link?

Code:

result = input * [x(1)*cos(x(2)), -x(1)*sin(x(2));
x(1)*sin(x(2)), x(1)*cos(x(2));
x(3), x(4)];
plot(result(:,1), result(:,2), 'go')

Funktion ohne Link?

ich weiss, ist viel verlangt, aber vielleicht nimmt sich jmd. die Zeit und testet selbst einmal und kann mir sagen ob ich was verbessern kann

MatLabNooB

Themenstarter

Forum-Guru


	Beiträge: 262

	Anmeldedatum: 27.03.09

	Wohnort: ---

	Version: ---

Verfasst am: 22.11.2010, 20:36 Titel:

push


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.