pbj0812의 코딩 일기

globe data download : https://www.ngdc.noaa.gov/mgg/topo/gltiles.htmlGet Data Tiles-Global Land One-km Base Elevation Project | NCEIwww.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개의 숫자로 저장되는데, 이 숫자의 자세한 의미를 잘..

clear all; close all; clc;[lon lat zeta] = textread('G:\2016_UM\glob_map.dat');%파일 읽기 lon, lat, zeta 한 줄 형태 a = 1321 LAT = reshape(lat,a,length(lat)/a); LON = reshape(lon,a,length(lon)/a); ZETA = reshape(zeta,a,length(zeta)/a);% 실제 형태 n*m 으로 변환 pcolor(LON,LAT,ZETA); shading flatx = 1:a; y = 1:length(lat)/a; x_v2 = linspace(LON(1),LON(end),620); y_v2 = linspace(LAT(1),LAT(end),807); [X,Y]=meshgri..

hd=colorbar('southoutside')% 칼라바 이름 지정 caxis([colormin colormax]);% 칼라바 범위 지정 yt=get(hd,'XTick');% 다른 변수에 칼라바 정보 저장 xlabel(hd,'Tsunami Height (cm)') set(hd,'XTickLabel',sprintf('%3.1f|',yt));% 칼라바를 실수 단위로 지정