Jump to content
PDS Geosciences Node Community
  • Recent Topics

    • March 13, 2025 - New MEX HRSC RDR Version 4 data loaded into ODE -    Added HRSC new radiometrically calibrated (RDR version 4) image data into ODE, with coverage from August 6, 2024 to December 29, 2024 (Orbit 26500). ODE product search page: https://ode.rsl.wustl.edu/mars/productsearch ODE map search page: https://ode.rsl.wustl.edu/mars/mapsearch ODE data holdings page: https://ode.rsl.wustl.edu/odeholdings/Mars_holdings.html
    • Hi Tejay,   I've confirmed your assessment that the BD1900_2 values reported are actually calculations for BD1930, instead. I've also tracked down the issue with BD1900R2 - in the SU file you referenced, the slopes for RC#### were calculated using R1815 and R2132, instead of R1850 and R2060.   Your calculations for BD1900_2 and BD1900R2 are correct, so I would recommend you keep using those. I'll reach out to the CRISM team regarding these issues you identified in the parameter calculations so we can get them corrected.   John
    • Hi Tejay,   Thanks for bringing this to our attention! There's definitely something odd going on with those two parameters in the SU file you referenced, and I agree with your initial assessment that the BD1900_2 values reported look like they're reporting BD1930, instead. I'll look into this in more detail, and report back here once I figure out what's going on.   John
    • Upon looking in this further, all BD1900_2 seems to be returning is BD1930:   The error between BD1900_2 as per the CRISM DPSIS documentation is this: But when I just return BD1930 instead:   Is anyone able to verify my suspicion. The next question is how is BD1900r2 really calculated as that one is going to be a lot hard to tease out?
    • Hi All, What I am seeking to do is to replicate/verify the calculation of the various bands contained within su164j MTRDR files and apply the same calculations to spectral data of minerals from a lab so I can compare like to like thereafter. I can replicate all of the MTRDR bands I have selected from the su164j file (with convincing accuracy) via the spectral data contained in the if164j file – except for BD1900_2 and BD1900r2 and thus the reason for my post. I am wondering if anyone else has tried to replicate the same and found an error with the specified formula? I have found an error in the CRISM DPSIS when trying to replicate other bands via cross referencing different sources… but this one I can’t resolve. I've detailed all my calculations: The formula From the CRISM DPSIS document:   The data From the file frt0000b385_07_su164j_mtr3.img and cell reference [580,595] I get BD1900_2 = -0.009848237. I am now going to try and replicate this result using the spectral data contained within frt0000b385_07_if164j_mtr3.img.   Kernel Width As I understand it, I need to retrieve the five closest wavelengths to 1850, 1930, 1985 and 2067, extract the reflectance values and calculate the median, this is as follows:   Calculating "a" and "b" values: The formula for this is: However, I could interpret this in 3 ways:| 1) Ideal - Stick strictly to the formula in the CRISM DPSIS document and use the wavelengths of 1850, 1930, 1985 and 2067; OR 2) Closest to Ideal - Use the wavelength closest to those ideal wavelengths, so for example instead of 1850 I would use 1848.95; OR 3) Median - Use the wavelength corresponding to where the median reflectance value come from in the case of R1850:5 that would be 1842.34. Calculating all 3 possibilities gives:   Calculating BD1900_2 Following the formula for BD1900_2 from there, with the 3 possibilities for the calculation of "a" and "b" gives:   Results discussion As per the file frt0000b385_07_su164j_mtr3.img at cell reference [580,595], the value of BD1900_2 is -0.009848237. Trying to calculate the same, but from first principals via the spectral data contained in the file frt0000b385_07_if164j_mtr3.img does not give the same result - not even the same sign... in fact for any cell there is a concerning difference between the BD1900_2 values in the su164j file and that calculated from the if164j. With other bands, I can pretty much reproduce nearly identical result (to 6-7 decimal places typically), which leads me to ask: - Is there a typo or something in the stated formula for BD1900_2 that is preventing me from replicating the same? Attached are my excel workings if someone can assist. Note: I am a Financial Math Masters student whose designated research project is identifying minerals on Mars via machine learning so if my lingo is not up to scratch, it's because just 5 weeks ago I didn't have any clue about spectral data! Thank you!! BD1900_2.xlsx
×
×
  • Create New...