Skip to content

GitLab

  • Projects
  • Groups
  • Snippets
  • Help
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
otb
otb
  • Project overview
    • Project overview
    • Details
    • Activity
    • Releases
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 273
    • Issues 273
    • List
    • Boards
    • Labels
    • Service Desk
    • Milestones
  • Merge Requests 8
    • Merge Requests 8
  • CI / CD
    • CI / CD
    • Pipelines
    • Jobs
    • Schedules
  • Operations
    • Operations
    • Incidents
    • Environments
  • Packages & Registries
    • Packages & Registries
    • Container Registry
  • Analytics
    • Analytics
    • CI / CD
    • Repository
    • Value Stream
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Members
    • Members
  • Collapse sidebar
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
  • Main Repositories
  • otbotb
  • Issues
  • #1739

Closed
Open
Opened Oct 10, 2018 by Sadewizzz@Sadewizzz

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.*

Edited Oct 30, 2018 by Manuel Grizonnet
Assignee
Assign to
None
Milestone
None
Assign milestone
Time tracking
None
Due date
None
Reference: orfeotoolbox/otb#1739