Mercurial > repos > galaxyp > msi_ion_images

changeset 5:a851b4e8fba7 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/msi_ion_images commit a7be47698f53eb4f00961192327d93e8989276a7
author galaxyp
date Mon, 11 Jun 2018 17:31:51 -0400
parents 729a8bf3ffa9
children 1a7cc5e14430
files msi_ion_images.xml test-data/Example_Continuous.ibd test-data/Example_Continuous.imzML test-data/Heatmaps_LM8_file16.pdf test-data/Heatmaps_analyze75.pdf test-data/Heatmaps_imzml.pdf test-data/Heatmaps_rdata.pdf
diffstat 7 files changed, 411 insertions(+), 450 deletions(-) [+]
line wrap: on
line diff
--- a/msi_ion_images.xml	Mon May 28 12:34:05 2018 -0400
+++ b/msi_ion_images.xml	Mon Jun 11 17:31:51 2018 -0400
@@ -1,6 +1,6 @@
-<tool id="mass_spectrometry_imaging_ion_images" name="MSI ion images" version="1.10.0.0">
+<tool id="mass_spectrometry_imaging_ion_images" name="MSI ion images" version="1.10.0.1">
     <description>
-        mass spectrometry imaging heatmaps
+        mass spectrometry imaging m/z heatmaps
     </description>
     <requirements>
         <requirement type="package" version="1.10.0">bioconductor-cardinal</requirement>
@@ -45,9 +45,9 @@
 
 ###################################### file properties in numbers ##############
 
-## Number of features (mz)
+## Number of features (m/z)
 maxfeatures = length(features(msidata))
-## Range mz
+## Range m/z
 minmz = round(min(mz(msidata)), digits=2)
 maxmz = round(max(mz(msidata)), digits=2)
 ## Number of spectra (pixels)
@@ -64,7 +64,7 @@
 medint = round(median(spectra(msidata)[]), digits=2)
 ## Number of intensities > 0
 npeaks= sum(spectra(msidata)[]>0)
-## Spectra multiplied with mz (potential number of peaks)
+## Spectra multiplied with m/z (potential number of peaks)
 numpeaks = ncol(spectra(msidata)[])*nrow(spectra(msidata)[])
 ## Percentage of intensities > 0
 percpeaks = round(npeaks/numpeaks*100, digits=2)
@@ -103,43 +103,34 @@
     peakpickinginfo=processinginfo@peakPicking
 }
 
-##################################### read and filter input masses ##############
+##################################### read and filter input m/z ##############
 
 input_list = read.delim("$massfile", header = FALSE, stringsAsFactors = FALSE)
 
-### in case input file had only one column with mz values but not names, duplicate mz values and use as names:
+### in case input file had only one column with m/z values but not names, duplicate m/z values and use as names:
+if (ncol(input_list) == 1){
+    input_list = cbind(input_list, input_list)}
 
-if (ncol(input_list) == 1)
-{
-    input_list = cbind(input_list, input_list)
-}
+### calculate how many input m/z are valid: 
+inputmasses = input_list[input_list[,1]>minmz & input_list[,1]<maxmz,]
+inputmz = as.numeric(inputmasses[,1])
+inputnames = as.character(inputmasses[,2])
 
-### calculate how many input masses are valid: 
-inputmasses = input_list[input_list[,1]>minmz & input_list[,1]<maxmz,]
-
-inputmz = inputmasses[,1]
-inputnames = inputmasses[,2]
 
-if (length(inputmz) == 1)
-{
-    countpixels = sum(spectra(msidata)[features(msidata, mz = inputmz), ] >0)
-    percentpixels = round(countpixels/pixelcount*100, digits=1)
-    valuesdataframe = cbind(inputmz, cbind(countpixels, percentpixels))
-    write.table(valuesdataframe, file="$pixel_count", quote = FALSE, row.names = FALSE, col.names=TRUE, sep = "\t")
-}else if (length(inputmz) >1) {
-    countpixels = rowSums(spectra(msidata)[features(msidata, mz=inputmz),] >0)
-    percentpixels = round(countpixels/pixelcount*100, digits=1)
-    valuesdataframe = cbind(inputmz, cbind(countpixels, percentpixels))
-    write.table(valuesdataframe, file="$pixel_count", quote = FALSE, row.names = FALSE, col.names=TRUE, sep = "\t")
-}else{
-    valuesdataframe = data.frame(0,0)
-    write.table(valuesdataframe, file="$pixel_count", quote = FALSE, row.names = FALSE, col.names=TRUE, sep = "\t")
-}
+############################## PDF #############################################
+################################################################################
 
-############################ summarize file properties in numbers ##############
+pdf("heatmaps.pdf", fonts = "Times", pointsize = 12)
+plot(0,type='n',axes=FALSE,ann=FALSE)
+#if not $filename:
+    #set $filename = $infile.display_name
+#end if
+title(main=paste("\nHeatmap images\n\n", "Filename:\n", "$filename"))
 
