Coordinates of upper left corner in sensor geometry should match ... But
you have :
In the first case (single tile) :
|Origin:[0.5,0.5]
||Spacing:[1,1]
|
In the second one (mosaic) :
|Origin:[1,-1]
Spacing:[1,-1]
|
Can you try removing the J2W files or use the skipcarto extended
filename option BEFORE generating the mosaic ? If this does not work, I
would like the gdalinfo output for both.
Thanks,
Julien
Le 23/08/2016 à 12:13, André Stumpf a écrit :
Sure...
Upper left tile image in:
|
otbcli_ReadImageInfo
-incrise/PROD_SPOT6_001/VOL_SPOT6_001_A/IMG_SPOT6_PMS_001_A/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.JP2
2016Aug2311:41:10:Application.logger (INFO)
Imagegeneral information:
Numberof bands :4
Nodata flags :Notfound
Startindex :[0,0]
Size:[12288,14336]
Origin:[0.5,0.5]
Spacing:[1,1]
Estimatedground spacing (inmeters):[1.77451,1.76453]
Imageacquisition information:
Sensor:SPOT 6
Imageidentification number:1285949101
Acquisitiontime :2015-05-02T04:47:00
Country:Nepal
Town:Lidiṅ
Imagefootprint coordinates:
Upperleft corner (latitude,longitude)=[28.3104,84.6953]
Upperright corner (latitude,longitude)=[28.2831,85.0394]
Lowerleft corner (latitude,longitude)=[28.0326,84.6948]
Lowerright corner (latitude,longitude)=[28.0037,85.0405]
ImagedefaultRGB composition:
[R,G,B]=[0,1,2]
Groundcontrol points information:
Numberof GCPs=0
GCPsprojection =
Outputparameters value:
indexx:0
indexy:0
sizex:12288
sizey:14336
spacingx:1
spacingy:1
originx:0.5
originy:0.5
estimatedgroundspacingx:1.774505496
estimatedgroundspacingy:1.76453495
numberbands:4
sensor:SPOT 6
id:1285949101
time:2015-05-02T04:47:00
ullat:28.31040001
ullon:84.69526672
urlat:28.28313637
urlon:85.03935242
lrlat:28.0037384
lrlon:85.04052734
lllat:28.03257179
lllon:84.69483185
town:Lidiṅ
country:Nepal
rgb.r:0
rgb.g:1
rgb.b:2
projectionref:
keyword:
gcp.count:0
gcp.proj:
gcp.ids:
gcp.info:
gcp.imcoord:
gcp.geocoord:
|
upper left tile orthorectification log:
|
otbcli_OrthoRectification
-io.in"crise/PROD_SPOT6_001/VOL_SPOT6_001_A/IMG_SPOT6_PMS_001_A/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.JP2"-io.out"crise/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1_ortho.tif?&gdal:co:TILED=YES"uint16
-map utm -map.utm.zone 45-outputs.spacingx 2-outputs.spacingy
-2-outputs.mode autosize -elev.dem SRTM30 -opt.ram 1000
2016Aug2311:32:04:Application.logger (INFO)Elevationmanagement:setting
defaultheight above ellipsoid to 0meters
2016Aug2311:32:04:Application.logger
(INFO)Elevationmanagement:usingDEM directory (SRTM30)
2016Aug2311:32:04:Application.logger (INFO)Generatingoutput withsize
=[11061,13789]
2016Aug2311:32:04:Application.logger (INFO)Generatingoutput withpixel
spacing =[2,-2]
2016Aug2311:32:04:Application.logger (INFO)Generatingoutput withorigin
=[273642,3.13359e+06]
2016Aug2311:32:04:Application.logger (INFO)Areaoutside input image
bounds will have a pixel value of [0,0,0,0]
2016Aug2311:32:04:Application.logger (INFO)Usinga deformation grid
witha physical spacing of 4
2016Aug2311:32:04:Application.logger (INFO)Usinga deformation grid of
size [5530,6894]
ossimImageElevationDatabase::createCell WARN:
Couldnotopen:/home/stumpf/Data/Nepal/SRTM30/egm96.grd
Removingfile frommap!
Writingcrise/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1_ortho.tif?&gdal:co:TILED=YES...:100%[**************************************************](2e+02seconds)
|
mosaiced image info:
|
otbcli_ReadImageInfo
-incrise/PROD_SPOT6_001/VOL_SPOT6_001_A/IMG_SPOT6_PMS_001_A/IMG_SPOT6_PMS_201505020447120_SEN_1285949101.tif
2016Aug2311:44:35:Application.logger (INFO)
Imagegeneral information:
Numberof bands :4
Nodata flags :Notfound
Startindex :[0,0]
Size:[19081,17533]
Origin:[1,-1]
Spacing:[1,-1]
Estimatedground spacing (inmeters):[1.76262,1.76281]
Imageacquisition information:
Sensor:NULL
Imageidentification number:NULL
Country:Nepal
Town:Lidiṅ
Imagefootprint coordinates:
Upperleft corner (latitude,longitude)=[28.3104,84.6953]
Upperright corner (latitude,longitude)=[28.2831,85.0394]
Lowerleft corner (latitude,longitude)=[28.0326,84.6948]
Lowerright corner (latitude,longitude)=[28.0037,85.0405]
ImagedefaultRGB composition:
[R,G,B]=[0,1,2]
Groundcontrol points information:
Numberof GCPs=0
GCPsprojection =
Outputparameters value:
indexx:0
indexy:0
sizex:19081
sizey:17533
spacingx:1
spacingy:-1
originx:1
originy:-1
estimatedgroundspacingx:1.762622833
estimatedgroundspacingy:1.762808681
numberbands:4
sensor:NULL
id:NULL
time:
ullat:28.31039238
ullon:84.69527435
urlat:28.28312683
urlon:85.03936005
lrlat:28.00373077
lrlon:85.0405426
lllat:28.03256226
lllon:84.69484711
town:Lidiṅ
country:Nepal
rgb.r:0
rgb.g:1
rgb.b:2
projectionref:
keyword:
gcp.count:0
gcp.proj:
gcp.ids:
gcp.info:
gcp.imcoord:
gcp.geocoord:
|
and the log ofthe mosiac orthorectification:
|
otbcli_OrthoRectification
-io.in"crise/PROD_SPOT6_001/VOL_SPOT6_001_A/IMG_SPOT6_PMS_001_A/IMG_SPOT6_PMS_201505020447120_SEN_1285949101.tif?&skipcarto=true"-io.out"crise/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_ortho.tif?&gdal:co:TILED=YES"uint16
-map utm -map.utm.zone 45-outputs.spacingx 2-outputs.spacingy
-2-outputs.mode autosize -elev.dem SRTM30 -opt.ram 1000
2016Aug2311:49:19:Application.logger (INFO)Elevationmanagement:setting
defaultheight above ellipsoid to 0meters
2016Aug2311:49:19:Application.logger
(INFO)Elevationmanagement:usingDEM directory (SRTM30)
2016Aug2311:49:19:Application.logger (INFO)Generatingoutput withsize
=[17174,17294]
2016Aug2311:49:19:Application.logger (INFO)Generatingoutput withpixel
spacing =[2,-2]
2016Aug2311:49:19:Application.logger (INFO)Generatingoutput withorigin
=[273534,3.13358e+06]
2016Aug2311:49:19:Application.logger (INFO)Areaoutside input image
bounds will have a pixel value of [0,0,0,0]
2016Aug2311:49:19:Application.logger (INFO)Usinga deformation grid
witha physical spacing of 4
2016Aug2311:49:19:Application.logger (INFO)Usinga deformation grid of
size [8587,8647]
ossimImageElevationDatabase::createCell WARN:
Couldnotopen:/home/stumpf/Data/Nepal/SRTM30/egm96.grd
Removingfile frommap!
Writingcrise/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_ortho.tif?&gdal:co:TILED=YES...:100%[**************************************************](4e+02seconds)
|
This is OTB 5.4 on Ubuntu 16.04. Thanks for having a look and let me
know if further infos are needed.
Cheers,
André
On Tuesday, August 23, 2016 at 9:44:32 AM UTC+2, Julien Michel wrote:
Hi, can you post the logs of the upper left tile for case 1 and
also the
complete logs for case 2 ?
A ReadImageInfo on both would also be helpful.
Thanks,
Julien
Le 19/08/2016 à 15:18, André Stumpf a écrit :
> ...just realized that the GIF doesn't seem to work after the
upload,
> so here two PNGs attached.
>
> On Friday, August 19, 2016 at 3:13:17 PM UTC+2, André Stumpf wrote:
>
> Hello,
>
> I noticed a small offset between orthoimges generated from
> individual tiles and mosaics. Not sure if this is a bug or a
> feature...
>
> I have a SPOT 6 image delivered in four separate tiles:
> |
>
> IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.JP2
> IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R2C1.JP2
> IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R2C2.JP2
> IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C2.JP2
> |
>
> ...and its respective DIMAPV2 file:
>
> |
> DIM_SPOT6_PMS_201505020447120_SEN_1285949101.XML
> RPC_SPOT6_PMS_201505020447120_SEN_1285949101.XML
> |
>
> One strategy would be an orthorectification of each tile
indivually...
>
> |
> otbcli_OrthoRectification
> -io.in
<http://io.in>"PROD_SPOT6_001/VOL_SPOT6_001_A/IMG_SPOT6_PMS_001_A/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.JP2"-io.out"PROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES"uint16
> -map utm -map.utm.zone 45-outputs.spacingx 2-outputs.spacingy
> -2-outputs.mode autosize -elev.dem ../SRTM30 -opt.ram 1000
> 2016Aug1912:44:01:Application.logger
> (INFO)Elevationmanagement:setting defaultheight above
ellipsoid to
> 0meters
> 2016Aug1912:44:02:Application.logger
> (INFO)Elevationmanagement:usingDEM directory (../SRTM30)
> 2016Aug1912:44:02:Application.logger (INFO)Generatingoutput
> withsize =[11061,13789]
> 2016Aug1912:44:02:Application.logger (INFO)Generatingoutput
> withpixel spacing =[2,-2]
> 2016Aug1912:44:02:Application.logger (INFO)Generatingoutput
> withorigin =[273642,3.13359e+06]
> 2016Aug1912:44:02:Application.logger (INFO)Areaoutside input
image
> bounds will have a pixel value of [0,0,0,0]
> 2016Aug1912:44:02:Application.logger (INFO)Usinga
deformation grid
> witha physical spacing of 4
> 2016Aug1912:44:02:Application.logger (INFO)Usinga
deformation grid
> of size [5530,6894]
>
WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:0%[WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:2%[*WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:4%[**WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:6%[***WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:8%[****WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:10%[*****WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:12%[******WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:14%[*******WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:16%[********WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:18%[*********WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:20%[**********WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:22%[***********WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:24%[************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:26%[*************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:28%[**************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:30%[***************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:32%[****************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:34%[*****************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:36%[******************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:38%[*******************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:40%[********************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:42%[*********************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:44%[**********************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:46%[***********************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:48%[************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:50%[*************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:52%[**************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:54%[***************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:56%[****************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:58%[*****************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:60%[******************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:62%[*******************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:64%[********************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:66%[*********************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:68%[**********************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:70%[***********************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:72%[************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:74%[*************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:76%[**************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:78%[***************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:80%[****************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:82%[****************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:84%[****************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:86%[****************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:88%[****************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:90%[****************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:92%[****************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:94%[****************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:96%[****************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:98%[****************************************WritingPROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_R1C1.tif?&gdal:co:TILED=YES...:100%[**************************************************](2e+02seconds)
>
>
> |
>
>
> ...the other (better?) strategy would be to fuse with
> otbcli_TileFusion:
> |IMG_SPOT6_PMS_201505020447120_SEN_1285949101.tif|
> |
> |...and then orthorectify the mosaic:
>
> |
> otbcli_OrthoRectification
> -io.in
<http://io.in>"IMG_SPOT6_PMS_201505020447120_SEN_1285949101.tif"-io.out"PROD_SPOT6_001/IMG_SPOT6_PMS_201505020447120_SEN_1285949101_ortho.tif?&gdal:co:TILED=YES"uint16
> -map utm -map.utm.zone 45-outputs.spacingx 2-outputs.spacingy
> -2-outputs.mode autosize -elev.dem ../SRTM30 -opt.ram 1000
> |
>
>
> In their overlap area both orthos have an sytematic offset of
> about one pixel (see attached GIF). Unfortunately I lack a
> reference to evaluate which of the orthos is more correct
but in
> any case I guess this should not happen since the input RPCs
are
> the same in both cases?
>
> Cheers,
> André
>
>
>
>
>
>
> --
> --
> Check the OTB FAQ at
> http://www.orfeo-toolbox.org/FAQ.html
<http://www.orfeo-toolbox.org/FAQ.html>
>
> You received this message because you are subscribed to the Google
> Groups "otb-users" group.
> To post to this group, send email to [email protected]
<javascript:>
> To unsubscribe from this group, send email to
> [email protected] <javascript:>
> For more options, visit this group at
> http://groups.google.com/group/otb-users?hl=en
<http://groups.google.com/group/otb-users?hl=en>
> ---
> You received this message because you are subscribed to the Google
> Groups "otb-users" group.
> To unsubscribe from this group and stop receiving emails from
it, send
> an email to [email protected] <javascript:>
> <mailto:[email protected] <javascript:>>.
> For more options, visit https://groups.google.com/d/optout
<https://groups.google.com/d/optout>.
--
Julien MICHEL
CNES - DCT/SI/AP - BPI 1219
18, avenue Edouard Belin
31401 Toulouse Cedex 09 - France
Tel: +33 561 282 894 - Fax: +33 561 283 109
--
--
Check the OTB FAQ at
http://www.orfeo-toolbox.org/FAQ.html
You received this message because you are subscribed to the Google
Groups "otb-users" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/otb-users?hl=en
---
You received this message because you are subscribed to the Google
Groups "otb-users" group.
To unsubscribe from this group and stop receiving emails from it, send
an email to [email protected]
<mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout.
--
Julien MICHEL
CNES - DCT/SI/AP - BPI 1219
18, avenue Edouard Belin
31401 Toulouse Cedex 09 - France
Tel: +33 561 282 894 - Fax: +33 561 283 109
--
--
Check the OTB FAQ at
http://www.orfeo-toolbox.org/FAQ.html
You received this message because you are subscribed to the Google
Groups "otb-users" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/otb-users?hl=en
---
You received this message because you are subscribed to the Google Groups "otb-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
