克里金插值的matlab实现

克里金插值法的计算机实现多种多样，这里是克里金插值的matlab实现。克里金法（Kriging）是依据协方差函数对随机过程/随机场进行空间建模和预测（插值）的回归算法。在特定的随机过程，例如固有平稳过程中，克里金法能够给出最优线性无偏估计（Best Linear Unbiased Prediction, BLUP），因此在地统计学中也被称为空间最优无偏估计器（spatial BLUP）。

function	initialization（）
% function initialization perform initialization

global para color data hdl

    
curDir=pwd;							% current directory
para.curDir=curDir;
cmd=[‘which（‘‘startkrig‘‘）‘];
startfilename=eval（cmd）;
indx=find（startfilename == ‘.‘）;
HDIR=startfilename（1:max（indx）-11）;
DirBackslash=length（find（HDIR == ‘\‘））;
DirColon=length（find（HDIR == ‘:‘））;
DirSlash=length（find（HDIR == ‘/‘））;
HDIRpath=[‘‘ HDIR ‘‘];

if DirBackslash >= 1			% WINDOWS
   DirMark=‘\‘;
   AddPathCmd=[‘addpath ‘  ‘‘‘‘ HDIRpath ‘\general‘‘‘ ...
                ‘ ‘‘‘ HDIRpath ‘\dataprep‘‘‘ ‘ ‘‘‘ HDIRpath ‘\variogram‘‘‘ ...
                ‘ ‘‘‘ HDIRpath ‘\krig‘‘‘ ‘ ‘‘‘ HDIRpath ‘\visualization‘‘‘ ...
                ‘ ‘‘‘ HDIRpath ‘\help‘‘‘ ‘ ‘‘‘ HDIRpath ‘\images‘‘‘ ...
                ‘ ‘‘‘ HDIRpath ‘\bin‘‘‘ ‘ ‘‘‘ HDIRpath ‘\data‘‘‘ ...
                ‘ ‘‘‘ HDIRpath ‘\output‘‘‘ ‘ ‘‘‘ HDIRpath ‘\‘‘‘  ‘ -begin‘];
   PlatForm=1;
   markersize=6;
elseif DirSlash >= 1			% UNIX/LINUX
   DirMark=‘/‘;
   AddPathCmd=[‘addpath ‘  ‘‘‘‘ HDIRpath ‘/general‘‘‘ ...
                ‘ ‘‘‘ HDIRpath ‘/dataprep‘‘‘ ‘ ‘‘‘ HDIRpath ‘/variogram‘‘‘ ...
                ‘ ‘‘‘ HDIRpath ‘/krig‘‘‘ ‘ ‘‘‘ HDIRpath ‘/visualization‘‘‘ ...
                ‘ ‘‘‘ HDIRpath ‘/help‘‘‘ ‘ ‘‘‘ HDIRpath ‘/images‘‘‘ ...
                ‘ ‘‘‘ HDIRpath ‘/bin‘‘‘ ‘ ‘‘‘ HDIRpath ‘/data‘‘‘ ...
                ‘ ‘‘‘ HDIRpath ‘/output‘‘‘ ‘ ‘‘‘ HDIRpath ‘/‘‘‘  ‘ -begin‘];
   PlatForm=2;  
   markersize=8;
end
eval（AddPathCmd）

%% Optimization Toolbox
OptimOption=0;
pp=path;
pindx=find（pp == DirMark）;
if pindx+5 <= length（pp）
   ppl=length（pindx）;
else
   ppl=length（pindx）-1;
end
for i=1:ppl
   StrIndx=pindx（i）+1:pindx（i）+5;
   if max（StrIndx） <= length（pp）
     DirStr=pp（StrIndx）;
     if DirStr == ‘optim‘
       OptimOption=1;
     end
   end
end

data.in.dim=2;              % default is 2D case

% set parameters
para.home_dir=HDIR;
para.platform=PlatForm;
para.optim=OptimOption;
para.status=0;
para.file_dir.data_conversion=HDIR;
para.file_dir.datafile=HDIR;
para.file_dir.data_format_file=HDIR;
para.file_dir.gridfile=HDIR;
para.file_dir.parafile=HDIR;
para.file_dir.batch_filename=HDIR;
para.file_dir.batch_log=HDIR;
para.file_dir.mat_file_in=HDIR;
para.file_dir.mat_file_out=HDIR;

