双线性插值用于Micaps

小小人

小小人 发表于 2014-1-20 16:44

                               
登录/注册后可看大图

嗯嗯，谢谢您提的建议。我在上面给出的X1,Y1是一个例子，它可以是任意四个经纬圈中的点，但是，这四个点必 ...

还有就是您说的已知四个站点插值一个站点，这个我还没有考虑过的。不过，我会努力的。

小小人

嗯嗯，谢谢大家的厚爱，我下去后有改进了下这个程序，是他可以直接在终端输入时显示坐标。源代码我把它贴出来。
program main
implicit none
integer::i,j
character*3 judge
real::x=101.87,y=25.73,x1(2),y1(2)
real::w(2,2),Vq(2,2),Vr(2),Vp,g,p
DATA(x1(i),i=1,2)/99.98,102.49/
DATA(y1(j),j=1,2)/25.02,27.52/
333  write(*,*)"Please enter the Value of Q"
do i=1,2
do j=1,2
write(*,"('Vq(',I1,',',I1,')=')")i,j
        read(*,*) Vq(i,j)
        end do
        end do
g=((x-x1(1))**2)+((x-x1(2))**2)
p=((y-y1(1))**2)+((y-y1(2))**2)
do i=1,2
w(1,i)=1-(((x-x1(i))**2)/g)
end do
do j=1,2
w(2,j)=1-(((y-y1(j))**2)/p)
end do

do i=1,2
Vr(i)=Vq(1,i)*w(1,1)+Vq(2,i)*w(1,2)
end do
Vp=Vr(1)*w(2,1)+Vr(2)*w(2,2)
write(*,*)  vp
!write(*,*)'是否需要继续？，是请按Y,其他键退出'
!  read(*,*) judge
  ! if(judge=="Y".or.judge=="y") then
   ! goto 333
   !endif

end
其中，被注释的那几句，我仿照了下论坛里面的，我想把它改为循环输入，

kc121

我不是很懂fortran，但我想

DATA(x1(i),i=1,2)/99.98,102.49/
DATA(y1(j),j=1,2)/25.02,27.52/

这里改一下，改成四个值如：
DATA(x1(i),i=1,4)/99,102,100,103/
DATA(y1(j),j=1,4)/25,27,24,26/
其他不变，这样是否能实现之前我说的任意四个站插值到一个站的问题呢？

小小人

嗯嗯，谢谢您的关注。这样的四个站点在我的这个程序里是不行，我在micaps里的插值，插的是由两条经纬线之间那个点，用图来表示就，它是在兰伯特投影下的四个点，所以在X方向上是在同一条经度，在y方向上是在同一条维度上。所以虽然有四个点，但是，他的值分别只有两个。

红河岸

LZ有木有不用输入这四个点，直接调用生成的几个网格点数据的方法？

lysx

贡献不够，先收藏了啊！

lysx

贡献不够，先收藏了啊！

qxtjls

呵呵，等待改进

daohua201314

谢谢楼主分享

2011wjj

好东东，值得分享啊！！！！！！！