gsdview.coor

Coordinate conversion

Classes and functions for coordinate conversion among the following reference systems:

• rc: image reference system (rc_s for products in slat range

  projection and rc_g for ptoducts in ground) (row, col) [pixel, pixel]

• ttau: times coordinates both in acquisition or zero Doppler geometry

  (slow time, fast time ) (azimuth, range) [s, s]

• xyz: cartesian coordinates in the terrestrian reference system (ECEF)

  (x, y, z) [m, m, m]

• llh: geodetic coordinates

  (lat, lon, h) [deg, deg, m] (lat, lon, h) [rad, rad, m]

There can be different method to perform conversions according to the input information availability:

• sttv (state vector mode): based on the information about the satellite position, velocity, acceleration, the Doppler Centroid ferequency, the wavelegth.
• grid (grid mode): based on the conversion grid given as input parameter; the conversion grid states the corresponces between (row, col) image points and (lat, lon, h) geodetic coordinate points.

src	dst	OK	function/method	mode	note
rc_s	ttau	OK	SlantRC2TTau.rc2ttau	direct
ttau	rc_s	OK	SlantRC2TTau.rc2ttau	direct
rc_g	ttau	NO	GroundRC2TTau.ttau2rc	rc_g -> rc_s -> ttau
ttau	rc_g	NO	GroundRC2TTau.ttau2rc	ttau -> rc_s -> rc_g
rc_s	xyz	OK	SlantRC2XYZ.rc2xyz	rc_s -> ttau -> xyz
xyz	rc_s	OK	SlantRC2XYZ.xyz2rc	xyz -> ttau -> rc_s
rc_g	xyz	NO	GroundRC2XYZ.rc2xyz	rc_g -> rc_s -> ttau -> xyz
xyz	rc_g	NO	GroundRC2XYZ.xyz2rc	xyz -> ttau -> rc_s -> rc_g
rc_s	llh	OK	X.rc2xyz + xyz2llh	rc_s -> ttau -> xyz -> llh
llh	rc_s	OK	llh2xyz + X.xyz2rc	llh -> xyz -> ttau -> rc_s
rc_g	llh	NO	GroundRC2LLH.rc2llh	rc_g -> rc_s -> ttau -> xyz -> llh
llh	rc_g	NO	GroundRC2LLH.llh2rc	llh -> xyz -> ttau -> rc_s -> rc_g
ttau	xyz	OK	TTau2XYZ.ttau2xyz	direct
xyz	ttau	OK	TTau2XYZ.xyz2ttau	direct
ttau	llh	OK	X.ttau2xyz + xyz2llh	ttau -> xyz -> llh
llh	ttau	OK	llh2xyz + X.xyz2ttau	llh -> xyz -> ttau
xyz	llh	OK	xyz2llh	direct
llh	xyz	OK	llh2xyz	direct
rc	llh	OK	GridRC2LLH.rc2llh	direct	grid based
llh	rc	OK	GridLLH2RC.llh2rc	direct	grid based

Functions

all(a[, axis, out])	Test whether all array elements along a given axis evaluate to True.
arange	arange([start,] stop[, step,], dtype=None)
asarray(a[, dtype, order])	Convert the input to an array.
atleast_1d(*arys)	Convert inputs to arrays with at least one dimension.
doppler(R, V, lambda_[, fd])	Doppler equation
ellipsoid(P[, ell])	Ellipsoid equation
isscalar(num)	Returns True if the type of num is a scalar type.
llh2xyz_deg(lat, lon, h[, ell])	Convert geocentric coordinates to cartesian
llh2xyz_rad(lat, lon, h[, ell])	Convert geocentric coordinates to cartesian
only_projected_projs()
osr_datums()	Return most significant datums
projections()	Return most significant projections
ravel(a[, order])	Return a flattened array.
size(a[, axis])	Return the number of elements along a given axis.
slantrange(R, tau)	Slant range equation
xyz2llh_deg(x, y, z[, ell])	Convert cartesian coordinates to geocentric
xyz2llh_rad(x, y, z[, ell])	Convert cartesian coordinates to geocentric

Classes

