read my data in wmsread

Question

Ebolo Nkongo 2023-8-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2012382-read-my-data-in-wmsread

评论： Voss 2023-8-25

clear;close all; clc
Tink_wms='https://nrt.cmems-du.eu/thredds/wms/global-analysis-forecast-bio-001-028-daily?request=GetCapabilities&service=WMS'
WMSinfo=wmsinfo(Tink_wms)
WMSinfo.LayerNames
WMSinfo.Layer(1).Abstract
%definir la variable
variable= WMSinfo.LayerNames(1)
%Define the parameters of the plot
date='2021-03-10T00:00:00.000Z';
lonlim=[-180:180];
latlim= [50:90];
style= 'boxfill/alg'
%read data corresponding to the parameters
[A,R]=wmsread(variable,'Time',date,'Stylename',style,'latlim',latlim,'lonlim',lonlim);
%plot data
figure(1)
gcolor [1 1 1 ] %grid in white
landcolor=[204 204 204]/255 %land in grey
axe=axems('egaazim',...
    ....,'Origin',[90 0 0],'MapLatLim',[50 90],...
    ....,'grid','on','Gcolor',gcolor)................%select the view and the grid
 axis off
geoshow(A,B) %plot ice thickness
land= shaperead ('landareas','UseGeoCoords',true) %plot coastlines
geoshow(land,'FaceColor',landcolor) %plot land
title('inpunt the title') %titre
date='2021-08-10T00:00:00.000z';

I don't know how to limit the study area only for the golf of guinea. lneed help for the modification of the coordinates

lonlim=[-180:180];
latlim= [50:90];
style= 'boxfill/alg'

I would like to represent the distribution of ph in the Gulf of Guinea. at line 17 of my code there is a bug. please I would like you to help me to correct this and also to limit my study area at the level of the golf of guinea thank you in advance

%read data corresponding to the parameters
[A,B]=wmsread(variable,'Time',date,'Stylename',style,'latlim',latlim,'lonlim',lonlim); %my problem starts at this line

%plot data
figure(2)
figure(1)
gcolor[1 1 1 ] %grid in white
landcolor=[204 204 204]/255 %land in grey
axe=axems('egaazim',...
    ....,'Origin',[90 0 0],'MapLatLim',[50 90],...
    ....,'grid','on','Gcolor',gcolor)................%select the view and the grid
 axis off
geoshow(A,B) %plot ice thickness
land= shaperead ('landareas','UseGeoCoords',true) %plot coastlines
geoshow(land,'FaceColor',landcolor) %plot land
title('inpunt the title') %titre

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

请先登录，再进行评论。

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

Answer 1

Voss 2023-8-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2012382-read-my-data-in-wmsread#answer_1293222

在 MATLAB Online 中打开

clear;close all; clc

Tink_wms='https://nrt.cmems-du.eu/thredds/wms/global-analysis-forecast-bio-001-028-daily?request=GetCapabilities&service=WMS'

Tink_wms = 'https://nrt.cmems-du.eu/thredds/wms/global-analysis-forecast-bio-001-028-daily?request=GetCapabilities&service=WMS'

WMSinfo=wmsinfo(Tink_wms)

WMSinfo =

WMSCapabilities Properties: ServerTitle: 'CMEMS DU' ServerURL: 'https://nrt.cmems-du.eu/thredds/wms/global-analysis-forecast-bio-001-028-daily?request=GetCapabilities&service=WMS' ServiceName: 'WMS' Version: '1.3.0' Abstract: 'Scientific Data' OnlineResource: '' ContactInformation: [1x1 struct] AccessConstraints: 'none' Fees: 'no conditions apply' KeywordList: {6x1 cell} ImageFormats: {5x1 cell} LayerNames: {12x1 cell} Layer: [12x1 WMSLayer] AccessDate: '24-Aug-2023' Methods

ph_idx = find(strcmpi(WMSinfo.LayerNames,'ph'));

WMSinfo.Layer(ph_idx).Abstract

ans = 'PH'

%Define the parameters of the plot

% date = '2021-03-10T00:00:00.000Z';

lonlim = [-10 10];

latlim = [-10 10];

style = 'boxfill/alg';

%read data corresponding to the parameters

[A,R] = wmsread(WMSinfo.Layer(ph_idx),'StyleName',style,'latlim',latlim,'lonlim',lonlim);

%plot data

figure

ax = axesm('eqaazim','MapLatLimit',latlim,'FlatLimit',[],'MapLonLimit',lonlim);

Warning: Ignoring the value provided for 'MapLonLimit'. If you are trying to use the 'MapLonLimit' property to initialize the equatorial aspect of an azimuthal projection, you should precede it with 'FLatLimit',[] in your call to AXESM.

geoshow(ax,A,R)

landcolor=[204 204 204]/255; %land in grey

land = shaperead ('landareas','UseGeoCoords',true); %plot coastlines

geoshow(land,'FaceColor',landcolor) %plot land