Opened 4 years ago

Last modified 3 years ago

#6293 closed defect

Multiple Pleiades/Airbus Tiles Don't Georegister Properly With gdalwarp — at Version 4

Reported by: edjansing Owned by: warmerdam
Priority: high Milestone: 2.0.3
Component: default Version: unspecified
Severity: major Keywords: gdalwarp, tiles, tfw, geotiff, georegistration
Cc: bishop

Description (last modified by edjansing)

Attempting to use gdalwarp to georegister several MSI tiles from Airbus' Pleiades sensor. A single area is provided in several tiles. In this example, here is the MSI directory listing with two tiles:

DIM_PHR1A_MS_201509141606268_SEN_1627093101-002.XML
ICON_PHR1A_MS_201509141606268_SEN_1627093101-002.JPG
IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TFW
IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TIF
IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TFW
IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TIF
INDEX.HTM
LIBRARY
LINEAGE
MASKS
PREVIEW_PHR1A_MS_201509141606268_SEN_1627093101-002.JPG
PREVIEW_PHR1A_MS_201509141606268_SEN_1627093101-002.KMZ
RPC_PHR1A_MS_201509141606268_SEN_1627093101-002.XML

The TFW files are unique between the two tiles. When I import the tiles into qgis, without any modification, it stitches them together properly, but they are not in the proper CRS (they are essentially referenced to pixel space, not to any CRS coordinate system).

When I use gdalwarp to put them in EPSG:4326, the first image is perfectly registered (the image with R1C1 as a suffix). The command used:

gdalwarp -rpc IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TIF IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1_projected.TIF

When I then gdalwarp the second image:

gdalwarp -rpc IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TIF IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1_projected.TIF

The second image does not get stitched at the bottom of the first image, but is placed at the top of the image, as if the tile starts at the same lat/lon coordinate (or at the same pixel location, whichever).

Is this normal behavior? Is there another procedure I should use? Can provide image data if helpful.

gdalwarp reports 2.0.0 as version gdalinfo on two original images:

$ gdalinfo IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TIF
Driver: GTiff/GeoTIFF
Files: IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TIF
       IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TIF.aux.xml
       ./DIM_PHR1A_MS_201509141606268_SEN_1627093101-002.XML
       ./RPC_PHR1A_MS_201509141606268_SEN_1627093101-002.XML
       IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TFW
