How to grid data with coordinates to create a spatial plot using geoshow

12 次查看（过去 30 天）

Ronnie 2017-8-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/354663-how-to-grid-data-with-coordinates-to-create-a-spatial-plot-using-geoshow

评论： TAPAS 2018-6-12

I have 3 vectors: data, lat, lon which I am trying to plot spatially for the continental US. Is there a function where I can organize my lat and lon vectors in the appropriate gridded format which geoshow will plot properly while ensuring the data vector is organized in the same fashion so the data points remain with their respective coordinates?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Chad Greene 2017-8-31

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/354663-how-to-grid-data-with-coordinates-to-create-a-spatial-plot-using-geoshow#answer_279878

在 MATLAB Online 中打开

Is it possible that your data vectors are in fact regular, but simply not in gridded format? To check, try

scatterm(lat,lon,20,data)

If the datapoints are regularly spaced, you should be able to use my xyz2grid function like this:

[LON,LAT,DATA] = xyz2grid(lon,lat,data);

where lowercase are your data vectors and upper case are gridded. Then you can do

pcolorm(LAT,LON,DATA)

13 个评论
显示 11更早的评论隐藏 11更早的评论

Ronnie 2017-9-6

Update: I got it to work by trimming down the data even further (getting rid of NaN values) and then xyz2grid worked; and I am able to show the data using geoshow. Thank you so much for your help!

Chad Greene 2017-9-7

Great, I'm glad you got it working!

请先登录，再进行评论。

更多回答（1 个）

KSSV 2017-8-31

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/354663-how-to-grid-data-with-coordinates-to-create-a-spatial-plot-using-geoshow#answer_279776

编辑：Chad Greene 2017-8-31

在 MATLAB Online 中打开

Let data be your nx3 array which has lon, lat and data in the first, second and column respectively.

    % Get longitude and latitude vectors 
    x = unique(data(:,1)) ;
    y = unique(data(:,2)) ;
    % dimensions of the data
    nx = length(x) ; 
    ny = length(y) ;
    % Frame matrix of grid 
    D = reshape(data(:,3),[ny,nx]) ;
    % flip  matrix to adjust for plot
    H = flipud(H) ;
    % Transpose the matrix 
    H = H' ;  % Check if is required
    surf(x,y,H) ;