克里金插值法MATLAB工具箱

克里金插值法MATLAB工具箱，亲测可以用，附相关论文。

function	initialization（）
% function initialization perform initialization

%%  Kriging Software Package  version 3.0   May 1 2004
%%  Copyright （c） 1999 2001 2004 property of Dezhang Chu and Woods Hole Oceanographic
%%  Institution.  All Rights Reserved.

%% This program is distributed in the hope that it will be useful but WITHOUT ANY
%% WARRANTY.

%% Contact Dr. Chu at dchu@whoi.edu with enhancements or suggestions for changes.

%% Modified on 11-16-2016 to remove an inconsistent definition for the new Matlab version

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=12;
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.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.datapre

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-02-01 03:06  EasyKrig_V3.0_Matlab2016b\
     目录           0  2016-11-17 00:42  EasyKrig_V3.0_Matlab2016b\bin\
     文件         453  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\bin\bessel0.m
     文件        2442  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\bin\degmins.m
     文件         845  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\bin\ll2xy.m
     文件        2934  2016-11-17 00:34  EasyKrig_V3.0_Matlab2016b\bin\mapax.m
     文件        1029  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\bin\mean_nan.m
     文件        1173  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\bin\median_nan.m
     文件        1095  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\bin\sort_mex.m
     文件        1028  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\bin\std_nan.m
     文件         965  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\bin\sum_nan.m
     文件        1462  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\bin\xy2ll.m
     文件        3127  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\color_info.dat
     目录           0  2013-10-18 01:52  EasyKrig_V3.0_Matlab2016b\data\
     目录           0  2013-11-01 21:52  EasyKrig_V3.0_Matlab2016b\data\GLOBEC_gridfile\
     文件       80059  2013-11-01 21:52  EasyKrig_V3.0_Matlab2016b\data\GLOBEC_gridfile\globec_grid_lat_lon.txt
     文件      135992  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\data\GLOBEC_gridfile\globec_grid_lon_lat.dat
     目录           0  2013-05-01 02:46  EasyKrig_V3.0_Matlab2016b\data\pressure\
     文件        3900  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\data\pressure\layer10.dat
     文件        3900  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\data\pressure\layer11.dat
     文件        3900  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\data\pressure\layer12.dat
     文件        3900  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\data\pressure\layer13.dat
     文件        3900  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\data\pressure\layer14.dat
     文件         135  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\data\pressure_batch_file.dat
     文件       18600  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\data\salinity.dat
     文件      161602  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\data\Temperature3d.dat
     文件         900  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\data\zooplankton.dat
     目录           0  2013-05-01 02:46  EasyKrig_V3.0_Matlab2016b\dataprep\
     文件        2865  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\dataprep\check_unitsfig.m
     文件         434  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\dataprep\check_unit_action.m
     文件        4982  2012-11-24 07:09  EasyKrig_V3.0_Matlab2016b\dataprep\datachk.m
............此处省略111个文件信息