Size is 9458, 15000
Coordinate System is `'
Origin = (-1.000000000000000,3.000000000000000)
Pixel Size = (4.000000000000000,-4.000000000000000)
Metadata:
  METADATATYPE=DIMAP
  TIFFTAG_DATETIME=2015:12:23 18:32:39
  TIFFTAG_IMAGEDESCRIPTION=B2 B1 B0 B3
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
  INTERLEAVE=PIXEL
RPC Metadata:
  HEIGHT_OFF=45
  HEIGHT_SCALE=95
  LAT_OFF=38.87154971192919
  LAT_SCALE=0.1416732123880315
  LINE_DEN_COEFF= 1 0.00026341144632701 -9.80822136398733e-05 -2.35953356209866e-06 -1.62170380697022e-05 -1.22640028446246e-08 1.55023313626127e-07 -1.5537111237091e-05 4.58036130936434e-05 -2.31198284873618e-05 3.2864533726195e-09 -3.06927464009735e-08 2.08793797077866e-06 -1.21475710839433e-08 1.54788420545811e-07 -4.44599257156653e-07 1.70368489248657e-09 1.33031932377729e-10 -1.54424096736091e-08 1.05222239987211e-10
  LINE_NUM_COEFF= -0.000119096471843013 0.00273796193211152 -1.00367656228961 -0.000859086493278233 0.000261713135774831 2.33979663689761e-07 -2.25402196362394e-06 -0.000340119441043266 2.19730836249524e-05 1.66642358408085e-08 -4.29909520634416e-08 4.95298240422465e-08 -2.62364432140933e-05 -6.28870839700389e-08 1.52651889380179e-05 8.08563520366722e-05 2.3195367391933e-05 1.1804689412366e-08 1.00010692110933e-07 1.98020806117544e-08
  LINE_OFF=7772.5
  LINE_SCALE=7771.5
  LONG_OFF=-77.13546672147527
  LONG_SCALE=0.1165198796426168
  SAMP_DEN_COEFF= 1 0.000472894792539169 0.000881419369988542 -7.53324742935084e-05 -3.46629852664066e-05 -4.87301506925132e-08 1.32614018054834e-08 -2.1872832721835e-07 1.49479836237352e-05 -6.59860144983547e-06 -1.48728353597154e-08 -3.52095449470832e-08 -2.14688694596795e-07 -3.3168465870712e-09 4.03572496745229e-08 3.97707771445365e-07 -1.15020506050241e-08 2.66109369887284e-09 1.05965377023676e-08 9.65442134520645e-10
  SAMP_NUM_COEFF= -0.000434608495681351 1.00210264963283 -0.000281300697070519 0.000316194516045025 -0.000879304997048266 7.52970914458278e-05 -5.20633271426872e-05 -3.73769088263009e-05 0.00138775409674169 -2.73849487440847e-08 5.19133691835542e-08 -1.24502711479507e-05 -2.1363347686686e-05 -6.59563969341986e-06 5.20138704377129e-05 -3.41383248102939e-05 -2.02958513833057e-09 -3.71155366535213e-08 -2.25289299184377e-07 -2.06504369348464e-09
  SAMP_OFF=5000.5
  SAMP_SCALE=4999.5
Corner Coordinates:
Upper Left  (  -1.0000000,   3.0000000)
Lower Left  (      -1.000,  -59997.000)
Upper Right (   37831.000,       3.000)
Lower Right (   37831.000,  -59997.000)
Center      (   18915.000,  -29997.000)
Band 1 Block=9458x1 Type=UInt16, ColorInterp=Red
  Min=0.000 Max=3924.000
  Minimum=0.000, Maximum=3924.000, Mean=241.479, StdDev=247.675
  Metadata:
    STATISTICS_MAXIMUM=3924
    STATISTICS_MEAN=241.47873181735
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=247.67541356986
Band 2 Block=9458x1 Type=UInt16, ColorInterp=Green
  Min=0.000 Max=3956.000
  Minimum=0.000, Maximum=3956.000, Mean=300.351, StdDev=253.536
  Metadata:
    STATISTICS_MAXIMUM=3956
    STATISTICS_MEAN=300.35078661846
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=253.53558755775
Band 3 Block=9458x1 Type=UInt16, ColorInterp=Blue
  Min=0.000 Max=4095.000
  Minimum=0.000, Maximum=4095.000, Mean=348.032, StdDev=266.569
  Metadata:
    STATISTICS_MAXIMUM=4095
    STATISTICS_MEAN=348.03234066915
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=266.56884586952
Band 4 Block=9458x1 Type=UInt16, ColorInterp=Undefined
$ gdalinfo IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TIF
Driver: GTiff/GeoTIFF
Files: IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TIF
       IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TIF.aux.xml
       ./DIM_PHR1A_MS_201509141606268_SEN_1627093101-002.XML
       ./RPC_PHR1A_MS_201509141606268_SEN_1627093101-002.XML
       IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TFW
Size is 9458, 544
Coordinate System is `'
Origin = (-1.000000000000000,-59997.000000000000000)
Pixel Size = (4.000000000000000,-4.000000000000000)
Metadata:
  METADATATYPE=DIMAP
  TIFFTAG_DATETIME=2015:12:23 18:32:39
  TIFFTAG_IMAGEDESCRIPTION=B2 B1 B0 B3
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
  INTERLEAVE=PIXEL
