Cannot process through Lyzenga 2006 method about opticalrs

[jkibele]

It's hard to say for sure what the problem is, but here's my best guess. I think that all of the training pixels are getting masked. That .compressed() bit only returns unmasked pixels. Try changing lyzdepth = de.lyzenga_depth_estimation() to lyzdepth = de.lyzenga_depth_estimation(n_std=3). If that works, you could try smaller values too. ...or if it doesn't you could try larger values.

If that does work and get you results, then I'm guessing that maybe you're working in temperate waters? There's something that's been called "over deduction" in the literature. I tried to address it in the data preprocessing section and the discussion of my depth estimation paper. It's difficult to explain clearly, but it's basically just that really dark bottom types cause the Lyzenga method to try to take the log of negative numbers. Where that happens, my code masks out those pixels. If too many pixels got masked, your error could be the result. I had to use a n_std of 2 for the Lyzenga method in my paper.

Please let me know if n_std fixes your problem. ...and please keep me posted if you get some comparable results. I'd love to know how the KNN method works vs. the Lyzenga method in other parts of the world.

[CyanBC]

Thank you for the direction. I'll keep pushing on it. My preliminary results for KNN for the waters off Okinawa (Motobu region) have been very, very good. My colleagues here have used Lyzenga processed through a different method that was not very convincing.

[jkibele]

That's great to hear. Please keep me posted.

Also, if you're interested in habitat mapping, please check out my thesis. I intend to get a couple more publications out that'll put that work into a more concise format, but I've been too busy to make much progress. There's a lot more in OpticalRS than depth estimation.