GRASS GIS v.surf.rst 명령으로 등고선으로 DEM 만들기

아래 그림과 같은 등고선도에서 DEM을 만듭니다.


우선은 고행 같은 디지타이즈. 5시간 정도 노력해 트레이스 완료.

GRASS GIS로 가져오기

GRASS GIS로 가져오는 것은 QGIS에서 수행하는 쉬운 방법입니다.

GRASS로 분석할 영역을 설정합니다. 쓸데없이 넓은 범위를 설정하면 작성하는 DEM의 데이터가 거대화하기 때문에 필요한 범위를 빠듯하게 설정하는 것이 포인트입니다.
QGIS의 GRASS 플러그인을 사용하면 마우스 드래그로 영역을 설정할 수 있습니다.
해상도는 DEM의 메쉬 해상도를 지정합니다. 여기에서는 0.5m 메쉬로 합니다.


GRASS 도구 → 파일 관리 → GRASS로 벡터 가져오기 → 로드된 벡터 가져오기

GRASS에서 데이터 열기

GRASS GIS를 열어 데이터를 가져올 수 있는지 확인합니다.

v.surf.rst 명령

여기서 v.surf.rst 명령을 사용합니다. 벡터 등고선을 래스터화할 필요도 없습니다.

##ベクタ等高線と標高フィールドを指定するだけ
v.surf.rst input=contour_Tokyo11@contour zcolumn=level elevation=contour_ras_rst --overwrite

입력 할 벡터 등고선 = contour_Tokyo11@contour

출력 DEM=contour_ras_rst

Z 좌표 = level (벡터 등고선의 고도 필드 값)

벡터 등고선의 고도 필드 값에 입력 실수가 있었기 때문에, 이런 에러가 나오고 있습니다.
이 규모의 데이터의 경우, 순조롭게 진행되어 10분 정도로 처리가 끝났습니다.


출력된 DEM 래스터. 좋은 출력입니다.


GeoTiff로 출력합니다.

## GeoTiffに出力
r.out.gdal input=contour_ras_rst@contour output=DEM_Tokyo11.tif format=GTiff

QGIS로 열어 CS 입체도로 한다

QGIS로 열어 음영도와 겹치면 이런 상태. 해자와 수로의 요철도 잘 표현되어 있습니다.


좀 더, 미지형을 관찰하고 싶기 때문에 CS 입체도를 작성합니다.
CS 입체도를 만드는 방법은 키타 코이치 씨가 소개하고 있기 때문에 할애합니다.
QGIS로 「CS 입체도」의 「입체도」를 만들어 보았다

조금 에로그로 느낌으로 완성되고 있습니다. 좀처럼 조정이 어렵고 도중에 좌절했습니다만, 미지형은 「이것이라도!!」라고 하는 만큼 잘 알 수 있습니다.


브라우저에서 글리그 움직이는 3D관성 도 만들어 보았습니다.