[gdal-dev] gdal & terrasar-x

[Vincent Schut]

I'm having problems warping terrasar data, both strip and spot. 
Everything seems to run smooth, but all resulting pixels are 0. Running 
gdalinfo -mm on the original terrasar xml header file shows a decent 
min/max, and opening the xml with openev shows the characteristic radar 
speckle image. However, somehow gdalwarp seems to convert all values to 0.
I'm running gdal from svn.

This is the gdalinfo -mm on the original:

Driver: TSX/TerraSAR-X Product
Files: TSX1_SAR__SSC______SL_S_SRA_20080203T094750_20080203T094751.xml
Size is 6432, 8123
Coordinate System is `'
GCP Projection =
GCP[  0]: Id=, Info=
           (3186,4061) -> 
(4.58577884873432,-58.7378260671197,88.2895278671792)
GCP[  1]: Id=, Info=
           (1,1) -> (4.62120875202606,-58.6758682205696,88.2895278671792)
GCP[  2]: Id=, Info=
           (6354,1) -> (4.6416261082315,-58.7806954917766,88.2895278671792)
GCP[  3]: Id=, Info=
           (1,8123) -> (4.52969375238776,-58.6937801528285,88.2895278671792)
GCP[  4]: Id=, Info=
           (6353,8123) -> 
(4.5500912421315,-58.7984603452316,88.2895278671792)
Metadata:
   OPERATIONAL_MODE=operational
   ORBIT_CYCLE=22
   ABSOLUTE_ORBIT=3542
   ORBIT_DIRECTION=DESCENDING
   IMAGING_MODE=SL
   PRODUCT_VARIANT=SSC
   IMAGE_TYPE=COMPLEX
   ROW_SPACING=6.06688654869458333E-09
   COL_SPACING=1.78571426658891141E-04
   AZIMUTH_LOOKS=1.00000000000000000E+00
   RANGE_LOOKS=1.00000000000000000E+00
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 8123.0)
Upper Right ( 6432.0,    0.0)
Lower Right ( 6432.0, 8123.0)
Center      ( 3216.0, 4061.5)
Band 1 Block=6432x1 Type=CInt16, ColorInterp=Undefined
     Computed Min/Max=-1329.000,1444.000
   Metadata:
     POLARIMETRIC_INTERP=VV

I then run gdalwarp:

$ gdalwarp 
TSX1_SAR__SSC______SL_S_SRA_20080203T094750_20080203T094751.xml 
TSX1_SAR__SSC______SL_S_SRA_20080203T094750_20080203T094751.tif
Creating output file that is 8144P x 8994L.
Processing input file 
TSX1_SAR__SSC______SL_S_SRA_20080203T094750_20080203T094751.xml.
0...10...20...30...40...50...60...70...80...90...100 - done.

$  gdalinfo -mm 
TSX1_SAR__SSC______SL_S_SRA_20080203T094750_20080203T094751.tif
Driver: GTiff/GeoTIFF
Files: TSX1_SAR__SSC______SL_S_SRA_20080203T094750_20080203T094751.tif
Size is 8144, 8994
Coordinate System is `'
Origin = (4.529703226723088,-58.675978519735985)
Pixel Size = (0.000013777380802,-0.000013777380802)
Image Structure Metadata:
   INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (   4.5297032, -58.6759785)
Lower Left  (   4.5297032, -58.7998923)
Upper Right (   4.6419062, -58.6759785)
Lower Right (   4.6419062, -58.7998923)
Center      (   4.5858047, -58.7379354)
Band 1 Block=8144x1 Type=CInt16, ColorInterp=Gray
     Computed Min/Max=0.000,0.000

See, I get georeferencing but lost all the data...

Cheers,
Vincent.