read my data in wmsread

1 visualización (últimos 30 días)
Victor Ebolo
Victor Ebolo el 24 de Ag. de 2023
Comentada: Voss el 25 de Ag. de 2023
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

Respuesta aceptada

Voss
Voss el 24 de Ag. de 2023
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
  2 comentarios
Victor Ebolo
Victor Ebolo el 25 de Ag. de 2023
thank you
Voss
Voss el 25 de Ag. de 2023
You're welcome!

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Geographic Plots en Help Center y File Exchange.

Etiquetas

Productos


Versión

R2019a

Community Treasure Hunt

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

Start Hunting!

Translated by