Coor	Base class for coordinate conversion.
CoorAffineTransform(geotransform, projection)
CoorGCP
CoorGeo(bottom_left_east_north, ...)
CoorGeoBase(bottom_left_east_north, ...)
CoorGeoCSK(bottom_left_east_north, ...[, ...])
CoorGrid(rows, cols, lats, lons[, hs])	Image to geodetic coordinate conversion based on [(r,c),(l,l,h)] grid.
CoorGridBSpline(rows, cols, lats, lons[, hs])
CoorGridTimes(t0, tau0, dt, dtau, rows, ...)	Image, geodetic and time coordinate conversion based on [(r,c),(l,l,h)] grid.
CoorGridTimesGround(t0, tau0, dt, dtau, ...)	Image, geodetic and time coordinate conversion based on [(r,c),(l,l,h)] grid in ground projection.
CoorGridTimesSlant(t0, tau0, dt, dtau, rows, ...)	Image, geodetic and time coordinate conversion based on [(r,c),(l,l,h)] grid in slant projection.
CoorPhModel(t0, tau0, dt, dtau, lambda_, t, p)	Image, geodetic, cartesian and time coordinate conversion based on physical model.
CoorPhModelGround(t0, tau0, dt, dtau, ...[, ...])	Image, geodetic, cartesian and time coordinate conversion based on physical model in ground projection.
CoorPhModelSlant(t0, tau0, dt, dtau, ...[, ...])	Image, geodetic, cartesian and time coordinate conversion based on physical model in slant projection.
GeoCart([ellipsoid_designator, projection_id])	Base class for conversion between geodetic and map projection coordinates.
GeoCartStere(centre_lat, centre_lon, ...[, ...])	Class for coordinate conversion between lat-long <-> Stereographic coordinates.
GeoCartUTM([utm_zone, north_emisphere, ...])	Class for coordinate conversion between lat-long <-> UTM coordinates.
GeoCartWGS84UPS([north_pole])	Class for coordinate conversion between WGS84 <-> UPS coordinates.
GridLLH2RC(rows, cols, lats, lons, hs)	Conversion (lat, lon, h) to (row, col) coordinates based on well know grid.
GridRC2LLH(rows, cols, lats, lons, hs)	Conversion (row, col) to (lat, lon, h) coordinates based on well know grid.
GroundProjLLH2RC(t0, tau0, dt, dtau, ...[, ...])	rc –> ttau –> xyz –> llh
GroundProjRC2LLH(t0, tau0, dt, dtau, ...[, ...])	rc –> ttau –> xyz –> llh
GroundProjRC2TTau(t0, tau0, dt, dtau, gr2sl, ...)	Image to geodetic coordinate conversion in slant projection.
GroundProjRC2XYZ(t0, tau0, dt, dtau, ...[, ...])	rc –> ttau –> xyz
GroundProjXYZ2RC(t0, tau0, dt, dtau, ...[, ...])	xyz –> ttau –> rc
SlantProjLLH2RC(t0, tau0, dt, dtau, lambda_, ...)	rc –> ttau –> xyz –> llh
SlantProjRC2LLH(t0, tau0, dt, dtau, lambda_, ...)	rc –> ttau –> xyz –> llh
SlantProjRC2TTau(t0, tau0, dt, dtau)	Image to geodetic coordinate conversion in slant projection.
SlantProjRC2XYZ(t0, tau0, dt, dtau, lambda_, ...)	rc –> ttau –> xyz
SlantProjXYZ2RC(t0, tau0, dt, dtau, lambda_, ...)	xyz –> ttau –> rc
StateVectors(t, p, v[, a, att, attnor])
TTau2XYZ(lambda_, sttv[, ell, solver, fd])
XYZ2TTau(lambda_, sttv[, solver, fd])
matrix	matrix(data, dtype=None, copy=True)

Exceptions

InvalidCoorMapperError
InvalidInterpolatorError
InvalidProjectionError

class gsdview.coor.Coor

Base class for coordinate conversion.

geoToImgPoint(lat, lon, h=0)

Coordinate conversion: (lat,lon,h) –> (row,col).

