Ausiàs R Posted March 2, 2020 Posted March 2, 2020 (edited) Good afternoon. I'm an undergraduate student of geomatics and remote sensing, and I was trying to use the CRISM images to classify minerals and chemical compounds on the surface of Mars. I collected some information about all the products to use, also I followed the workflow proposed here and I installed CAT on ENVI. I was trying to use a an image of the Jezero Crater data product: frt00005c5e_07_if166j_mtr3, which is projected, atmosferically and geometrically corrected. The data is on I/F which as is explained here, on the FAQs means: Q: What does I/F stand for? A: CRISM is a reflectance spectrometer, and I/F is how reflectance is represented algebraically: I is the energy (actually radiance) observed by the CRISM instrument, and F is the energy (actually solar irradiance) incident at the top of the Martian atmosphere. I/F is a ratio of energies (radiance/irradiance), with some additional scaling so the ratio is unitless. Acording to my knowledge radiance/irradiance=reflectance (unitless), maybe I'm wrong. The case is that I tried to classify the minerals on the image of the Jezero Crater using this library from Viviano et al using the "CRISM RATIOED I/F CORRECTED" spectra but I get wrong or absurd results. I noticed that the spectra of my images is clearly different from the spectra that I find on the Viviano library ratioed i/f. My Spectra always follow this pattern, and never reaches a reflectance of 0.5. Nevertheless on the Viviano library I found spectra like this one of epidote: Why the spectra of my image is so different? Is my data really on reflectance? How can I get real reflectances? Edited March 3, 2020 by Ausiàs R error
Ray Arvidson Posted March 3, 2020 Posted March 3, 2020 Hi There: I am the manager of the PDS Geosciences Node and also a CRISM Science Team Member. The I/F data you are using come from the spectral radiance for each band and pixel location divided by the solar spectral irradiance divided by pi. The is a ratio of spectral radiances and the pi converts solar spectral irradiance to spectral radiance. The solar value used is specific for the heliocentric distance when the scene was acquired, i.e., Mars distance relative to he sun. Your I/F spectra, put through the volcano scan correction using CAT, look good. Here is what you are missing. The next step typically is to then divide the volcano scan corrected I/F spectra by a spectrum from the same corrected data for an area that looks like it is spectrally neutral. This removes any instrument artifacts and tends to de-emphasize the ubiquitous basaltic sands and regolith and thus emphasizes what might be unique in terms of mineral signatures for the area of interest. It kind is like in the lab when you use reflectance spectrometers you ratio your data relative to some standard surface. Remember there are no pure outcrops of given minerals on Mars, rather wind mixes things a lot. I hope this note helps. Please ask more questions if needed. Ray Arvidson
Ausiàs R Posted March 4, 2020 Author Posted March 4, 2020 (edited) Hi Ray, First of all I want to thank for your reply, it is really apreciated. You said that: On 3/3/2020 at 4:37 PM, Ray Arvidson said: The next step typically is to then divide the volcano scan corrected I/F spectra by a spectrum from the same corrected data for an area that looks like it is spectrally neutral. I understand that the area that looks like it and is spectrally neutral means an area that exhibes a constant reflectivity or transmissivity over the electromagnetic spectrum. But, where can I find this kind of surface/area? Is the data that I need on the ddr files of the image? I saw that there are 14 bands on the ddr images. On this document (pages 76-77) it is explained what those bands represents. I'm not pretty sure which band I shoud use, or if any. Once again, thanks for your answer, Ausiàs. Edited March 4, 2020 by Ausiàs R Grammar
Ray Arvidson Posted March 9, 2020 Posted March 9, 2020 You need to find an area in the corrected spectral data cube that looks spectrally neutral and divide your spectrum of interest by the spectrum that looks spectrally neutral. No DDR involved. Ray
Ausiàs R Posted April 1, 2020 Author Posted April 1, 2020 I'm sorry, I don't know what that area can be or how can I get it. Nevertheless thanks for your help. Ausiàs.
PRIYANKA Posted November 24, 2020 Posted November 24, 2020 On 3/3/2020 at 9:07 PM, Ray Arvidson said: Hi There: I am the manager of the PDS Geosciences Node and also a CRISM Science Team Member. The I/F data you are using come from the spectral radiance for each band and pixel location divided by the solar spectral irradiance divided by pi. The is a ratio of spectral radiances and the pi converts solar spectral irradiance to spectral radiance. The solar value used is specific for the heliocentric distance when the scene was acquired, i.e., Mars distance relative to he sun. Your I/F spectra, put through the volcano scan correction using CAT, look good. Here is what you are missing. The next step typically is to then divide the volcano scan corrected I/F spectra by a spectrum from the same corrected data for an area that looks like it is spectrally neutral. This removes any instrument artifacts and tends to de-emphasize the ubiquitous basaltic sands and regolith and thus emphasizes what might be unique in terms of mineral signatures for the area of interest. It kind is like in the lab when you use reflectance spectrometers you ratio your data relative to some standard surface. Remember there are no pure outcrops of given minerals on Mars, rather wind mixes things a lot. I hope this note helps. Please ask more questions if needed. Ray Arvidson hello sir, here you have mentioned sir to divide the volcano scan corrected I/F spectra by a spectrum from the same corrected data for an area that looks like it is spectrally neutral area.? what is process of choosing spectrally neutral area in a corrected i/f image
Ausiàs R Posted November 24, 2020 Author Posted November 24, 2020 Hi, Take a look to these pdfs (more precisely to the last one): https://pds-geosciences.wustl.edu/workshops/CRISM_Workshop_Mar09.htm Ausiàs.
PRIYANKA Posted December 7, 2020 Posted December 7, 2020 On 11/24/2020 at 3:16 PM, Ausiàs R said: Hi, Take a look to these pdfs (more precisely to the last one): https://pds-geosciences.wustl.edu/workshops/CRISM_Workshop_Mar09.htm Ausiàs. thanks a lot
PRIYANKA Posted January 3, 2021 Posted January 3, 2021 as it is mentioned here that choosing denominator will be similar size spectrally neural dark area nearby numerator. what is the reason of choosing area nearby numerator ? why we are not taking any dark material throughout the images? and how we are getting stretch limit for each of summary parameter?
Ray Arvidson Posted January 4, 2021 Posted January 4, 2021 See attached paper, Section 5.3, for what we do for thresholds for parameter maps. Choosing a neutral spectrum is a case of trial and error. I do not do this, rather I just look at the volcano scanned spectra for features such as the 2.3 um metal-OH absorption. Ray Arvidson Viviano-Beck_et_al-2014-Journal_of_Geophysical_Research__Planets.pdf
Ray Arvidson Posted June 28, 2021 Posted June 28, 2021 Hi There: I am the Geosciences Node Manager and also a CRISM Science Team Member. I can answer your questions. I/F is defined as the radiance on sensor divided by the solar radiance at the top of the atmosphere that would have occurred during the observation. The TRDRs are presented with these units. It is a form of reflectance but one that has the atmosphere and surface terms included. The Viviano-Beck spectral parameters operate on these data sets. They do not operate on he MTRDRs, which have been processed to remove atmospheric effects and normalized to surface reflectance, with S and L bands joined. This document explains the processing and is a must for reading if you wish to use CRISM data: https://pds-geosciences.wustl.edu/mro/mro-m-crism-2-edr-v1/mrocr_0001/document/crism_dpsis.pdf Ray Arvidson
TARUN SINGH PANU Posted June 29, 2021 Posted June 29, 2021 On 6/28/2021 at 8:51 PM, Ray Arvidson said: Hi There: I am the Geosciences Node Manager and also a CRISM Science Team Member. I can answer your questions. I/F is defined as the radiance on sensor divided by the solar radiance at the top of the atmosphere that would have occurred during the observation. The TRDRs are presented with these units. It is a form of reflectance but one that has the atmosphere and surface terms included. The Viviano-Beck spectral parameters operate on these data sets. They do not operate on he MTRDRs, which have been processed to remove atmospheric effects and normalized to surface reflectance, with S and L bands joined. This document explains the processing and is a must for reading if you wish to use CRISM data: https://pds-geosciences.wustl.edu/mro/mro-m-crism-2-edr-v1/mrocr_0001/document/crism_dpsis.pdf Ray Arvidson
TARUN SINGH PANU Posted August 30, 2021 Posted August 30, 2021 Respected Sir/Madam I am Mtech student in RS GIS and working on the CRISM data for my major project. I want MTRDR data for this project for the year of 2020 but the data updated in the PDS Geosciences Node is till the May 21, 2012. So can you please update me about when the directory will get updated with the latest MTRDR data and can you please guide me regarding the processing for getting it out from the raw EDR data. I will be thankful to you.
Jennifer Ward Posted September 3, 2021 Posted September 3, 2021 Hello. I spoke with members of the CRISM team, and unfortunately there are no plans to produce MTRDR products for CRISM observations acquired during 2020. There will only be 2 or 3 more deliveries from the team, and these deliveries will not include any MTRDR data past May 21, 2012. These deliveries will just fill in the gaps prior to that date. Hope this helps, Jennifer Ward PDS Geosciences Node
Recommended Posts
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 accountSign in
Already have an account? Sign in here.
Sign In Now