Orthorectification Issue with OTB 6.6.0
Description
Hello,
I am trying to orthorectify a Pleiades images with OTB library and i have some difficult to do this on OTB 6.6.0. My output is smallest than expected (31 000 ko rather than 3 Go). However on OTB 6.4.0 the process works fine.
My full process bellow.
Steps to reproduce
SET DEM=DEM
SET OUT=OUTFOLDER
SET GEOID=GEOID
SET EPSG=EPSG
SET PAN1=PANCHROMATIQUE_TIFF_1
SET PAN2=PANCHROMATIQUE_TIFF_2
SET MS=MULTISPECTRAL_TIFF
SET IMGREF=IMAGE_REFERENCE (that is a .vrt)
@echo Fusion pan + ms
otbcli_tilefusion -il %PAN1%"?&skipcarto=true" %PAN2%"?&skipcarto=true" -cols 1 -rows 2 -out %OUT%\PANCHRO.tif"?&gdal:co:COMPRESS=LZW&gdal:co:TILED=YES&gdal:co:PREDICTOR=2&gdal:of:GTiff" uint16
@echo etape2: Pansharp 50cm
otbcli_BundleToPerfectSensor -inp %OUT%\PANCHRO.tif"?&skipcarto=true" -inxs %MS%"?&skipcarto=true" -out %OUT%\PXS.tif"?&gdal:co:COMPRESS=LZW&gdal:co:TILED=YES&gdal:co:PREDICTOR=2&gdal:of:GTiff" uint16 -elev.dem %DEM% -mode phr -ram 15000
@echo Create Envelope
otbcli_ImageEnvelope -in %OUT%\PXS.tif"?&skipcarto=true" -proj "EPSG:%EPSG%" -out %OUT%\ENVELOP_PXS.shp
@echo SIFT
otbcli_HomologousPointsExtraction -in1 %OUT%\PXS.tif -in2 %IMGREF% -algorithm sift -mode geobins -mode.geobins.binstep 1000 -mode.geobins.binsize 1000 -mfilter 1 -precision 20 -2wgs84 1 -threshold 0.6 -backmatching 1 -out %OUT%\HOMOLOGOUS_POINTS_SIFT.txt -outvector %OUT%\HOMOLOGOUS_POINTS_SIFT.shp -elev.dem %DEM% -elev.geoid %GEOID%
@echo Refine sensor model
otbcli_RefineSensorModel -ingeom %OUT%\PXS.geom -outgeom %OUT%\PXS_REFINE.geom -inpoints %OUT%\HOMOLOGOUS_POINTS_SIFT.txt -map epsg -map.epsg.code %EPSG% -elev.dem %DEM% -elev.geoid %GEOID% -outstat %OUT%\PXS_GEOM_REFINE_STATS.txt -outvector %OUT%\PXS_GEOM_REFINE_VECTOR.shp
@echo Orthorectification avec gcp
otbcli_OrthoRectification -io.in %OUT%\PXS.tif"?&skipcarto=true&geom=%OUT%\PXS_REFINE.geom" -io.out %OUT%\PXS_ORTHO.tif"?&gdal:co:COMPRESS=LZW&gdal:co:TILED=YES&gdal:co:PREDICTOR=2&gdal:co:BIGTIFF=YES&gdal:of:GTiff" uint16 -elev.dem %DEM% -elev.geoid %GEOID% -map epsg -map.epsg.code %EPSG% -opt.gridspacing 20 -outputs.spacingx 0.5 -outputs.spacingy -0.5 -opt.ram 15000
Configuration information
Windows 10 x64 OTB 6.6.0 OTB 6.4.0 Python 2.7.*