imgToGeoPoint(row, col)

Coordinate conversion: (row,col) –> (lat,lon,h).

llh2rc()

rc2llh()

class gsdview.coor.CoorAffineTransform(geotransform, projection)

geoToImgPoint(lat, lon, h=0)

Coordinate conversion: (lat,lon) –> (row,col).

imgToGeoPoint(row, col, fd=0.0)

Coordinate conversion: (row,col) –> (lat,lon).

class gsdview.coor.CoorGCP

class gsdview.coor.CoorGeo(bottom_left_east_north, line_spacing, column_spacing, ellipsoid_designator, projection_id, map_projection_zone, map_projection_orientation)

llh2rc(lat, lon, h)

rc2llh(row, col)

class gsdview.coor.CoorGeoBase(bottom_left_east_north, line_spacing, column_spacing)

col2east(col)

east2col(east)

eastnorth2rc(east, north)

llh2rc()

north2row(north)

rc2eastnorth(row, col)

rc2llh()

row2north(row)

class gsdview.coor.CoorGeoCSK(bottom_left_east_north, line_spacing, column_spacing, ellipsoid_designator, projection_id, map_projection_zone, map_projection_orientation, map_projection_centre=None, map_projection_scale_factor=None, map_projection_false_easting=None, map_projection_false_northing=None)

class gsdview.coor.CoorGrid(rows, cols, lats, lons, hs=None)

Image to geodetic coordinate conversion based on [(r,c),(l,l,h)] grid.

A conversion grid is a set of image points whose geodetic coordinates are known. Given the conversion grid for initialization, then it is possible to convert (r,c)<->(l,l,h).

Parameters :

• rows: image rows of the conversion grid
• cols: image columns of the conversion grid
• lats: latitudes [deg/rad]
• lons: longitudes [deg/rad]
• hs: heights [m]

cols

img cols of the conversion grid

geoToImgPoint(lat, lon, h=0)

Coordinate conversion: (lat,lon) –> (row,col).

hs

heigths of the conversion grid [m]

imgToGeoPoint(row, col, fd=0.0)

Coordinate conversion: (row,col) –> (lat,lon,h).

lats

latitudes of the conversion grid (deg/rad)

llh2rc(lat, lon, h, mode='deg')

Returns (row,col) image coodinates given (lat,lon,h) geodetic coordinates.

Parameters :

• lat: latitude of the point to convert (deg/rad)
• lon: longitude of the point to convert (deg/rad)
• h: height of the point to convert [m] (optional)
• mode: no use

lons

longitudes of the conversion grid (deg/rad)

rc2llh(row, col, mode='deg', fd=0.0)

Returns (lat,lon,h) geodetic coordinates given (row,col) image coodinates.

Parameters :

• row: image row of the point to convert
• col: image column of the point to convert
• mode: no use

rows

img rows of the conversion grid

class gsdview.coor.CoorGridBSpline(rows, cols, lats, lons, hs=None)

geoToImgPoint(lat, lon, h=0)

Coordinate conversion: (lat,lon) –> (row,col).

imgToGeoPoint(row, col, fd=0.0)

Coordinate conversion: (row,col) –> (lat,lon,h).

llh2rc(lat, lon, h, mode='deg')

Returns (row,col) image coodinates given (lat,lon,h) geodetic coordinates.

Parameters :

• lat: latitude of the point to convert (deg/rad)
• lon: longitude of the point to convert (deg/rad)
• h: height of the point to convert [m] (optional)
• mode: no use

rc2llh(row, col, mode='deg', fd=0.0)

Returns (lat,lon,h) geodetic coordinates given (row,col) image coodinates.

Parameters :

• row: image row of the point to convert
• col: image column of the point to convert
• mode: no use

class gsdview.coor.CoorGridTimes(t0, tau0, dt, dtau, rows, cols, lats, lons, hs=None, ell='WGS84')

Image, geodetic and time coordinate conversion based on [(r,c),(l,l,h)] grid.

A conversion grid is a set of image points whose geodetic coordinates are known.

Parameters :

