MATLAB - toridai_suikou @ Wiki

緯度経度座標からUTM及び第Ⅲ座標系への変換
data=load('Shinjiko_lon_lat.dat')
lat=data(:,1)/180*pi
lon=data(:,2)/180*pi
lon=data(:,2)/180*pi
[N,E,Zone]=ell2utm(lat,lon,a,e2)
lon1=(132+10/60)/180*pi
lat1=(36)/180*pi
[N1,E1,Zone]=ell2utm(lat1,lon1,a,e2)
N=N-N1;
E=E-E1;

x,yデータからspline補間をもとめ、yに対応したxを求める方法の例
x = sort(rand(10,1));
y = exp(x);
pp = spline(x,y);

y0 = 1.5;

x0=fzero(@(x) ppval(x,pp)-y0,.1)
x0 =

   0.4055

雲量データのようにとびとびデータを補完する方法の例
雲量の観測値があるデータのみ取り出したデータを作成（時間順？）
例，ファイル名= unryo.dat
data=load('unryo.dat')
補完したい範囲の時間軸の大きさを確認（例，10000)
xx=[1:10000];
yy=interp1(data(:,1),data(:,2),xx);
yy=yy'

同水位データがあり，水位変化後のデータだけを利用して補完する方法の例
data=load('sample.dat');
[m,n]=size(data)
data2=data;
for i=2:m

if(data(i,2)==data(i-1,2)) 
   data2(i,2)=NaN;
end;

end
i=find(~isnan(data2(:,2)))
data3=data2(i,:);
xx=[1:m];
yy=zeros(m,2);
yy(:,1)=xx;
yy(:,2)=interp1(data3(:,1),data3(:,2),xx);