Jump to content
PDS Geosciences Node Community

Mineral lab reference color calculation for browse product


Nathan Weinstein

Recommended Posts

Hello,

For my thesis completion:

I’m trying to calculate the color of a lab reference mineral with the summary parameters formulations as Viviano et al. (2014) did in this figure (9):

 

image.thumb.png.ca74d3d4eca5609c8ba190fdfba5e648.png

 

For that, I used the relevant summary parameters formulations described in Viviano et al. (2014), table 2, and section 4. But first, to check that my calculations are correctly implemented (using Matlab), I tried to calculate them on an IF_TER cube to reproduce the same browse product as I would get if I opened it by the pre-calculated summary parameters bands in the SU_TER cube (as was shown originally in the data processing pipeline by Seelos et al. 2011).

Unfortunately, it wasn’t a matching result. What could be the problem?

 

I appreciate any help you can provide.

Nathan

Link to comment
Share on other sites

Hello Ray,

Thank you so much for your reply.

Attaching my Matlab `browse_product()` code (just change from `.txt` to `.m`) and a spectral library. Please run it and just follow the GUI instructions. I recommend using the default settings for the first time by hitting ‘Enter’ (and clicking on ‘OK’ in stretch parameters choose) through the process.

The core function of the summary parameters calculations is inside the main code under the function name `calc_summary_parameters()`, which calls `calc_summary_parameter()` with its nested helper functions for each parameter name.

I try to reproduce the work of Wray et al. (2010) with the FRT00009326 scene and bassanite mineral.

Also, I’ve uploaded these files to my GitHub repo.

 

Thanks a lot,

Nathan

Sulfates and Sulfite.txt browse_product.txt

Link to comment
Share on other sites

Hello Feng,

Thank you for your reply. Yes, of course.

As mentioned, I used the FRT00009326 scene.

When I load the custom browse product of the pre-calculated summary parameters of SINDEX2, D2300, BD1900_2 (stretch: [0.003, 0.020], [0.010, 0.038], [0.028, 0.058], respectively), by loading from the SU_TER cube – I get this (as should be):

SU_TER_bp.thumb.png.61f8dcb065ab642f72cd134ac3406f33.png

 

But when I am calculating these summary parameters with IF_TER cube by myself with the relevant algebra, I get different values in each parameter matrix (band), which get me this:

IF_TER_bp.thumb.png.454f3677ef0dae68c0ccda6cba43fede.png

 

Thanks,

Nathan

Link to comment
Share on other sites

Hi Nathan,

 

I am unable to run your Matlab code. I got some errors. One error is:

image.png.42ff3ebd32622ca0154c49ce88a24a75.png

I use CAT function to compute these summary parameters from IF_TER and then compare to the SU_TER. I also got the similar results like yours. I wonder if it is caused by calculation rounding error.

 

Thanks,

Feng

Link to comment
Share on other sites

Hi Nathan,

 

We emailed the question to the data producer one week ago but didn't get his response. On the other hand, we think your result and code are correct. We did get the same result as yours through different method. There is some difference between SU TER and summary parameters from IF TER. I suggest you can continue your way for the further steps.

 

Thanks,


Feng

Link to comment
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now
×
×
  • Create New...