RPC Metadata:
  HEIGHT_OFF=45
  HEIGHT_SCALE=95
  LAT_OFF=38.87154971192919
  LAT_SCALE=0.1416732123880315
  LINE_DEN_COEFF= 1 0.00026341144632701 -9.80822136398733e-05 -2.35953356209866e-06 -1.62170380697022e-05 -1.22640028446246e-08 1.55023313626127e-07 -1.5537111237091e-05 4.58036130936434e-05 -2.31198284873618e-05 3.2864533726195e-09 -3.06927464009735e-08 2.08793797077866e-06 -1.21475710839433e-08 1.54788420545811e-07 -4.44599257156653e-07 1.70368489248657e-09 1.33031932377729e-10 -1.54424096736091e-08 1.05222239987211e-10
  LINE_NUM_COEFF= -0.000119096471843013 0.00273796193211152 -1.00367656228961 -0.000859086493278233 0.000261713135774831 2.33979663689761e-07 -2.25402196362394e-06 -0.000340119441043266 2.19730836249524e-05 1.66642358408085e-08 -4.29909520634416e-08 4.95298240422465e-08 -2.62364432140933e-05 -6.28870839700389e-08 1.52651889380179e-05 8.08563520366722e-05 2.3195367391933e-05 1.1804689412366e-08 1.00010692110933e-07 1.98020806117544e-08
  LINE_OFF=7772.5
  LINE_SCALE=7771.5
  LONG_OFF=-77.13546672147527
  LONG_SCALE=0.1165198796426168
  SAMP_DEN_COEFF= 1 0.000472894792539169 0.000881419369988542 -7.53324742935084e-05 -3.46629852664066e-05 -4.87301506925132e-08 1.32614018054834e-08 -2.1872832721835e-07 1.49479836237352e-05 -6.59860144983547e-06 -1.48728353597154e-08 -3.52095449470832e-08 -2.14688694596795e-07 -3.3168465870712e-09 4.03572496745229e-08 3.97707771445365e-07 -1.15020506050241e-08 2.66109369887284e-09 1.05965377023676e-08 9.65442134520645e-10
  SAMP_NUM_COEFF= -0.000434608495681351 1.00210264963283 -0.000281300697070519 0.000316194516045025 -0.000879304997048266 7.52970914458278e-05 -5.20633271426872e-05 -3.73769088263009e-05 0.00138775409674169 -2.73849487440847e-08 5.19133691835542e-08 -1.24502711479507e-05 -2.1363347686686e-05 -6.59563969341986e-06 5.20138704377129e-05 -3.41383248102939e-05 -2.02958513833057e-09 -3.71155366535213e-08 -2.25289299184377e-07 -2.06504369348464e-09
  SAMP_OFF=5000.5
  SAMP_SCALE=4999.5
Corner Coordinates:
Upper Left  (      -1.000,  -59997.000)
Lower Left  (      -1.000,  -62173.000)
Upper Right (   37831.000,  -59997.000)
Lower Right (   37831.000,  -62173.000)
Center      (   18915.000,  -61085.000)
Band 1 Block=9458x1 Type=UInt16, ColorInterp=Red
  Min=0.000 Max=3131.000
  Minimum=0.000, Maximum=3131.000, Mean=303.059, StdDev=222.224
  Metadata:
    STATISTICS_MAXIMUM=3131
    STATISTICS_MEAN=303.05929724396
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=222.22368384095
Band 2 Block=9458x1 Type=UInt16, ColorInterp=Green
  Min=0.000 Max=3071.000
  Minimum=0.000, Maximum=3071.000, Mean=389.381, StdDev=199.449
  Metadata:
    STATISTICS_MAXIMUM=3071
    STATISTICS_MEAN=389.38147159371
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=199.44912482837
Band 3 Block=9458x1 Type=UInt16, ColorInterp=Blue
  Min=0.000 Max=3020.000
  Minimum=0.000, Maximum=3020.000, Mean=453.803, StdDev=187.787
  Metadata:
    STATISTICS_MAXIMUM=3020
    STATISTICS_MEAN=453.8029269754
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=187.78652664007
Band 4 Block=9458x1 Type=UInt16, ColorInterp=Undefined

TFW files:

$ cat IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TFW
4.00000000000000
0
0
-4.00000000000000
1.00000000000000
1.00000000000000
$ cat IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TFW
4.00000000000000
0
0
-4.00000000000000
1.00000000000000
-59999.00000000000000

Change History (6)

Changed 4 years ago by edjansing

Properly stitched, not in CRS

Changed 4 years ago by edjansing

Improperly stitched, but in EPSG:4326 CRS

comment:1 Changed 4 years ago by Even Rouault

Two things :

  • You should pass the 2 images in the same gdalwarp command "gdalwarp -rpc in1.tif in2.tif out.tif". If you do in several times, it will not increase the raster dimension & georeferenced extent of the first output.
  • You could try removing the .TFW files as I'm wondering if they don't confuse gdalwarp

comment:2 Changed 4 years ago by edjansing

Tried the following:

  • Passing both images into the same gdalwarp command. No improvement.
  • Passing both images into the same gdalwarp command, without the TFW files. No improvement.
  • Processing each image separately, without the TFW files. No improvement.

comment:3 Changed 4 years ago by edjansing

Used:

gdalwarp -rpc -dstnodata 0 -r cubic first.tif second.tif output.tif

comment:4 Changed 4 years ago by edjansing

Description: modified (diff)
Note: See TracTickets for help on using tickets.