1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56
|
% Compute metric scale factors.
%
% [pm,pn] = divand_metric(lon,lat)
%
% Compute metric scale factors pm and pn based on
% longitude lon and latitude lat. The variables pm and pn
% represent the inverse of the local resolution in meters using
% the mean Earth radius.
function [pm,pn] = divand_metric(lon,lat)
sz = size(lon);
i = 2:sz(1)-1;
j = 2:sz(2)-1;
dx = distance(lat(i-1,:),lon(i-1,:),lat(i+1,:),lon(i+1,:))/2;
dx = cat(1,dx(1,:),dx,dx(end,:));
dy = distance(lat(:,j-1),lon(:,j-1),lat(:,j+1),lon(:,j+1))/2;
dy = cat(2,dy(:,1),dy,dy(:,end));
dx = real(dx);
dy = real(dy);
dx = deg2m(dx);
dy = deg2m(dy);
pm = 1./dx;
pn = 1./dy;
end
function dy = deg2m(dlat)
% Mean radius (http://en.wikipedia.org/wiki/Earth_radius)
R = 6371.009e3;
dy = dlat*(2*pi*R)/360;
end
% Copyright (C) 2014 Alexander Barth <a.barth@ulg.ac.be>
%
% This program is free software; you can redistribute it and/or modify it under
% the terms of the GNU General Public License as published by the Free Software
% Foundation; either version 2 of the License, or (at your option) any later
% version.
%
% This program is distributed in the hope that it will be useful, but WITHOUT
% ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
% FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
% details.
%
% You should have received a copy of the GNU General Public License along with
% this program; if not, see <http://www.gnu.org/licenses/>.
|