<div dir="ltr"><span style="font-size:12.8px">Hi Jonathan,</span><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Jonathan Greenberg-4 wrote</span><span class="gmail-im" style="font-size:12.8px"><br>> USGS, in their infinite wisdom, decided to embed what is effectively a<br>> 1-bit, 16-band byte-interleaved-by-pixel image into a 16-bit, single-band<br>> geotiff for the Landsat quality mask.  I could write some crazy<br>> translation<br>> program that "reverse engineers" the data, but I was wondering if there is<br>> some way to "hack" or translate the file by simply redefining the "header"<br>> (change the number type to 1-bit, number of layers to 16, interleave to<br>> BIP).  Any suggestions for how to do this?</span><br><div class="gmail_extra"><br></div><div class="gmail_extra"><a href="https://github.com/mapbox/rio-l8qa">https://github.com/mapbox/rio-l8qa</a> might be what you're looking for; provides a python API and CLI for extracting arrays out of the Landsat QA band, both the pre-Collection and collection format. It's at version 0.1 so caveat emptor but could be useful. Let me know how it works for you.<br></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature">Matt Perry</div></div>
<br></div></div></div>