• t0: azimuth first time [s]
• tau0: range first time [s]
• dt: time line interval [s]
• dtau: column time interval [s]
• rows: image rows of the conversion grid
• cols: image columns of the conversion grid
• lats: latitudes (deg/rad)
• lons: longitudes (deg/rad)
• hs: heights [m]
• ell: ellipsoid descriptor ‘WGS84’

dt

Line Time Interval [s]

dtau

Column Time Interval [s]

ell

rc2ttau()

Returns (t,tau) time coordinates given (row,col) image coodinates.

Parameters :

• row: image row of the point to convert
• col: image column of the point to convert

t0

Azimuth First Time [s]

tau0

Range First Time [s]

ttau2rc()

Returns (row,col) image coordinates given (t,tau) time coordinates.

Parameters :

• t: azimuth time [s]
• tau: range time [s]

class gsdview.coor.CoorGridTimesGround(t0, tau0, dt, dtau, rows, cols, lats, lons, hs=None, ell='WGS84')

Image, geodetic and time coordinate conversion based on [(r,c),(l,l,h)] grid in ground projection.

A conversion grid is a set of image points whose geodetic coordinates are known.

Parameters :

• t0: azimuth first time [s]
• tau0: range first time [s]
• dt: time line interval [s]
• dtau: column time interval [s]
• rows: image rows of the conversion grid
• cols: image columns of the conversion grid
• lats: latitudes (deg/rad)
• lons: longitudes (deg/rad)
• hs: heights [m]
• ell: ellipsoid descriptor ‘WGS84’

rc2ttau(row, col)

Returns (t,tau) time coordinates given (row,col) image coodinates.

Parameters :

• row: image row of the point to convert
• col: image column of the point to convert

ttau2rc(t, tau)

Returns (row,col) image coordinates given (t,tau) time coordinates.

Parameters :

• t: azimuth time [s]
• tau: range time [s]

class gsdview.coor.CoorGridTimesSlant(t0, tau0, dt, dtau, rows, cols, lats, lons, hs=None, ell='WGS84')

Image, geodetic and time coordinate conversion based on [(r,c),(l,l,h)] grid in slant projection.

A conversion grid is a set of image points whose geodetic coordinates are known.

Parameters :

• t0: azimuth first time [s]
• tau0: range first time [s]
• dt: time line interval [s]
• dtau: column time interval [s]
• rows: image rows of the conversion grid
• cols: image columns of the conversion grid
• lats: latitudes (deg/rad)
• lons: longitudes (deg/rad)
• hs: heights [m]
• ell: ellipsoid descriptor ‘WGS84’

rc2ttau(row, col)

Returns (t,tau) time coordinates given (row,col) image coodinates.

Parameters :

• row: image row of the point to convert
• col: image column of the point to convert

ttau2rc(t, tau)

Returns (row,col) image coordinates given (t,tau) time coordinates.

Parameters :

• t: azimuth time [s]
• tau: range time [s]

class gsdview.coor.CoorPhModel(t0, tau0, dt, dtau, lambda_, t, p, v=None, a=None, att=None, attnor=None, ell='WGS84', Pllh_guess=None)

Image, geodetic, cartesian and time coordinate conversion based on physical model.

The physical model is based on satellite motion information.

Parameters :

• t0: azimuth first time [s]
• tau0: range first time [s]
• dt: time line interval [s]
• dtau: column time interval [s]
• lambda_: wavelength [m]
• t: 1D array-like state vector times [s]
• p: 3D array-like ECEF satellite positions [m]
• v: 3D array-like ECEF satellite velocities [m/s]
• a: 3D array-like ECEF satellite accelerations [m/s2]
• att: array-like satellite attidutes
• attnor: array-like attnor
• ell: ellipsoid descriptor ‘WGS84’
• Pllh_guess: initial llh point for the iterative algorithm

a

ECEF Satellite Acceleration [m/s2]

dt

Line Time Interval [s]

dtau

Column Time Interval [s]

ell

geoToImgPoint(lat, lon, h=0)

Coordinate conversion: (lat,lon) –> (row, col).

imgToGeoPoint(row, col, fd=0.0)

Coordinate conversion: (row, col) –> (lat,lon).

