Compare the difference of two Images.

Cited from https://rosettacode.org/wiki/Percentage_difference_between_images#Python.

 

I chose PYTHON script. But the script has some problem. I need to change “import Image” to “From PIL import Image”.

 

The output will give the pixel difference in percentage.

 

But I don’t understand the output. A more clear way is got from http://stackoverflow.com/questions/29229535/measure-similarity-of-two-images-in-java-or-imagemagick. In this method, the different pixel number will be output.

This used ImageMagick in Linux. The command to compare a.png and b.png is:

compare -metric ae a.png b.png null:

The output is the number of different pixel. If you want the percentage. You need to use the following command to get the X and Y for the image and use X*Y to calculate the total pixel. Then use the above result divided by the below result to get the percentage.

identify a.png

Advertisements

Use Jump Host for RSYNC

I have computer A, storage server B and HPCC.

Previously, I use A to SSH B and RSYNC files from HPCC to B. One day, B can not PING or SSH or RSYNC or SCP to HPCC. But B can PING or SSH or RSYNC or SCP to A, and A can also communicate with HPCC. So I am trying to use A as a jump host to help B to download files from HPCC.

 

Got help form https://www.freeture.ch/?p=815.

The first thing is make sure the SSH KEY works between A and B, and between A and HPCC.

In ~/.ssh directory on B, build a new file config. The content is below.

Host remote_server.com
HostName remote_server.com
User remote_server_user_id
ProxyCommand ssh jump_host_user_id@jump_host.com nc remote_server.com 22

Change VMD startup script.

I use iMac. The version of VMD is 1.9.3. There is a .vmdrc file at /Applications/VMD 1.9.3.app/Contents/vmd/

You can copy this file to your home directory ~ and then modify it.

What I want to modify is turn on all four lights, remove the axes and change the background color to gray.

The following modification is based on the original file from VMD.

  1. line 25: light 2 on
  2. line 26: light 3 on
  3. line 29: axes location off
  4. add color Display Background gray somewhere

Change VMD Movie Maker Default Setting.

I use it every week to make gif animation and show them in the group meeting. But the default setting is not proper for me. I use an iMac. The default working directory is lkjadfljaslkfjlasdfjlsadfdd……. (which is the system variable $TEMPDIR). The default animation format is MPEG-1. The default movement of molecules is Rock and Roll.

 

What I want is GIF, rotation along y axis and a specific working directory. If you need to generate tons of animation and sometimes VMD crashed, the repeatable setting will make me crazy.

 

Let’s work on this. The file is /Applications/VMD 1.9.3.app/Contents/vmd/plugins/noarch/tcl/vmdmovie1.9/vmdmovie.tcl

The first thing I change is the working directory in Line 63, delete frames and write to your specific directory. And delete line 130 to 147, this part will find a working directory according to different system.

 

The second thing is to change the movement of molecules. In line 67, change rockandroll to rotation.

 

The third thing is to change the animation format. In line 80, change ppmtompeg to imgif.

 

That’s all.

 

Enjoy~~

Calculating hydrogen bonds life time.

I was asked to do such an analysis that is calculating average lifetime of hydrogen bonds formed between certain atom and surrounding water molecules . And I found that AMBERTOOLS had such a tool called “lifetime” but the calculation seems wrong when you set up certain values of “fuzz”. And I think it’s also easy to do this by myself, so I used following method to do the calculation.

  1. You need a trajectory file. I choose a 5ns trajectory, with 1ps per frame.
  2. Find surrounding water molecules of certain atom.
  3. Calculate hydrogen bonds formed between these surrounding water molecules and the target atom.
  4. Calculate the average lifetime of all hydrogen bonds formed between the target atom and surrounding water molecules in the 5ns.

The following is the script for Step 2 and 3. O1 on RES is the target atom.

rm -f *_temp.dat
for res in RES; do
 for atom in O1; do
  rm -f ${res}_${atom}.dat
cpptraj ../*.prmtop << EOF
trajin *.dcd
mask "(:${res}@${atom}<:3.6)&:WAT" maskout mask_temp.dat
run
exit
EOF
  water_res=`cat mask_temp.dat | awk '{print $4}' | sort -n | uniq -c | sed /Res/d | awk '{print $2}' | xargs echo`

   for water in $water_res; do
cpptraj ../*prmtop << EOF
trajin *.dcd
hbond A acceptormask :${res}@${atom} donormask :${water}@O= dist 3.5 angle 130.0 out hb_temp.dat
hbond B acceptormask :${water}@O= donormask :${res}@${atom} dist 3.5 angle 130.0 out hb_temp.dat
run
C = A[UU] + B[UU]
writedata ${res}_${atom}_temp.dat C
quit
EOF
     cat ${res}_${atom}_temp.dat >> ${res}_${atom}.dat
    done

 done
done

rm -f *_temp.dat

 

The following is the script for Step 4. All result will output file “lifetime.dat”. The “fuzz” here is the same with that in AMBERTOOLS.

rm -f lifetime.dat
((fuzz=0))
for res in RES; do
 for atom in O1; do
  echo ${res}_${atom}
  count=(`cat ${res}_${atom}.dat | awk '{print $2}' | cut -f1 -d"." | uniq -c | awk '{print $1}'`)
  hb=(`cat ${res}_${atom}.dat | awk '{print $2}' | cut -f1 -d"." | uniq -c | awk '{print $2}'`)
  for ((i=0;i<${#count[*]};i++)); do
   if [ ${hb[$i]} == 0 ] && [ ${count[$i]} == $fuzz ];then hb[$i]=1; fi
   if [ ${hb[$i]} == C ];then hb[$i]=0; fi
   ((temp[$i] = ( ${count[$i]} * ${hb[$i]} )))
  done
  hb_num=`echo ${temp[*]} | tr ' ' '\n' | uniq -c | awk '{print $2}' | grep -vw 0 | wc -l`
  hb_sum=`echo ${temp[*]} | tr ' ' '\n' | awk '{s+=$1} END {print s}'`
  hb_avg=`echo "$hb_sum / $hb_num * 1.0" | bc -l | cut -c -5`
  echo ${res}_${atom} $hb_avg >> lifetime.dat
  unset temp count hb
 done
done