高手帮看下这个问题

dahaoshanhe

有两个文件,
第一个文件如下,有三列,第一列为x坐标,第2列为y坐标,第3列为高程,用z表示.
文件2也有三列,第一列为x坐标,第2列为y坐标,第3列为待求高程.
利用文件1与文件2中的xy坐标得到文件2中的第3列,
规则为:如果文件2中的xy与文件1中的xy距离最近,也就是两点间距离最小
则把文件1中的第3列z的值赋给文件2中.
两点间距离公式:   distance=sqrt(x1-x2)^2+(y1-y2)^2

1. 文件1
   
2. x y   z
   
3. 1 2   34
   
4. 3 5   56
   
5. 6 7   20
   
6. 9 4   17
   
7. 
   
8. 文件2
   
9. x y    z
   
10. 1 3
    
11. 2 4  
    
12. 6 6
    
13. 
    
14. 
    
15. 输出结果:
    
16. x y    z
    
17. 1 3    34
    
18. 2 4    56
    
19. 6 6    20

复制代码

blackold

回复 1# dahaoshanhe

1. awk 'NR==FNR{a[$1" "$2]=$3;next}{d=-1;for(i in a){split(i,t);td=(t[1]-$1)^2+(t[2]-$2)^2;if(td<d||d<0){$3=a[i];d=td}}}1' file1.txt file2.txt
   

复制代码

dahaoshanhe

感谢黑哥 高手  佩服~~~谢啦   研读下你的思路~

dahaoshanhe

回复 2# blackold


    td=(t[1]-$1)^2+(t[2]-$2)^2;if(td<d||d<0)  这个没有看太明白  ,怎么没有用到sqrt这个呢?

blackold

回复 4# dahaoshanhe


    没有必要。

    a>b  等价于 sqrt(a) > sqrt(b)

hxre_cu

#!/bin/sh

x1=(`awk '{print $1}' file1`)
y1=(`awk '{print $2}' file1`)
z1=(`awk '{print $3}' file1`)

x2=(`awk '{print $1}' file2`)
y2=(`awk '{print $2}' file2`)

mindist=9999.999
echo x y z

for ((i=0;i<${#x2[*]};i=i+1))
do
        for ((j=0;j< ${#x1[*]};j=j+1))
        do
                x=$[${x2[$i]} - ${x1[$j]}]
                y=$[${y2[$i]} - ${y1[$j]}]
                z=$(expr "scale=3;sqrt ($[$x**2 + $y**2])"|bc)
                [[ $z > $mindist ]]   || z2=${z1[$j]}
                [[ $z > $mindist ]]   || mindist=$z
        done
        mindist=9999.999
        echo ${x2[$i]} ${y2[$i]} $z2
done

hxre_cu

回复 2# blackold

厉害，awk确实强大