-properties = c("Number of mz features",
-               "Range of mz values [Da]",
+############################# I) numbers ####################################
+
+properties = c("Number of m/z features",
+               "Range of m/z values [Da]",
                "Number of pixels", 
                "Range of x coordinates", 
                "Range of y coordinates",
@@ -153,7 +144,7 @@
                "Baseline reduction",
                "Peak picking",
                "Centroided", 
-                paste0("# valid masses in \n", "$massfile.display_name"))
+                paste0("# valid m/z in \n", "$massfile.display_name"))
 
 values = c(paste0(maxfeatures), 
            paste0(minmz, " - ", maxmz), 
@@ -174,38 +165,53 @@
 
 property_df = data.frame(properties, values)
 
-
-############################## PDF #############################################
-
-pdf("heatmaps.pdf", fonts = "Times", pointsize = 12)
-plot(0,type='n',axes=FALSE,ann=FALSE)
-#if not $filename:
-    #set $filename = $infile.display_name
-#end if
-title(main=paste("\nHeatmap images\n\n", "Filename:\n", "$filename"))
-
-############################# I) numbers ####################################
-
 grid.table(property_df, rows= NULL)
 
 ############################# II) images ####################################
 
-### only plot images when file has peaks and valid input mz: 
+### only plot images when file has peaks and valid input m/z: 
+
+if (npeaks > 0){
+    if (length(inputmz) != 0){
+      for (mass in 1:length(inputmz)){
+
+###standard image
+
+#if str($image_cond.image_type) == "standard_image":
+    print("standard image")
+
+    print(image(msidata, mz=inputmz[mass],plusminus = $plusminus_dalton, contrast.enhance = "$image_contrast", 
+    smooth.image = "$image_smoothing", strip=$strip, colorkey=$colorkey,
+    main= paste0(inputnames[mass], " (", round(inputmz[mass], digits = 2)," ± ", $plusminus_dalton, " Da)")))}
+
+
+###lattice image
 
-if (npeaks > 0)
-{
-    if (length(inputmz) != 0)
-    {
-      for (mass in 1:length(inputmz))
-      {
-      print(image(msidata, mz=inputmz[mass], strip = strip.custom(bg="lightgrey", par.strip.text=list(col="black", cex=.9)),
-              lattice=TRUE, plusminus = $plusminus_dalton, contrast.enhance = "$image_contrast", smooth.image = "$image_smoothing", 
-              main= paste0(mass, ") ", inputnames[mass], " (", round(inputmz[mass], digits = 2)," ± ", $plusminus_dalton, " Da)")))
-      }
-    } else {print("The input masses were invalid")}
+#elif str($image_cond.image_type) == "lattice_image":
+    print("lattice image")
+
+    #if str($strip) =="TRUE": 
+
+    print(image(msidata, mz=inputmz[mass], strip = strip.custom(bg="lightgrey", par.strip.text=list(col="black", cex=.9)),
+    lattice=TRUE, plusminus = $plusminus_dalton, contrast.enhance = "$image_contrast", smooth.image = "$image_smoothing", 
+    colorkey=$colorkey,
+    main= paste0(inputnames[mass], " (", round(inputmz[mass], digits = 2)," ± ", $plusminus_dalton, " Da)")))}
+
+    #elif str($strip) =="FALSE": 
+
+    print(image(msidata, mz=inputmz[mass], strip = $strip,
+    lattice=TRUE, plusminus = $plusminus_dalton, contrast.enhance = "$image_contrast", smooth.image = "$image_smoothing", 
+    colorkey=$colorkey,
+    main= paste0(inputnames[mass], " (", round(inputmz[mass], digits = 2)," ± ", $plusminus_dalton, " Da)")))}
+
+    #end if
+
+#end if 
+
+    } else {print("The input m/z were invalid")}
     dev.off()
 }else{
-  print("inputfile has no intensities > 0")
+    print("inputfile has no intensities > 0")
 dev.off()
 }
     ]]></configfile>
@@ -213,9 +219,9 @@
     <inputs>
         <param name="infile" type="data" format="imzml,rdata,analyze75" label="Inputfile as imzML, Analyze7.5 or Cardinal MSImageSet saved as RData"
             help="Upload composite datatype imzml (ibd+imzML) or analyze75 (hdr+img+t2m) or regular upload .RData (Cardinal MSImageSet)"/>
-        <param name="filename" type="text" value="" label="Title" help="will appear in the quality report. If nothing given it will take the dataset name"/>
-        <param name="massfile" type="data" format="tabular" label="Tabular file with masses and names"
-            help="first column mass (m/z), second column mass name, tab separated file"/>
+        <param name="filename" type="text" value="" label="Title" help="will appear in the pdf output. If nothing given it will take the dataset name"/>
+        <param name="massfile" type="data" format="tabular" label="Tabular file with m/z and names"
+            help="first column m/z, second column m/z name, tab separated file"/>
         <param name="image_contrast" type="select" label="Select a contrast enhancement function for the heatmap images" help="The 'histogram' equalization method flatterns the distribution of intensities. The hotspot 'suppression' method uses thresholding to reduce the intensities of hotspots">
             <option value="none" selected="True">none</option>
             <option value="suppression">suppression</option>
@@ -226,11 +232,20 @@
             <option value="gaussian">gaussian</option>
             <option value="adaptive">adaptive</option>
         </param>
-        <param name="plusminus_dalton" value="0.25" type="float" label="Mass range" help="plusminus mass window in Dalton"/>
+        <param name="plusminus_dalton" value="0.25" type="float" label="M/z range" help="plusminus m/z window in Dalton"/>
+        <param name="strip" type="boolean" checked="True" display="radio" truevalue="TRUE" falsevalue="FALSE" label="Display m/z value in plot"/>
+        <param name="colorkey" type="boolean" checked="True" display="radio"  truevalue="TRUE" falsevalue="FALSE" label="Display colorkey in plot"/>
+        <conditional name="image_cond">
+            <param name="image_type" type="select" label="Select the image type">
+                <option value="standard_image" selected="True">standard</option>
+                <option value="lattice_image">lattice</option>
+            </param>
+        <when value="standard_image"/>
+        <when value="lattice_image"/>
+        </conditional>
     </inputs>
     <outputs>