lambda_

Wavelength [Hz]

llh2rc()

Returns (row,col) image coodinates given (lat,lon,h) geodetic coordinates.

Parameters :

• lat: latitude of the point to convert (deg/rad)
• lon: longitude of the point to convert (deg/rad)
• h: height of the point to convert [m] (optional)
• mode: degree (deg) or radians (rad)
• Pllh_guess: (lat,lon,h) guess point to start iterative algorithm

p

ECEF Satellite Positions [m]

rc2llh()

Returns (lat,lon,h) geodetic coordinates given (row,col) image coodinates.

Parameters :

• row: image row of the point to convert
• col: image column of the point to convert
• mode: degree (deg) or radians (rad)
• Pllh_guess: (lat,lon,h) guess point to start iterative algorithm

t

State Vector Times [s]

t0

Azimuth First Time [s]

tau0

Range First Time [s]

v

ECEF Satellite Velocity [m/s]

class gsdview.coor.CoorPhModelGround(t0, tau0, dt, dtau, lambda_, t, p, v, a, att=None, attnor=None, ell='WGS84', gr2sl=None, sl2gr=None, range_ref_m=None, gr_surface=None, light_speed=None, Pllh_guess=None)

Image, geodetic, cartesian and time coordinate conversion based on physical model in ground projection.

The physical model is based on satellite motion information.

Parameters :

• t0: azimuth first time [s]
• tau0: range first time [s]
• dt: time line interval [s]
• dtau: column time interval [s]
• lambda_: wavelength [m]
• t: 1D array-like state vector times [s]
• p: 3D array-like ECEF satellite positions [m]
• v: 3D array-like ECEF satellite velocities [m/s]
• a: 3D array-like ECEF satellite accelerations [m/s2]
• att: array-like satellite attidutes
• attnor: array-like attnor
• ell: ellipsoid descriptor ‘WGS84’
• gr2sl: Ground to Slant Polynomial Coefficients
• sl2gr: Slant to Ground Polynomial Coefficients
• range_ref_m: Ground Projection Polynomial Reference Range
• light_speed: the speed of light
• gr_surface: Ground Projection Reference Surface
• Pllh_guess: initial llh point for the iterative algorithm

llh2rc(lat, lon, h, mode='deg', Prc_guess=None)

Returns (row,col) image coodinates given (lat,lon,h) geodetic coordinates.

Parameters :

• lat: latitude of the point to convert (deg/rad)
• lon: longitude of the point to convert (deg/rad)
• h: height of the point to convert [m] (optional)
• mode: degree (deg) or radians (rad)
• Pllh_guess: (lat,lon,h) guess point to start iterative algorithm

rc2llh(row, col, mode='deg', Pllh_guess=None)

Returns (lat,lon,h) geodetic coordinates given (row,col) image coodinates.

Parameters :

• row: image row of the point to convert
• col: image column of the point to convert
• mode: degree (deg) or radians (rad)
• Pllh_guess: (lat,lon,h) guess point to start iterative algorithm

class gsdview.coor.CoorPhModelSlant(t0, tau0, dt, dtau, lambda_, t, p, v=None, a=None, att=None, attnor=None, ell='WGS84', Pllh_guess=None)

Image, geodetic, cartesian and time coordinate conversion based on physical model in slant projection.

The physical model is based on satellite motion information.

Parameters :

• t0: azimuth first time [s]
• tau0: range first time [s]
• dt: time line interval [s]
• dtau: column time interval [s]
• lambda_: wavelength [m]
• t: 1D array-like state vector times [s]
• p: 3D array-like ECEF satellite positions [m]
• v: 3D array-like ECEF satellite velocities [m/s]
• a: 3D array-like ECEF satellite accelerations [m/s2]
• att: array-like satellite attidutes
• attnor: array-like attnor
• ell: ellipsoid descriptor ‘WGS84’
• Pllh_guess : initial llh point for the iterative algorithm

imgToGeoPoint(row, col, fd=0.0)

Coordinate conversion: (row, col) –> (lat,lon).

Override standard method

llh2rc(lat, lon, h=0, mode='deg', Prc_guess=None, fd=0.0)

