On 6/14/2012 12:27 PM Ben Temperton said...
Hi there,
I am working with mass spectroscopy data in the mzXML format that looks like
this:
<mzXML>
<msRun>
<scan num="1">...</scan>
<scan num="2">...</scan>
<scan num="3">...</scan>
<scan num="4">...</scan>
.....
</msRun>
<index>
<offset id="1">160409990</offset>
<offset id="2">160442725</offset>
<offset id="3">160474927</offset>
<offset id="4">160497386</offset>
....
</index>
</mzXML>
Where the offset element contains the byte offset
byte offset to what in what?
of the scan element that shares the id. I am trying to write a python script to
remove scan elements and their respective offset, but I can't figure out how I
re-calculate the byte offset for each remaining element once the elements have
been removed.
Removing the reference will render the location as inaccessible so you
may not need to recalculate the remaining offsets. I'm assuming these
are offset in some other file and that the lengths of the data contained
at that location is known or knowable. But, if you don't point to that
location you won't try to get the data at the offset.
If you're trying to minimize the size on disk you'll probably want to
write a new file and leave the original alone. Initialize a buffer,
then for each valid offset id add the content and bump the offset by the
content length (or whatever is appropriate)
Emile
My plan was to write the file out, the read it back in again and search through the file for a
particular string (e.g. '<scan num="1">') and then use the tell() method to
return the current byte location in the file. However, I'm not sure how I would implement this.
Any ideas?
Many thanks,
Ben
--
http://mail.python.org/mailman/listinfo/python-list