| 両方とも前のリビジョン 前のリビジョン 次のリビジョン | 前のリビジョン |
| mitosearchデータ検証 [2022/06/03 04:52] – 118.240.79.152 | 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 mapwater.py | |
| import sys | |
| from staticmap import StaticMap | |
| map = StaticMap(50, 50, url_template='http://a.tile.openstreetmap.org/{z}/{x}/{y}.png') | |
| 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 == 170 and g == 211 and b == 223): | |
| water = 1 | |
| |
| print(water) | |
| |
| [yoshitake@m50v251n3 yoshitake]$ pwd | |
| /home/yoshitake/yoshitake | |
| ``` | |
| |
| ``` | |
| more ~/mitosearch/Mitosearch/data/fish/lat-long-date.txt |awk -F'\t' '$2~"^[0-9]"{split($2,arr," "); if(arr[2]=="N"){k=arr[1]}else{k=-arr[1]}; if(arr[4]=="E"){k2=arr[3]}else{k2=-arr[3]}; print $1"\t"k"\t"k2}'|awk '{print "echo "$1"; docker run -v $PWD:$PWD -w $PWD -i --rm python3-staticmap python mapwater.py "$3" "$2}'|xargs -I{} bash -c "{}" > mapwater.result | |
| |
| cat mapwater.result|paste - - > mapwater.result.txt | |
| ``` | |
| |
| ## ANEMOMNEデータベース | |
| |
| ``` | |
| wget --http-user=username --http-passwd=pass -r -np --wait=1 https://db.anemone.bio/dist/MiFish/ --no-check-certificate | |
| find ../db.anemone.bio/|grep fastq.xz|while read i; do ln -s $i `echo $i|sed 's%../db.anemone.bio/dist/MiFish/ANEMONE/%%; s%/%___%g'`; done | |
| find ../db.anemone.bio/|grep sample.tsv.xz|while read i; do ln -s $i `echo $i|sed 's%../db.anemone.bio/dist/MiFish/ANEMONE/%%; s%/%___%g'`; done | |
| ``` | |