Databases
StackOverflow2013
Postsdbo
POphp uncompress & unpack gives array containing a few (9 in 10309) undefined values
 

Note that there are some explanatory texts on larger screens.

plurals
  1. POphp uncompress & unpack gives array containing a few (9 in 10309) undefined values
    primarykey
    data
    text
    <p>I have a binary data string that i uncompress and unpack into an array using php with the following code (complete code for this php page is included at the bottom of this question):</p> <pre><code>while ($row = mysql_fetch_array($result, MYSQL_NUM)) { $glycopeptide[$hits] = $row[1]; echo $row[4]; // $row[4] contains the binaryString $mz = base64_decode($row[4]); $unc_mz = gzuncompress($mz); $max = strlen($unc_mz); $counter = 0; for ($i = 0; $i &lt; $max; $i = $i+4) { $temp = substr($unc_mz,$i,4); $temp = unpack("f",$temp); $mz_array[$counter] = $temp[1]; $counter++; } $hits++; } </code></pre> <p>I find that both the mz_array (the above code, X-coords) has 9 undefined values (all at the end) but that the int_array (similar code, Y-coords) also has 9 undefined values (distributed throughout the whole array (not grouped or at the beginning/end).</p> <p>Here is an example of a small testing chunk that I added to my page:</p> <p>Test code:</p> <pre><code>for ($i = 0; $i &lt; $counter; $i++) { echo $i; echo " - "; echo $mz_array[$i]; echo " - "; echo $int_array[$i]; echo "&lt;br/&gt;"; } </code></pre> <p>A "selection" of the output (Note the missing values):</p> <pre><code>671 - 274.20001220703 - 429 672 - 274.39999389648 - 673 - 274.60000610352 - 1098 -- skipping a few lines -- 10299 - 2199.8000488281 - 0 10300 - 2200 - 0 10301 - - 0 10302 - - 0 </code></pre> <p>The weirdest part is that if i manually enter the strings in the complete/original code (see bottom of page) I get the undefined values while if I manually enter the string that is returned by doing "echo $row[4]" (contains the binaryString) in the following code it yields no undefined values. </p> <pre><code>&lt;?php $string = " /* Copy the string in the spoiler (on this page) here */ "; $int = base64_decode($string); $unc_int = gzuncompress($int); $max = strlen($unc_int); $counter = 0; $max_int = 0; for ($i = 0; $i &lt; $max; $i = $i + 4) { $temp= substr($unc_int,$i,4); $temp = unpack("f",$temp); $int_array[$counter] = $temp[1]; echo $counter; echo " -- "; echo $int_array[$counter]; echo "&lt;br/&gt;"; $counter++; } ?&gt; </code></pre> <p>Does anyone have any suggestions or idea's why this might occur?</p> <p>PS: Can someone add the tag gzuncompress? (I don't have the reputation for it).</p> <p><strong>EDIT 1</strong></p> <p>I have included an example binaryString (WARNING: Huge!)</p> <p>Y-Coords (retrieved by the commented echo $row[4] in the code):</p> <blockquote class="spoiler"> <p> eJztXQuQFMUZblQEQfEQUQSJYwRBSuLG+AAFHdmeRYOPS2lRgEgNvgpLhUuQkooYB5GXBLIaHwRQRwMpCqOuiciJhU4KjUoMnm98r2KUaKkXreBpSCXf3z3Lzu7Na3fn9nZhvqop2N2enu6///7fPcdYjBgxYsSIESNGjBhVgvFGkqmLOWPjk77tlNs5a35AY9kdnBlPo+3GMZ5tm87VmDULfXL/PnOwzkGfP0P7RLj2MUrAaNB0l2OtJuDzVu+188XJnKmnYZ0a5DpZE7Fu3TVm/BOft5TZZ4yqwvwf1u+ne89ea7wM/Pke5vvI3jHfGGVA56xxEnTWKeCT/h588i5k3cEp1vyuxtQN+H/an5+ya9GuBf2t6zi+s2ZAHh/KC+V7BTDXoK/15Y03uxD3zo1urmoKND4efXb17tPqgTZz0OYy/+eaB2uMHYP1UKOjVc1gDuy30TUuz3fWP82NJuy1nZyZx4GXDgO9M/acXgPd9w+mvXFrx66P+TLGNBnXBbXBB4tWwldYhvFcGmI8GbRJ1xGPhFjvPQXqvljDbZjv6jpan1pHfxcfjHyqm0Hra8vnraa+sGF68prmz4aeGtM3Y4wr3fd8I35nrxfZTU/g88YK/NbORAJ6owvG31Y49mwT/JJFXPrjYQBZSrGaAj97UMTrfFYJ4wkD9GUYtcuLBHU7Z/rzoO3dtT1OgQ0Y42neNqyFuSjrOcuamNMo7KPzecHvR45MsYtwKZfAFh7GBS9lv0LbH+BzH+67v7KvcBEDdNrarV0hb77Ddy9VWT9AFmRfwJqtqR+9pF8MGq/iwl7PId1bY62DQcMGB+0hF5W1th8WIO/a7a0TsN+maDK+M9qdnynWao3H76fVAb87oN8K2/s3XNDKfE5jzc0aW7sRNuYhkKM7bFpNt23g/XhhfBp6RtDKoW/UXmjzkaPNefg8j2Jj9UWXmgN4z7gnT0N1IuTRZ5yle4Avdzi+HwK+/0md+eNmsLwxTsT8z8zHyMvGDOzlbujnag//cmcyOtq11c8aNDCK5Ugb1/hvsmw5ps7GOl1dphycHs5fNS7EM26prjzJ/htz+rJ6tnLrnZC/74ePczadhX0/LZp4VWY59OcB4IdmXqALaf+JPNUAj2dUQBvql3KbhqKx6ceMZVNnjGWT1qeY/jZ00Z/xTMP2bbzkBPEudLA6ksaHfxVvf02HbaV3wfy2JJn+MeyvP+K5C8D3m5NMfRb3rfaYH+k+l/yseRXGuAT3XSHvI5tXjPk2Rz8Yi/kwl3GtiGyrxl6wff6FPr8oHK91v23D2dAX4/NDhW1Itxvbq8DPB/jwixM9bNodBT6Yg/EO5OHoRD7t6dgr/fnekSvNBM+RYvw6fHt1n0Iakk9GMXurL/e5uz3onpYbQeMr2vvbyrd5uTgtrbXn+xJgvJeUcbGPQuiAr5IF7VoHYF4Dg2M0jSTXBvHd/NaZoDiFscuWC6Ycz/LFGN8MXA/ityXyO30bZwmsQevJ2O8j/HM3Tlh3oI/BlcWtyJ43+uT7UFdCXn5ox418oDwlec28OYQtQDo/xF4XebkWjCdj10kkOGu8VGON+4I3t+WfYS7P047sZv3X0j/M3pm/t1oQ9LflnwIby/g95kA0acuPIzEUdDrdzn9UIo8z0pcjWzmsHWCcgfFcX2J8akNSxsx8oA5Bvydx4Tu7Pvd59PEkxvspZy3naywN3m7dCv7ejPGbEawR2ZLTi/oZxlnDZMgJkmMjbX1DvtNheb/WmMKlv/t3XJuSrvI2vQJy7gKN6eR7ebQhmGtBgwcwH/CtiAGG3Ldh0KhpgbVd5C8KOVkEfW7p+fPGQaDbqTaPusgTiq1YS6LNkxnvS1uJ9q1ym3ctnToTzx4V0s7IYWR7X6OxK+Z3anvfWdnIhV1Idp2i2nQfHPwsYY+XQg/y8buF039+0EdoMhaV4w1ar2uwt6bl+b4joHbDM36bfy7VNLa9jL09H3SdmV8f62aNKX/F5zvKH0ua6g5IbiX9+yBZKPiyWBaEhKihXBD+XtE2QIZnn0OfOzEmj1qXzFSs3zUdEDMjWkWZ24AMURt59PmSMIDOJx9O6CGv/TKssO6oBTaigfFSjILNDh7zCOJb2FBsZlHbIBscMsu4qjAm7WVvks3LniiKoUIO6W/lbUJXvIXf3uC79TDZQtZk7/h0tWEZLvy9PskSkJ+Jh2HvPIbx4uq+RmPTntbY2tUayxyF9XlX0kLUKOX2Nvxz6yJc2DeZE7EmvfIy2lwhbbrmp3DvX/D90iSzYAfp6xwxIuhoYfdh3T7ulmLKLLTVvW2THIxL0McrnE2cN5a1cFxDxrI+08ey759MsSOvSrF+A1OCD4WtSfo0rGz9hks9c7h/e5oDa/G2Lwgq5KjgoWbZhnyg3XH7jsS3oDOrTFdZgxw1EPAnjH7h47sm/Avj0/Y6nGpjjOv89YJxkFx74g/rStgW4zThPwqb3PZvyRYQ+ykn20h3UezGpz6/ZJyA/slugX2jki15LN+9jk60wjfIwm+opA6QbBi23Z+X/KCei/23gbMJy6DPj9Rk7p3WPoFxXxfen6Z8p+BPoiPuV8gf/yI8v6p4dsPlcm1EnsNhEyqQH8p4/Das0EYlH5T0soW92jJQkzbWljEiZiTqPF3Gnk3j+39wUZ/GJsJ3/B7/P1zOWe2LZ9+ej40akEGkW9iupOv6GZANxjP+sW2KW6rw+dkS2Y7sXWNT0bgg4422wu+Mn0OWLi3y36B3lCMwzy2gVXeMays+H4r//8iOk9l7lmL+nrLiJsyR4s+P+Iy7aO9TvaTwN5w8Rr4/aJ9QME7Ib2sSrpVcxDTYH3D9yvaN34QuwHibbsD6XSBjOSq+Ux5Fnzm+xzq0jJO5cZXiTI9hTX+J/iBLlc+4jNGvt+XJpfjtfFmjri9Ae/hz1ufyHr91MJs5S3+sse7Pgld+gX0HW8wgfgI9s/BxzPu5sAmdvkbD0WjDg2UX5S0S+2A9IK/3651iIwal2IT78JxzMKdN6LeFs22fYX/tj/5ecPQFHaAO8K+Nng39p8ImNsGzZj/015WLOAjRxNhfxhUaxqDf2fgNulZZhd9hY1ivc+HbsQ/Qrq+7/BGgfTYP+ncVnrPa234ScWdHrfMiPItkNbsBazTX5g/7zID6DL77BD45fNmh92oiLpLokWKH9kyxlgupzsKO6Tn4rAV+TBb6Wx0bTa5jb4A1GHvgXtBrVZ5elD8OistRTkbEzmZUl84UAyYbTVngbc9SnZdCtnkFdX1U6yLyRBMD9m0ftCPbJODcHcVGhV1eho2gfC5jweXaU43TNVlnEXDesGpYxvN6drytZwPyM+pwO14dkDc3P+TC34kyhk+1kyLvDTmpv4PxbgB/fUm2mXxGGvaguQKyGvrWSuZlkrEQ93UpLZ/ihyx0YgOHfBuKPfA1l/mQsHX+k3HfTFzQJ41vQk6+XVrswhVe+sABw+LijEu1c/ICdyVlLqKMOnfrP1jjRzSRpzA+8e5Dgc2oTnD4i21jAvep2sBZE/pW4ds2L9TYbOg38x061wl9e70mfEvqr61NE3tf1CDWE0LwRcWAjSXORDvWJTMJe3AO6KqXsC+qBOtJ2Fs97HO7FcZwd/cJ+4h9hcvK8V5S+DND18FWvhYyfwz5jnL/ZWAbmWTv5WL7F0m7K8qYPNValnPmRvhOcyLm8anBtarWG3bdWEC8vntaY4m+WLvvku31FPlhAbEaNyxaAD79U2VyUf0dF7kyMcfiOcDPMkLUzkUOet5jZT53Cvj1tnC2u/B3Ovr80vakjAV67FeRD/TaP3RPB9URUHyFfFFnzWcxWg9MgQc0lrlDY9vWaOzygSk26qgUa4Wvag6A/n8VY3fUATGKITnz31SnAfnBXnZ5RiXxQzdaboVsfBrzuU/GA8z+mqwndmlLdST6Jv+zAxQjJ5+4VDlrLQddMuh7scMXgm2qQ06oFJd9ycFzsKWF/2nLA1qTxBl2/QjFn0E/GqNoE5ADscjPyN3XySAbky0N9nuqBeNF8Ab8TLU7F7QnnzO9n50LyMnifljv4aB9z8IYhb4IvgDFtVSbF6gOkHyDWYVzI70s4oW5OjOKe1VaSxwh9FOohiVft2eMA00eLD1HquzEfbAvzfmFPEl+LsVwnXLXJN9+RXs+SJC9+0PH3oIPpE5x8SVesnMAJ3Dv8/1uc30e+w1rSblD8W4Ar3WA/UcxXVFDGJE9VQzSr1mi/bjgfFCtgOSQyBfcHa1NRXJh1KAUW/6Kxhou0dhJQ1Iy/nmxw6bsBNA7JJbPxXwvzNv8+i7w0HHwc5fh+zcxvtW23O7dPmaUOUjWtivYX+qxmoiHUg2bONNZj2cfw2CqowbOBcpd9L4pLmuqHTCSnZBDPs+n1uEWHvgeihj+oPOLYWuWIsEWHvjuMzp31kD78psq5Ko7EcKHKYF/jeGQa4/KXEnY97yJ+2wfU/gNz9TOGU7rBY2NOj3F+p2RYme/Ax97oSbrK+eBR5Z1zN6mvJf+YuEZCFFLU4ENoc+C3/A1L7lu0twMPZPqpBqpWgbldGGLCru9I+X7pqSM77rUwUaFpj6arDXzkmPQb8ZNPjquTnFg7xRTeqVE7YzhUfsQBehcVnHMkfRLhs54Jxy+AXxFytlMWEpnKBx2TG8X/yGELKD+ldWytkPkCTvIB6l5wAYzbsy/S5Tywgr5BfTeAkdtZLabXRcUUY1UtpWzxCBZy+CM9YvzpS45H/0hLmu0XGKGon58e8fvP7ULaPC4d33EbmTBVwei7SEydyj86Ffb62x1Yek6JweRO6mWzRcAqjvSH8d+erG9z+ELt7j7BvT1iZ0zc8TqyZ9TH/X3e1zHNj9Zc7mkGBECezHovZKegMwXccMw9lud6nfDeVbEicH2WZdSztXsSUhwz/dpVRPZD+xYVdD7JNqSsq7rLp8xU13s3NrwDWsRdG6q4pqZHIh3cFHdZ4FsqAY/Xcs9zxd1JuhcsMg5OM9jefC119ke5bXy6+zCgPJPmW5a4RnXHXYNW1TA+mTn82jfCxcG9PcDyH4e3gm8QTKn2M6i3M/RkG9t8CmWyZiXsN8c8knUSt1Ze7y8N0DUH5R5ZqSusML++xt7SCyazqQYr/Hy63RcYK0riqnWGyZGkPMf7VIfRqin9yxXE2SLmp1vQ5cD8ovS52rMPFNjVgTn8mN0LoTsqpH3uO/poNxeZH5UB8D6MZfngMvV998G5BMTdl1gBfrGOoTLOqw9PDbZcIQm/uZIqPe3Ua6wjt6rGRpkY+8NdnZnoMz30e2poPNzIl/rst+MBSFyV27IFMbYqR+j1T67sJci009jjbdo8t3ZPnk48e5RF79KvLOkBurFax3qlZxlNxbl/SBL0/SuruMr08H1DuVvVa79ixEjRowYMWJUDeZHnFln188ZphgxYsSIESNGjBgxYsSIESNGjBgxooB473oF783vaDT11Vj2njJz7jFixOhwiL/PuDE+Mxej5vF/owjA3A== </p> </blockquote> <p>Complete code:</p> <pre><code>&lt;?php require 'phplot/phplot.php'; $type = $_GET['type']; $gp = $_GET['gp']; $site = $_GET['site']; $prec = $_GET['prec']; $link = mysql_connect("localhost","reader","") or die (mysql_error()); mysql_select_db('leidenGlycoPeptide') or die (); $query = sprintf("select precursor.mzValue, glycoPeptide.protein, binaryDataArray.arrayLength, binaryDataArray.encodedLength, binaryDataArray.arrayData, precursor.chargeState, run.pepMass, run.PepSeq from glycoPeptide, spectrum, binaryDataArray, run, precursor where run.glycoPeptide = glycoPeptide.id AND spectrum.run = run.id AND precursor.run = run.id AND binaryDataArray.spectrum = spectrum.id AND precursor.id = spectrum.precursor AND spectrum.spectrum like 'm/z' AND precursor.mzValue like '%s' and glycoPeptide.protein like '%s' and run.glycoSite like '%s' and run.glycoType like '%s' ORDER by glycoPeptide.protein, spectrum.spectrum",(string)$prec, (string)$gp, (string)$site, (string)$type); $result = mysql_query($query); $hits = 0; while ($row = mysql_fetch_array($result, MYSQL_NUM)) { $charge = $row[5]; $pepmass = $row[6]; $pepseq = $row[7]; $glycopeptide[$hits] = $row[1]; /* Manually entering string here also gives undefined values */ /* $mz = " I was not able to include the mz string due to message size limit "; */ $mz = base64_decode($row[4]); $unc_mz = gzuncompress($mz); $max = strlen($unc_mz); $counter = 0; for ($i = 0; $i &lt; $max; $i = $i+4) { $temp = substr($unc_mz,$i,4); $temp = unpack("f",$temp); $mz_array[$counter] = $temp[1]; $counter++; } $hits++; } $query = sprintf("select precursor.mzValue, glycoPeptide.protein, binaryDataArray.arrayLength, binaryDataArray.encodedLength, binaryDataArray.arrayData from glycoPeptide, spectrum, binaryDataArray, run, precursor where run.glycoPeptide = glycoPeptide.id AND spectrum.run = run.id AND precursor.run = run.id AND binaryDataArray.spectrum = spectrum.id AND precursor.id = spectrum.precursor AND spectrum.spectrum like 'intensity' AND precursor.mzValue like '%s' and glycoPeptide.protein like '%s' and run.glycoSite like '%s' and run.glycoType like '%s' ORDER by glycoPeptide.protein, spectrum.spectrum",(string)$prec, (string)$gp, (string)$site, (string)$type); $result = mysql_query($query); while ($row = mysql_fetch_array($result, MYSQL_NUM)) { /* Manually entering string here also gives undefined values */ /* $int = " Copy the string from above in here "; */ $int = base64_decode($row[4]); /* The result from this is the above binaryString */ /* echo $row[4]; */ $unc_int = gzuncompress($int); $max = strlen($unc_int); $counter = 0; $max_int = 0; for ($i = 0; $i &lt; $max; $i = $i + 4) { $temp= substr($unc_int,$i,4); $temp = unpack("f",$temp); $int_array[$counter] = $temp[1]; $counter++; if ($temp[1] &gt; $max_int) { $max_int = $temp[1]; $counter++; } } } /* The following chunk is just to test the arrays */ for ($i = 0; $i &lt; $counter; $i++) { echo $i; echo " - "; echo $mz_array[$i]; echo " - "; echo $int_array[$i]; echo "&lt;br/&gt;"; } for ($i = 0; $i &lt; $counter; $i++) { $plot_data[$i] = array('',$mz_array[$i],$int_array[$i]); } // Plot the regular spectrum $width = 1024; $height = 768; $plot = new PHPlot($width,$height); $plot-&gt;SetMarginsPixels(NULL,NULL,NULL,35); $plot-&gt;SetPrintImage(False); $plot-&gt;SetPlotType('thinbarline'); //$plot-&gt;SetXTitle('m/z Values'); $plot-&gt;SetXTickAnchor('400'); $plot-&gt;SetDataColors('red'); $plot-&gt;SetXTickIncrement('200'); $plot-&gt;SetXDataLabelPos('none'); $plot-&gt;SetYTitle('Intensity'); $plot-&gt;SetYTickAnchor('0'); //Might need to define this dynamically with nested if/else loops $plot-&gt;SetYTickIncrement('100000'); $plot-&gt;SetDataType('data-data'); $plot-&gt;SetDataValues($plot_data); $plot-&gt;SetTitle('Fragmentation Spectrum'); //$plot-&gt;DrawGraph(); mysql_close($link); ?&gt; &lt;!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/html4/loose.dtd"&gt; &lt;html&gt; &lt;head&gt; &lt;title&gt;Query result page&lt;/title&gt; &lt;script src="jquery-1.9.1.min.js"&gt;&lt;/script&gt; &lt;/head&gt; &lt;script&gt; var gp = '&lt;?php echo htmlspecialchars($_GET['gp']); ?&gt;'; $(document).ready(function() { $('.button').click(function() { window.open('http://www.uniprot.org/uniprot/?query='+gp+'+AND+organism:human&amp;sort=score'); }); $('.XY').click(function() { var mz_array = &lt;?php echo json_encode($mz_array) ?&gt;; var int_array = &lt;?php echo json_encode($int_array) ?&gt;; var table = "&lt;table border=\"1\"&gt;" +"&lt;tr&gt;" +"&lt;th&gt;m/z&lt;/th&gt;" +"&lt;th&gt;intensity&lt;t/h&gt;" +"&lt;/tr&gt;"; var max = &lt;?php echo $counter ?&gt;; for (var i = 0; i &lt; max; i++) { table += "&lt;tr&gt;" +"&lt;td&gt;"+mz_array[i]+"&lt;/td&gt;" +"&lt;td&gt;"+int_array[i]+"&lt;/td&gt;" +"&lt;/tr&gt;"; } table += "&lt;/table&gt;"; var disp = window.open(); $(disp.document.body).html(table); }); }); &lt;/script&gt; &lt;body&gt; &lt;p&gt;The displayed spectrum belongs to &lt;?php echo $gp ?&gt; with a precursor [M+H] of &lt;?php echo (($prec*$charge)-($charge+1)); ?&gt;.&lt;br&gt; The peptide belonging to this glycopeptide has a mass of &lt;?php echo $pepmass ?&gt; and sequence: &lt;?php echo $pepseq ?&gt;.&lt;br&gt; &lt;button class="button"&gt;Uniprot search&lt;/button&gt; &lt;button class="XY"&gt;Display XY data&lt;/button&gt;&lt;/p&gt; &lt;img src="&lt;?php echo $plot-&gt;EncodeImage();?&gt;" alt="Plot Image"&gt; &lt;/body&gt; &lt;/html&gt; </code></pre>
    singulars
    1. PO
      singulars
      1. PTAnswer
    1. This table or related slice is empty.
    1. PTQuestion
    1. USBas Jansen
    1. USBas Jansen
    plurals
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. PO
      singulars
      1. PTAnswer
    2. PO
      singulars
      1. PTAnswer
    1. VO
      singulars
      1. POphp uncompress & unpack gives array containing a few (9 in 10309) undefined values
      1. This table or related slice is empty.
      1. VTUpMod
    2. VO
      singulars
      1. POphp uncompress & unpack gives array containing a few (9 in 10309) undefined values
      1. This table or related slice is empty.
      1. VTUpMod
    3. VO
      singulars
      1. POphp uncompress & unpack gives array containing a few (9 in 10309) undefined values
      1. This table or related slice is empty.
      1. VTUpMod
    1. COI have just added a example binaryString that I use and that gives me 9 undefined values (randomly throughout the array).
      singulars
      1. POphp uncompress & unpack gives array containing a few (9 in 10309) undefined values
      1. USBas Jansen
    2. COunable to reproduce from the base64 dump. Getting all `0` float values in `$mz_array`
      singulars
      1. POphp uncompress & unpack gives array containing a few (9 in 10309) undefined values
      1. USSparKot
    3. COStrange, I will double check tomorrow morning when I get back in the office
      singulars
      1. POphp uncompress & unpack gives array containing a few (9 in 10309) undefined values
      1. USBas Jansen
 

Querying!

 
Guidance
A row detail

Detail views are divided into sections. All the information in the data section comes from columns in the selected row. The other sections display data from other, related rows.

Related data can be related in a to-one or a to-many fashion. Captions of data related in a to-many fashion link to a list view showing a filtered view of the table.

Try moving around until you find a non-empty to-many entry and click on the label to get to one. You can move back to the root by clicking on the database name in the header.

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload