우리나라의 기준원점 경도,위도,높이를 사용하여서 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에 저장합니다.
학교정보
2주간 다운받은 학생의 학교정보를 보여줍니다.(5P 소요)
저작권 정보
위 정보 및 게시물 내용의 진실성에 대하여 해피캠퍼스는 보증하지 아니하며, 해당 정보 및 게시물 저작권과 기타 법적 책임은 자료 등록자에게 있습니다. 위 정보 및 게시물 내용의 불법적 이용, 무단 전재·배포는 금지되어 있습니다.저작권침해, 명예훼손 등 분쟁요소 발견시 고객센터의 저작권침해 신고센터를 이용해 주시기 바랍니다.
문서 요약정보
TEL : 02) 890-3333 | FAX : 02) 890-3334
*gps*우리나라의 기준원점 경도,위도,높이를 사용하여서 wgs84 타원체에 맞는 좌표를 구하고 변환하는 소스입니다. matlab(맵틀랩) 소스도 첨부합니다.
