GMT

보간하지 않은 그리드 파일 생성 [xyz2grd]

작성자
heukchoo
작성일
2021-01-17 20:22
조회
557
 

그리드 파일 생성 툴 중 surface는 빈 공간을 보간하여 그리드 파일을 생성한다고 하였다.

그리드 파일 영역 내에 데이터 공백이 크지 않을 때는 surface를 통한 보간이 유용하다.

그러나 데이터 공백이 클 경우 surface는 큰 공백을 억지로 보간하여 그리드 파일을 제작한다.


이것을 해결하기 위해서는 보간없이 그리드 파일을 생성해야 한다.

xyz2grd는 설정한 격자내에서 데이터를 평균하여 격자에 하나의 대표값을 갖도록 한다.
데이터가 없는 격자는 보간을 하지 않고 NaN 값으로 채운다.

 
gmt surface data.txt -Gdata.grd -R128/129/35/36 -I0.2m -T0.35 gmt xyz2grd data.txt -Gdata.grd -R128/129/35/36 -I2m
 

위의 예시는 128/129/35/36 범위의 고도 데이터를 임의로 가운데 데이터를 삭제한 후 surface와 xyz2grd로 각각 그리드 파일을 생성한 결과이다.
(기존 데이터의 해상도가 1m 이상인 경우도 있어 xyz2grd의 -I 옵션은 2m으로 설정하였다.)

surface는 가운데 빈 공간을 보간하여 z값을 채워 넣었지만 공백이 너무 커서 왜곡된 정보로 표현된다.
이러한 경우 보간보다는 공백으로 유지하는 것이 좋을 것이다.

xyz2grd는 surface와 달리 가운데 공백을 그대로 표현하고 있다.
회색은 GMT에서 NaN 값을 나타내는 기본색이다.
또한, xyz2grd는 보간을 하지 않기 때문에 기존 데이터의 해상도 이상을 표현할 수 없다.






Synopsis

xyz2grd [ table ] -Ggrdfile -Iincrement -Rregion [ -A[f|l|m|n|r|s|u|z] ] [ -Dxname/yname/zname/scale/offset/invalid/title/remark] [ -S[zfile] ] [ -V[level] ] [ -Z[flags] ] [ -bibinary ] [ -dinodata ] [ -fflags ] [ -hheaders ] [ -iflags ] [ -r ] [ -:[i|o] ]


필수 옵션은 -G, -I, -R로서 사용방법은 surface와 동일하다.
다만 -I 옵션으로 설정한 해상도가 입력 데이터의 해상도 보다 높을 경우 바둑판식 점으로만 표현된다.

그 외에 옵션은 특별한 경우에만 사용한다.

gmt xyz2grd data.txt -Gdata.grd -I2m -R128/129/35/36
 

 

 
전체 0