c c comment ___ GEOCODING ___ c GEO_OUT_LAM Outdata/lam.raw GEO_OUT_PHI Outdata/phi.raw
If you want to obtain the latitude/longitude of the pixels in an interferogram that was created, but you do not have a DEM in radarcoordinates available, you will have to create one. This means that you will have to edit the products result file and create a SLANT2H section, see Chapter 32 for a description of this section in the products result file. The section will look something like:
******************************************************************* *_Start_slant2h: ******************************************************************* Method: schwabisch Data_output_file: Outdata/dummy_height.raw Data_output_format: real4 First_line (w.r.t. original_master): 1001 Last_line (w.r.t. original_master): 2105 First_pixel (w.r.t. original_master): 501 Last_pixel (w.r.t. original_master): 700 Multilookfactor_azimuth_direction: 10 Multilookfactor_range_direction: 2 Ellipsoid (name,a,b): WGS84 6.37814e+06 6.35675e+06 ******************************************************************* * End_slant2h:_NORMAL *******************************************************************
(Also set the pcf to 1 on top of the products result file.) If you have an external DEM you can compute the required file (Outdata/dummy_height.raw) using step COMPREFDEM. The dimensions and multilooking can be copied from the interferogram section.
If your area is flat, you may want to use a dummy file filled with zeros. You can create such a file with appropriate dimensions using Matlab for example. Alternatively, the much faster Unix way would be along these lines. First compute the height (number of lines) of the dummy file:
echo "(2105-1001+1)/10" | bc -l 110.5
then the width (number of pixels):
echo "(700-501+1)/2" | bc -l 100
(ie., the file should be 110 lines by 100 pixels of 4 byte). Now create the file using dd:
dd if=/dev/zero of=Outdata/dummy_height.raw count=110 bs=400