*gps*우리나라의 기준원점 경도,위도,높이를 사용하여서 wgs84 타원체에 맞는 좌표를 구하고 변환하는 소스입니다. matlab(맵틀랩) 소스도 첨부합니다.

소개글

우리나라의 기준원점 경도,위도,높이를 사용하여서 wgs84 타원체에 맞는 좌표를 구하고
다시 그좌표를 경도 위도 높이로 바꾸는 프로그램입니다.
matlab 소스도 함께 첨부합니다.

목차

% 우리나라의 기준원점 경도,위도,높이를 사용하여서 wgs84 타원체에 맞는 좌표를 구하고
% 다시 그좌표를 경도 위도 높이로 바꾸는 프로그램입니다.

본문내용

clc;
format long g % 출력 값을 정수로 하기 위해 g 를 입력합니다.
a = 6378137; % wgs 84 타원체의 장반경을 변수 a에 저장합니다.
f = 1/298.257223563; % wgs 84 타원체의 평편도를 변수 f에 저장합니다.
b = a-a*f; % b를 구합니다.
e = 2*f - f^2; % 식을 이용해 이심율의 제곱 값을 E에 저장합니다.


% 위도 경도를 라디안으로 고치기
lati = 0.650530638; % 위도 37,16,21.576 * pi * 180를 곱해서 구한값입니다.
long = 2.217547372; % 경도 127,3,21.979 * pi * 180를 곱해서 구한값입니다.
n = a/((1-e*(sin(lati))^2)^(1/2));
h = (158*30.48)/100000 ;% 표고 값을 km 단위로 환산하여 입력합니다.

% 위도, 경도, 표고를 이용해 얻은 X, Y, Z 값입니다.
fprintf(경도,위도,높이값을 이용하여 얻은 3차원공간의 좌표 값입니다.`)
x = (n+h)*cos(lati)*cos(long)
y = (n+h)*cos(lati)*sin(long)
z = (n*(1-e)+h)*sin(lati)

% 좌료를 사용해서 경도 위도 높이를 구하는 공식에 각 값을 넣어 경도 위도 높이를 구합니다.
fprintf(각 좌표값을 다시 경도,위도,높이로 바꾼 값입니다.`)
p = (x^2 + y^2)^(1/2); % p값을 변수 p에 저장합니다.
th = atan((z*a)/(p*b)); % 세타 값을 변수 theta 에 저장합니다
ep = (a^2 - b^2)/(b^2); % e#/squar값을 변수 Eprime 에 저장합니다.

selati = atan((z+ep*b*(sin(th))^3)/(p-e*a*cos(th)^3))
% 지금까지 구한 값으로 구한 위도() 값을 Piw에 저장합니다.

참고 자료

없음

압축파일 내 파일목록

gps.hwp
lalong.m