Page Comparison

h2. Name

*geo* -- Geographical coordinates handling {alias:geo}

h2. Syntax - Converting between coordinate systems
{panel:|borderStyle=solid|borderWidth=1|bgColor=#ffffaa}

*geo* \[*cvt*] \[*xy*] \[*lla*]

{panel}

h2. Syntax - Displaying geographical lines
{panel:|borderStyle=solid|borderWidth=1|bgColor=#ffffaa}

*geo dra* \[*lon*] \[*lat*] \[*cut*] \[*fra*] \[*col*] \[*set*]

{panel}

h2. Syntax - Annotate geographical axes
{panel:|borderStyle=solid|borderWidth=1|bgColor=#ffffaa}

*geo ann* 
{panel}

h2. Syntax - Setting options for geographical coordinates
{panel:|borderStyle=solid|borderWidth=1|bgColor=#ffffaa}

*geo set* 
{panel}

h2. Syntax - Calculate distances between the points in a dataset
{panel:|borderStyle=solid|borderWidth=1|bgColor=#ffffaa}

*geo line point* 
{indent}
Will add a Field data column called *Point Distance* to active data and the content will be the distance to the previous point.{indent}
{panel}


h2. Description

Handling and display of geographical coordinates. 

The *geo* command has a syntax for coordinate conversion and display of geographical coordinates. The *geo* command is applied for instance in the menu script *geodesy.tcl* that comes up from *Tools > Geographical Management* . In order to learn the command and practical use one can look into the script. Look in the +Geocap installation directory+ under the +Geodesy+ module.

h2. Geodesywidget

The *geo* command uses coordinate information from the menu. This is handled through a special +widget+ called *geodesywidget* which is designed to be used in tcl scripts. One can think of *geodesywidget* as a complex tcl variable which contains the coordinate and datum information seen in the panel.

The conversion menu in Geodesy.tcl is organized in two parts: 1) *upper* and 2) *lower* . The conversion is from *upper* to *lower* or the other way. *Active data* is always part of the conversion as the input dataset.

Due to the *geodesywidget* the script code for conversion and display of geographical coordinates is relatively small.

Below is shown how the tcl code reads the *upper* and *lower* geodesywidget into two global variables called *_cu_* and *_cl_* . Those variables can be used in conversion control later on.

{noformat:|borderWidth=1|bgColor=#eeeeee}
        # Setting upper and lower geodesywidget into a tcl variables
	set ::cu [geodesywidget]
	set ::cl [geodesywidget]
{noformat}

The statements above is written into the menu part where *upper* and *lower* is coded. A change in a geodesywidget will then immediately update the global variables.

h2. Arguments
*geo cvt* \[*low*] *coord*
{indent}
Convert *active* data from *upper* to *lower*. Or the other way if *low* is used. The result is in *active*.

h3. Example showing coordinate conversion of a dataset
*Active* data will be converted from coordinates shown in upper to coordinates shown in lower geodesywidget. For this should work correct it is the users responsibility to secure that the coordinates in the *upper* geodesywidget is set correct to the coordinate system of *active*. Be aware that when a conversion command object is activated upon a dataset in the project, that dataset is make *active* before conversion starts.
 
{noformat:|borderWidth=1|bgColor=#eeeeee}
proc GcGeodesy::convertActiveDataFromUpperToLower { } {
	geodesywidget setasupper $::cu
	geodesywidget setaslower $::cl
	geo cvt coord
}
{noformat}

{indent}

*geo xy* _xcoordinate_ _ycoordinate_
{indent}
Convert a coordinate pair from *xy* to *Lon Lat* or from *Lon Lat* to *xy*.

h3. Example showing coordinate conversion of one single point (coordinate pair)

{noformat:|borderWidth=1|bgColor=#eeeeee}
proc GcGeodesy::convertSinglePoint { } {
	geodesywidget setasupper $::cu
	if {$GcGeodesy::convertOptionSingle == "From XY to LongLat"} {
		set x $GcGeodesy::easting
		set y $GcGeodesy::northing
		geo xy $x $y
		set GcGeodesy::longitude [gvar gen longitude]
		set GcGeodesy::latitude [gvar gen latitude]
	} elseif {$GcGeodesy::convertOptionSingle == "From LongLat to XY"} {
		set x $GcGeodesy::longitude
		set y $GcGeodesy::latitude
		geo lla $x $y
        # the result is saved in tcl variables
		set GcGeodesy::easting [gvar gen xcoord]
		set GcGeodesy::northing [gvar gen ycoord]
	}
}
{noformat}
{indent}

*geo xy | lla coord*
{indent}
Convert *active* data from *xy* to *Lon Lat* or from *Lon Lat* to *xy*.

h3. Example showing coordinate conversion of active data

{noformat:|borderWidth=1|bgColor=#eeeeee}
proc GcGeodesy::convertCoordinates { } {
	geodesywidget setasupper $::cu
	if {$GcGeodesy::convertOptionCoord == "From XY to LongLat"} {
		geo xy coord
	} elseif {$GcGeodesy::convertOptionCoord == "From LongLat to XY"} {
		geo lla coord
	}
}
{noformat}
{indent}

*geo dra* \[*lon*] \[*lat*] \[*inc*] \[*cut*] \[*fra*] \[*col*]
{indent}
Draw geographical coordinates
{indent}
{indent}
*lon | lat* *inc* _start stop increment_ - specify the start, stop and increment values for geographical lines.
{indent}

{indent}
*cut* - specify is the geographical lines is to be cut against the grid frame.
{indent}

{indent}
*fra* - specify is the grid frame is to be drawn.
{indent}

{indent}
*col whi | red | yel ...* - specify col of the geographical lines.
{indent}

{indent}
h3. Trimmed example from geodesy.tcl showing how to draw geographical coordinate lines

{noformat:|borderWidth=1|bgColor=#eeeeee}
proc GcGeodesy::detailedGeoDisplay { } {
	geodesywidget setasupper $::cu
	set cut " "
	if {$GcGeodesy::framecut == 1} {set cut cut}
	if {$GcGeodesy::applevelGrid == 1} {geo set geolev $GcGeodesy::levelGrid}
	if {$GcGeodesy::sticky_grid == 1} {
		geo set sticky
		set cut " "
	}
	set GcGeodesy::coordsys [geodesywidget getcoordinatesystem $::cu]
	eval geo dra $cut lon "$GcGeodesy::longitudeStart" "$GcGeodesy::longitudeStop" inc "$GcGeodesy::longitudeInc" lat "$GcGeodesy::latitudeStart" "$GcGeodesy::latitudeStop" inc "$GcGeodesy::latitudeInc" col $GcGeodesy::colorGeogrid
	geo set sticky off
	geo set geolev off
}
{noformat}
{indent}

*geo ann* _updown_ _annofreq_
{indent}
Draw geographical annotation
{indent}
{indent}
_updown_ - may contain *NE | SE | NESW* to describe where annotation is placed.
{indent}

{indent}
_annofreq_ - is a number from 1 to 6 to describe the annotation frequency.
{indent}

To see how to draw geographical annotation look into the procedure +GcGeodesy::annotateGeoGrid+ in geodesy.tcl .

*geo set* \[*oldgeo | sticky | geolev*] 
{indent}
Set options for drawing geographical lines.
{indent}
{indent}
*oldgeo* \[*off*] - use old (previous) position for geographical lines.
{indent}

{indent}
*sticky* \[*off*] - specify that the geographical lines will be glued to *sticky\_surface*.
{indent}

{indent}
*geolev* _value_ - specify that the geographical lines will be displayed at height level _value_.
{indent}


h2. Display examples of map with geographical coordinates

The first example is from the workflow *Seismic and reservoir models - Course B1*.

!att:_img_ref^geoexample.png|border=1!
_Cartographical display from course workflow_

The next examples are nice displays with cartographic standard.

{section}
{column:width=50%}
!att:_img_ref^atlantis_scale.png|border=1!
_Geographical coordinates from Atlantis demo project_
{column}
{column:width=50%}
!att:_img_gen^geocoord2.png|border=1!
_Cartographical display in Polar Stereographic projection_
{column}
{section}


h2. See also

[axe - Axes display]
{pagebreak}