jw2zj.m

吐血奉献雷达数据处理中GPS转大地坐标系源程序

function [x y] = Jw2zj(jd,wd);%%经纬度转直角坐标函数

A = 6378245.0;
B = 6356863.0;
E = 0.081813369;
EB = 2.718281828;


L1_map = 13;
L2_map = 57;
BB0_map = 117;
LL0_map = 35;

%L1_map = 13;
%L2_map = 57;
%BB0_map = 116.5;
%LL0_map = 14;

% L1_map = 18;
% L2_map = 62;
% BB0_map = 116;
% LL0_map = 19;

% F0_map = 88;
% G0_map = 27;
% FB0_map = 117;
% FL0_map = 34.1;

XOFF_map = 4000000;
YOFF_map = 0;

[alpha K Ps] = projection_parameter(L1_map,L2_map,LL0_map,A,E);

wd = wd * pi/180;
u_wd = log( tan(pi/4.0 + wd/2.0) ) + (E/2.0) * ( log((1-E*sin(wd))/(1+E*sin(wd))) );

a = alpha * (jd - BB0_map) * pi/180.0;

r = exp( log(K) - alpha*u_wd);

y = ( Ps - r * cos(a) ) + YOFF_map;
x = ( r * sin(a) ) + XOFF_map;

Endstring = '地理坐标转直角坐标完成!';