0

あるファイルのlatとlongの場所を、別のファイルのいくつかの名前付きフィールドに解決しようとしています。

このようなファイルが1つあります。

                 f1--f2--f3--------f4--------                f5---
 R               20175155 41273951N078593973W                18012              
 R               20175156 41274168N078593975W                18000              
 R               20175157 41274387N078593976W                17999              
 R               20175158 41274603N078593977W                18024              
 R               20175159 41274823N078593978W                18087

各文字は特定の場所にあるため、文字に基づいてフィールドを定義する必要があります。

f1 char 18-21; f2 char 22-25; f3 char 26-35; f4 char 36-45; f5文字62-66。

f3、f4、f5に対応するフィールド11、12、および13を持つ別のはるかに大きなcsvファイルがあります。

awk -F',' '{print $11, $12, $13}'
41.46703821 -078.98476926 519.21
41.46763555 -078.98477791 524.13
41.46824123 -078.98479015 526.67
41.46884129 -078.98480615 528.66
41.46943371 -078.98478482 530.50

ファイル2のフィールド11&&12;でファイル1のフィールド1&&2に最も近いものを見つける必要があります。

最も近い一致が見つかったら、ファイル1のフィールド1、2、3、4、5をファイル2のフィールド16、17、18、19、20に挿入する必要があります。

ご覧のとおり、形式は少し異なります。ファイル1は次のように分類されます。

ファイル1

f3-------f4--------

DDMMSSdd DDDMMSSdd 

41273951N078593973W

ファイル2

f11-------- f12---------

DD dddddddd  DDD dddddddd

41.46703821 -078.98476926

Nはf3が正の数であることを意味し、Wはf4が負の数であることを意味します。

私はファイル1をsedで変更しました、うまく機能するばかげた1つのライナー.. （より良い方法???）

cat $file1 |sed 's/.\{17\}//' |sed 's/\(.\{4\}\)\(.\{4\}\)\(.\{9\}\)\(.\)\(.\{9\}\)\(.\)\(.\{16\}\)\(.\{5\}\)/\1,\2,\3,\4,\5,\6,\8/'|sed 's/\(.\{10\}\)\(.\{3\}\)\(.\{2\}\)\(.\{2\}\)\(.\{2\}\)\(.\{3\}\)\(.\{3\}\)\(.\{2\}\)\(.*\)/\1\2,\3,\4.\5\6\7,\8\9/'|sed 's/\(.\{31\}\)\(.\{2\}\)\(.*\)/\1,\2.\3/'

2017,5155、41,27,39.51、N、078,59,39.73、W、18012
2017,5156、41,27,41.68、N、078,59,39.75、W、18000
2017,5157、41,27、 43.87、N、078,59,39.76、W、17999
2017,5158、41,27,46.03、N、078,59,39.77、W、18024
2017,5159、41,27,48.23、N、078,59、 39.78、W、18087

次に、形式を変換する必要があります。（これを解決しました（以下を参照）-問題-数値の四捨五入が多すぎます。小数点以下6桁以上にする必要があります。）

awk -F',' '{for (i=1;i<=NF;i++) {if (i <= 2) printf ($i","); else if (i == 3&&$6 == "S") printf("-"$3+($4/60)+($5/3600)","); else if (i == 3&&$6 == "N") printf($3+($4/60)+($5/3600)","); else if (i == 7&&$10 == "W") printf("-"$7+($8/60)+($9/3600)","); else if (i == 7&&$10 == "E") printf($7+($8/60)+($9/3600)","); if (i == 11) printf ($i"\n")}}'

2017,5155,41.461、-78.9944,18012
2017,5156,41.4616、-78.9944,18000
2017,5157,41.4622、-78.9944,17999
2017,5158,41.4628、-78.9944,18024
2017,5159,41.4634、-78.9944,18087

それが私がいるところです。

これを解決しました *この数式から小数点以下6桁以上の数値形式を取得する必要があります。*

printf（$ 3 +（$ 4/60）+（$ 5/3600）） 「％。8f」を追加

printf（ "％。8f"、$ 3 +（$ 4/60）+（$ 5/3600））

次の問題は、フィールドファイル1f3とf4をファイル2f11とf12の最も近い一致に一致させることです。

何か案は？

次に、フィールド間の距離を計算する必要があります。

Excelでは、式は次のようになります。

=ATAN2(COS(lat1)*SIN(lat2)-SIN(lat1)*COS(lat2)*COS(lon2-lon1), SIN(lon2-lon1)*COS(lat2))

その計算には何を使用できますか？

*更新---私は一致する場所の短い距離を見ています。ピタゴラスの定理のような単純なものを最も近い一致に適用することを考えていました。たぶん、小数点以下の桁数を減らすことさえできます。それは何倍も速くなければなりません。多分このようなもの.. *

x = (lon2-lon1) * Math.cos((lat1+lat2)/2);

y = (lat2-lat1);

d = Math.sqrt(x*x + y*y) * R;

次に、最終ファイルが更新された後、精度を高めるために必要な重い計算を行うことができます。

ありがとう

4

1 に答える 1

0

最も近い一致を実行した後、距離の計算を行うことはできません。最も近いとは、距離の値の比較によって定義されます。Awkは、必要な数式を評価できます（大円距離のように見えますか？）。この章を見て、必要なものを確認してください。

大きな問題は、最も近い一致を見つけることです。ファイル1の1行を取得し、追加の列を使用してファイル2の行を出力するawkスクリプトを記述します。その列は、距離の式に従ってポイントのペア間の距離を計算したものです。そのファイルを数値で並べ替えると（sort -n）、最も近いものが一番上になります。次に、ファイル1の各行をループし、awkスクリプトを呼び出しhead -n1、最も近い一致を引き出して、必要な形式で出力するスクリプトが必要です。

これはすべてbashとawkで可能ですが、Pythonでははるかに単純なスクリプトになります。どちらを好むかによります。

于 2012-07-08T09:34:52.733 に答える