UTM to Geographical Projection

GMT프로그램 중에 mapproject라는 프로그램을 이용하면 projection을 쉽게 바꿀 수 있다.

awk ‘BEGIN {FS =”,”} {if (NR>1) {print $10, $11}}’ M_T10101.csv | mapproject -Ju1j/1:1 -F -C -I | less

M_T10101.csv파일은 첫번째 라인은 주석(if (NR>1)),

칼럼은 “,”로 구별되기 때문에 FS=”,”를 넣어주고,

10번째 컬럼과 11번째 컬럼에 utmx, utmy가 기록되어 있으므로 이것만 출력({print $10, $11}).

그 다음에 mapproject를 이용해 geographical projection해주면 끝(mapproject -Ju1j/1:1 -F -C -I).

일단 기본적으로 알고 있어야 하는건 utm zone이 어떻게 되느냐.

그건 여기에서 찾아보던가.

http://www.apsalin.com/utm-zone-finder.aspx

아니면 여기에서 찾아보던가..

http://www.dmap.co.uk/utmworld.htm

더 자세한 내용은 메뉴얼페이지를.

http://www.soest.hawaii.edu/gmt/gmt/doc/gmt/html/man/mapproject.html

ps) utm zone 1S하는게 어떤 프로그램은 utm zone 1의 남반구를 뜻하기도 한다. (왜냐면 utm zone이라는게 사실 위도로 구분짓지는 않으니까)

ps2) 그러고 보니까 잃어버린 핸드폰에 soest에서 운영하는 gmt 서버 사진 들어있었는데 말이지 -,.-

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: