GRASS를 사용하여 DEM 스무딩?


45

SRTM DEM이 있는데 음영 처리 된 구호를 만들고 싶습니다. GRASS에서 음영 처리 된 구호를 만들었으며 결과는 매우 좋지만, 면적이 거의 평평하고 DEM의 해상도가 90m이므로 약간 거칠습니다.

내가 원하는 것은 부드러운 음영 처리를 생성하기 위해 DEM을 더 매끄럽게 만드는 것입니다. 그렇게하는 알고리즘이나 보간법이 있습니까?

아이디어를 얻으려면 음영 처리 된 구호가 있습니다.이 작은 융기를 평평하게하고 싶습니다.

영상


답변:


43

어떻게 존 스티븐슨의 약 r.denoise의에서 GRASS 애드온은 위키 :

r.denoise denoises (부드러운 / 디스 플) 지형 데이터, Xianfang Sun의 노이즈 제거 알고리즘을 사용하여 레이더 데이터 (SRTM 포함)에서 파생 된 특정 DEM . 원본 데이터를 최소한으로 변경하여 날카로운 모서리를 유지하고 노이즈를 제거하도록 설계되었습니다.

mdenoise

Esri ASCII Grid 파일을 사용하는 것이 더 일반적인 방법이라는 이 웹 사이트 에서 더 많은 내용을 읽었습니다 . 의 위치 mdenoise(에서 다운로드 썬의 웹 사이트는 ) PATH 변수에 있어야합니다 (예를 들어, Windows 사용자 : 드롭 MDenoise.exe에서 bin당신의 OSGeo4w 또는 FWTools와 폴더에 설치). 예를 들어, 다음 쉘 명령을 사용하여 ASCII 그리드 파일을 처리 할 수 ​​있습니다.

# gdal_translate -of AAIGrid my_dem.tif my_dem.asc      # convert to .asc
mdenoise -i my_dem.asc -n 5 -t 0.99 -o my_dem_DN.asc    # denoise
# gdal_translate -of GTiff my_dem_DN.asc my_dem_DN.tif  # convert back to .tif

Denoise는 GNU 라이센스하에 있습니다. 여기를 참조 하십시오


파이썬 / gdal 포티지?
Hugolpz

1
항상 Python에서 모든 명령을 실행할 수 있습니다.
Mike T

R을 필요로하지 않는 파이썬과 동등한 것을 의미합니다.
Hugolpz

1
이것은 R이 필요하지 않습니다. Python os.system('mdenoise ...')에서 실행 파일을 실행하는 데 사용 하십시오.
Mike T

명확히하기 위해 : r. 여기는 래스터를 나타냅니다. 명령은 R이 아닌 GRASS GIS에서 온 것입니다.
bugmenot123

11

US National Parks Service의 수석지도 제작자 인 Tom Patterson 은 DEM 데이터를 사용하여 아름다운 음영을 완화하는 방법에 대한 훌륭한 자습서 를 제공합니다. 그의 작업 과정 중 일부는 Natural Scene Designer 및 Adobe Photoshop을 사용하는 것입니다.

내 워크 플로의 경우 음영 처리를 렌더링하기 전에 GDAL을 사용하여 DEM의 크기를 다시 샘플링하고 싶습니다. 파일 크기는 말할 것도없고 디테일과 노이즈의 양을 줄이는 데 도움이됩니다. 나는이 Github에서에 튜토리얼 이 작업을 수행하는 방법을 보여줍니다.

기본적으로 프로세스는 원본 파일 크기보다 작은 출력 DEM의 새 너비 및 / 또는 높이 (픽셀)를 지정하는 것입니다. 예를 들어

gdalwarp -ts 3000 0 -r bilinear kings_canyon_2228.tif kings_canyon_2228_rs.tif

DEM의 너비를 원래 너비 3800 픽셀에서 3000 픽셀로 줄입니다. 높이를 0으로 설정하면 GDAL이 원본의 종횡비를 기반으로 새 파일의 최적 높이를 결정할 수 있습니다.


7

김프 플러그인 웨이블릿 노이즈 감소 도구를 사용했는데 좋은 결과와 빠른 결과를 얻었습니다.

김프 플러그인 웨이블릿 노이즈 감소

여기에 이미지 설명을 입력하십시오

여기에 이미지 설명을 입력하십시오


래스터 (힐 셰이드)를 부드럽게했지만 파일 (tiff)에서 공간 참조가 손실되었습니다. 따라서 qgis 프로젝트에 넣을 때 다른 레이어와 올바르게 오버레이되지 않았습니다. 이 일이 당신과 함께 일어 났습니까?
Diogo Caribé

지리적 위치를 다시 주입하는 방법이 있습니다.
Hugolpz

예, QGIS에서 래스터 메뉴의 언덕 모양을 내보내고 PNG로 추출합니다.
SamTux

6

DEM에서 간단한 필터 (예 : 5x5 이동 창 평균 필터)를 실행할 수 있습니다.


3

ArcMap에서 나는 공간 분석가> 초점 통계 도구를 사용하여 등고선을 만들기 전에 레이더를 부드럽게하고 싶습니다 ... 이것은 언덕 음영 처리에 도움이 될 수 있습니다. 그것은 내가 믿는 위의 게시물과 비슷한 기능입니다.


죄송합니다. GRASS 부분을 읽으십시오 ... 글쎄 비슷한 작업이있을 것입니다 ...
aug_aug
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.