Close

Service Interruption

The NSIDC Web site and data services are currently having intermittent problems and may be unavailable. We are working to restore these services as soon as possible and apologize for any inconvenience this may cause. Please contact NSIDC User Services for assistance.

Processing Radar Data with MATLAB

Below is an example of how to view binary data with MATLAB. The steps presented here, in addition to further processing steps, can be found on Example of how to analyze RES data from Uwashington RES system using Matlab.

>> read410  (read in the DAT file)
Enter file name: 0032006.DAT
Bike Wheel Diameter in meters? (0.5)
        Using Bike Wheel Diameter of 0.500000 meters

How Many Magnets on Bike Wheel? (1)
        Using 1.000000 magnet on Bike Wheel

Trigger pos (fraction of total record) (0.1)
Read 629 records.
>> save 0032006.dat
>> whos
  Name             Size           Bytes  Class

  Averages         1x629           5032  double array
  Bike_Count       1x629           5032  double array
  Hpos             1x629           5032  double array
  Pressure         1x629           5032  double array
  TEK              1x1                8  double array
  Xinc             1x629           5032  double array
  Yinc             1x629           5032  double array
  ans              1x1                8  double array
  data          1000x629        5032000  double array
  dday             1x629           5032  double array
  time          1000x1             8000  double array

Grand total is 634405 elements using 5075240 bytes

RELEVANT ARRAYS:

Averages: number of stacks in each record
Bike_Count: number of magnet clicks
Hpos: calculated array giving horizontal displacement (m) for the profile
Pressure: digitized pressure record relative to first measurement in bits. Needs to be multiplied by -0.1078 meters/bit
Xinc: Volts/bit for converting values (not normally useful)
Yinc: Seconds/sample - used for calculating associated time array and for calculating filter parameters.
data: matrix of the records. Each column is a complete waveform. In this case there are 494 records of 1000 samples length.
dday: decimal day of year (time) each record was acquired.
time: calculated travel time array for plotting and calculating purposes.

>> plot(time,data);axis([-1e-6 19e-6 0 1024]);