Surface de la parcelle avec 3 vecteurs Matlab

J'ai besoin d'être en mesure de faire une surface de la parcelle à l'aide de données à partir de 3 vecteurs. J'ai trouvé des informations similaires, mais aucune méthode semble fonctionner avec mes données. Mon X et Y colonnes sont régulièrement espacés, mais pas dans l'ordre croissant. J'ai essayé différentes méthodes, mais aucune d'entre elles semblent me donner ce que je veux, qui est une simple surface linking à proximité des points ensemble. J'ai essayé le suivant:

[X Y]=meshgrid(x,y);
Z=griddata(x,y,z, X,Y);
surf(X,Y,Z);

Ce n'est pas exactement ce que je veux, parce qu'il crée une surface à z=0 et en fait ressemble plus à un volume de la parcelle que juste une surface. Il fonctionne aussi très lentement sur mon ordinateur (probablement à partir de la création de tous les gridpoints). Si je pouvais obtenir quelque chose qui n'a pas besoin d'autant de mémoire qu'elle serait l'idéal (à mon vecteurs ont environ 20k valeurs de chacun), mais ce n'est pas une nécessité.

***Edit: j'ai aussi essayé d'utiliser le scatteredInterpolant méthode à la ici,mais la fonction ne semble pas être reconnu par MATLAB et j'obtiens cette erreur: 

Undefined function 'scatteredInterpolant' for input arguments of type 'double'.

Aussi, voici une image de mon problème:
Surface de la parcelle avec 3 vecteurs Matlab

Vous pouvez voir que nous ne pouvons pas voir sous la surface, il y a quelques plan z=0, la bloque.

Si vous avez quelque chose pour moi, toute aide est très appréciée.

Merci d'avance.

**Edit 2: j'ai ajouté de l'échantillon vecteurs, ils sont mon x,y et z les valeurs de la gauche vers la droite.

***Edit 3: Voici une image de la triangulation que je reçois. Comme vous pouvez le voir, certains points sont ignorés pour une raison quelconque, ce qui donne à ces long et bizarre à la recherche triangles bleus.

Surface de la parcelle avec 3 vecteurs Matlab

Mike

  • Il aide si vous ne Z(~Z)=NaN? Ou si vous utilisez un autre method comme décrit dans doc griddata? Sinon, s'il vous plaît poster des exemples de code et des vecteurs de parfaitement reproduire le problème.
  • Quand je fais Z(~Z)=NaN-je obtenir un flatter les bas mais je ne peux toujours pas voir à travers elle, elle le ferme et forme un volume. Mes vecteurs sont assez grands, je vais essayer de trouver un moyen de les lier.
  • Êtes-vous sûr que vos données ne contiennent pas les valeurs faibles qui sont tracées? Quel est le minimum de z où x et y sont entre 0,8 et 1,2? Aussi, n'oubliez pas d'essayer les différentes griddata options.
  • Oui, dans le cas d'un nuage de points, nous pouvons clairement voir qu'il n'y a pas de points en vertu de la surface. J'ai essayé toutes les méthodes dans griddata et ils ont tous se ressembler. Est-il possible de le faire sans l'aide d'une grille? Il n'y a pas d'options qui semble juste à divers points et tente de faire le lien avec une surface si elles sont suffisamment proches l'un de l'autre? La grille de rendre les bords carrés et il n'est pas vraiment souhaitable/ La méthode de delaunay presque atteint, mais il semble y avoir une erreur avec la méthode de triangulation, qui donne d'étranges liens, les couleurs, et le même plan z=0 problème.
