search

Applied-GeoSolutions / gips

Geospatial Image Processing System
GNU General Public License v3.0
17 stars 5 forks source link

L8 Pan product #495

Open bhbraswell opened 5 years ago

bhbraswell commented 5 years ago

I'd like to export L8 band 8, 15m panchromatic as its own thing. Because it is 15m. I think I can do this if no one else is working on it. But, I know it is currently handled as part of the Asset. Any input on that is welcome. I think atmospheric correction is probably not an option and that is probably OK.

ircwaves commented 5 years ago

Looks like a panchromatic product would be a good addition, and applicable to both LC8 and LE7 (although their bandwidths are very different). We don't have it in queue yet over here. If you wanted just TOA, it could be as simple as:

                    out_fn = os.path.join(self.path, self.basename + '_pan.tif')
                    vsi_pan_path = [os.path.join('/vsitar/' + asset_obj.filename, f)
                        for f in md['filenames'] if f.endswith('B8.tif')][0]
                    os.symlink(out_fn, vsi_pan_path)
                    imgout = gippy.GeoImage(out_fn)

If you want to atmospherically correct, then you would need to extract the image, and compute an ESUN value for LC8's panchromatic band. And then we could get back into discussion about the appropriate ESUN values to use, keep it :100:%! And maybe even talk about how badly we need to get VRT products in the mix, and deprecate symlinking :open_mouth: