sg_xy2ll.m

一些制作正交曲线网格的matlab源程序

function [lon, lat] = sg_xy2ll(x, y)

% sg_ll2xy -- Stub for m_map-like "m_ll2xy" function.
%  [lon, lat] = sg_xy2ll(x, y) returns the (lon, lat)
%   geographic positions, in degrees, that correspond
%   to the given (x, y).
 
% Copyright (C) 1999 Dr. Charles R. Denham, ZYDECO.
%  All Rights Reserved.
%   Disclosure without explicit written consent from the
%    copyright owner does not constitute publication.
 
% Version of 27-Dec-1999 13:37:03.
% Updated    25-Oct-2000 10:53:22.

[theProjection, theCenter] = sg_proj;
if isempty(theCenter)
	theCenter = [0 0 0];
end

plon = theCenter(1);
plat = theCenter(2);
pang = theCenter(3);

switch lower(theProjection)
case 'mercator'
	[lon, lat] = sg_mercator(x, y, 1, plon, plat, pang);
case 'transverse mercator'
	[lon, lat] = sg_mercator(x, y, 1, 0, 0, 90);
case 'stereographic'
	[lon, lat] = sg_stereographic(x, y, 1, plon, plat, pang);
case 'lambert equal area'
	[lon, lat] = sg_lambert_equal_area(x, y, 1, plon, plat, pang);
otherwise
	disp([' ## No such projection: ' theProjection])
end