https://github.com/briochemc/matlab_plotting_routines
Some MATLAB plotting functions for OCIM-modelling-like research
https://github.com/briochemc/matlab_plotting_routines
Last synced: 3 months ago
JSON representation
Some MATLAB plotting functions for OCIM-modelling-like research
- Host: GitHub
- URL: https://github.com/briochemc/matlab_plotting_routines
- Owner: briochemc
- Created: 2017-10-24T19:24:34.000Z (over 7 years ago)
- Default Branch: master
- Last Pushed: 2021-05-21T03:57:45.000Z (about 4 years ago)
- Last Synced: 2024-10-12T10:30:03.152Z (7 months ago)
- Language: MATLAB
- Homepage:
- Size: 31.3 KB
- Stars: 2
- Watchers: 2
- Forks: 0
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
Awesome Lists containing this project
README
# MATLAB plotting routines
## Installation
You don't really need to install this for the functions to work, although some of them call on another, so it's probably best to follow these instructions to properly install this on your laptop.
1. Clone this repository somewhere on your computer by typing:
```bash
git clone https://github.com/briochemc/MATLAB_plotting_routines.git
```
> ***Note***: I suggest creating a directory at the root of your file system called `MATLAB_utilities`.
> For example, this can be done from the command-line interface (CLI) by typing `mkdir ~/MATLAB_utilities`.
> You can then fill this directory with all the functions that you regularly use in MATLAB.
1. Then, you need to add the newly added functions of this repo to the MATLAB path.
This is done by typing, in MATLAB:
```matlab
addpath(genpath('~/MATLAB_utilities'))
savepath
```
After that, any function in your `~/MATLAB_utilities` directory should be available to you whenever you start MATLAB.
## To Do
- Update documentation for
- [ ] `yzplot_v4`
- Add functions
- [x] `inpaint_nans3`
## Examples
### Latitude–depth plot:
```matlab
% load the ao MAT file (you might have to change that)
load ~/Projects/AWESOME-OCIM/data/ao.mat
% dummy tracer
tracer_3d = (cosd(ao.LAT) + sind(ao.LON) + cos(ao.DEPTH / 1000)) ;
tracer_3d(find(~ao.OCN)) = nan ;
tracer_2d = squeeze(tracer_3d(:,100,:)) ;
% plot options
% Originally made for the origina OCIM grid,
% but I could adapt it to the AO...
opt.grid.yt = ao.lat ;
opt.grid.zt = ao.depth ;
opt.grid.dzt = ao.height ;
opt.clevs = -2:0.1:2 ;
opt.dash_clevs = -2:0.5:2 ;
opt.black_clevs = 0:1:2 ;
opt.white_clevs = -2:1:0 ;
opt.yTicks = -6000:1000:0 ;
opt.yTickLabels = cellstr(num2str((6:-1:0)')) ;
opt.xTicks = -90:30:90 ;
opt.xTickLabels = {'SP';'60S';'30S';'EQ';'30N';'60N';'NP'} ;
opt.ylabel = 'Depth' ;
opt.xlabel = 'Latitude' ;
opt.grid.dyt = (ao.lat(2) - ao.lat(1)) * ones(1,length(ao.lat)) ;
opt.grid.yv = ao.lat + opt.grid.dyt / 2 ;
opt.grid.dzt = ao.height ;
opt.grid.zw = ao.depth - ao.height / 2
% nice plot
yzplot(tracer_2d', opt) ;
colorbar
```
### joint PDF
```matlab
figure(1)
N = 2000 ; % number of data points / 2
obs = 5 + 2 * randn(N, 1) ; % fake observations data
model = 1.1 * obs + 0.6 * randn(N, 1) ; % fake model data
obs = obs + 0.2 * randn(N, 1) ; % add noise to obs
v = ones(N, 1) ; % fake volume boxes (can use AO volumes)
jpdf = subplot(1,2,2)
clevs = 0:10:100
jointPDF(obs, model, v, 'clevs', clevs/100) ; % the plotting function
axis square
xlim([-2, 12])
ylim([-2, 12])
xlabel('observations')
ylabel('model')
title('jointPDF')
% colorbar
subplot('position', [.93 .3 0.02 .4]) % Set the correct sizes
percentiles = 0:0.1:100 ;
contourf([0 1], percentiles, kron([1 1], percentiles'), clevs(1:end-1), 'linestyle', 'none')
%caxis([0 100]) ;
set(gca, 'YAxisLocation', 'right')
set(gca, 'xtick', [])
set(gca, 'ytick', clevs)
ylabel('percentile')
cAx = caxis ;
axes(jpdf) ; caxis(cAx/100) ;
subplot(1,2,1)
scatter(obs, model)
axis square
xlim([-2, 12])
ylim([-2, 12])
xlabel('observations')
ylabel('model')
title('scatter')
```
which should look like:
