Grid search and grid plot with colours

10 次查看（过去 30 天）

显示更早的评论

Anitha Limann 2021-12-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1613925-grid-search-and-grid-plot-with-colours

评论： Anitha Limann 2021-12-20

采纳的回答： Voss

在 MATLAB Online 中打开

Hello,

I would like to do a grid search by assining zeros to each grid point first. And then I would like to calculate required values for each grid point and finally I am required to plot the grid with relevant values. here my grid represent error values for each grid point. So I want to plot this grid using some colour scheme method.

Attached here is the type of graph I am looking for.

Could someone please help me with this?

Below is the code I wrote for this. However I am not sure how much of assigning zeros, for loop and then plotting is correct.

I tried to use plot and quiver fuctions but different length x and y axis appeared to be an error in that case. For now i used surf function, though i am not sure if that is correct. But it gives me a z axis too, which is not what i want.

Please let me know if someone could help me with this.

% Obserevd Z location km/Myr ## 176.7 W, 19.2 S, 64 mm/yr (km/Myr) velocity; 
Zlon=-176.324; 
Zlat=-18.395; 
Zobsw=64; 
R=6371; %km
% evaluate error at each grid point
AClon =-177.5:0.01:-176.7;
AClat =-20.5:0.01:-19;
[X,Y]=meshgrid(AClon,AClat);
AC=zeros(151,81);
for i=[1:151]
    for j=[1:81]
        arclen = distance(Zlon,Zlat,AClon(j),AClat(i)); % arc distant between 1.Z location and 2.each grid point.
        w =(Zobsw/(R*sin(theta*pi/180)))*180/pi; % deg/Myr
        vl = plate_vel(Zlon,Zlat,AClon(j),AClat(i),w); % i found this plate-vel funtion online. please use any set of data as an example here.
        misfit=(sum((Zobsw-vl).^2)./(numel(vl)));
        AC(i,j)=misfit;
    end
end
surf(X,Y,AC)

18 个评论
显示 16更早的评论隐藏 16更早的评论

Anitha Limann 2021-12-19

在 MATLAB Online 中打开

sure. Thank you. posted below is the code.

clear;clc;
% Obserevd Z location km/Myr ## 176.7 W, 19.2 S, 64 mm/yr (km/Myr) velocity; 
Zlon=-176.324; 
Zlat=-18.395; 
Zomega=64; 
R=6371; %km
% evaluate error at each grid point
AClon =-176.8:0.01:-176.5;
AClat =-19.5:0.01:-19;
[X,Y]=meshgrid(AClon,AClat);
AC=zeros(51,31);
for i=[1:51]
    for j=[1:31]
        theta = arcdist(Zlon,Zlat,AClon(j),AClat(i)); % arc distant between 1.Z location and 2.each grid point.
        AComega =(Zomega/(R*sin(theta*pi/180)))*180/pi; % deg/Myr
        [pvel,azi] = plate_motion(Zlon,Zlat,AClon(j),AClat(i),AComega);
        misfit=(sum((Zomega-pvel).^2));
        AC(i,j)=misfit;
    end
end
a1 = median(AC(:));
a2 = median(min(AC(:)):a1);
a3 = median(a1:max(AC(:)));
AC_group = zeros(size(AC));
AC_group(AC <=a2)=1; % 4 groups. not sure if i need to make groups at all.
AC_group(AC >a2 & AC <=a1)=2;
AC_group(AC >a1 & AC <=a3)=3;
AC_group(AC >a3)=4;
% Both scatter and surf were used, since i am not sure how to the required type of a map. 
figure(1)
scatter3(X(:),Y(:),AC_group(:),[],AC_group(:),'Marker','.');
set(gca(),'View',[0 90]);
set(gcf(),'Colormap',winter); % color is [blue; red] <-> AC_group is [1; 2]
figure(2)
surf(X,Y,AC);
% surf(X,Y,zeros(size(X)),AC,'EdgeColor','none'); 
% set(gca(),'View',[0 90]);
% colormap('winter');

Voss 2021-12-20

No problem! Don't hesitate to ask another question if you suspect your calculations are off.

Anitha Limann 2021-12-20

Sure will do. Thank you!

请先登录，再进行评论。

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

采纳的回答

Voss 2021-12-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1613925-grid-search-and-grid-plot-with-colours#answer_858510

在 MATLAB Online 中打开

surf(X,Y,zeros(size(X)),AC,'EdgeColor','none','FaceColor','interp'); 
set(gca(),'View',[0 90]);
set(gcf(),'Colormap',flip(gray(),1));

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

请先登录，再进行评论。

类别

MATLAB Graphics 2-D and 3-D Plots Surfaces, Volumes, and Polygons Surface and Mesh Plots

在 Help Center 和 File Exchange 中查找有关 Surface and Mesh Plots 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by

Grid search and grid plot with colours

18 个评论
显示 16更早的评论隐藏 16更早的评论

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

Grid search and grid plot with colours

18 个评论 显示 16更早的评论隐藏 16更早的评论

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

WeChat

18 个评论
显示 16更早的评论隐藏 16更早的评论

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