Returns (row,col) image coodinates given (lat,lon,h) geodetic coordinates.

Parameters :

• lat: latitude of the point to convert (deg/rad)
• lon: longitude of the point to convert (deg/rad)
• h: height of the point to convert [m] (optional)
• mode: degree (deg) or radians (rad)
• Pllh_guess: (lat,lon,h) guess point to start iterative algorithm

rc2llh(row, col, mode='deg', Pllh_guess=None, fd=0.0)

Returns (lat,lon,h) geodetic coordinates given (row,col) image coodinates.

Parameters :

• row: image row of the point to convert
• col: image column of the point to convert
• mode: degree (deg) or radians (rad)
• Pllh_guess: (lat,lon,h) guess point to start iterative algorithm
• fd : Doppler Centroid value

class gsdview.coor.GeoCart(ellipsoid_designator=None, projection_id=None)

Base class for conversion between geodetic and map projection coordinates.

The physical model is based on satellite motion information.

Parameters :

• ellipsoid_designator: ellipsoid descriptor (‘WGS84’)
• projection_id: identifier of the projection (‘UTM’, ‘UPs’)

cart2geo()

Returns geodetic coodinates given cartographic coordinates.

ellipsoid_designator

ellipsoid designator defines geodetic system

geo2cart()

Returns cartographic coodinates given geodetic coordinates.

projection_id

cartographic projection id

class gsdview.coor.GeoCartStere(centre_lat, centre_lon, false_east, false_north, scale_factor=1.0, ellipsoid_designator='WGS84')

Class for coordinate conversion between lat-long <-> Stereographic coordinates.

Parameters :

• lat_0: Latitude at projection center
• lon_0: Longitude at projection center
• x_0: False Easting
• y_0: False Northing
• ellipsoid_designator: name of datum to use (in osr_datums()) [default ‘WGS84’]

cart2geo(east, north, height=0.0)

Returns geodetic coodinates given UPS cartographic coordinates.

Parameters :

• east: east coordinate [m]
• north: north coordinate [m]
• heigth: height elevation [m], optional

Output :

• (lat,lon,h): - lat: latitude - lon: longitude - h: height

geo2cart(lat, lon, h=0.0)

Returns Stereographic coordinates given geodetic coodinates.

Parameters:

• lat: latidute [decimal degrees]
• lon: longitude [decimal degrees]
• heigth: [m], optional

Output :

• east: east coordinate [m]
• north: north coordinate [m]
• height: [m]

class gsdview.coor.GeoCartUTM(utm_zone=None, north_emisphere=1, ellipsoid_designator='WGS84')

Class for coordinate conversion between lat-long <-> UTM coordinates.

Parameters :

• utm_zone: UTM zone (optional, in case it is known)
• north_emisphere: int - ‘1’ indicates the north emisphere - ‘0’ indicates the south emisphere - ‘ellipsoid_designator’ name of datum to use (in osr_datums()) [default ‘WGS84’]

Note

UTM projection is applied for latitudes <= 84deg and >= -80deg.

cart2geo(east, north, height=0.0, utm_zone=None, north_emisphere=1)

Returns geodetic coodinates given UTM cartographic coordinates.

Parameters :

• east: UTM east coordinate [m]
• north: UTM north coordinate [m]
• heigth: UTM height elevation [m], optional
• utm_zone: UTM zone (it can be set when init the class)
• north_emisphere: int - ‘1’ indicates the northern emisphere - ‘0’ indicates the southern emisphere

Output :

• (lat,lon,h): - lat: latitude [decimal degrees] - lon: longitude [decimal degrees] - h: height [m]

geo2cart(lat, lon, h=0.0)

Returns UTM cartographic coordinates given geodetic coodinates.

Parameters:

• lat: latidute [decimal degrees]
• lon: longitude [decimal degrees]
• heigth: [m], optional

Output :

• (east, north, h, utm_zone) - east: UTM east coordinate [m] - north: UTM north coordinate [m] - height: [m] - utm_zone: UTM zone

utm_zone

UTM projection zone

