Changes between Initial Version and Version 4 of Ticket #6293
- Timestamp:
- Jan 6, 2016, 12:04:50 PM (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Ticket #6293 – Description
initial v4 1 1 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: 2 2 3 {{{DIM_PHR1A_MS_201509141606268_SEN_1627093101-002.XML 3 {{{ 4 DIM_PHR1A_MS_201509141606268_SEN_1627093101-002.XML 4 5 ICON_PHR1A_MS_201509141606268_SEN_1627093101-002.JPG 5 6 IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TFW … … 13 14 PREVIEW_PHR1A_MS_201509141606268_SEN_1627093101-002.JPG 14 15 PREVIEW_PHR1A_MS_201509141606268_SEN_1627093101-002.KMZ 15 RPC_PHR1A_MS_201509141606268_SEN_1627093101-002.XML}}} 16 RPC_PHR1A_MS_201509141606268_SEN_1627093101-002.XML 17 }}} 16 18 17 19 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). … … 19 21 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: 20 22 21 {{{gdalwarp -rpc IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TIF IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1_projected.TIF}}} 23 {{{ 24 gdalwarp -rpc IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TIF IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1_projected.TIF 25 }}} 22 26 23 27 When I then gdalwarp the second image: 24 28 25 {{{gdalwarp -rpc IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TIF IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1_projected.TIF}}} 29 {{{ 30 gdalwarp -rpc IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TIF IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1_projected.TIF 31 }}} 26 32 27 33 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). … … 32 38 gdalinfo on two original images: 33 39 34 {{{$ gdalinfo IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TIF 40 {{{ 41 $ gdalinfo IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TIF 35 42 Driver: GTiff/GeoTIFF 36 43 Files: IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TIF … … 96 103 STATISTICS_MINIMUM=0 97 104 STATISTICS_STDDEV=266.56884586952 98 Band 4 Block=9458x1 Type=UInt16, ColorInterp=Undefined}}} 105 Band 4 Block=9458x1 Type=UInt16, ColorInterp=Undefined 106 }}} 99 107 100 {{{$ gdalinfo IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TIF 108 {{{ 109 $ gdalinfo IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TIF 101 110 Driver: GTiff/GeoTIFF 102 111 Files: IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TIF … … 162 171 STATISTICS_MINIMUM=0 163 172 STATISTICS_STDDEV=187.78652664007 164 Band 4 Block=9458x1 Type=UInt16, ColorInterp=Undefined}}} 173 Band 4 Block=9458x1 Type=UInt16, ColorInterp=Undefined 174 }}} 165 175 166 176 TFW files: 167 177 168 {{{$ cat IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TFW 178 {{{ 179 $ cat IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R1C1.TFW 169 180 4.00000000000000 170 181 0 … … 172 183 -4.00000000000000 173 184 1.00000000000000 174 1.00000000000000}}} 185 1.00000000000000 186 }}} 175 187 176 {{{$ cat IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TFW 188 {{{ 189 $ cat IMG_PHR1A_MS_201509141606268_SEN_1627093101-002_R2C1.TFW 177 190 4.00000000000000 178 191 0 … … 180 193 -4.00000000000000 181 194 1.00000000000000 182 -59999.00000000000000}}} 195 -59999.00000000000000 196 }}}