Finding the closest coordinate from a surface plot based on a X, Y location
86 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Jorge Luis Paredes Estacio
le 16 Août 2024 à 16:48
Commenté : Voss
le 20 Août 2024 à 14:26
Hello, I want to extrapolate a point (longitude,latitude) from the coordinates of a surface file (attached here as "slab_strike") as it is empty (Strike=NaN) when using interp2 as the points are outside the boundary. Despite I used the option "nearest", it is empty anyway.
% Coordinates of the points:
lat_GMM= -17.8990;
lon_GMM=-73.5295;
% The surface plot
load slab_strike % Loading the slab strike
Slab_strike.x=x;
Slab_strike.y=y;
Slab_strike.z=z;
Strike = interp2(Slab_strike.x,Slab_strike.y,Slab_strike.z,lon_GMM,lat_GMM)
As Strike=NaN, there is a way I can choose the closest point value from the surface avoiding any NaN value and select the closest non-NaN value.
I would appreciate the help
0 commentaires
Réponse acceptée
Voss
le 16 Août 2024 à 18:17
Déplacé(e) : Torsten
le 16 Août 2024 à 21:11
load slab_strike
is_ok = ~isnan(z);
[X,Y] = meshgrid(x,y);
X = X(is_ok);
Y = Y(is_ok);
Z = z(is_ok);
lon_GMM = -73.5295;
lat_GMM = -17.8990;
[~,idx] = min((X-lon_GMM).^2+(Y-lat_GMM).^2);
figure()
hold on
surf(x,y,z,'EdgeColor','none')
h1 = plot(lon_GMM,lat_GMM,'.r');
h2 = plot3(X(idx),Y(idx),Z(idx),'.k');
legend([h1 h2],{'Requested Point','Nearest Non-NaN z'})
xlim(lon_GMM+[-5 5])
ylim(lat_GMM+[-5 5])
6 commentaires
Jorge Luis Paredes Estacio
le 20 Août 2024 à 13:02
Modifié(e) : Jorge Luis Paredes Estacio
le 20 Août 2024 à 13:02
Plus de réponses (2)
Torsten
le 16 Août 2024 à 17:08
Modifié(e) : Torsten
le 16 Août 2024 à 17:22
lat_GMM and lon_GMM are not within the rectangle in which data for z are supplied.
In this case, interp2 returns NaN because it does not extrapolate.
Choose lat_GMM and lon_GMM in the limits for x and y where values are given for z.
Else you could use
[~,idx] = min(abs(x-lon_GMM))
[~,idy] = min(abs(y-lat_GMM))
Strike = z(idy,idx)
But it seems that your z-matrix contains NaN values.
6 commentaires
Jorge Luis Paredes Estacio
le 16 Août 2024 à 17:50
Modifié(e) : Jorge Luis Paredes Estacio
le 16 Août 2024 à 17:52
Torsten
le 16 Août 2024 à 17:56
Then try whether z is NaN in the 8 points surrounding (x(idx),y(idy)).
If this also doesn't work, test z for NaN in the 16 following points and so on.
Matt J
le 16 Août 2024 à 17:15
Modifié(e) : Matt J
le 16 Août 2024 à 17:19
F=griddedInterpolant({Slab_strike.x,Slab_strike.y},Slab_strike.z,'linear','nearest');
Strike = F(lon_GMM,lat_GMM);
Voir également
Catégories
En savoir plus sur Geographic Plots dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!