class gsdview.coor.GeoCartWGS84UPS(north_pole=1)

Class for coordinate conversion between WGS84 <-> UPS coordinates.

Parameters :

• north_pole: int - ‘1’ indicates the north pole (default) - ‘0’ indicates the south pole

Note

UPS projection is applied for latitudes >= 84deg or <= -80deg.

cart2geo(east, north, height=0.0, north_pole=None)

Returns WGS84 geodetic coodinates given UPS cartographic coordinates.

Parameters :

• east: UPS east coordinate [m]
• north: UPS north coordinate [m]
• heigth: UPS height elevation [m], optional
• north_pole: int - ‘1’ indicates the north pole - ‘0’ indicates the south pole

Output :

• (lat,lon,h) WGS84: - lat: latitude - lon: longitude - h: height

geo2cart(lat, lon, h=0.0)

Returns UPS cartographic coordinates given WGS84 geodetic coodinates.

Parameters :

• lat: latidute [decimal degrees]
• lon: longitude [decimal degrees]
• heigth: [m], optional

Output :

• (east, north, h, north pole) - east: UPS east coordinate [m] - north: UPS north coordinate [m] - height: [m] - north pole: int

  • ‘1’ indicates the north pole
  • ‘0’ indicates the south pole

north_pole

North pole flag

class gsdview.coor.GridLLH2RC(rows, cols, lats, lons, hs)

Conversion (lat, lon, h) to (row, col) coordinates based on well know grid.

Parameters :

• rows: array of int. Rows of the well known grid.
• cols: array of int. Columns of the well known grid.
• lats: array of float. Latitudes of the well known grid.
• lons: array of float. Longitudes of the well known grid.
• hs: array of float. Heights of the well known grid.

llh2rc(Pll)

Conversion (row, col) to (lat, lon, h) coordinates based on well know grid.

Parameters :

• Pllh: array of float. Point [lat, lon, h] to be converted.

class gsdview.coor.GridRC2LLH(rows, cols, lats, lons, hs)

Conversion (row, col) to (lat, lon, h) coordinates based on well know grid.

Parameters :

• rows: array of int. Rows of the well known grid.
• cols: array of int. Columns of the well known grid.
• lats: array of float. Latitudes of the well known grid.
• lons: array of float. Longitudes of the well known grid.
• hs: array of float. Heights of the well known grid.

rc2llh(Prc)

Conversion (row, col) to (lat, lon, h) coordinates based on well know grid.

Parameters :

• Prc: array of int. Point [row, col] to be converted.

class gsdview.coor.GroundProjLLH2RC(t0, tau0, dt, dtau, lambda_, sttv, ell, gr2sl, sl2gr, range_ref_m, gr_surface, light_speed, solver=None)

rc –> ttau –> xyz –> llh

llh2rc_deg(Pllh, Prc_guess=None)

llh2rc_rad(Pllh, Prc_guess=None)

class gsdview.coor.GroundProjRC2LLH(t0, tau0, dt, dtau, lambda_, sttv, ell, gr2sl, sl2gr, range_ref_m, gr_surface, light_speed, solver=None, fd=0.0)

rc –> ttau –> xyz –> llh

rc2llh_deg(Prc, Pllh_guess)

rc2llh_rad(Prc, Pllh_guess)

class gsdview.coor.GroundProjRC2TTau(t0, tau0, dt, dtau, gr2sl, sl2gr, range_ref_m, gr_surface, light_speed)

Image to geodetic coordinate conversion in slant projection.

Parameters :

• t0: azimuth first time
• tau0: range first time
• dt: time line interval
• dtau: time column interval

c2tau(col)

tau2c(tau)

class gsdview.coor.GroundProjRC2XYZ(t0, tau0, dt, dtau, lambda_, sttv, ell, gr2sl, sl2gr, range_ref_m, gr_surface, light_speed, solver=None, fd=0.0)

rc –> ttau –> xyz

rc2xyz(Prc, Pxyz_guess=None)

class gsdview.coor.GroundProjXYZ2RC(t0, tau0, dt, dtau, lambda_, sttv, gr2sl, sl2gr, range_ref_m, gr_surface, light_speed, solver=None)

