안녕하세요? 이번 글은 ESA(European Space Agency, 유럽 우주국)이 제공하는 전세계 10m급 토지피복지도 데이터, WorldCover 다운로드 방법을 정리해 보겠습니다. 현재 WorldCover은 Sentinel-1호와 2호 위성 데이터에 기반하여 2020년 전세계 토지피복지도 산출물을 10m급 해상도로 제공하고 있습니다. 공식 홈페이지는 다음과 같습니다.
일단 데이터를 받아보겠습니다. 공식 홈페이지에서 아래와 같이 'EXPLORE(탐색)' 버튼을 클릭합니다.
데이터 다운로드를 위해 일단 회원가입 및 로그인을 진행합니다. 로그인 후 왼쪽 'WorldCover-Map'을 체크하고 우리나라 영역으로 확대해 봅니다.
아래와 같이 보일 텐데요, 여기서 'DOWNLOAD' 버튼을 클릭합니다.
Extent 목록 중 'Administrative borders'를 선택해 봅니다. 남한 경계와 중첩되는 데이터만 받아보려고 합니다.
Borders를 'Countries'로 지정한 후 지도상에서 남한 영역(REBULIC OF KOREA)을 클릭해 줍니다.
아래와 같이 3도 간격으로 설정된 4개 지역이 선택됐는데요, 전부 체크하고 'NEXT(다음)' 버튼을 클릭합니다.
'DOWNLOAD(다운로드)' 버튼을 클릭합니다.
압축 파일을 만드는데 약간의 시간이 소요됩니다.
다운로드 후 압축을 해제하면 다음과 같은 구조로 되어 있습니다.
각각의 폴더에서 '_Map.tif' 파일만 QGIS에서 열어보겠습니다.
아래와 같이 4개 레이어를 추가해 봤습니다. 이 파일들을 하나로 병합해 보겠습니다.
'공간 처리 툴박스 > GDAL > 래스터 기타 > 병합'을 클릭합니다.
병합 실행 창은 다음과 같습니다.
입력 레이어에서 다음과 같이 4개 파일을 체크해 줍니다.
참고로 ESA WorldCover2020 공식 웹사이트에서 Product User Manual(산출물 사용자 매뉴얼) PDF 파일을 내려받을 수 있습니다. 지도 코드는 다음과 같습니다.
병합 창에서 '첫 번째 레이어에서 유사 색상표 가져오기'를 체크하고 산출 데이터형은 기본값인 Float32에서 'Byte'로 변경해 주겠습니다. Byte는 부호가 없는 정수형(8-bit unsigned integer)으로 픽셀 값이 0-255 범위를 갖습니다.
병합 결과는 다음과 같습니다.
ESA World Cover 2020 웹사이트로 접속하시면 60도 간격으로 된 데이터셋을 한번에 내려받으실 수 있습니다. 여기까지 ESA가 제공하는 전세계 10m급 토지피복지도 다운로드 및 전처리 방법을 간략히 정리해 봤습니다.