-        <data format="pdf" name="plots" from_work_dir="heatmaps.pdf" label = "${tool.name} ${on_string}"/>
-        <data format="tabular" name="pixel_count" label="Number of peaks (intensity > 0) per mz"/>
+        <data format="pdf" name="plots" from_work_dir="heatmaps.pdf" label = "$infile.display_name heatmaps"/>
     </outputs>
     <tests>
         <test>
@@ -242,8 +257,10 @@
             <param name="plusminus_dalton" value="0.25"/>
             <param name="filename" value="Testfile_imzml"/>
             <param name="image_contrast" value="histogram"/>
+            <param name="strip" value="True"/>
+            <param name="colorkey" value="True"/>
+            <param name="image_type" value="lattice_image"/>
             <output name="plots" file="Heatmaps_imzml.pdf" compare="sim_size" delta="20000"/>
-            <output name="pixel_count" file="tabular_imzml.tabular"/>
         </test>
         <test>
             <param name="infile" value="" ftype="analyze75">
@@ -255,24 +272,28 @@
             <param name="plusminus_dalton" value="0.5"/>
             <param name="filename" value="Testfile_analyze75"/>
             <param name="image_smoothing" value="gaussian"/>
+            <param name="strip" value="False"/>
+            <param name="colorkey" value="True"/>
             <output name="plots" file="Heatmaps_analyze75.pdf" compare="sim_size" delta="20000"/>
-            <output name="pixel_count" file="tabular_analyze75.tabular"/>
         </test>
         <test>
             <param name="infile" value="preprocessed.rdata" ftype="rdata"/>
             <param name="massfile" value="inputpeptides.tabular" ftype="tabular"/>
             <param name="plusminus_dalton" value="0.5"/>
+            <param name="strip" value="True"/>
+            <param name="colorkey" value="True"/>
+            <param name="image_type" value="lattice_image"/>
             <param name="filename" value="Testfile_rdata"/>
             <output name="plots" file="Heatmaps_rdata.pdf" compare="sim_size" delta="20000"/>
-            <output name="pixel_count" file="tabular_rdata.tabular"/>
         </test>
         <test>
             <param name="infile" value="empty_spectra.rdata" ftype="rdata"/>
             <param name="massfile" value="inputpeptides2.tabular" ftype="tabular"/>
             <param name="plusminus_dalton" value="0.5"/>
+            <param name="strip" value="True"/>
+            <param name="colorkey" value="False"/>
             <param name="filename" value="Testfile_rdata"/>
             <output name="plots" file="Heatmaps_LM8_file16.pdf" compare="sim_size" delta="20000"/>
-            <output name="pixel_count" file="tabular_LM8file16.tabular"/>
         </test>
     </tests>
     <help><![CDATA[
@@ -280,32 +301,32 @@
 
 Cardinal is an R package that implements statistical & computational tools for analyzing mass spectrometry imaging datasets. `More information on Cardinal <http://cardinalmsi.org//>`_
 
-This tool uses the Cardinal image function to plot the intensity distribution of interesting masses of mass-spectrometry imaging data. 
+This tool uses the Cardinal image function to plot the intensity distribution of interesting m/z of mass spectrometry imaging data. 
 Input data: 
 
-3 types of mass-spectrometry imaging data can be used:
+3 types of mass spectrometry imaging data can be used:
 
 - imzml file (upload imzml and ibd file via the "composite" function) `Introduction to the imzml format <https://ms-imaging.org/wp/imzml/>`_
 - Analyze7.5 (upload hdr, img and t2m file via the "composite" function)
 - Cardinal "MSImageSet" data (with variable name "msidata", saved as .RData)
 
-Tabular file with masses:
+Tabular file with m/z:
 
 - tab separated file (.tabular), datatype in Galaxy must be tabular otherwise file will not appear in selection window (if Galaxy auto-detection was wrong, datatype can be changed by pressing button with the pen (edit attributes))
-- first column must contain masses (separate point numbers by point, not comma)
-- optionally a second column with names for the masses can be provided
+- first column must contain m/z (separate point numbers by point, not comma)
+- optionally a second column with names for the m/z can be provided
 - no empty fields or letters are allowed in the first column
 
 Output:
 
 - Pdf with the heatmap images
-- Tabular with masses that were in the mass range and their occurence over all pixels (absolute and in %)
 
 Troubleshooting:
 
 - no heatmaps are plotted when tabular file doesn't fulfill the criteria described above
 - no heatmaps are plotted when the input mass spectrometry imaging file has no intensities > 0
-- out of thetabular file only masses with > 1.5-2% pixel coverage can be used with the contrast enhance and image smoothing functions, as both crash when a mass has not enough intensity values
+- the contrast enhance and image smoothing functions require a certain number of m/z with intensities > 0 (empirical value > 2% of spectra)
+- the standard image function should work for all files while the lattice function works not on every file (nicely)
 
     ]]>
     </help>