xyz –> ttau –> rc

xyz2rc(Pxyz, Prc_guess=None)

exception gsdview.coor.InvalidCoorMapperError

exception gsdview.coor.InvalidInterpolatorError

exception gsdview.coor.InvalidProjectionError

class gsdview.coor.SlantProjLLH2RC(t0, tau0, dt, dtau, lambda_, sttv, ell='WGS84', solver=None)

rc –> ttau –> xyz –> llh

llh2rc_deg(Pllh, Prc_guess=None)

llh2rc_rad(Pllh, Prc_guess=None)

class gsdview.coor.SlantProjRC2LLH(t0, tau0, dt, dtau, lambda_, sttv, ell='WGS84', solver=None, fd=0.0, Pllh_guess=None)

rc –> ttau –> xyz –> llh

rc2llh_deg(Prc, Pllh_guess, fd=0.0)

rc2llh_rad(Prc, Pllh_guess, fd=0.0)

class gsdview.coor.SlantProjRC2TTau(t0, tau0, dt, dtau)

Image to geodetic coordinate conversion in slant projection.

Parameters :

• t0: azimuth first time
• tau0: range first time
• dt: time line interval
• dtau: time column interval

c2tau(col)

fs_azimuth

fs_range

r2t(row)

rc2ttau(Prc)

t2r(t)

tau2c(tau)

ttau2rc(Pttau)

class gsdview.coor.SlantProjRC2XYZ(t0, tau0, dt, dtau, lambda_, sttv, ell='WGS84', solver=None, fd=0.0)

rc –> ttau –> xyz

rc2xyz(Prc, Pxyz_guess=None, fd=0.0)

class gsdview.coor.SlantProjXYZ2RC(t0, tau0, dt, dtau, lambda_, sttv, solver=None)

xyz –> ttau –> rc

xyz2rc(Pxyz, Prc_guess=None)

class gsdview.coor.StateVectors(t, p, v, a=None, att=None, attnor=None)

acceleration(t)

attitude(t)

attitude_normal(t)

position(t)

velocity(t)

class gsdview.coor.TTau2XYZ(lambda_, sttv, ell='WGS84', solver=None, fd=0.0)

ellipsoid_descriptor

fd

lambda_

solver

ttau2xyz(Pttau, Pxyz_guess=None)

class gsdview.coor.XYZ2TTau(lambda_, sttv, solver=None, fd=0.0)

fd

lambda_

solver

xyz2ttau(Pxyz, Pttau_guess=None)

gsdview.coor.doppler(R, V, lambda_, fd=0)

Doppler equation

Parameters :

• R: satellite target distace vector
• V: satellite target relative velocity vector
• lambda_: radar wavelenght
• fd: Doppler frequency (default=0)

gsdview.coor.ellipsoid(P, ell='WGS84')

Ellipsoid equation

Parameters :

• P: cartesian coordnates (3D) of the point
• ell: a vactor containing the three semi-axis of the ellipsoid

gsdview.coor.llh2xyz_deg(lat, lon, h, ell='WGS84')

Convert geocentric coordinates to cartesian

gsdview.coor.llh2xyz_rad(lat, lon, h, ell='WGS84')

Convert geocentric coordinates to cartesian

OGP Surveying and Positioning Guidance Note number 7, part 2 – August 2008 To facilitate improvement, this document is subject to revision. The current version is available at www.epsg.org.

gsdview.coor.only_projected_projs()

gsdview.coor.osr_datums()

Return most significant datums

gsdview.coor.projections()

Return most significant projections

gsdview.coor.slantrange(R, tau)

Slant range equation

Parameters :

• R: satellite target distace vector
• tau: two way slant rage time

gsdview.coor.xyz2llh_deg(x, y, z, ell='WGS84')

Convert cartesian coordinates to geocentric

gsdview.coor.xyz2llh_rad(x, y, z, ell='WGS84')

Convert cartesian coordinates to geocentric

OGP Surveying and Positioning Guidance Note number 7, part 2 – August 2008 To facilitate improvement, this document is subject to revision. The current version is available at www.epsg.org.