Myname is Andra Ureche and I am studying in Austria my Master in
AerospaceEngineering.Currently I am doing my thesis together with Schiebel in
finding aflight test area in Austria.Ihave created a code in Matlab which is
searching this area in a data file whichcovers Austria and it is in Lambert
EPSG 31287 coordinate system (N,E, andelevation). The next step is export the
solutions from Matlab (polygons definedby 4 points) into a .kml file for
visualization in Google Earth. Here theproblem arrived… I need to convert my
solutions to WGS84. Ihave made this Matlab file to convert one point from
Lambert to WGS84 but I ammissing 3 degrees and I am not sure what is done
wrong. I have attached thefile to this email. Canyou please help me with this
conversion? Thankyou in advance. Kindregards, Andra Ureche
a=6378137.0;
f=1/298.257223563;
e=sqrt(2*f-f^2);
phiF=(47+30/60)*(pi/180);
lambdaF=(13+20/60)*(pi/180);
phi1=49*(pi/180);
phi2=46*(pi/180);
Ef=400000;
Nf=400000;
E=576892.512074;
N=292949.267204;
%phi=0.49741884;
%lambda=-1.67551608;
m1=cos(phi1)/(1-e^2*(sin(phi1))^2)^0.5;
m2=cos(phi2)/(1-e^2*(sin(phi2))^2)^0.5;
%t=(tan((pi/4)-(phi/2))/((1-e*sin(phi))/(1+e*sin(phi)))^(e/2));
tF=(tan((pi/4)-(phiF/2))/((1-e*sin(phiF))/(1+e*sin(phiF)))^(e/2));
t1=(tan((pi/4)-(phi1/2))/((1-e*sin(phi1))/(1+e*sin(phi1)))^(e/2));
t2=(tan((pi/4)-(phi2/2))/((1-e*sin(phi2))/(1+e*sin(phi2)))^(e/2));
n=(log(m1)-log(m2))/(log(t1)-log(t2));
F=m1/(n*(t1^n));
%r=a*F*(t^n);
rF=a*F*(tF^n);
%theta=n*(lambda-lambdaF);
%E=Ef+r*sin(theta);
%N=Nf+rF-r*cos(theta);
rprim=((E-Ef)^2+(rF-(N-Nf))^2)^0.5;
tprim=(rprim/(a*F))^(1/n);
thetaprim=atan((E-Ef)/(rF-(N-Nf)));
phi_initial=(pi/2-2*atan(tprim))*(180/pi);
phi_2iteration=(pi/2-2*atan(tprim*((1-e*sin(phi_initial)/(1+e*sin(phi_initial)))^(e/2))))*(180/pi);
phi_3iteration=(pi/2-2*atan(tprim*((1-e*sin(phi_2iteration)/(1+e*sin(phi_2iteration)))^(e/2))))*(180/pi);
phi_4iter=(pi/2-2*atan(tprim*((1-e*sin(phi_3iteration)/(1+e*sin(phi_3iteration)))^(e/2))))*(180/pi);
lamba_initial=((thetaprim/n)+lambdaF)*(180/pi);
------------------------------------------------------------------------------
Find and fix application performance issues faster with Applications Manager
Applications Manager provides deep performance insights into multiple tiers of
your business applications. It resolves application problems quickly and
reduces your MTTR. Get your free trial!
https://ad.doubleclick.net/ddm/clk/302982198;130105516;z
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users
