I just wrote the decoding code based on the Level-0 products of Sentinel-1A, according to the library documents downloaded from the website: earth.esa.int, named “A6_S1-IF-ASD-PL-0007_SAR Space Packet Protocol Data Unit.pdf”. I have 3 questions to ask for help:
I tested several IW mode data(TOPS), and I found that statistical property of echo I decoded is abnormal (shown in Figure below), only focusing on several values(0, 13, 128, 245). I am wondering what’s the matter, whether the decoding of BAQ or statistical approach?
I did the range compression using the echo I decoded, based on the first burst. The parameters I used are: sampling rate is 6.434523812571428e7 Hz, bandwidth is 56.5MHz, pulse width is 5.240481033595628e-5, but the compression result is incorrect. I don’t know where the problem is, is it caused by compression method or echo decoding method ?
Is there a processing tool based on the Sentinel-A raw data ? Covert raw data to complex echo data.
The results are pretty strange indeed. I can confirm that several other people have used successfully this datasert to develop their own decoder. In other words, it is not a compression issue but a decoding issue, I am afraid.
I provided in addition to RAW data a version decoded such that you can verify pixel by pixel the result of your decoding.
The FDBAQ decompression is pretty complex. I suggest you carefully review the implementation of the section 4 of the S-1 SAR packet protocol data unit (S1-IF-ASD-PL-0007).
The simple fact of using correctly the tables provided in section 5.2 makes such that the Svalue is a float number and not an integer. This is another indication of something wrong in the decoding.
Thanks for your reply. @nuno.miranda. Maybe I didn’t express my statistical approach cleraly. The interger value shown in the picture was the result of quantization among 0~255. The decoding result is float type.
You said that you provided a version decoded. How can I get it ? Could you give me a sample result? I want to learn about it and verify my result. My email is elivesyy@gmail.com.
All the information you need is on this document (available in the library)
Sentinel-1 Level-0 Data Decoding Package (COPE-GSEG-EOPG-TN-15-0005, 28/01/2015 ). Section 5 and 6 desribes the test data (RAW and decoded)
The document I was referrring to was in the document library. Since the Jan 2015.
Wiht your last remark, I understood that you have some difficulties in finding. After checking, I realise that actually it has been removed. I can’t tell when or why.
I asked to put it back in the library and it is there (as free and open as one could wish.) and here it is:
Hi all
I’m trying to decode L0 by myself.
I get some abnormal values of ‘PulseStartFrequency’ parameter on some IW products. Like ~4.8GHz for IW1&IW2 and ~7GHz for IW3
Centre frequency is said to be 5.405 GHz, but what is operational range? Could’n find it anythere
This is my IDL code.
It is the exact formula as what is written in the document section 3.2.5.7
;+
; :Description:
; Describe the of the Tx Pulse Starting Frequency in [MHz]
;
; :Params:
; txpsfCode : code of the TXPSF. The code shall be signed already
; txprr : value of TX Pulse Ramp Rate in [MHz / us]
; :Author: nunomiranda
;-
function decode_txpsf, txpsfCode, txprr, HERTZ=HERTZ
fref =37.5347222D ;[MHz]
txpsf = (txprr / 4D / fref) + txpsfCode * fref / 2D^14
if keyword_set(HERTZ) then txpsf*=1e6
return, txpsf ;
end
Also the S-1 level-1 product SLC or GRD report the values decoded by the processor. For example:
I don’t believe that actually you need to add the radar frequency to it.
The txPulseLength, txPulseStartFrequency, txPulseRampRate are all what you need to construct the chirp and range compress the data
N
I downloaded the sample data decoded you referred. I found that there were several TIFF format files. I tried to read them in the Matlab, using the function ‘imread(‘xxx.tiff’)’, but an error happened: ‘Error using rtifc. Unsupported sample format 5’. How to read this type of TIFF file and get the pixel value ?
Yes, I think so. Most of current software are focusing on the processing based on the image product, not the raw echo data. It is important for many researchers who major imaging formation algorithm to use the decoded raw data. I look forward to use such decoding-software very much.
At ESA, RAW data have always been considered as users product. We always made them available starting from ERS and we continue for S-1. THis is not the case for other missions. However, we never made available RAW decoding software.
We have supported and will continue supporting anyone who would like to develop a “decoder” and we fully support if this is made available in public domain.
My personal opinion is that people willing to process S-1 RAW data (or any other mission) needs to learn about the sensor specificities.
If you want to process S-1 TOPS data, you have to learn about S-1. Being able to decode the RAW data is part of the learning curve. I can witness it>
The issue you are encountering is related to the matlab TIFF support that doesn’t understand complex format pixels. The standard tiff library and also GDAL seem able to properly understand this data (including IDL+ENVI).
Dear Nuno Miranda,
sorry to use this post relating to another product but maybe you can help me with my request.
For my PhD project I need to work with Level 0/RAW ERS data, collecting all possible info recorded during the acquisition step and used in the focusing stage. I have the permission from ESA to download them through Eolisa. Unfortunately the files have extension .E1 or .E2 and and not .dat (as requested in commercial software as ENVI that reads ERS-1 Level-0 with extension dat*.001).
Do you know some available software to open this type of data? Moreover, do you know if some document for decoding is available?
The native format of ERS was the CEOS format that was a folder based product containing several file like the leader file, the dat file (dat.*001).
When came ASAR into play ESA have started proposing to the users to get ERS in native or in ENvisat format. The two formats were maitained in parallel for many, many years. In the last years, ESA have put in place a huge activity aiming at consolidating the ERS archive: rapatriating the old magnetic tape from every stations in the world, dumping their content (when possible) and keeping, perform an analysis to segregate good from bad data in order to store the data with todays standard ensuring that users can still have access to this data in the future. The idea is to allow users to work on long time series ERS1 + ERS2 +ASAR +S-1.
As part of this activity, it could very well be that the decision was taken to no support the native CEOS format. I recommend to drop an email to eohelp@esa.int for a clear confirmation.
If this format is not supported by ENvi/SARSCAPe you need to go back to the Envi people to update the s.w accordingly. Another solution is that you get directly the SLC or PRI from ESA.
However, I recommend to ask confirmation to eohelp in first place.
Kind regards,
Nuno Miranda
Dear Nuno,
thank you for your reply and I’m sorry to respond to you so late. Unfortunately ENVI is not able to deal with radar raw data and SARscape is an expensive tool and not useful for my project. I need to decode the RAW data to collect all the info associated to it. I read your response about S1 RAW data and I would like to learn how it is possible to decode it. Do you know some open source code to do that?
Dear Nuno,
the document Sentinel 1 SAR Instrument Calibration and Characterisation Plan (S1-PL-ASD-PL-0001) is not available for download.
Will I need it too to decode level 0 data?
Regards,
Igor
Hi Nuno,
the FTP account in this document doesn’t work anymore. I want to debug the huffman decoder in my fdbaq decompressor. Is there any chance an example file could be posted again?