克里金插值程序国外一个基金支持下开发的程序，赞

EasyKrig_V3.0_Matlab2012a.zip。是大师对之前的那个版本的最新修改，修改了很多的错误，而且计算效率有所提高，是我在网上搜索很久才找到的最好的、最有效的克里金插值程序。适合二维和三维的数据、大量的数据点的插值！

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  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\
     文件         453  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\bessel0.m
     文件        2442  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\degmins.m
     文件         845  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\ll2xy.m
     文件        2645  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\mapax.m
     文件        1029  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\mean_nan.m
     文件        1173  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\median_nan.m
     文件        1095  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\sort_mex.m
     文件        1028  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\std_nan.m
     文件         965  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\sum_nan.m
     文件        1462  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\bin\xy2ll.m
     文件        3127  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\color_info.dat
     目录           0  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\
     目录           0  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\GLOBEC_gridfile\
     文件      135992  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\GLOBEC_gridfile\globec_grid.dat
     目录           0  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\pressure\
     文件        3900  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\pressure\layer10.dat
     文件        3900  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\pressure\layer11.dat
     文件        3900  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\pressure\layer12.dat
     文件        3900  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\pressure\layer13.dat
     文件        3900  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\pressure\layer14.dat
     文件         135  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\pressure_batch_file.dat
     文件       18600  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\salinity.dat
     文件      161602  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\Temperature3d.dat
     文件         900  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\data\zooplankton.dat
     目录           0  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\dataprep\
     文件         434  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\dataprep\check_unit_action.m
     文件        2865  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\dataprep\check_unitsfig.m
     文件        4982  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\dataprep\datachk.m
     文件        7852  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\dataprep\dataprep3d.m
     文件       22444  2012-11-23 16:09  EasyKrig_V3.0_Matlab2012a\dataprep\dataprep3dfig.m
............此处省略99个文件信息