PDBLoader.js 6.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217
  1. /**
  2. * @author alteredq / http://alteredqualia.com/
  3. * @author Mugen87 / https://github.com/Mugen87
  4. */
  5. THREE.PDBLoader = function ( manager ) {
  6. this.manager = ( manager !== undefined ) ? manager : THREE.DefaultLoadingManager;
  7. };
  8. THREE.PDBLoader.prototype = {
  9. constructor: THREE.PDBLoader,
  10. load: function ( url, onLoad, onProgress, onError ) {
  11. var scope = this;
  12. var loader = new THREE.FileLoader( scope.manager );
  13. loader.setPath( scope.path );
  14. loader.load( url, function ( text ) {
  15. onLoad( scope.parse( text ) );
  16. }, onProgress, onError );
  17. },
  18. setPath: function ( value ) {
  19. this.path = value;
  20. return this;
  21. },
  22. // Based on CanvasMol PDB parser
  23. parse: function ( text ) {
  24. function trim( text ) {
  25. return text.replace( /^\s\s*/, '' ).replace( /\s\s*$/, '' );
  26. }
  27. function capitalize( text ) {
  28. return text.charAt( 0 ).toUpperCase() + text.substr( 1 ).toLowerCase();
  29. }
  30. function hash( s, e ) {
  31. return 's' + Math.min( s, e ) + 'e' + Math.max( s, e );
  32. }
  33. function parseBond( start, length ) {
  34. var eatom = parseInt( lines[ i ].substr( start, length ) );
  35. if ( eatom ) {
  36. var h = hash( satom, eatom );
  37. if ( bhash[ h ] === undefined ) {
  38. bonds.push( [ satom - 1, eatom - 1, 1 ] );
  39. bhash[ h ] = bonds.length - 1;
  40. } else {
  41. // doesn't really work as almost all PDBs
  42. // have just normal bonds appearing multiple
  43. // times instead of being double/triple bonds
  44. // bonds[bhash[h]][2] += 1;
  45. }
  46. }
  47. }
  48. function buildGeometry() {
  49. var build = {
  50. geometryAtoms: new THREE.BufferGeometry(),
  51. geometryBonds: new THREE.BufferGeometry(),
  52. json: {
  53. atoms: atoms,
  54. bonds: bonds
  55. }
  56. };
  57. var geometryAtoms = build.geometryAtoms;
  58. var geometryBonds = build.geometryBonds;
  59. var i, l;
  60. var verticesAtoms = [];
  61. var colorsAtoms = [];
  62. var verticesBonds = [];
  63. // atoms
  64. for ( i = 0, l = atoms.length; i < l; i ++ ) {
  65. var atom = atoms[ i ];
  66. var x = atom[ 0 ];
  67. var y = atom[ 1 ];
  68. var z = atom[ 2 ];
  69. verticesAtoms.push( x, y, z );
  70. var r = atom[ 3 ][ 0 ] / 255;
  71. var g = atom[ 3 ][ 1 ] / 255;
  72. var b = atom[ 3 ][ 2 ] / 255;
  73. colorsAtoms.push( r, g, b );
  74. }
  75. // bonds
  76. for ( i = 0, l = bonds.length; i < l; i ++ ) {
  77. var bond = bonds[ i ];
  78. var start = bond[ 0 ];
  79. var end = bond[ 1 ];
  80. verticesBonds.push( verticesAtoms[ ( start * 3 ) + 0 ] );
  81. verticesBonds.push( verticesAtoms[ ( start * 3 ) + 1 ] );
  82. verticesBonds.push( verticesAtoms[ ( start * 3 ) + 2 ] );
  83. verticesBonds.push( verticesAtoms[ ( end * 3 ) + 0 ] );
  84. verticesBonds.push( verticesAtoms[ ( end * 3 ) + 1 ] );
  85. verticesBonds.push( verticesAtoms[ ( end * 3 ) + 2 ] );
  86. }
  87. // build geometry
  88. geometryAtoms.addAttribute( 'position', new THREE.Float32BufferAttribute( verticesAtoms, 3 ) );
  89. geometryAtoms.addAttribute( 'color', new THREE.Float32BufferAttribute( colorsAtoms, 3 ) );
  90. geometryBonds.addAttribute( 'position', new THREE.Float32BufferAttribute( verticesBonds, 3 ) );
  91. return build;
  92. }
  93. var CPK = { h: [ 255, 255, 255 ], he: [ 217, 255, 255 ], li: [ 204, 128, 255 ], be: [ 194, 255, 0 ], b: [ 255, 181, 181 ], c: [ 144, 144, 144 ], n: [ 48, 80, 248 ], o: [ 255, 13, 13 ], f: [ 144, 224, 80 ], ne: [ 179, 227, 245 ], na: [ 171, 92, 242 ], mg: [ 138, 255, 0 ], al: [ 191, 166, 166 ], si: [ 240, 200, 160 ], p: [ 255, 128, 0 ], s: [ 255, 255, 48 ], cl: [ 31, 240, 31 ], ar: [ 128, 209, 227 ], k: [ 143, 64, 212 ], ca: [ 61, 255, 0 ], sc: [ 230, 230, 230 ], ti: [ 191, 194, 199 ], v: [ 166, 166, 171 ], cr: [ 138, 153, 199 ], mn: [ 156, 122, 199 ], fe: [ 224, 102, 51 ], co: [ 240, 144, 160 ], ni: [ 80, 208, 80 ], cu: [ 200, 128, 51 ], zn: [ 125, 128, 176 ], ga: [ 194, 143, 143 ], ge: [ 102, 143, 143 ], as: [ 189, 128, 227 ], se: [ 255, 161, 0 ], br: [ 166, 41, 41 ], kr: [ 92, 184, 209 ], rb: [ 112, 46, 176 ], sr: [ 0, 255, 0 ], y: [ 148, 255, 255 ], zr: [ 148, 224, 224 ], nb: [ 115, 194, 201 ], mo: [ 84, 181, 181 ], tc: [ 59, 158, 158 ], ru: [ 36, 143, 143 ], rh: [ 10, 125, 140 ], pd: [ 0, 105, 133 ], ag: [ 192, 192, 192 ], cd: [ 255, 217, 143 ], in: [ 166, 117, 115 ], sn: [ 102, 128, 128 ], sb: [ 158, 99, 181 ], te: [ 212, 122, 0 ], i: [ 148, 0, 148 ], xe: [ 66, 158, 176 ], cs: [ 87, 23, 143 ], ba: [ 0, 201, 0 ], la: [ 112, 212, 255 ], ce: [ 255, 255, 199 ], pr: [ 217, 255, 199 ], nd: [ 199, 255, 199 ], pm: [ 163, 255, 199 ], sm: [ 143, 255, 199 ], eu: [ 97, 255, 199 ], gd: [ 69, 255, 199 ], tb: [ 48, 255, 199 ], dy: [ 31, 255, 199 ], ho: [ 0, 255, 156 ], er: [ 0, 230, 117 ], tm: [ 0, 212, 82 ], yb: [ 0, 191, 56 ], lu: [ 0, 171, 36 ], hf: [ 77, 194, 255 ], ta: [ 77, 166, 255 ], w: [ 33, 148, 214 ], re: [ 38, 125, 171 ], os: [ 38, 102, 150 ], ir: [ 23, 84, 135 ], pt: [ 208, 208, 224 ], au: [ 255, 209, 35 ], hg: [ 184, 184, 208 ], tl: [ 166, 84, 77 ], pb: [ 87, 89, 97 ], bi: [ 158, 79, 181 ], po: [ 171, 92, 0 ], at: [ 117, 79, 69 ], rn: [ 66, 130, 150 ], fr: [ 66, 0, 102 ], ra: [ 0, 125, 0 ], ac: [ 112, 171, 250 ], th: [ 0, 186, 255 ], pa: [ 0, 161, 255 ], u: [ 0, 143, 255 ], np: [ 0, 128, 255 ], pu: [ 0, 107, 255 ], am: [ 84, 92, 242 ], cm: [ 120, 92, 227 ], bk: [ 138, 79, 227 ], cf: [ 161, 54, 212 ], es: [ 179, 31, 212 ], fm: [ 179, 31, 186 ], md: [ 179, 13, 166 ], no: [ 189, 13, 135 ], lr: [ 199, 0, 102 ], rf: [ 204, 0, 89 ], db: [ 209, 0, 79 ], sg: [ 217, 0, 69 ], bh: [ 224, 0, 56 ], hs: [ 230, 0, 46 ], mt: [ 235, 0, 38 ], ds: [ 235, 0, 38 ], rg: [ 235, 0, 38 ], cn: [ 235, 0, 38 ], uut: [ 235, 0, 38 ], uuq: [ 235, 0, 38 ], uup: [ 235, 0, 38 ], uuh: [ 235, 0, 38 ], uus: [ 235, 0, 38 ], uuo: [ 235, 0, 38 ] };
  94. var atoms = [];
  95. var bonds = [];
  96. var histogram = {};
  97. var bhash = {};
  98. var x, y, z, index, e;
  99. // parse
  100. var lines = text.split( '\n' );
  101. for ( var i = 0, l = lines.length; i < l; i ++ ) {
  102. if ( lines[ i ].substr( 0, 4 ) === 'ATOM' || lines[ i ].substr( 0, 6 ) === 'HETATM' ) {
  103. x = parseFloat( lines[ i ].substr( 30, 7 ) );
  104. y = parseFloat( lines[ i ].substr( 38, 7 ) );
  105. z = parseFloat( lines[ i ].substr( 46, 7 ) );
  106. index = parseInt( lines[ i ].substr( 6, 5 ) ) - 1;
  107. e = trim( lines[ i ].substr( 76, 2 ) ).toLowerCase();
  108. if ( e === '' ) {
  109. e = trim( lines[ i ].substr( 12, 2 ) ).toLowerCase();
  110. }
  111. atoms[ index ] = [ x, y, z, CPK[ e ], capitalize( e ) ];
  112. if ( histogram[ e ] === undefined ) {
  113. histogram[ e ] = 1;
  114. } else {
  115. histogram[ e ] += 1;
  116. }
  117. } else if ( lines[ i ].substr( 0, 6 ) === 'CONECT' ) {
  118. var satom = parseInt( lines[ i ].substr( 6, 5 ) );
  119. parseBond( 11, 5 );
  120. parseBond( 16, 5 );
  121. parseBond( 21, 5 );
  122. parseBond( 26, 5 );
  123. }
  124. }
  125. // build and return geometry
  126. return buildGeometry();
  127. }
  128. };