InformationsquelleAutor Mike Nelson | 2013-11-21
  1. 1

    Que les méthodes classiques ne semblent pas arriver, je vous suggère de le faire manuellement.

    1. Créer un Z matrice complète de NaN valeurs. La taille de la matrice doit être la personne à charge sur votre x et y valeurs.
    2. En boucle sur tous les présentes x,y, paires et de mettre leur (moyen?) z valeur dans la bonne position de votre Z de la matrice.
    3. En boucle sur toutes les NaN et d'interpoler les valeurs de leur valeur. Peut-être en utilisant filter2.
    4. Utilisation surf le tracé de la surface résultante
    • Merci de répondre. Je n'ai pas vraiment tout suivi mais, que voulez-vous dire "mis à leur valeur de z dans la position de droite"? Mon vecteurs x,y,z sont déjà commandés de telle sorte que (x(i),y(i),z(i)) forme le point que je veux graphique (nuage de points fonctionne très bien). Aussi comment peut-on interpoler des valeurs NaN si ils n'ont pas de valeurs par définition? Merci encore il est très apprécié.
    • Ce que j'ai essayé de dire: La griddata approche génère une Z matrice composée de prédiction, de valeurs et d'interpolations. Essayez de créer cette matrice vous-même en quelque sorte. Peut-être en commençant seulement avec survenant valeurs et ensuite de faire votre propre interpolations.
    InformationsquelleAutor Dennis Jaheruddin
  2. 1

    Si vous avez des points qui sont décrits par des vecteurs, et vous voulez tracer eux, vous pouvez toujours utiliser une Delauny triangulation. La fonction dans matlab est appelé Tri=delauny(X,Y,Z). Les données générées par cette fonction peut être décrite avec soit trimesh(Tri,X,Y,Z) ou trisurf(Tri,X,Y,Z). Gardez à l'esprit trisurf est uniquement pour les données 3D. Si vous souhaitez ajuster la transparence de parcelles dans votre graphique utiliser le alpha réglage.

    J'espère que cette aide

    • Après avoir essayé trimesh sur mes données, j'obtiens quelque chose qui ressemble à ceci. On peut voir les lignes qui se croisent sous la surface. Le trisurf semble bon mais j'ai toujours ce plan z=0 à venir, je suppose, de la manière que je vois dans trimesh.
    • Désolé pour ma lenteur de réponse, j'étais absent pour le week-end. Si je vous comprends bien vous souhaitez vous débarrasser de l'triangles dans le plan z=0. Cela peut être fait facilement par l'identification des triangles qui décrivent ce domaine. Le trisurf commande a l'option de coloration spécifique de triangles avec un vecteur (appelons la comme C dans l'aide), trisurf(Tri,X,Y,Z,C). Avec C, vous pouvez faire le tour et d'identifier les triangles. Une fois identifiés, de faire un nouveau tableau Tri2, qui est le même que le Tri sans offenser les triangles. Alors ne trisurf(Tri2,X,Y,Z). C'est un peu maladroit, mais il est difficile de prévoir les triangles
    • L'autre, plus de façon générique, serait d'écrire un algorithme qui teste la composante z de chaque sommet de chaque triangle. Si tous les sommets ont composante z = 0, il décrit le plan z=0 et peut être supprimée.
    InformationsquelleAutor Martin
  3. 1

    Pour moi, il semble que vous juste besoin de trier vos données avant de traçage.

    Voici un exemple qui, je crois, est similaire à votre cas (puisque je ne pouvais pas télécharger vos données).

    x = [2 1 4 3 -1 -3 -4 -2];
y = [1 2 3 4 -1 -2 -3 -4];
z = 32 - x.*x - y.*y;

[X1 Y1] = meshgrid(x,y);
Z1 = 32 - X1.*X1 -Y1.*Y1; 
surf(X1,Y1,Z1)

aux = sort([x;y],2);
x = aux(1,:);
y = aux(2,:);

[X2 Y2] = meshgrid(x,y);

Z2 = 32 - X.*X - Y.*Y;
figure()
surf(X2,Y2,Z2)

    La première figure qui en résulte est très problématique surface:Surface de la parcelle avec 3 vecteurs Matlab

    La deuxième figure contient la surface désirée:Surface de la parcelle avec 3 vecteurs Matlab

    InformationsquelleAutor xuva