%反距离插值方法,使用ArcGIS很容易就办到,但是具体算法最近才学会
%师兄的程序,我要像师兄看齐~增大阅读面,多做数据,并且善于分析,我要多出成果!!!
function Z=IDW2015(data,station,lon,lat,e,num)
�ta数据 n3*n4(时间*站点数)
%station data 对应坐标的经纬度 n4*2 第一列纬度 第二列经度
%lon lat 需要插值格点的经纬度 数据格式 横着放(横着放,横着放~),1*n
%e参数 可设为2
%按变量选着 num为离插值点近的点 可设为4~12 ,也可按固定半径来选取点进行插值(没设计)
%e=2和num=12为arcgis的默认值
%例子
% DATA=IDW2015(P_m,latlon,lon',lat',2,12);
%
[n1 n2]=size(lon);
[n3 n4]=size(data);
% if station(1,2)>station(1,1) %%%%%%%%只适用于特定区域;目的:保证精度为x
x=
%师兄的程序,我要像师兄看齐~增大阅读面,多做数据,并且善于分析,我要多出成果!!!
function Z=IDW2015(data,station,lon,lat,e,num)
�ta数据 n3*n4(时间*站点数)
%station
%lon lat 需要插值格点的经纬度 数据格式 横着放(横着放,横着放~),1*n
%e参数 可设为2
%按变量选着 num为离插值点近的点 可设为4~12 ,也可按固定半径来选取点进行插值(没设计)
%e=2和num=12为arcgis的默认值
%例子
% DATA=IDW2015(P_m,latlon,lon',lat',2,12);
%
[n1 n2]=size(lon);
[n3 n4]=size(data);
% if station(1,2)>station(1,1) %%%%%%%%只适用于特定区域;目的:保证精度为x