para.dataprep.filename=‘‘;
para.dataprep.ext_prog=0;
para.dataprep.dat_conv_fname=‘‘;
para.dataprep.xy_switch=0;
para.krig.data_format_file=[];
para.status.dataprepfig=0;
para.status.dataprep=0;
para.status.variogramfig=0;
para.status.variogram=0;
para.status.krigingfig=0;
para.status.kriging=0;
para.status.dispkrigfig=0;
para.status.dispkrig=0;

hdl.status.dataprepfig=0;
hdl.status.variogramfig=0;
hdl.status.krigingfig=0;
hdl.status.dispkrigfig=0;

para.dataprep.ytox=1;
para.dataprep.ztox=1;
para.dataprep.ext_prog=0;
para.dataprep.filter_type=2;				% default filter = mean
para.dataprep.reduct_fac=1;
para.dataprep.fil

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2004-05-13 15:52  EasyKrig_V3.0\
     目录           0  2004-05-13 15:52  EasyKrig_V3.0\bin\
     文件         453  2000-09-05 18:18  EasyKrig_V3.0\bin\bessel0.m
     文件        2442  2000-03-07 12:42  EasyKrig_V3.0\bin\degmins.m
     文件         845  2000-07-21 12:44  EasyKrig_V3.0\bin\ll2xy.m
     文件        2645  2003-04-30 17:03  EasyKrig_V3.0\bin\mapax.m
     文件        1029  2001-10-02 10:00  EasyKrig_V3.0\bin\mean_nan.m
     文件        1173  1999-09-22 14:41  EasyKrig_V3.0\bin\median_nan.m
     文件        1095  2000-05-03 10:09  EasyKrig_V3.0\bin\sort_mex.m
     文件        1028  1999-09-22 14:42  EasyKrig_V3.0\bin\std_nan.m
     文件         965  1999-09-22 14:43  EasyKrig_V3.0\bin\sum_nan.m
     文件        1462  2000-07-21 12:45  EasyKrig_V3.0\bin\xy2ll.m
     文件        3127  2004-07-16 18:54  EasyKrig_V3.0\color_info.dat
     目录           0  2004-05-13 15:52  EasyKrig_V3.0\data\
     目录           0  2004-05-13 15:52  EasyKrig_V3.0\dataprep\
     文件        2865  2004-04-03 12:26  EasyKrig_V3.0\dataprep\check_unitsfig.m
     文件         434  2004-04-03 12:28  EasyKrig_V3.0\dataprep\check_unit_action.m
     文件        4982  2004-05-21 13:39  EasyKrig_V3.0\dataprep\datachk.m
     文件        7851  2004-06-17 17:48  EasyKrig_V3.0\dataprep\dataprep3d.m
     文件       22442  2004-07-20 09:22  EasyKrig_V3.0\dataprep\dataprep3dfig.m
     文件        2181  2004-04-03 12:31  EasyKrig_V3.0\dataprep\datareduction.m
     文件        1257  2004-04-03 12:31  EasyKrig_V3.0\dataprep\datatransform.m
     文件        6414  2004-08-02 13:31  EasyKrig_V3.0\dataprep\loaddatfile.m
     文件         519  2004-04-03 12:28  EasyKrig_V3.0\dataprep\loadvar_opt.m
     文件        1938  2004-07-12 13:34  EasyKrig_V3.0\dataprep\load_data_format_info.m
     文件        2418  2004-07-30 11:03  EasyKrig_V3.0\dataprep\save_data_format_info.m
     目录           0  2004-05-13 15:52  EasyKrig_V3.0\data\GLOBEC_gridfile\
     文件      135992  2004-06-24 14:07  EasyKrig_V3.0\data\GLOBEC_gridfile\globec_grid.dat
     目录           0  2004-05-13 15:53  EasyKrig_V3.0\data\pressure\
     文件        3900  2000-02-03 13:51  EasyKrig_V3.0\data\pressure\layer10.dat
     文件        3900  2000-02-03 13:51  EasyKrig_V3.0\data\pressure\layer11.dat
............此处省略98个文件信息