VectorDataReprojection¶
Reproject a vector data using support image projection reference, or a user specified map projection
Description¶
Reproject vector data using a support image as projection reference or to a user given map projection. If given an image keywordlist can be added to the reprojected vectordata.
Parameters¶
Input data¶
Input vector data -in.vd filename [dtype] Mandatory
The input vector data to reproject
Use image keywords list -in.kwl image
Optional input image to fill vector data with image metadata.
Output data¶
Output vector data -out.vd filename [dtype] Mandatory
The reprojected vector data
Output Projection choice -out.proj [image|user] Default value: image
- Use image projection ref
Vector data will be reprojected in image projection ref. - User defined projection
Use image projection ref options¶
Image used to get projection map -out.proj.image.in image Mandatory
Projection map will be found using image metadata
User defined projection options¶
Map Projection -out.proj.user.map [utm|lambert2|lambert93|wgs|epsg] Default value: utm
Defines the map projection to be used.
- Universal Trans-Mercator (UTM)
A system of transverse mercator projections dividing the surface of Earth between 80S and 84N latitude. - Lambert II Etendu
This is a Lambert Conformal Conic projection mainly used in France. - Lambert93
This is a Lambert 93 projection mainly used in France. - WGS 84
This is a Geographical projection - EPSG Code
This code is a generic way of identifying map projections, and allows specifying a large amount of them. See www.spatialreference.org to find which EPSG code is associated to your projection;
Universal Trans-Mercator (UTM) options¶
Zone number -out.proj.user.map.utm.zone int Default value: 31
The zone number ranges from 1 to 60 and allows defining the transverse mercator projection (along with the hemisphere)
Northern Hemisphere -out.proj.user.map.utm.northhem bool Default value: false
The transverse mercator projections are defined by their zone number as well as the hemisphere. Activate this parameter if your image is in the northern hemisphere.
EPSG Code options¶
EPSG Code -out.proj.user.map.epsg.code int Default value: 4326
See www.spatialreference.org to find which EPSG code is associated to your projection
Elevation management¶
This group of parameters allows managing elevation values.
DEM directory -elev.dem directory
This parameter allows selecting a directory containing Digital Elevation Model files. Note that this directory should contain only DEM files. Unexpected behaviour might occurs if other images are found in this directory. Input DEM tiles should be in a raster format supported by GDAL.
Geoid File -elev.geoid filename [dtype]
Use a geoid grid to get the height above the ellipsoid in case there is no DEM available, no coverage for some points or pixels with no_data in the DEM tiles. A version of the geoid can be found on the OTB website(https://gitlab.orfeo-toolbox.org/orfeotoolbox/otb-data/blob/master/Input/DEM/egm96.grd).
Default elevation -elev.default float Default value: 0
This parameter allows setting the default height above ellipsoid when there is no DEM available, no coverage for some points or pixels with no_data in the DEM tiles, and no geoid file has been set. This is also used by some application as an average elevation value.
Examples¶
From the command-line:
otbcli_VectorDataReprojection -in.vd VectorData_QB1.shp -out.proj image -out.proj.image.in ROI_QB_MUL_1.tif -out.vd reprojected_vd.shp
From Python:
import otbApplication
app = otbApplication.Registry.CreateApplication("VectorDataReprojection")
app.SetParameterString("in.vd", "VectorData_QB1.shp")
app.SetParameterString("out.proj","image")
app.SetParameterString("out.proj.image.in", "ROI_QB_MUL_1.tif")
app.SetParameterString("out.vd", "reprojected_vd.shp")
app.ExecuteAndWriteOutput()