Commit 3910fcf1 authored by Pierre Cazenave's avatar Pierre Cazenave

Tidy up some old example script and files which are not used anywhere.

parent eecebbec
This diff is collapsed.
function [out_east,out_north] = my_project(in_east,in_north,direction)
% Sample user-defined projection and inverse projection of (lon,lat) to (x,y)
% Copy to my_project (not a member of the toolbox) and modify to suite you
%
% function [out_east,out_north] = my_project(in_east,in_north,direction)
%
% DESCRIPTION:
% Define projections between geographical and Euclidean coordinates
%
% INPUT:
% in_east = 1D vector containing longitude (forward) x (reverse)
% in_north = 1D vector containing latitude (forward) y (reverse)
% direction = ['forward' ; 'inverse']
%
% OUTPUT:
% (lon,lat) or (x,y) depending on choice of forward or reverse projection
%
% EXAMPLE USAGE
% [lon,lat] = my_project(x,y,'reverse')
%
% Author(s):
% Geoff Cowles (University of Massachusetts Dartmouth)
%
% Revision history
%
%==============================================================================
%subname = 'my_project';
%fprintf('\n')
%fprintf(['begin : ' subname '\n'])
%------------------------------------------------------------------------------
% Parse input arguments
%------------------------------------------------------------------------------
ProjectDirection = 'forward';
if(direction == 'forward')
ProjectDirection = 'forward';
lon = in_east;
lat = in_north;
else
ProjectDirection = 'inverse';
x = in_east;
y = in_north;
end;
%------------------------------------------------------------------------------
% Perform the projection: USER DEFINED
% Example: project/inverse project to state plane 1802
%------------------------------------------------------------------------------
%if(ProjectDirection == 'forward')
% fprintf('Projecting from (lon,lat) to (x,y)\n');
% [x,y] = sp_proj('1802','forward',lon,lat,'m');
%
%else
% fprintf('Inverse Projecting from (x,y) to (lon,lat)\n')
% [lon,lat] = sp_proj('1802','inverse',x,y,'m');
%end;
%------------------------------------------------------------------------------
% Skagit, UTM, Zone 10 (see http://www.dmap.co.uk/utmworld.htm)
%------------------------------------------------------------------------------
m_proj('UTM','longitude',[-124,-122],'latitude',[47,49],'zone',10,'hemisphere','north','ellipsoid','wgs84')
%m_proj get
%[x,y] = m_ll2xy(-122.530820 , 48.363114);
%fprintf('x %f y %f\n',x,y-1e7);
%fprintf('should be 534752, 5356766.\n')
deltay = 1e7;
if(ProjectDirection == 'forward')
% fprintf('Projecting from (lon,lat) to (x,y)\n');
[x,y]=m_ll2xy(lon,lat);
y = y - deltay; %why?
else
% fprintf('Inverse Projecting from (x,y) to (lon,lat)\n')
[lon,lat]=m_xy2ll(x,y+deltay);
end;
% set the output
if(ProjectDirection == 'forward')
out_east = x;
out_north = y;
else
out_east = lon;
out_north = lat;
end;
nRivers = 1
River# 1 name tstRiver #Nodes 3
838 844 845
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment