read my data in wmsread
2 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
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 Kommentare
Akzeptierte Antwort
Voss
am 24 Aug. 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)
ph_idx = find(strcmpi(WMSinfo.LayerNames,'ph'));
WMSinfo.Layer(ph_idx).Abstract
%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);
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 Kommentare
Weitere Antworten (0)
Siehe auch
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!