[MATLAB] globe data 읽기, 영역 자르기

globe data download : https://www.ngdc.noaa.gov/mgg/topo/gltiles.html

Get Data Tiles-Global Land One-km Base Elevation Project | NCEI

www.ngdc.noaa.gov

globedems([32.1631 43.1645],[120.9934 132.3886]);

% 위 명령어를 치면 영역의 번호가 표출된다 ex) h10g

[map,maplegend] = globedem('h10g',1,[32.1631 43.1645],[120.9934 132.3886]);

% map에는 영역에 해당하는 zeta data가 입력되고, maplegend에는 그에 따른 레퍼런스가 3개의 숫자로 저장되는데, 이 숫자의 자세한 의미를 잘 모름.

% globedem('전에 나온 영역 번호',스케일 숫자(1이 기준),[위도 첫 값, 위도 끝 값],[경도 첫 값, 경도 끝 값]);

% 위 명령어를 입력하게 되면, 헤더 파일을 넣으라는 GUI 방식의 입력창이 표출되게 된다.

% 헤더 파일은 아래의 경로를 통해 구한다.

http://www.ngdc.noaa.gov/mgg/topo/elev/esri/hdr

Index of /mgg/topo/elev/esri/hdr

www.ngdc.noaa.gov

% 이후 globe data의 위치를 묻는 GUI 입력창이 표출된다.

% 넣는다. 끝.

x = linspace(120.9934, 132.3886,1369);
y = linspace(32.1631, 43.1645,1321);
m_pcolor(x,y,map);
 

% 그림 확인