ConvertCartoToGeoPoint

Convert cartographic coordinates to geographic ones.

Description

This application computes the geographic coordinates from cartographic ones. User has to give the X and Y coordinate and the cartographic projection (see mapproj parameter for details).

This application has several output images and supports “multi-writing”. Instead of computing and writing each image independently, the streamed image blocks are written in a synchronous way for each output. The output images will be computed strip by strip, using the available RAM to compute the strip size, and a user defined streaming mode can be specified using the streaming extended filenames (type, mode and value). Note that multi-writing can be disabled using the multi-write extended filename option: &multiwrite=false, in this case the output images will be written one by one. Note that multi-writing is not supported for MPI writers.

Parameters

Input cartographic coordinates

X cartographic coordinates -carto.x float Mandatory
X cartographic coordinates in the projection defined by mapproj parameter

Y cartographic coordinates -carto.y float Mandatory
Y cartographic coordinates in the projection defined by mapproj parameter

Map Projection -mapproj [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 -mapproj.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 -mapproj.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 -mapproj.epsg.code int Default value: 4326
See www.spatialreference.org to find which EPSG code is associated to your projection

Output long -long float Mandatory
Point longitude coordinates in decimal degrees.

Output lat -lat float Mandatory
Point latitude coordinates in decimal degrees.

Examples

From the command-line:

otbcli_ConvertCartoToGeoPoint -carto.x 367074.625 -carto.y 4835740 -mapproj utm -mapproj.utm.northhem true -mapproj.utm.zone 31

From Python:

import otbApplication

app = otbApplication.Registry.CreateApplication("ConvertCartoToGeoPoint")

app.SetParameterFloat("carto.x", 367074.625)
app.SetParameterFloat("carto.y", 4835740)
app.SetParameterString("mapproj","utm")
app.SetParameterString("mapproj.utm.northhem","true")
app.SetParameterInt("mapproj.utm.zone", 31)

app.ExecuteAndWriteOutput()