このページの2つのバージョン間の差分を表示します。
| 両方とも前のリビジョン 前のリビジョン 次のリビジョン | 前のリビジョン |
| mitosearchデータ検証 [2022/03/17 10:57] – 133.11.222.89 | mitosearchデータ検証 [Unknown date] (現在) – 削除 - 外部編集 (Unknown date) 127.0.0.1 |
|---|
| ## 東京湾の魚の月平均表 | |
| |
| ``` | |
| for i in `more data/fish/lat-long-date.txt |grep N|grep E|awk '$2>35&&$2<35.69 && $4>139.6&&$4<140.2'|cut -f 1`; do cat db_fish/$i.input; done|awk -F'\t' '$0~"^20"{split($1,arr,"-"); m=arr[2]; n[m]++} $0!~"^20"{d[$1][m]+=$2} END{for(i in d){max=0; for(j in d[i]){d[i][j]=d[i][j]/n[m]; if(d[i][j]>max){max=d[i][j]}}; print max"\t"i"\t"d[i]["01"]"\t"d[i]["02"]"\t"d[i]["03"]"\t"d[i]["04"]"\t"d[i]["05"]"\t"d[i]["06"]"\t"d[i]["07"]"\t"d[i]["08"]"\t"d[i]["09"]"\t"d[i]["10"]"\t"d[i]["11"]"\t"d[i]["12"]}}'|sort -k1,1nr|awk -F'\t' '{ORS=""; print $1"\t"$2; for(i=3;i<=14;i++){print "\t"$i/$1*100}; print "\n"}' > /tmp/aa | |
| ``` | |
| |
| ## 陸か水かどうか | |
| |
| ``` | |
| [yoshitake@m50v251n3 yoshitake]$ more testmap.py | |
| import sys | |
| from staticmap import StaticMap | |
| map = StaticMap(50, 50) | |
| img = map.render(zoom=17, center=[float(sys.argv[1]), float(sys.argv[2])]) | |
| pixelSizeTuple = img.size | |
| water = 0 | |
| for i in range(pixelSizeTuple[0]): | |
| for j in range(pixelSizeTuple[1]): | |
| r,g,b = img.getpixel((i,j)) | |
| if r == 150 and g == 200 and b == 219: | |
| water = 1 | |
| |
| print(water) | |
| |
| [yoshitake@m50v251n3 yoshitake]$ pwd | |
| /home/yoshitake/yoshitake | |
| ``` | |
| |
| ``` | |
| more ~/mitosearch/Mitosearch/data/fish/lat-long-date.txt |grep N|grep E|awk '{print "echo "$1"; docker run -v $PWD:$PWD -w $PWD -i --rm python3-staticmap python testmap.py "$4" "$2}'|xargs -I{} bash -c "{}" > testmap.result | |
| ``` | |