Madagascar2011 - 06 - OTB geometry processing
1 like928 views
The document discusses the geometric pre-processing capabilities of the Orfeo Toolbox (OTB), including orthorectification of images to transform them from satellite to geographic coordinates using geometric models, superimposing images for tasks like pan-sharpening, and manual registration of images through selection of homologous points. It provides examples of orthorectifying and superimposing images using both the graphical user interface and command line tools in OTB.
Madagascar2011 - 06 - OTB geometry processing
- 1. Orfeo Toolbox Geometric pre-processing Stéphane MAY stephane.may@cnes.fr orfeo-toolbox.org 1
- 2. Geometric corrections Goal Orthorectify images ➢ Transform images from satellite referential to terrestrial referential Cartographic or geographic projection ➢ Use of the geometrical model Impact of the geometrical model precision Superpose images ➢ Superpose images for pan-sharpening, change detection OTB orthorectification works for : ➢ Spot 5, Ikonos, Quickbird, GeoEye, WorldView orfeo-toolbox.org 2
- 3. Geometric corrections Menu Geometry > Reproject Image (1/2) orfeo-toolbox.org 3
- 4. Geometric corrections Menu Geometry > Reproject Image (2/2) orfeo-toolbox.org 4
- 5. Use case 1 : Orthorectify an image without DEM Menu File > Open ./img_kalideos_reunion/2003-05-04/j_1a/1.0/IMAGERY.TIF Menu Geometry > Reproject image UTM / Linear User defined Parameters : ➢ Longitude : 55.612 Latitude : -21.290 Spacing X : 10, Spacing Y : -10 (m) Size X : 520, Size Y : 500 File > Save dataset orfeo-toolbox.org 5
- 6. Use case 1 : Orthorectify an image without DEM orfeo-toolbox.org 6
- 7. Use case 2 : Orthorectify an image with DEM Menu File > Open ./img_kalideos_reunion/2003-05-04/j_1a/1.0/IMAGERY.TIF Menu Geometry > Reproject image UTM / Linear User defined Parameters : ➢ Longitude : 55.612 Latitude : -21.290 Spacing X : 10, Spacing Y : -10 (m) Size X : 520, Size Y : 500 DEM : enter the directory where .hgt files are File > Save dataset Visualization > Viewer > Add 2 images and compare orfeo-toolbox.org 7
- 8. Use case 2b : Orthorectify an image Use the command line tools gdalinfo image.tif otbObtainUTMZoneFromGeoPoint -lat -21.290 -lon 55.612 otbOrthoRectification-cli -in ./img_kalideos_reunion/2003-05- 04/j_1a/1.0/IMAGERY.TIF -out output.tif -ul 978771 7638749 -rpc 10 -size 500 520 -spacing 10 -10 -mapProj UTM 39 S -interp LINEAR otbConvertGeoToCartoPoint -lat -21.290 -lon 55.612 -mapProj UTM 39 S otbConvertCartoToGeoPoint -x 978771 -y 7638749 -mapProj UTM 39 S orfeo-toolbox.org 8
- 9. Monteverdi – Superimpose two images Menu Geometry > Superimpose two images Combination of geometrical models of the two images (note : works only with images with geometrical model) ➢ Image to be superimposed : use of the direct geometrical model to build Longitude/Latitude projection ➢ Combination with the inverse model of the reference image. ➢ Reprojection orfeo-toolbox.org 9
- 10. Monteverdi – Superimpose two images Menu Geometry > Superimpose two images orfeo-toolbox.org 10
- 11. Use case 3 : Superimpose an image orfeo-toolbox.org 11
- 12. Use case 3b : Superimpose an image Use the command line tool otbSuperimpose-cli -inR img_geom1.tif -inM ../data_formation/img_kalideos_reunion/2004-08- 19/j_1a/1.0/IMAGERY.TIF -out output2.tif orfeo-toolbox.org 12
- 13. Monteverdi – Manual registration Menu Geometry > Homologous points extraction orfeo-toolbox.org 13
- 14. Use case 4 : Manual registration Menu File > Open Last 2 ortho images Menu Geometry > Homologous points extraction Select homologous points File > Save dataset Visualization > Viewer > Add 2 images (2003 and 2004 registered) and compare orfeo-toolbox.org 14
- 15. Orfeo Toolbox - Registration Use command line tools otbRigidTransformResample-cli -in output2.tif -out output2_tr.tif -tx 150 -ty -90 otbFineRegistration-cli -ref output.tif -sec output2.tif -w output2.tif -out disparity.tif -wo output2_reg.tif -er 10 10 -mr 9 9 -co 150 -90 -m CC -spa 0.5 -ssr 10 10 orfeo-toolbox.org 15
- 16. Monteverdi – Build a simple sensor model Menu Geometry > GCP to sensor model orfeo-toolbox.org 16
- 17. Use case 5 : P+XS fusion Menu File > Open ./img_kalideos_reunion/2003-05-04/thr_1a/1.0/IMAGERY.TIF ./img_kalideos_reunion/2003-05-04/j_1a/1.0/IMAGERY.TIF Menu Geometry Orthorectification of PAN image over a selected ROI (Region of Interest) Menu Geometry > Superimpose 2 images Image to reproject : Reader1 (XS) Reference image : Projection>OutputImage Option : DEM : enter the directory where .hgt files are Filtering > Pan-sharpening > Simple RSC pan sharpening File > Save dataset Filtering > Pan-sharpening > Bayesian fusion File > Save dataset orfeo-toolbox.org Visualization > Viewer > Add 2 images (2 fusions) and compare 17
- 18. Monteverdi Thank you for your attention ! orfeo-toolbox.org 18