Invitation:
Prior to the 2015 NCK days icebreaker, OpenEarth.nl organizes a hands-on sprint session to explain how to access and visualize some great open access core-description and grain size datasets of the Dutch coast:
...
Flyer: NCK days 2015 sprintsession north sea grain sizes.pdf
View File | ||||
---|---|---|---|---|
|
Course information:
Access information for the data can be found at Dataset documentation.
...
Code Block |
---|
set from=http://opendap.deltares.nl/thredds/fileServer/opendap/tno/ncp/ set into=d:\opendap\tno\ bitsadmin /transfer Job001 /download /priority normal %from%slib_juli2007.nc %into%slib_juli2007.nc bitsadmin /transfer Job001 /download /priority normal %from%grind_fbr2007.nc %into%grind_fbr2007.nc bitsadmin /transfer Job001 /download /priority normal %from%dz10_juli2007.nc %into%dz10_juli2007.nc bitsadmin /transfer Job001 /download /priority normal %from%dz50_juli2007.nc %into%dz50_juli2007.nc bitsadmin /transfer Job001 /download /priority normal %from%dz90_juli2007.nc %into%dz90_juli2007.nc |
Yoou You can access subsets these files with the Matlab netCDF package called ncread and with the python package netcdf4-python, see our Matlab example and
python example.
TNO grain size maps
A Matlab simple example how to make a publication quality figure is given here
Code Block |
---|
%% laod openearthtools for handy tools such as nc2struct, pcolorcorcen, colorbarwithvtext, tickmap, axislat, findfont oetsettings %% load data (local cache), D contains data, M contains all metadata [D,M] = nc2struct('slib_juli2007.nc') %% load data (via opendap) L = nc2struct('http://opendap.deltares.nl/thredds/dodsC/opendap/deltares/landboundaries/northsea.nc') %% plot pcolorcorcen(D.longitude,D.latitude,D.mud) colorbarwithvtext([M.mud.long_name,'[',M.mud.units,']']) % add text inside colorbar to maximize use of paper space. automatically use metdata from netCDF file tickmap('ll') % add lat/lon as last tickmarks, no need for space, that dilute efficient use of paper space (google: Edward Tufte) grid on axislat % fix aspect ratio plot(L.lon,L.lat,'k') set(findfont,'fontsize',8) title({M.nc_global.institution,'NCP map mud [version jul 2007]'}) % automatically use metdata from netCDF file as title print2screensize('slib_juli2007') |
...
We used this image as Fig. 8 in our 2011 paper Mechanisms controlling the intra-annual mesoscale variability of SST and SPM in the southern North Sea.
Sedimentatlas waddensea
A second Matlab example is a dataset used for Fig 9 in our 2008 paper Modeling large-scale cohesive sediment transport affected by small-scale biological activity. Here's an exmaple how to work with these data, find more script in
Code Block |
---|
%oetsettings
%% load data
F = 'http://opendap.deltares.nl/thredds/dodsC/opendap/rijkswaterstaat/sedimentatlas_waddenzee/korrel.nc';
D.lon = ncread(F,'lon');
D.lat = ncread(F,'lat');
D.cumphi = ncread(F,'cumphi');
D.diameter = ncread(F,'diameter');
[D.x,D.y] = convertCoordinates(D.lon,D.lat,'CS1.code',4326,'CS2.code',28992);% wgs84 to RD
fraction = 10; % fraction to work with
%% add coastal data
L = nc2struct('http://opendap.deltares.nl/thredds/dodsC/opendap/deltares/landboundaries/northsea.nc')
[L.x,L.y] = convertCoordinates(L.lon,L.lat,'CS1.code',4326,'CS2.code',28992);% wgs84 to RD
%% plot data
caxis ([0 100])
scatter (D.x,D.y,20,D.cumphi(fraction,:),'.')
colorbarwithvtext(['D < ',num2str(D.diameter(fraction)),' \mum [%]'])
axis equal
grid on
tickmap ('xy')
axis(axis)
hold on
plot(L.x,L.y,'k')
box on
print2screensize(['sedimentatlas_waddenzee_fraction_',num2str(D.diameter(fraction)),'mm'])
%% plot in Google Earth
KMLscatter(D.lat ,D.lon ,D.cumphi(10,:),'fileName',['sedimentatlas_waddenzee.kml'],...
'CBcolorTitle',['D < ',num2str(D.diameter(fraction)),' \mum [%]'],...
'CBcolorbarlocation',{'W'});
|
...