Binary file test-data/Example_Continuous.ibd has changed
--- a/test-data/Example_Continuous.imzML	Mon May 28 12:34:05 2018 -0400
+++ b/test-data/Example_Continuous.imzML	Mon Jun 11 17:31:51 2018 -0400
@@ -1,373 +1,313 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
-<mzML xmlns="http://psi.hupo.org/ms/mzml" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://psi.hupo.org/ms/mzml http://psidev.info/files/ms/mzML/xsd/mzML1.1.0_idx.xsd" version="1.1">
-  <cvList count="3">
-    <cv id="MS" fullName="Proteomics Standards Initiative Mass Spectrometry Ontology" version="1.3.1" URI="http://psidev.info/ms/mzML/psi-ms.obo"/>
-    <cv id="UO" fullName="Unit Ontology" version="1.15" URI="http://obo.cvs.sourceforge.net/obo/obo/ontology/phenotype/unit.obo"/>
-    <cv id="IMS" fullName="Imaging MS Ontology" version="0.9.1" URI="http://www.maldi-msi.org/download/imzml/imagingMS.obo"/>
-  </cvList>
-  <fileDescription>
-    <fileContent>
-      <cvParam cvRef="MS" accession="MS:1000579" name="MS1 spectrum" value=""/>
-      <cvParam cvRef="MS" accession="MS:1000128" name="profile spectrum" value=""/>
-      <cvParam cvRef="IMS" accession="IMS:1000080" name="universally unique identifier" value="{554A27FA-79D2-4766-9A2C-862E6D78B1F3}"/>
-      <cvParam cvRef="IMS" accession="IMS:1000091" name="ibd SHA-1" value="A5BE532D25997B71BE6D20C76561DDC4D5307DDD"/>
-      <cvParam cvRef="IMS" accession="IMS:1000030" name="continuous" value=""/>
-    </fileContent>
-    <sourceFileList count="1">
-      <sourceFile id="sf1" name="Example.raw" location="C:\Users\Thorsten Schramm\Documents\Promotion\imzML\Website\files\Beispiel-Dateien\Example images\">
-        <cvParam cvRef="MS" accession="MS:1000563" name="Thermo RAW file" value=""/>
-        <cvParam cvRef="MS" accession="MS:1000768" name="Thermo nativeID format" value=""/>
-        <cvParam cvRef="MS" accession="MS:1000569" name="SHA-1" value="7623BE263B25FF99FDF017154B86FAB742D4BB0B"/>
-      </sourceFile>
-    </sourceFileList>
-    <contact>
-      <cvParam cvRef="MS" accession="MS:1000586" name="contact name" value="Thorsten Schramm"/>
-      <cvParam cvRef="MS" accession="MS:1000590" name="contact organization" value="Institut für Anorganische und Analytische Chemie"/>
-      <cvParam cvRef="MS" accession="MS:1000587" name="contact address" value="Schubertstraße 60, Haus 16, Gießen, Germany"/>
-      <cvParam cvRef="MS" accession="MS:1000589" name="contact email" value="thorsten.schramm@anorg.chemie.uni-.giessen.de"/>
-    </contact>
-  </fileDescription>
-  <referenceableParamGroupList count="4">
-    <referenceableParamGroup id="mzArray">
-      <cvParam cvRef="MS" accession="MS:1000576" name="no compression" value=""/>
-      <cvParam cvRef="MS" accession="MS:1000514" name="m/z array" value="" unitCvRef="MS" unitAccession="MS:1000040" unitName="m/z"/>
-      <cvParam cvRef="IMS" accession="IMS:1000101" name="external data" value="true"/>
-      <cvParam cvRef="MS" accession="MS:1000521" name="32-bit float" value=""/>
-    </referenceableParamGroup>
-    <referenceableParamGroup id="intensityArray">
-      <cvParam cvRef="MS" accession="MS:1000576" name="no compression" value=""/>
-      <cvParam cvRef="MS" accession="MS:1000515" name="intensity array" value="" unitCvRef="MS" unitAccession="MS:1000131" unitName="number of counts"/>
-      <cvParam cvRef="IMS" accession="IMS:1000101" name="external data" value="true"/>
-      <cvParam cvRef="MS" accession="MS:1000521" name="32-bit float" value=""/>
-    </referenceableParamGroup>
-    <referenceableParamGroup id="scan1">
-      <cvParam cvRef="MS" accession="MS:1000093" name="increasing m/z scan" value=""/>
-      <cvParam cvRef="MS" accession="MS:1000095" name="linear" value=""/>
-      <cvParam cvRef="MS" accession="MS:1000512" name="filter string" value="ITMS - p NSI Full ms [100,00-800,00]"/>
-    </referenceableParamGroup>
-    <referenceableParamGroup id="spectrum1">
-      <cvParam cvRef="MS" accession="MS:1000579" name="MS1 spectrum" value=""/>
-      <cvParam cvRef="MS" accession="MS:1000511" name="ms level" value="0"/>
-      <cvParam cvRef="MS" accession="MS:1000128" name="profile spectrum" value=""/>
-      <cvParam cvRef="MS" accession="MS:1000129" name="negative scan" value=""/>
-    </referenceableParamGroup>
-  </referenceableParamGroupList>
-  <sampleList count="1">
-    <sample id="sample1" name="Sample1">
-      <cvParam cvRef="MS" accession="MS:1000001" name="sample number" value="1"/>
-    </sample>
-  </sampleList>
-  <softwareList count="2">
-    <software id="Xcalibur" version="2.2">
-      <cvParam cvRef="MS" accession="MS:1000532" name="Xcalibur" value=""/>
-    </software>
-    <software id="TMC" version="1.1 beta">
-      <cvParam cvRef="MS" accession="MS:1000799" name="custom unreleased software tool" value=""/>
-    </software>
-  </softwareList>
-  <scanSettingsList count="1">
-    <scanSettings id="scansettings1">
-      <cvParam cvRef="IMS" accession="IMS:1000401" name="top down" value=""/>
-      <cvParam cvRef="IMS" accession="IMS:1000413" name="flyback" value=""/>
-      <cvParam cvRef="IMS" accession="IMS:1000480" name="horizontal line scan" value=""/>
-      <cvParam cvRef="IMS" accession="IMS:1000491" name="linescan left right" value=""/>
-      <cvParam cvRef="IMS" accession="IMS:1000042" name="max count of pixel x" value="3"/>
-      <cvParam cvRef="IMS" accession="IMS:1000043" name="max count of pixel y" value="3"/>
-      <cvParam cvRef="IMS" accession="IMS:1000044" name="max dimension x" value="300" unitCvRef="UO" unitAccession="UO:0000017" unitName="micrometer"/>
-      <cvParam cvRef="IMS" accession="IMS:1000045" name="max dimension y" value="300" unitCvRef="UO" unitAccession="UO:0000017" unitName="micrometer"/>
-      <cvParam cvRef="IMS" accession="IMS:1000046" name="pixel size x" value="100" unitCvRef="UO" unitAccession="UO:0000017" unitName="micrometer"/>
-      <cvParam cvRef="IMS" accession="IMS:1000047" name="pixel size y" value="100" unitCvRef="UO" unitAccession="UO:0000017" unitName="micrometer"/>
-      <cvParam cvRef="MS" accession="MS:1000836" name="dried dropplet" value=""/>
-      <cvParam cvRef="MS" accession="MS:1000835" name="matrix solution concentration" value="10"/>
-      <cvParam cvRef="MS" accession="MS:1000834" name="matrix solution" value="DHB"/>
-    </scanSettings>
-  </scanSettingsList>
-  <instrumentConfigurationList count="1">
-    <instrumentConfiguration id="LTQFTUltra0">
-      <cvParam cvRef="MS" accession="MS:1000557" name="LTQ FT Ultra" value=""/>
-      <cvParam cvRef="MS" accession="MS:1000529" name="instrument serial number" value="none"/>
-      <componentList count="3">
-        <source order="1">
-          <cvParam cvRef="MS" accession="MS:1000073" name="electrospray ionization" value=""/>
-          <cvParam cvRef="MS" accession="MS:1000485" name="nanospray inlet" value=""/>
-          <cvParam cvRef="MS" accession="MS:1000843" name="wavelength" value="337"/>
-          <cvParam cvRef="MS" accession="MS:1000844" name="focus diameter x" value="10"/>
-          <cvParam cvRef="MS" accession="MS:1000845" name="focus diameter y" value="10"/>
-          <cvParam cvRef="MS" accession="MS:1000846" name="pulse energy" value="10"/>
-          <cvParam cvRef="MS" accession="MS:1000847" name="pulse duration" value="10"/>
-          <cvParam cvRef="MS" accession="MS:1000848" name="attenuation" value="50"/>
-          <cvParam cvRef="MS" accession="MS:1000850" name="gas laser" value=""/>
-          <cvParam cvRef="IMS" accession="IMS:1000202" name="target material" value="Conductive Glas"/>
-        </source>
-        <analyzer order="2">
-          <cvParam cvRef="MS" accession="MS:1000264" name="ion trap" value=""/>
-          <cvParam cvRef="MS" accession="MS:1000014" name="accuracy" value="0" unitCvRef="MS" unitAccession="MS:1000040" unitName="m/z"/>
-        </analyzer>
-        <detector order="3">
-          <cvParam cvRef="MS" accession="MS:1000253" name="electron multiplier" value=""/>
-          <cvParam cvRef="MS" accession="MS:1000120" name="transient recorder" value=""/>
-        </detector>
-      </componentList>
-      <softwareRef ref="Xcalibur"/>
-    </instrumentConfiguration>
-  </instrumentConfigurationList>
-  <dataProcessingList count="2">
-    <dataProcessing id="XcaliburProcessing">
-      <processingMethod order="1" softwareRef="Xcalibur">
-        <cvParam cvRef="MS" accession="MS:1000594" name="low intensity data point removal" value=""/>
-      </processingMethod>
-    </dataProcessing>
-    <dataProcessing id="TMCConversion">
-      <processingMethod order="2" softwareRef="TMC">
-        <cvParam cvRef="MS" accession="MS:1000544" name="Conversion to mzML" value=""/>
-      </processingMethod>
-    </dataProcessing>
-  </dataProcessingList>
-  <run defaultInstrumentConfigurationRef="LTQFTUltra0" defaultSourceFileRef="sf1" id="Experiment01" sampleRef="sample1" startTimeStamp="2009-08-11T15:59:44">
-    <spectrumList count="9" defaultDataProcessingRef="XcaliburProcessing">
-      <spectrum id="Scan=1" defaultArrayLength="0" index="0">
-        <referenceableParamGroupRef ref="spectrum1"/>
-        <scanList count="1">
-          <cvParam cvRef="MS" accession="MS:1000795" name="no combination" value=""/>
-          <scan instrumentConfigurationRef="LTQFTUltra0">
-            <referenceableParamGroupRef ref="scan1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="1"/>
-          </scan>
-        </scanList>
-        <binaryDataArrayList count="2">
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="mzArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="intensityArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="33612"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-        </binaryDataArrayList>
-      </spectrum>
-      <spectrum id="Scan=2" defaultArrayLength="0" index="1">
-        <referenceableParamGroupRef ref="spectrum1"/>
-        <scanList count="1">
-          <cvParam cvRef="MS" accession="MS:1000795" name="no combination" value=""/>
-          <scan instrumentConfigurationRef="LTQFTUltra0">
-            <referenceableParamGroupRef ref="scan1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="2"/>
-            <cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="1"/>
-          </scan>
-        </scanList>
-        <binaryDataArrayList count="2">
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="mzArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="intensityArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="67208"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-        </binaryDataArrayList>
-      </spectrum>
-      <spectrum id="Scan=3" defaultArrayLength="0" index="2">
-        <referenceableParamGroupRef ref="spectrum1"/>
-        <scanList count="1">
-          <cvParam cvRef="MS" accession="MS:1000795" name="no combination" value=""/>
-          <scan instrumentConfigurationRef="LTQFTUltra0">
-            <referenceableParamGroupRef ref="scan1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="3"/>
-            <cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="1"/>
-          </scan>
-        </scanList>
-        <binaryDataArrayList count="2">
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="mzArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="intensityArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="100804"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-        </binaryDataArrayList>
-      </spectrum>
-      <spectrum id="Scan=4" defaultArrayLength="0" index="3">
-        <referenceableParamGroupRef ref="spectrum1"/>
-        <scanList count="1">
-          <cvParam cvRef="MS" accession="MS:1000795" name="no combination" value=""/>
-          <scan instrumentConfigurationRef="LTQFTUltra0">
-            <referenceableParamGroupRef ref="scan1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="2"/>
-          </scan>
-        </scanList>
-        <binaryDataArrayList count="2">
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="mzArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="intensityArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="134400"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-        </binaryDataArrayList>
-      </spectrum>
-      <spectrum id="Scan=5" defaultArrayLength="0" index="4">
-        <referenceableParamGroupRef ref="spectrum1"/>
-        <scanList count="1">
-          <cvParam cvRef="MS" accession="MS:1000795" name="no combination" value=""/>
-          <scan instrumentConfigurationRef="LTQFTUltra0">
-            <referenceableParamGroupRef ref="scan1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="2"/>
-            <cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="2"/>
-          </scan>
-        </scanList>
-        <binaryDataArrayList count="2">
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="mzArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="intensityArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="167996"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-        </binaryDataArrayList>
-      </spectrum>
-      <spectrum id="Scan=6" defaultArrayLength="0" index="5">
-        <referenceableParamGroupRef ref="spectrum1"/>
-        <scanList count="1">
-          <cvParam cvRef="MS" accession="MS:1000795" name="no combination" value=""/>
-          <scan instrumentConfigurationRef="LTQFTUltra0">
-            <referenceableParamGroupRef ref="scan1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="3"/>
-            <cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="2"/>
-          </scan>
-        </scanList>
-        <binaryDataArrayList count="2">
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="mzArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="intensityArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="201592"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-        </binaryDataArrayList>
-      </spectrum>
-      <spectrum id="Scan=7" defaultArrayLength="0" index="6">
-        <referenceableParamGroupRef ref="spectrum1"/>
-        <scanList count="1">
-          <cvParam cvRef="MS" accession="MS:1000795" name="no combination" value=""/>
-          <scan instrumentConfigurationRef="LTQFTUltra0">
-            <referenceableParamGroupRef ref="scan1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="3"/>
-          </scan>
-        </scanList>
-        <binaryDataArrayList count="2">
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="mzArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="intensityArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="235188"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-        </binaryDataArrayList>
-      </spectrum>
-      <spectrum id="Scan=8" defaultArrayLength="0" index="7">
-        <referenceableParamGroupRef ref="spectrum1"/>
-        <scanList count="1">
-          <cvParam cvRef="MS" accession="MS:1000795" name="no combination" value=""/>
-          <scan instrumentConfigurationRef="LTQFTUltra0">
-            <referenceableParamGroupRef ref="scan1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="2"/>
-            <cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="3"/>
-          </scan>
-        </scanList>
-        <binaryDataArrayList count="2">
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="mzArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="intensityArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="268784"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-        </binaryDataArrayList>
-      </spectrum>
-      <spectrum id="Scan=9" defaultArrayLength="0" index="8">
-        <referenceableParamGroupRef ref="spectrum1"/>
-        <scanList count="1">
-          <cvParam cvRef="MS" accession="MS:1000795" name="no combination" value=""/>
-          <scan instrumentConfigurationRef="LTQFTUltra0">
-            <referenceableParamGroupRef ref="scan1"/>
-            <cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="3"/>
-            <cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="3"/>
-          </scan>
-        </scanList>
-        <binaryDataArrayList count="2">
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="mzArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-          <binaryDataArray encodedLength="0">
-            <referenceableParamGroupRef ref="intensityArray"/>
-            <cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="8399"/>
-            <cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="302380"/>
-            <cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="33596"/>
-            <binary/>
-          </binaryDataArray>
-        </binaryDataArrayList>
-      </spectrum>
-    </spectrumList>
-  </run>
-</mzML>
\ No newline at end of file
+<?xml version="1.0"?>
+<mzML version="1.1" xmlns="http://psi.hupo.org/ms/mzml" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://psi.hupo.org/ms/mzml http://psidev.info/files/ms/mzML/xsd/mzML1.1.0_idx.xsd">
+	<cvList count="3">
+		<cv id="MS" fullName="Proteomics Standards Initiative Mass Spectrometry Ontology" version="1.3.1" URI="http://psidev.info/ms/mzML/psi-ms.obo" />
+		<cv id="UO" fullName="Unit Ontology" version="1.15" URI="http://obo.cvs.sourceforge.net/obo/obo/ontology/phenotype/unit.obo" />
+		<cv id="IMS" fullName="Imaging MS Ontology" version="0.9.1" URI="http://www.maldi-msi.org/download/imzml/imagingMS.obo" />
+	</cvList>
+	<fileDescription>
+		<fileContent>
+			<cvParam cvRef="MS" accession="MS:1000579" name="MS1 spectrum" value="" />
+			<cvParam cvRef="MS" accession="MS:1000128" name="profile spectrum" value="" />
+			<cvParam cvRef="IMS" accession="IMS:1000080" name="universally unique identifier" value="{51BB7C6F-9974-4626-B35F-5B65547BAE6B}" />
+			<cvParam cvRef="IMS" accession="IMS:1000091" name="ibd SHA-1" value="26C26A63DCF7CD1AD768A67AEDB7131C5D88C397" />
+			<cvParam cvRef="IMS" accession="IMS:1000030" name="continuous" value="" />
+		</fileContent>
+	</fileDescription>
+	<referenceableParamGroupList count="4">
+		<referenceableParamGroup id="spectrum1">
+			<cvParam cvRef="MS" accession="MS:1000579" name="MS1 spectrum" value="" />
+			<cvParam cvRef="MS" accession="MS:1000511" name="ms level" value="0" />
+			<cvParam cvRef="MS" accession="MS:1000128" name="profile spectrum" value="" />
+		</referenceableParamGroup>
+		<referenceableParamGroup id="scan1">
+			<cvParam cvRef="MS" accession="MS:1000093" name="increasing m/z scan" value="" />
+			<cvParam cvRef="MS" accession="MS:1000095" name="linear" value="" />
+		</referenceableParamGroup>
+		<referenceableParamGroup id="mzArray">
+			<cvParam cvRef="MS" accession="MS:1000576" name="no compression" value="" />
+			<cvParam cvRef="MS" accession="MS:1000514" name="m/z array" value="" unitCvRef="MS" unitAccession="MS:1000040" unitName="m/z" />
+			<cvParam cvRef="IMS" accession="IMS:1000101" name="external data" value="true" />
+			<cvParam cvRef="MS" accession="MS:1000521" name="32-bit float" value="32-bit float" />
+		</referenceableParamGroup>
+		<referenceableParamGroup id="intensityArray">
+			<cvParam cvRef="MS" accession="MS:1000576" name="no compression" value="" />
+			<cvParam cvRef="MS" accession="MS:1000515" name="intensity array" value="" unitCvRef="MS" unitAccession="MS:1000131" unitName="number of counts" />
+			<cvParam cvRef="IMS" accession="IMS:1000101" name="external data" value="true" />
+			<cvParam cvRef="MS" accession="MS:1000521" name="32-bit float" value="32-bit float" />
+		</referenceableParamGroup>
+	</referenceableParamGroupList>
+	<sampleList count="1">
+		<sample id="sample1" name="Sample1">
+			<cvParam cvRef="MS" accession="MS:1000001" name="sample number" value="1" />
+		</sample>
+	</sampleList>
+	<softwareList count="1">
+		<software id="Cardinal" version="1.12.0">
+			<cvParam cvRef="MS" accession="MS:1000799" name="custom unreleased software tool" value="" />
+		</software>
+	</softwareList>
+	<scanSettingsList count="1">
+		<scanSettings id="scansettings1">
+			<cvParam cvRef="IMS" accession="IMS:1000042" name="max count of pixel x" value="3" />
+			<cvParam cvRef="IMS" accession="IMS:1000043" name="max count of pixel y" value="3" />
+		</scanSettings>
+	</scanSettingsList>
+	<instrumentConfigurationList count="1">
+		<instrumentConfiguration id="IC1" />
+	</instrumentConfigurationList>
+	<dataProcessingList count="1">
+		<dataProcessing id="CardinalWriteImzML">
+			<processingMethod order="1" softwareRef="Cardinal">
+				<cvParam cvRef="MS" accession="MS:1000544" name="Conversion to mzML" value="" />
+			</processingMethod>
+		</dataProcessing>
+	</dataProcessingList>
+	<run defaultInstrumentConfigurationRef="IC1" id="Experiment01" sampleRef="sample1">
+		<spectrumList count="9" defaultDataProcessingRef="CardinalWriteImzML">
+			<spectrum id="Spectrum=1" defaultArrayLength="0" index="1">
+				<referenceableParamGroupRef ref="spectrum1" />
+				<scanList count="1">
+					<cvParam cvRef="MS" accession="MS:1000795" name="no combination" value="" />
+					<scan instrumentConfigurationRef="IC1">
+						<referenceableParamGroupRef ref="scan1" />
+						<cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="1" />
+						<cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="1" />
+					</scan>
+				</scanList>
+				<binaryDataArrayList count="2">
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="mzArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="intensityArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="4812" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+				</binaryDataArrayList>
+			</spectrum>
+			<spectrum id="Spectrum=2" defaultArrayLength="0" index="2">
+				<referenceableParamGroupRef ref="spectrum1" />
+				<scanList count="1">
+					<cvParam cvRef="MS" accession="MS:1000795" name="no combination" value="" />
+					<scan instrumentConfigurationRef="IC1">
+						<referenceableParamGroupRef ref="scan1" />
+						<cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="2" />
+						<cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="1" />
+					</scan>
+				</scanList>
+				<binaryDataArrayList count="2">
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="mzArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="intensityArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="9608" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+				</binaryDataArrayList>
+			</spectrum>
+			<spectrum id="Spectrum=3" defaultArrayLength="0" index="3">
+				<referenceableParamGroupRef ref="spectrum1" />
+				<scanList count="1">
+					<cvParam cvRef="MS" accession="MS:1000795" name="no combination" value="" />
+					<scan instrumentConfigurationRef="IC1">
+						<referenceableParamGroupRef ref="scan1" />
+						<cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="3" />
+						<cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="1" />
+					</scan>
+				</scanList>
+				<binaryDataArrayList count="2">
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="mzArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="intensityArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="14404" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+				</binaryDataArrayList>
+			</spectrum>
+			<spectrum id="Spectrum=4" defaultArrayLength="0" index="4">
+				<referenceableParamGroupRef ref="spectrum1" />
+				<scanList count="1">
+					<cvParam cvRef="MS" accession="MS:1000795" name="no combination" value="" />
+					<scan instrumentConfigurationRef="IC1">
+						<referenceableParamGroupRef ref="scan1" />
+						<cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="1" />
+						<cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="2" />
+					</scan>
+				</scanList>
+				<binaryDataArrayList count="2">
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="mzArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="intensityArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="19200" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+				</binaryDataArrayList>
+			</spectrum>
+			<spectrum id="Spectrum=5" defaultArrayLength="0" index="5">
+				<referenceableParamGroupRef ref="spectrum1" />
+				<scanList count="1">
+					<cvParam cvRef="MS" accession="MS:1000795" name="no combination" value="" />
+					<scan instrumentConfigurationRef="IC1">
+						<referenceableParamGroupRef ref="scan1" />
+						<cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="2" />
+						<cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="2" />
+					</scan>
+				</scanList>
+				<binaryDataArrayList count="2">
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="mzArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="intensityArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="23996" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+				</binaryDataArrayList>
+			</spectrum>
+			<spectrum id="Spectrum=6" defaultArrayLength="0" index="6">
+				<referenceableParamGroupRef ref="spectrum1" />
+				<scanList count="1">
+					<cvParam cvRef="MS" accession="MS:1000795" name="no combination" value="" />
+					<scan instrumentConfigurationRef="IC1">
+						<referenceableParamGroupRef ref="scan1" />
+						<cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="3" />
+						<cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="2" />
+					</scan>
+				</scanList>
+				<binaryDataArrayList count="2">
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="mzArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="intensityArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="28792" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+				</binaryDataArrayList>
+			</spectrum>
+			<spectrum id="Spectrum=7" defaultArrayLength="0" index="7">
+				<referenceableParamGroupRef ref="spectrum1" />
+				<scanList count="1">
+					<cvParam cvRef="MS" accession="MS:1000795" name="no combination" value="" />
+					<scan instrumentConfigurationRef="IC1">
+						<referenceableParamGroupRef ref="scan1" />
+						<cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="1" />
+						<cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="3" />
+					</scan>
+				</scanList>
+				<binaryDataArrayList count="2">
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="mzArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="intensityArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="33588" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+				</binaryDataArrayList>
+			</spectrum>
+			<spectrum id="Spectrum=8" defaultArrayLength="0" index="8">
+				<referenceableParamGroupRef ref="spectrum1" />
+				<scanList count="1">
+					<cvParam cvRef="MS" accession="MS:1000795" name="no combination" value="" />
+					<scan instrumentConfigurationRef="IC1">
+						<referenceableParamGroupRef ref="scan1" />
+						<cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="2" />
+						<cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="3" />
+					</scan>
+				</scanList>
+				<binaryDataArrayList count="2">
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="mzArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="intensityArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="38384" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+				</binaryDataArrayList>
+			</spectrum>
+			<spectrum id="Spectrum=9" defaultArrayLength="0" index="9">
+				<referenceableParamGroupRef ref="spectrum1" />
+				<scanList count="1">
+					<cvParam cvRef="MS" accession="MS:1000795" name="no combination" value="" />
+					<scan instrumentConfigurationRef="IC1">
+						<referenceableParamGroupRef ref="scan1" />
+						<cvParam cvRef="IMS" accession="IMS:1000050" name="position x" value="3" />
+						<cvParam cvRef="IMS" accession="IMS:1000051" name="position y" value="3" />
+					</scan>
+				</scanList>
+				<binaryDataArrayList count="2">
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="mzArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="16" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+					<binaryDataArray encodedLength="0">
+						<referenceableParamGroupRef ref="intensityArray" />
+						<cvParam cvRef="IMS" accession="IMS:1000102" name="external offset" value="43180" />
+						<cvParam cvRef="IMS" accession="IMS:1000103" name="external array length" value="1199" />
+						<cvParam cvRef="IMS" accession="IMS:1000104" name="external encoded length" value="4796" />
+						<binary />
+					</binaryDataArray>
+				</binaryDataArrayList>
+			</spectrum>
+		</spectrumList>
+	</run>
+</mzML>
Binary file test-data/Heatmaps_LM8_file16.pdf has changed
Binary file test-data/Heatmaps_analyze75.pdf has changed
Binary file test-data/Heatmaps_imzml.pdf has changed
Binary file test-data/Heatmaps_rdata.pdf has changed