分享我写FORTRAN程序（反距离权重法插值）的一些步骤和心得

max020729

卤煮 反距离插值不应该是距离的-2次方么 我看你程序怎么还是2次方呢

zheng

楼主你计算的距离是用度数表示的，而不是空间距离（m或者km），这样是不是不太妥当

咕噜

max020729 发表于 2017-1-2 20:37
卤煮 反距离插值不应该是距离的-2次方么 我看你程序怎么还是2次方呢

谢谢发现问题，稍后改正！

咕噜

zheng 发表于 2017-1-3 12:50
楼主你计算的距离是用度数表示的，而不是空间距离（m或者km），这样是不是不太妥当

确实没有考虑好，谢谢指正！

zljln00

此处必须有掌声！加油

草木尽欲言哉

路过，学习一下~~~~~

taoyiwei

感谢分享！！！！！！！！！！！！

xdzhao2013yk

没有改进版的么？点在点上，点在线上不用管吗？

HAHAHAH

楼主，你好，最近刚好要用反距离插值法把格点数据插值到站点上，但输出数据不对，自己找不出问题，可以帮我指点下吗？
代码如下：
program id_interp
implicit none

!声明变量
real lat,lon,d,dtotal,temp(141,81,6)
integer i,j,k,irec
integer,parameter::nx=141,ny=81
real :: disdum(141,81,6)
real :: valstn(6)
!声明结束

!开始程序
open(11,file='e:\test\1981_1982.grd',form='unformatted',access='direct',recl=141*81*1)
open(12,file='e:\test\c52446.txt',form='formatted')
read(*,*) lon,lat
irec=0
do k=1,6
  irec=irec+1
  read(11,rec=irec) ((temp(i,j,k),i=1,141),j=1,81)
enddo

do k=1,6 !! 需要给出时间循环个数
dtotal=0.
disdum(:,:,:)=0.
do j=1,81
  do i=1,141
    if (abs((i-1)*0.5+70-lon)<=0.5.and.abs((j-1)*0.5+15-lat)<=0.5) then
        d=1/(((i-1)*0.5+70-lon)**2+((j-1)*0.5+15-lat)**2)
           disdum(i,j,k)=d
            dtotal=dtotal+d
    endif
  enddo
enddo
  
do j=1,81
  do i=1,141
    if(disdum(i,j,k)>0.) then
      valstn(k)=temp(i,j,k)*(disdum(i,j,k)/dtotal)     
    endif
  enddo
enddo

enddo !! 时间循环结束

write(12,*) (valstn(k),k=1,6)

end
!程序结束

万分感谢！！！

xdzhao2013yk

谢谢，打扰了