www.pudn.com > ld.rar > Jw2zj.m, change:2008-12-23,size:800b


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 = '地理坐标转直角坐标完成!';