read my data in wmsread

3 次查看（过去 30 天）
Victor Ebolo
Victor Ebolo 2023-8-24
评论： 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

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

采纳的回答

Voss
Voss 2023-8-24
Ran in:
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

更多回答（0 个）

类别

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

标签

产品


版本

R2019a

Community Treasure Hunt

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

Start Hunting!

Translated by