Textur-mapping in MATLAB

Beachten Sie, dass die Funktion trisurf, dass Sie ursprünglich versucht, gibt ein handle zu einem patch Objekt. Wenn man sich die 'FaceColor' - Eigenschaft für patch-Objekte, können Sie sehen, dass es keine 'texturemap' option. Diese option ist nur gültig für die 'FaceColor' - Eigenschaft von Oberfläche Objekte. Sie müssen also einen Weg finden, um zu Plotten Ihre dreieckigen Fläche als Oberfläche - Objekt anstelle eines patch Objekt. Hier sind zwei Möglichkeiten, dies zu nähern:

, Wenn die Daten in ein einheitliches raster...

Wenn die Koordinaten Ihrer Oberfläche Daten repräsentieren ein einheitliches raster, so dass z ist eine rechteckige Menge der Punkte, die Spanne aus xmin zu xmax in der x-Achse und ymin zu ymax in der y-Achse können Sie zeichnen Sie es mit surfen statt trisurf:

Z = ...  % N-by-M matrix of data
x = linspace(xmin, xmax, size(Z, 2));  % x-coordinates for columns of Z
y = linspace(ymin, ymax, size(Z, 1));  % y-coordinates for rows of Z
[X, Y] = meshgrid(x, y);               % Create meshes for x and y
C = imread('image1.jpg');              % Load RGB image
h = surf(X, Y, Z, flipdim(C, 1), ...   % Plot surface (flips rows of C, if needed)
         'FaceColor', 'texturemap', ...
         'EdgeColor', 'none');
axis equal

Zur illustration der Ergebnisse der oben genannten code, ich initialisiert die Daten als Z = peaks;, verwendet die integrierten sample-Bild 'peppers.png', und legen Sie die x und y Werte erstrecken sich von 1 bis 16. Dies führte in den folgenden textur-mapped-Oberfläche:

, Wenn Ihre Daten nicht-gleichmäßiger Abstand...

Wenn Sie Ihre Daten nicht regelmäßig verteilt sind, können Sie eine Reihe von regelmäßig angeordneten X und Y Koordinaten (wie ich oben mit meshgrid) und verwenden Sie anschließend die Funktionen griddata oder TriScatteredInterp zur Interpolation in einem regelmäßigen Gitter von Z Werte aus Ihrer unregelmäßigen Reihe von z Werte. Ich erörtern, wie diese beiden Funktionen in meine Antwort auf eine andere Frage ALSO. Hier ist eine verfeinerte version der code, den Sie mit gebucht TriScatteredInterp (Hinweis: da der R2013a scatteredInterpolant ist die empfohlene alternative):

x = ...  % Scattered x data
y = ...  % Scattered y data
z = ...  % Scattered z data
xmin = min(x);
xmax = max(x);
ymin = min(y);
ymax = max(y);
F = TriScatteredInterp(x(:), y(:), z(:));  % Create interpolant
N = 50;  % Number of y values in uniform grid
M = 50;  % Number of x values in uniform grid
xu = linspace(xmin, xmax, M);         % Uniform x-coordinates
yu = linspace(ymin, ymax, N);         % Uniform y-coordinates
[X, Y] = meshgrid(xu, yu);            % Create meshes for xu and yu
Z = F(X, Y);                          % Evaluate interpolant (N-by-M matrix)
C = imread('image1.jpg');             % Load RGB image
h = surf(X, Y, Z, flipdim(C, 1), ...  % Plot surface
         'FaceColor', 'texturemap', ...
         'EdgeColor', 'none');
axis equal

In diesem Fall müssen Sie ersten wählen Sie die Werte von N und M für die Größe der matrix Z. Um zu veranschaulichen, die Ergebnisse der oben genannten code, ich initialisiert die Daten für x, y, und z wie folgt und verwendet die eingebaute sample-Bild 'peppers.png':

x = rand(1, 100)-0.5;  % 100 random values in the range -0.5 to 0.5
y = rand(1, 100)-0.5;  % 100 random values in the range -0.5 to 0.5
z = exp(-(x.^2+y.^2)./0.125);  % Values from a 2-D Gaussian distribution

Diese führte in den folgenden textur-mapped-Oberfläche:

Feststellen, dass es scharfe Kanten in der Nähe der Ecken der Oberfläche. Dies sind die Orte, an denen waren es zu wenige Punkte für TriScatteredInterp ausreichend fit mit einer interpolierten Oberfläche. Die Z Werte an diesen stellen sind daher nan, was in der Fläche Punkt nicht gezeichnet werden.