GMT

경사를 명암으로 표현하는 그리드 파일 [grdgradient]

작성자
heukchoo
작성일
2021-01-17 20:26
조회
469
 

Synopsis

grdgradient in_grdfile -Gout_grdfile [ -Aazim[/azim2] ] [ -D[a][c][o][n] ] [ -E[m|s|p]azim/elev[+aambient][+ddiffuse][+pspecular][+sshine] ] [ -Lflag ] [ -N[e|t][amp][+ssigma][+ooffset] ] [ -Rregion ] [ -Sslopefile ] [ -V[level] ] [ -fg ] [ -nflags ]

 




앞서 소개한 grdimage에서 -I 옵션을 이용하여 z값 경사도에 따른 명암을 표현할 수 있음을 보여주었다.

-I 옵션에 입력되는 파일은 grdgradient 툴로 제작한다.
 

xy공간에 빛을 비추었을때 z값의 경사도에 따라 그늘지고 밝아지는 부분이 생길것이다.
이를 빛 강도 값 혹은 명암값(-1~+1 범위)의 그리드 파일로 제작하는 툴이 grdgradient 이다.

예시자료는 128/129/35/36의 고도자료를 사용하였다.

 
gmt grdgradient data.grd -Gintens.grd -A45 -N1
 

필수옵션은 -G이다.
-G 옵션은 제작할 intensity 그리드 파일명을 지정한다.
z값으로 제작된 그리드 파일도 필수로 입력해 주어야 한다.


-A 옵션은 빛이 비춰지는 방향을 지정한다.
'-A90'과 같이 양수 값으로 지정하면 북쪽을 0도로 하여 시계방향으로 인식한다. 예는 동쪽에서 빛을 비추는 경우다.
음수 값은 북쪽을 0으로 하여 반시계 방향으로 인식한다.
'-A90/-45'와 같이 빛을 두 방향에서 비추는 설정도 가능하다.






경사도의 조정


-N 옵션은 경사도를 조정하여 경사에 따른 명암 표현의 정도를 바꾼다.
참고로 경사도를 조정하는 방법은 이해하기 어려울 수 있지만 굳이 이해할 필요는 없다.
예시를 참고하여 여러 조건으로 시험해 본 뒤, 원하는 그림이 나오는 설정을 사용하면된다.


'-N1', '-Ne1', '-Nt1'과 같이 첨자에 따라 다음의 공식으로 경사도를 조정한다.
-N gn = amp * (g - offset)/max(abs(g - offset))
-Ne gn = amp * (1.0 - exp(sqrt(2) * (g - offset)/ sigma))
-Nt gn = (2 * amp / PI) * atan( (g - offset)/ sigma)
 

g는 z값으로 계산된 경사도, gn은 조정된 경사도 이다.
'-N1'의 1은 amp 이다.
amp 값이 클 수록 명암 표현이 극명해 진다.
amp 값은 0보다 커야한다.


offset은 밝기 조정의 효과가 있다.
'-N1+o0.2'와 같이 첨자 +o를 사용하여 지정할 수 있다. (-1~+1 범위)
offset을 지정하지 않은 경우 평균 g를 사용한다.


sigma는 첨자 e나 t를 사용한 경우에 사용하며, 첨자 +s를 사용하여 지정할 수 있다.
sigma는 경사도의 편차를 조정하며 0~1의 값을 지정하여 값이 클 수록 경사도의 편차가 완화된다.
값을 지정하지 않거나 0을 지정하면 g-offset을 sigma로 지정한다.


다음의 각 조건별 결과를 토대로 원하는 조건을 찾도록 한다.
-A45 -N1 -A45 -Ne1 -A45 -Nt1
-A45 -Ne0.5 -A45 -Ne1 -A45 -Ne1.5
-A45 -N1+o-0.2 -A45 -N1+o0 -A45 -N1+o0.2
-A45 -Ne1+s0 -A45 -Ne1+s0.2 -A45 -Ne1+s1
 

 
전체 0