pbj0812의 코딩 일기

clear all; clc; close all; wait_time=input('write update time(sec) : ') num_file=textread('G:\2016_대지진\m_file\00__gui\test_empty_data\tsunami.txt'); length_file=size(num_file); length_file=length_file(1); a1=num_file(end,1); a2=num_file(end,2); aaa=(a1~=a2); bbb=(length_file

출처 - matlab을 이용한 해양자료처리(2006, 유학렬)

ex)r = 150; % pixels per inch set(gcf, 'PaperUnits', 'inches', 'PaperPosition', [0 0 1500 700]/r); print(gcf,'-dpng',sprintf('-r%d',r), ['G:\2016_swan\WAVE_KO_2015\result\',num2str(spot),'_month.png']);

ncdisp('~.nc','변수명')

datestr(739529 + unix_time/86400,'yyyy-mm--dd HH:MM:SS') ans = 2015-05-03 01:00:00 숫자를 변형하고 싶을때는 datestr(739529 + unix_time/86400 + 원하는 숫자,'yyyy-mm--dd HH:MM:SS') 여기서 +1을 하게 되면 하루가 증가한다. 따라서 1시간을 더하고 싶으면 +1/24를 한다.

앞에 \를 붙인다. m_vec을 만들고 싶으면 m\_vec 입력

버튼을 지웠는데 그에 해당하는 function이 m-file 안에 남아있나 확인. cd 등의 명령어를 통하여 다른 위치로 path 가 이동했는지 확인(방지하기 위해 addpath 등의 명령어를 통해 gui가 위치한 path 저장하기)

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));% 칼라바를 실수 단위로 지정