subobject] = general_region_subset(object,xcoords,ycoords)[subobject] = general_region_subset(object,xcoords,ycoords) subsets object into a polygon region specified by xcoords and ycoords (vertices of the polygon); where object is a struct created by any of the _build functions in ocean_data_tools (e.g. argo, cruise, hycom, mercator, woa, wod).
xcoords and ycoords are vectors of coordinates. Rows or columns are fine. -180/180 or 0/360 notation should match that of object
subobject is a struct which is structurally identical to object but contains only data within the polygon region
% Get variable information:
argo_dir = '/Users/lnferris/Documents/GitHub/ocean_data_tools/data/argo/*profiles*.nc';
listing = dir(argo_dir);
ncdisp([listing(1).folder '/' listing(1).name]) % Peek at netCDF header info to inform choice of variable_list.
% Load Argo data from west of New Zealand:
region = [-60.0 -50.0 150.0 160.0]; % Search region [-90 90 -180 180]
start_date = '01-Nov-2015 00:00:00';
end_date = '01-Jan-2017 00:00:00';
variable_list = {'TEMP_ADJUSTED','PSAL_ADJUSTED'};
[argo,matching_files] = argo_build(argo_dir,region,start_date,end_date,variable_list);
% Choose a region for subsetting the uniform struct:
bathymetry_dir = '/Users/lnferris/Documents/data/bathymetry/topo_20.1.nc';
general_map(argo,bathymetry_dir)
[xcoords,ycoords] = region_select(); % click desired region on the figure
% Subset the struct and remap:
[subargo] = general_region_subset(argo,xcoords,ycoords);
general_map(subargo,bathymetry_dir)
