Suite

Matrice de rotation dans l'EPSG

Matrice de rotation dans l'EPSG


J'utilise Proj4j pour un WebGIS et je dois créer un nouvel ESPG. J'ai essayé sur QGIS de convertir les coordonnées d'un système local (en Grèce) en un système de coordonnées mondial.

Tout d'abord, je dois faire une transformation 2D : 2 translations (Tx ans Ty) et une rotation (Rz) d'environ 0,58931 radians autour de "l'axe z" (dans le plan du repère local) entre les coordonnées locales et les coordonnées grecques (EGSA87).

Je pense donc que la matrice de rotation est :[cos(Rz),-sin(Rz);sin(Rz),cos(Rz)]et sans le cos et le péché que j'ai[0.8313242,-0.5557878;0.5557878,0.8313242].

Ensuite, je peux réutiliser les paramètres de la définition EPSG:2100 :

+proj=tmerc +lat_0=0 +lon_0=24 +k=0.9996 +x_0=500000 +y_0=0 +ellps=GRS80 +towgs84=-199.87,74.79,246.62,0,0,0,0 +units=m + no_defs

Je pense que je peux insérer les traductions (Tx et Ty) dans les x_0 et y_0, mais je ne sais pas comment insérer la matrice de rotation. Comment insérer une transformation avec une matrice dans un nouvel EPSG ?

J'ai trouvé sur Internet que le+xform=La commande semble bonne, mais je ne sais pas comment l'utiliser. J'ai essayé:

+proj=tmerc +lat_0=0 +lon_0=24 +k=0.9996 +x_0=-869388.894 +y_0=-4260031.458 +xform=[0.8313242,-0.5557878,0.5557878,0.8313242] +ellps=GRS80 +towgs84=-199.87,74.79 ,246.62,0,0,0,0 +units=m +no_defs

Mais apparemment ce n'est pas une bonne solution. Alors, est-ce que quelqu'un sait comment je peux utiliser cette commande dans mon nouvel EPSG ou où je peux trouver la documentation proj4 ?


Vous pouvez configurer une projection de mercator oblique locale pour obtenir la rotation, comme je l'ai décrit dans ma réponse ici :

Vous utilisez un système de coordonnées personnalisé dans ArcGIS Desktop ?

Votre paramètre de rotation+alphaserait de 33,765 degrés. Vous devez exprimer le signe à utiliser.

Les paramètres Proj.4 sont répertoriés ici : http://trac.osgeo.org/proj/wiki/GenParms, mais xform n'est pas inclus.


Matrice de rotation dans l'EPSG - Systèmes d'Information Géographique

Les paramètres Bursa-Wolf ont deux objectifs :

    Repli pour les décalages de référence
    Il existe différentes méthodes pour transformer les coordonnées d'un datum géodésique en un autre datum, et les paramètres Bursa-Wolf sont utilisés avec certaines d'entre elles. Cependant, différents ensembles de paramètres peuvent exister pour la même paire de données (source, cible), il n'est donc souvent pas suffisant de connaître ces données. Le couple ( source , cible ) de SIR est souvent nécessaire, parfois accompagné de l'étendue géographique des coordonnées à transformer.

Apache SIS recherche les méthodes de décalage de données (y compris les paramètres Bursa-Wolf) dans la base de données EPSG lorsqu'une CoordinateOperation ou une MathTransform est demandée pour une paire de CRS. C'est ce qu'on appelle l'approche de la liaison tardive. Si une méthode de décalage de données est trouvée dans la base de données, elle aura la priorité sur toute instance BursaWolfParameters associée à ce DefaultGeodeticDatum . Ce n'est que si aucune méthode de décalage de datum n'est trouvée dans la base de données, que les BursaWolfParameters associés à la datum peuvent être utilisés comme solution de secours.


Spécification des paramètres de l'étalonnage du montage

Les paramètres de l'étalonnage du montage peuvent être spécifiés à l'aide d'une chaîne selon la syntaxe EBNF (Extended Backus-Naur Form) suivante :

Les paramètres de montage suivants peuvent être spécifiés :

décalage horaire: C'est le décalage horaire du système de temps de trajectoire par rapport au système de scanner laser (c'est-à-dire que cette valeur doit être soustraite d'un horodatage de trajectoire afin d'obtenir le même moment dans le système de temps du scanner)

système de balayage: Une séquence de 3 directions hors de l'ensemble est utilisé pour spécifier le attitude approximative du système de balayage à inclinaison zéro - note : PAS celui du système de scanner (!) - par rapport au système corporel (voir figure 1). La séquence de direction est limitée aux 24 systèmes possibles pour droitiers. Sur la base de cette entrée, la matrice de rotation (mentionné dans la section précédente) est déterminé. Par exemple, dans le cas de la figure 1, les 3 axes du système pointer vers le bas, l'arrière et la gauche (ici, il n'est pas précisé lequel des axes est X,Y,Z). Comme nous ne prenons en charge que les systèmes pour droitiers, la configuration de cet exemple correspond à l'une des 3 séquences possibles "D-B-L", "B-L-D", "L-D-B".

Par exemple, "SCANNERSYS(D-B-L)" renvoie .

rotation de montage: La rotation du système dans le système

décalage de montage: La traduction (shift) du système dans le système

rotation d'inclinaison: La rotation du système dans le système

changement d'inclinaison: La traduction (shift) du système dans le système

Ainsi, du point de vue de ces éléments de transformation, les systèmes mentionnés ci-dessus jouent le rôle de système global ou de système local.

La section suivante décrit comment spécifier une rotation et/ou un décalage à l'aide de la syntaxe fournie.


6.3 Prise en charge des coordonnées locales

Spatial fournit un niveau de prise en charge des systèmes de coordonnées locaux. Les systèmes de coordonnées locales sont souvent utilisés dans les systèmes de CAO, et ils peuvent également être utilisés dans les enquêtes locales où la relation entre le site étudié et le reste du monde n'est pas importante.

Plusieurs systèmes de coordonnées locaux sont prédéfinis et inclus avec Spatial dans la table SDO_COORD_REF_SYS (décrite dans la Section 6.7.9). Ces systèmes de coordonnées locaux fournis, dont les noms commencent par Non-Earth , définissent des systèmes de coordonnées cartésiens non terrestres basés sur différentes unités de mesure ( Meter , Millimeter , Inch , etc.).

Dans la version actuelle, vous ne pouvez pas effectuer de transformation de système de coordonnées entre des systèmes de coordonnées locaux et terrestres et lors de la transformation d'une géométrie ou d'une couche de géométries entre des systèmes de coordonnées locaux, vous pouvez uniquement convertir les coordonnées d'un système de coordonnées local d'une unité de mesure en un autre (par exemple, pouces à millimètres). Cependant, vous pouvez effectuer toutes les autres opérations spatiales (par exemple, en utilisant SDO_RELATE, SDO_WITHIN_DISTANCE et d'autres opérateurs) avec des systèmes de coordonnées locaux.


Constantes IPJ_ORIENT¶

pas d'orientation particulière - vue en plan. Toutes les vues des cartes créées avant la v5.1.3 renverront cette valeur.

gxapi. IPJ_ORIENT_DEFAULT = 0 IPJ_ORIENT_PLAN

Une vue en plan avec une élévation de référence et une rotation facultative.

gxapi. IPJ_ORIENT_PLAN = 1 IPJ_ORIENT_SECTION ¶

A un azimut et un swing. La vue en coupe projette tous les objets tracés HORIZONTALEMENT sur le plan de visualisation afin de préserver les élévations, même si la coupe a une oscillation.

gxapi. IPJ_ORIENT_SECTION = 2 IPJ_ORIENT_SECTION_NORMAL ¶

Identique à IPJ_ORIENT_SECTION , mais la projection est perpendiculaire à la section, pas horizontale, donc les élévations ne sont pas conservées sur les sections inclinées.

gxapi. IPJ_ORIENT_SECTION_NORMAL = 5 IPJ_ORIENT_DEPTH_SECTION ¶

Cette section simple n'a pas d'azimut ou de balancement défini, seule la profondeur est importante, et elle est sortie en tant que paramètre Y, augmentant vers le bas. Utilisé (par exemple) pour les journaux de bande dans Wholeplot.

gxapi. IPJ_ORIENT_DEPTH_SECTION = 3 IPJ_ORIENT_3D ¶

Une orientation de rotation/mise à l'échelle/translation 3D

gxapi. IPJ_ORIENT_3D = 4 IPJ_ORIENT_3D_MATRIX ¶

gxapi. IPJ_ORIENT_3D_MATRIX = 7 IPJ_ORIENT_SECTION_CROOKED ¶

Il s'agit d'une section verticale qui suit un chemin sinueux, comme une rivière ou un cheminement d'arpentage. L'emplacement de la section horizontale est la distance le long du chemin, tandis que l'axe vertical donne l'élévation.

gxapi. IPJ_ORIENT_SECTION_CROOKED = 6


Matrice de rotation dans l'EPSG - Systèmes d'Information Géographique

Défini dans le module de référencement sis

Résumé de la méthode

Méthodes héritées de la classe Object

Détails de la méthode

ParseAxisDirection

  • Directions cardinales comme "nord" et "est" .
  • Les directions intercardinales "nord-est" et "sud-sud-est" , en utilisant soit '-' , '_' ou des espaces comme séparateur entre les points cardinaux.
  • Directions à partir d'un pôle comme "Sud à 180 degrés" et "Sud à 90° Est" , en utilisant le symbole "deg" ou "°". Notez que le méridien n'est pas nécessairement relatif à Greenwich (voir directionAlongMeridian(…) pour plus d'informations).

DirectionAlongMeridian

Angle

  • L'angle de l'EST au NORD est de 90°
  • L'angle du SUD à l'OUEST est de -90°
  • L'angle entre le « Nord le long de 90° Est » et le « Nord le long de 0° » est de 90°.

Tous les angles sont des approximations car cette méthode ne prend pas en compte la forme ellipsoïdale ou géoïde de la Terre.

ÉchangeEtEchelleAxes

RemplacerAxes

RemplacerUnitéLinéaire

Cette méthode pratique est équivalente au code suivant :

ReplaceAngularUnit

Cette méthode pratique est équivalente au code suivant :

GetAxisDirections

GetEpsgCode

La mise en œuvre actuelle utilise une liste codée en dur de systèmes de coordonnées connus, elle n'analyse pas encore la base de données EPSG (cela peut changer dans la future version d'Apache SIS). La liste actuelle des systèmes de coordonnées connus est donnée ci-dessous.

Systèmes de coordonnées connus (CS)
EPSG type CS Directions des axes Unité horizontale
6424 Ellipsoïdale est Nord degré
6422 Ellipsoïdale Nord est degré
6425 Ellipsoïdale est Nord diplômés
6403 Ellipsoïdale Nord est diplômés
6429 Ellipsoïdale est Nord radian
6428 Ellipsoïdale Nord est radian
6426 Ellipsoïdale est Nord en haut degré
6423 Ellipsoïdale Nord est en haut degré
6427 Ellipsoïdale est Nord en haut diplômés
6421 Ellipsoïdale Nord est en haut diplômés
6431 Ellipsoïdale est Nord en haut radian
6430 Ellipsoïdale Nord est en haut radian
4400 cartésien est Nord mètre
4500 cartésien Nord est mètre
4491 cartésien Ouest Nord mètre
4501 cartésien Nord Ouest mètre
6503 cartésien Ouest Sud mètre
6501 cartésien Sud Ouest mètre
1039 cartésien est Nord le pied
1029 cartésien Nord est le pied
4403 cartésien est Nord Le pied de Clarke
4502 cartésien Nord est Le pied de Clarke
4497 cartésien est Nord Pied d'enquête américain

GetEpsgCode

La mise en œuvre actuelle utilise une liste codée en dur de systèmes de coordonnées connus, elle n'analyse pas encore la base de données EPSG (cela peut changer dans la future version d'Apache SIS). La liste actuelle des systèmes de coordonnées connus est documentée ci-dessus.


Contenu

L'origine des coordonnées du WGS 84 est censée être située au centre de masse de la Terre, l'incertitude serait inférieure à 2 cm . [5]

Le méridien WGS 84 de longitude zéro est le méridien de référence IERS, [6] 5,3 secondes d'arc ou 102 mètres (335 pieds) à l'est du méridien de Greenwich à la latitude de l'observatoire royal. [7] [8]

La surface de référence WGS 84 est un sphéroïde aplati de rayon équatorial a = 6 378 137 m à l'équateur et d'aplatissement f = 1/ 298,257 223 563 . La valeur affinée de la constante gravitationnelle WGS 84 (masse de l'atmosphère terrestre incluse) est GM = 3 986 004 ,418 × 10 8 m 3 /s 2 . La vitesse angulaire de la Terre est définie comme étant ω = 72,921 15 × 10 −6 rad/s . [9]

Cela conduit à plusieurs paramètres calculés tels que le demi-petit axe polaire b qui est égal à une × (1 − F) = 6 356 752 .3142 m , et la première excentricité au carré, e 2 = 6.694 379 990 14 × 10 −3 . [9]

Actuellement, WGS 84 utilise le Earth Gravitational Model 2008. [10] Ce géoïde définit la surface nominale du niveau de la mer au moyen d'une série d'harmoniques sphériques de degré 2160. [11] C'est une amélioration par rapport au modèle EGM96 de 1996, qui était à son tour une amélioration par rapport au géoïde WGS 84 original, appelé EGM84. Les écarts du géoïde EGM96 par rapport à l'ellipsoïde de référence WGS 84 vont d'environ -105 m à environ +85 m . [12]

Le WGS 84 utilise actuellement le World Magnetic Model 2020. [13] La prochaine mise à jour régulière (WMM2025) aura lieu en décembre 2024.

Les efforts visant à compléter les divers systèmes nationaux d'arpentage ont commencé au 19ème siècle avec F.R. Le célèbre livre de Helmert Mathematische und Physikalische Theorien der Physikalischen Geodäsie (Théories mathématiques et physiques de la géodésie physique). L'Autriche et l'Allemagne ont fondé le Zentralbüro für die Internationale Erdmessung (Bureau central de géodésie internationale), et une série d'ellipsoïdes globaux de la Terre ont été dérivés (par exemple, Helmert 1906, Hayford 1910/1924).

Un système géodésique unifié pour le monde entier est devenu indispensable dans les années 1950 pour plusieurs raisons :

  • Science spatiale internationale et début de l'astronautique.
  • Le manque d'informations géodésiques intercontinentales.
  • L'incapacité des grands systèmes géodésiques, tels que le système de référence européen (ED50), le système de référence nord-américain (NAD) et le système de référence de Tokyo (TD), à fournir une base de données géographiques mondiale
  • Besoin de cartes mondiales pour la navigation, l'aviation et la géographie.
  • La préparation à la guerre froide occidentale nécessitait un système de référence géospatial standardisé à l'échelle de l'OTAN, conformément à l'Accord de normalisation de l'OTAN

À la fin des années 1950, le ministère de la Défense des États-Unis, en collaboration avec des scientifiques d'autres institutions et pays, a commencé à développer le système mondial nécessaire auquel les données géodésiques pourraient être référées et la compatibilité établie entre les coordonnées de sites d'intérêt très éloignés. Les efforts de l'armée américaine, de la marine et de l'armée de l'air ont été combinés pour aboutir au DoD World Geodetic System 1960 (WGS 60). Le terme données tel qu'il est utilisé ici, fait référence à une surface lisse définie arbitrairement comme une élévation nulle, cohérente avec un ensemble de mesures de distances entre différentes stations et des différences d'élévation, le tout réduit à une grille de latitudes, de longitudes et d'élévations. Les méthodes d'arpentage du patrimoine ont trouvé des différences d'élévation par rapport à une horizontale locale déterminée par le niveau à bulle, le fil à plomb ou un appareil équivalent qui dépend du champ de gravité local (voir géodésie physique). En conséquence, les altitudes dans les données sont référencées au géoïde, une surface qui n'est pas facilement trouvée en utilisant la géodésie par satellite. Cette dernière méthode d'observation est plus adaptée à la cartographie globale. Par conséquent, une motivation et un problème substantiel dans le WGS et les travaux similaires est de rassembler des données qui ont été non seulement faites séparément, pour différentes régions, mais de ré-référencer les élévations à un modèle ellipsoïde plutôt qu'au géoïde.

Dans la réalisation du WGS 60, une combinaison de données gravimétriques de surface disponibles, de données astro-géodétiques et de résultats des levés HIRAN [14] et canadiens SHORAN ont été utilisés pour définir un ellipsoïde le mieux adapté et une orientation centrée sur la Terre pour chacun des datum initialement sélectionnés. (Chaque donnée est relativement orientée par rapport à différentes parties du géoïde par les méthodes astro-géodétiques déjà décrites.) La seule contribution des données satellitaires au développement de WGS 60 était une valeur pour l'aplatissement de l'ellipsoïde qui a été obtenu à partir du mouvement nodal d'un satellite.

Avant le WGS 60, l'armée américaine et l'armée de l'air américaine avaient chacune développé un système mondial en utilisant différentes approches de la méthode d'orientation de référence gravimétrique. Pour déterminer leurs paramètres d'orientation gravimétrique, l'Air Force a utilisé la moyenne des différences entre les déflexions gravimétriques et astro-géodétiques et les hauteurs du géoïde (ondulations) à des stations spécifiquement sélectionnées dans les zones des principaux systèmes de référence. L'armée a effectué un ajustement pour minimiser la différence entre les géoïdes astro-géodétiques et gravimétriques. En faisant correspondre les géoïdes astro-géodétiques relatifs des datums sélectionnés avec un géoïde gravimétrique centré sur la terre, les datums sélectionnés ont été réduits à une orientation centrée sur la terre. Étant donné que les systèmes de l'armée et de l'armée de l'air s'accordaient remarquablement bien pour les zones NAD, ED et TD, ils ont été consolidés et sont devenus WGS 60.

Les améliorations apportées au système mondial comprenaient l'astrogéoïde d'Irene Fischer et le datum astronautique de Mercure. En janvier 1966, un comité du système géodésique mondial composé de représentants de l'armée, de la marine et de l'air des États-Unis a été chargé de développer un WGS amélioré, nécessaire pour répondre aux exigences de cartographie, de cartographie et de géodésie. Des observations de gravité de surface supplémentaires, des résultats de l'extension des réseaux de triangulation et de trilatération, et de grandes quantités de données Doppler et optiques satellitaires étaient devenus disponibles depuis le développement du WGS 60. En utilisant les données supplémentaires et des techniques améliorées, le WGS 66 a été produit qui a répondu aux besoins du DoD. pendant environ cinq ans après sa mise en œuvre en 1967. Les paramètres déterminants de l'ellipsoïde WGS 66 étaient l'aplatissement (1/298,25 déterminé à partir des données satellitaires) et le demi-grand axe (6 378 145 m déterminé à partir d'une combinaison de satellite Doppler et Les données). Un champ d'anomalies de gravité dans l'air libre moyen mondial de 5° × 5° a fourni les données de base pour la production du géoïde gravimétrique WGS 66. En outre, un géoïde référencé à l'ellipsoïde WGS 66 a été dérivé des données astrogéodétiques disponibles pour fournir une représentation détaillée de zones terrestres limitées.

Après un effort considérable sur une période d'environ trois ans, le système géodésique mondial du ministère de la Défense 1972 a été achevé. Des données satellitaires, gravimétriques de surface et astrogéodétiques sélectionnées disponibles jusqu'en 1972 à partir de sources DoD et non DoD ont été utilisées dans une solution WGS unifiée (un ajustement des moindres carrés à grande échelle). Les résultats de l'ajustement consistaient en des corrections des coordonnées initiales de la station et des coefficients du champ gravitationnel.

La plus grande collection de données jamais utilisée aux fins du WGS a été rassemblée, traitée et appliquée au développement du WGS 72. Des données satellitaires optiques et électroniques ont été utilisées. Les données satellitaires électroniques se composaient, en partie, de données Doppler fournies par la marine américaine et de stations de suivi par satellite non-DoD coopérantes établies à l'appui du système de navigation par satellite (NNSS) de la marine. Les données Doppler étaient également disponibles à partir des nombreux sites établis par GEOCEIVERS en 1971 et 1972. Les données Doppler étaient la principale source de données pour le WGS 72 (voir image). Des données satellitaires électroniques supplémentaires ont été fournies par le réseau équatorial SECOR (Sequential Collation of Range) complété par l'armée américaine en 1970. Les données satellitaires optiques du Worldwide Geometric Satellite Triangulation Program ont été fournies par le système de caméra BC-4 (voir image). Les données de l'observatoire d'astrophysique Smithsonian ont également été utilisées, notamment une caméra (Baker-Nunn) et des télémétries laser.

Le champ de gravité de surface utilisé dans la solution WGS unifiée consistait en un ensemble de 410 anomalies de gravité moyenne en air libre de 10° × 10° de surface égale déterminées uniquement à partir de données terrestres. Ce champ gravimétrique comprend des valeurs moyennes d'anomalies compilées directement à partir des données gravimétriques observées lorsque celles-ci étaient disponibles en quantité suffisante. La valeur pour les zones de données d'observation éparses ou inexistantes a été développée à partir d'approximations de gravité géophysiquement compatibles en utilisant des techniques de corrélation gravité-géophysique. Environ 45 pour cent des 410 valeurs moyennes d'anomalies gravimétriques à l'air libre ont été déterminées directement à partir des données gravimétriques observées.

Les données astrogéodétiques dans leur forme de base consistent en la déviation des composantes verticales rapportées aux différents systèmes géodésiques nationaux. Ces valeurs de déflexion ont été intégrées dans les cartes du géoïde astrogéodétique en référence à ces référentiels nationaux. Les hauteurs du géoïde ont contribué à la solution WGS unifiée en fournissant des données supplémentaires et plus détaillées pour les zones terrestres. Des données d'enquête au sol conventionnelles ont été incluses dans la solution pour imposer un ajustement cohérent des coordonnées des sites d'observation voisins des systèmes BC-4, SECOR, Doppler et Baker-Nunn. De plus, huit cheminements précis de longue ligne géodimétriques ont été inclus dans le but de contrôler l'échelle de la solution.

La solution WGS unifiée, comme indiqué ci-dessus, était une solution pour les positions géodésiques et les paramètres associés du champ gravitationnel basé sur une combinaison optimale de données disponibles. Les paramètres de l'ellipsoïde WGS 72, les décalages de référence et d'autres constantes associées ont été dérivés séparément. Pour la solution unifiée, une matrice d'équation normale a été formée sur la base de chacun des ensembles de données mentionnés. Ensuite, les matrices d'équations normales individuelles ont été combinées et la matrice résultante résolue pour obtenir les positions et les paramètres.

La valeur du demi-grand axe (a) de l'ellipsoïde WGS 72 est de 6 378 135 m. L'adoption d'une valeur a inférieure de 10 mètres à celle de l'ellipsoïde WGS 66 était basée sur plusieurs calculs et indicateurs, notamment une combinaison de données gravimétriques satellitaires et de surface pour la détermination de la position et du champ gravitationnel. Des ensembles de coordonnées de stations dérivées du satellite et de déviation gravimétrique des données de hauteur verticale et du géoïde ont été utilisés pour déterminer les décalages de référence locaux à géocentriques, les paramètres de rotation de référence, un paramètre d'échelle de référence et une valeur pour le demi-grand axe de l'ellipsoïde WGS. Huit solutions ont été faites avec les différents ensembles de données d'entrée, à la fois d'un point de vue d'enquête et aussi en raison du nombre limité d'inconnues qui pourraient être résolues dans une solution individuelle en raison des limitations informatiques. Des stations sélectionnées de poursuite par satellite Doppler et d'orientation de systèmes géodésiques ont été incluses dans les diverses solutions. Sur la base de ces résultats et d'autres études connexes réalisées par le Comité, une valeur a de 6 378 135 m et un aplatissement de 1/298,26 ont été adoptés.

Dans le développement de décalages de datum locaux vers WGS 72, les résultats de différentes disciplines géodésiques ont été étudiés, analysés et comparés. Ces décalages adoptés étaient basés principalement sur un grand nombre de coordonnées de stations Doppler TRANET et GEOCEIVER qui étaient disponibles dans le monde entier. Ces coordonnées avaient été déterminées à l'aide de la méthode de positionnement ponctuel Doppler.


Classe GXIPJ¶

Azimut (degrés CW) du nord géographique à partir du nord de la grille à un emplacement.

  • argument1 ( geosoft.gxapi.float_ref ) – Origine du plan X
  • arg2 ( geosoft.gxapi.float_ref ) – Origine du plan Y
  • arg3 ( geosoft.gxapi.float_ref ) – Plan Origine Z
  • arg4 ( geosoft.gxapi.float_ref ) – Azimut plan (section) ou Rotation (plan)
  • arg5 ( geosoft.gxapi.float_ref ) – Plan Swing (section)
  • argument1 (flotter) – Min X de la surface
  • arg2 (flotter) – Min Y de surface
  • arg3 (flotter) – Max X de la surface
  • arg4 (flotter) – Max Y de surface
  • arg5 ( geosoft.gxapi.float_ref ) – Angle de tangage (entre -360 et 360)
  • arg6 ( geosoft.gxapi.float_ref ) – Angle de lacet (entre -360 et 360)
  • arg7 ( geosoft.gxapi.float_ref ) – Angles de roulis (entre -360 et 360)
  • arg8 ( geosoft.gxapi.float_ref ) – Décalage X du plan
  • arg9 ( geosoft.gxapi.float_ref ) – décalage Y du plan
  • arg10 ( geosoft.gxapi.float_ref ) – décalage Z du plan
  • arg11 ( geosoft.gxapi.float_ref ) – Échelle X
  • arg12 ( geosoft.gxapi.float_ref ) – Échelle Y
  • arg13 ( geosoft.gxapi.float_ref ) – Échelle Z
  • argument1 ( geosoft.gxapi.GXIPJ ) – objet IPJ pour les valeurs de sortie
  • arg2 (flotter) – Min X de la surface (en coordonnées de grille)
  • arg3 (flotter) – Min Y de surface
  • arg4 (flotter) – Max X de la surface
  • arg5 (flotter) – Max Y de surface
  • arg6 ( geosoft.gxapi.float_ref ) – Angle d'inclinaison (entre -360 et 360) (dans les coordonnées de vue)
  • arg7 ( geosoft.gxapi.float_ref ) – Angle de lacet (entre -360 et 360)
  • arg8 ( geosoft.gxapi.float_ref ) – Angles de roulis (entre -360 et 360)
  • arg9 ( geosoft.gxapi.float_ref ) – Décalage X du plan (dans les coordonnées de la vue)
  • arg10 ( geosoft.gxapi.float_ref ) – Décalage Y du plan
  • arg11 ( geosoft.gxapi.float_ref ) – décalage Z du plan
  • arg12 ( geosoft.gxapi.float_ref ) – Échelle X (dans les coordonnées de la vue)
  • arg13 ( geosoft.gxapi.float_ref ) – Échelle Y
  • arg14 ( geosoft.gxapi.float_ref ) – Échelle Z

get_units ( (float_ref)arg1, (str_ref)arg2 ) &rarr Aucun :¶

Vérifiez si l'instance de geosoft.gxapi.GXIPJ est nulle (non définie)

Retour:True s'il s'agit d'une instance nulle de geosoft.gxapi.GXIPJ , False sinon.
Type de retour :bool`
make_geographic ( ) &rarr Aucun :¶

  • argument1 (flotter) – longitude minimale
  • arg2 (flotter) – latitude minimale
  • arg3 (flotter) – longitude maximale
  • arg4 (flotter) – latitude maximale
  • argument1 ( geosoft.gxapi.GXIPJ ) – nouvel IPJ
  • arg2 (flotter) – résolution des données dans l'IPJ d'origine
  • arg3 (flotter) – X minimum de cadre de délimitation dans le nouvel IPJ
  • arg4 (flotter) – Y minimum
  • arg5 (flotter) – X maximum
  • arg6 (flotter) – Y maximum
  • arg7 ( geosoft.gxapi.float_ref ) – résolution minimale des données dans le nouvel IPJ,
  • arg8 ( geosoft.gxapi.float_ref ) – résolution maximale des données dans le nouvel IPJ
  • arg9 ( geosoft.gxapi.float_ref ) – résolution diagonale des données dans le nouvel IPJ

Une instance nulle (non définie) de geosoft.gxapi.GXIPJ

Retour:Un geosoft.gxapi.GXIPJ nul
Type de retour : geosoft.gxapi.GXIPJ
orientations_sont_les_mêmes ( (GXIPJ)arg1 ) &rarr int :¶

orientations_are_the_same_within_a_small_tolerance ( (GXIPJ)arg1 ) &rarr int :¶
Paramètres:argument1 ( geosoft.gxapi.GXIPJ ) – IPJ 2
Retour:0 - Non 1 - Oui
Type de retour :entier
  • argument1 (entier) – Constantes IPJ_TYPE
  • arg2 (str) – chaîne 1
  • arg3 (str) – chaîne 2
  • arg4 (str) – chaîne 3
  • argument1 ( geosoft.gxapi.GXIPJ ) – IPJ reprojeté en entrée (n'a pas besoin d'inclure une orientation). En sortie contient le même type d'orientation que l'IPJ initial, ajusté pour être au même endroit.
  • arg2 ( geosoft.gxapi.float_ref ) – X origine de la grille (entrée valeur initiale, sortie nouvelle valeur)
  • arg3 ( geosoft.gxapi.float_ref ) – Y origine de la grille (entrée valeur initiale, sortie nouvelle valeur)
  • arg4 ( geosoft.gxapi.float_ref ) – X taille de cellule de la grille (entrée valeur initiale, sortie nouvelle valeur)
  • arg5 ( geosoft.gxapi.float_ref ) – Taille de cellule Y de la grille (valeur initiale d'entrée, nouvelle valeur de sortie)
  • arg6 ( geosoft.gxapi.float_ref ) – Rotation de la grille (degrés CCW) (entrée valeur initiale, sortie nouvelle valeur)
  • argument1 (flotter) – Ligne 0 Élément 0
  • arg2 (flotter) – Ligne 0 Élément 1
  • arg3 (flotter) – Ligne 0 Élément 2
  • arg4 (flotter) – Ligne 0 Élément 3
  • arg5 (flotter) – Ligne 1 Élément 0
  • arg6 (flotter) – Ligne 1 Élément 1
  • arg7 (flotter) – Ligne 1 Élément 2
  • arg8 (flotter) – Ligne 1 Élément 3
  • arg9 (flotter) – Ligne 2 Élément 0
  • arg10 (flotter) – Ligne 2 Élément 1
  • arg11 (flotter) – Ligne 2 Élément 2
  • arg12 (flotter) – Ligne 2 Élément 3
  • arg13 (flotter) – Ligne 3 Élément 0
  • arg14 (flotter) – Ligne 3 Élément 1
  • arg15 (flotter) – Ligne 3 Élément 2
  • arg16 (flotter) – Ligne 3 Élément 3
  • argument1 (flotter) – X emplacement de l'origine de la vue
  • arg2 (flotter) – Y emplacement de l'origine de la vue
  • arg3 (flotter) – Emplacement Z de l'origine de la vue
  • arg4 (flotter) – Rotation en X
  • arg5 (flotter) – Rotation en Y
  • arg6 (flotter) – Rotation en Z
  • arg7 (flotter) – Mise à l'échelle en X
  • arg8 (flotter) – Mise à l'échelle en Y
  • arg9 (flotter) – Mise à l'échelle en Z
  • argument1 (flotter) – X emplacement de l'origine de la vue
  • arg2 (flotter) – Y emplacement de l'origine de la vue
  • arg3 (flotter) – Z emplacement de l'origine de la vue
  • arg4 (flotter) – Rotation en X
  • arg5 (flotter) – Rotation en Y
  • arg6 (flotter) – Rotation en Z
  • arg7 (flotter) – Mise à l'échelle en X
  • arg8 (flotter) – Mise à l'échelle en Y
  • arg9 (flotter) – Mise à l'échelle en Z
  • arg10 (entier) – Constantes IPJ_3D_ROTATE
  • arg11 (entier) – Constantes IPJ_3D_FLAG
  • argument1 (flotter) – X emplacement de l'origine de la vue
  • arg2 (flotter) – Y emplacement de l'origine de la vue
  • arg3 (flotter) – Emplacement Z de l'origine de la vue
  • arg4 (flotter) – Axe X Composant X
  • arg5 (flotter) – composante Y de l'axe X
  • arg6 (flotter) – Composante Z de l'axe X
  • arg7 (flotter) – Composant X de l'axe Y
  • arg8 (flotter) – Composante Y de l'axe Y
  • arg9 (flotter) – Composante Z de l'axe Y
  • arg10 (flotter) – Mise à l'échelle en X
  • arg11 (flotter) – Mise à l'échelle en Y
  • arg12 (flotter) – Mise à l'échelle en Z
  • argument1 ( geosoft.gxapi.GXVV ) – Emplacements de la section X (par exemple, distance le long de la courbe)
  • arg2 ( geosoft.gxapi.GXVV ) – Vrai X
  • arg3 ( geosoft.gxapi.GXVV ) – Vrai Y
  • arg4 (entier) – Utiliser l'axe Y logarithmique (généralement pour les profils de données) 0 : Non, 1 : Oui
Paramètres:argument1 (str) – ESRI chaîne de projection au format prj
Retour:Rien
Type de retour :Rien
  • argument1 (str) – “nom de projection” ou PCS_NAME de ipj_pcs.csv (référence / projection) ou numéro de code du système de coordonnées EPSG ou “<file.prj>” nom de fichier de projection ou “<file.wrp>” nom de fichier de distorsion
  • arg2 (str) – “nom du datum”[, grand axe, ellipticité, premier méridien] ou DATUM de datum.csv ou numéro de code de datum EPSG
  • arg3 (str) – “nom de la méthode”, paramètres (P1 à P8) ou “nom de la projection”[,”nom de la méthode”,”Unités”,P1,P2. ] ou TRANSFORM à partir du numéro de code de la méthode de transformation transform.csv ou EPSG
  • arg4 (str) – “nom de l'unité”, conversion en mètres ou UNIT_LENGTH à partir de units.csv
  • arg5 (str) – “nom de la transformation locale”[,dX,dY,dZ,rX,rY,rZ,Scale] ou DATUM_TRF de datumtrf.csv ou AREA_OF_USE de ldatum.csv ou numéro de code de transformation de datum local EPSG
  • argument1 (str) – “nom de projection” ou PCS_NAME de ipj_pcs.csv (référence / projection) ou numéro de code du système de coordonnées EPSG ou “<file.prj>” nom de fichier de projection ou “<file.wrp>” nom de fichier de distorsion
  • arg2 (str) – “nom du datum”[, grand axe, ellipticité, premier méridien] ou DATUM de datum.csv ou numéro de code de datum EPSG
  • arg3 (str) – “nom de la méthode”, paramètres (P1 à P8) ou “nom de la projection”[,”nom de la méthode”,”Unités”,P1,P2. ] ou TRANSFORM à partir du numéro de code de la méthode de transformation transform.csv ou EPSG
  • arg4 (str) – “nom de l'unité”, conversion en mètres ou UNIT_LENGTH à partir de units.csv
  • arg5 (str) – “nom de la transformation locale”[,dX,dY,dZ,rX,rY,rZ,Scale] ou DATUM_TRF de datumtrf.csv ou AREA_OF_USE de ldatum.csv ou numéro de code de transformation de datum local EPSG

0 - erreur lors de la définition d'IPJ, entrée IPJ inchangée. 1 - succès : IPJ défini à l'aide des valeurs d'entrée.

set_method_parm ( (int)arg1, (flotter)arg2 ) &rarr Aucun :¶


Coordonner les modes opératoires à mettre en œuvre

Il s'agit d'une tâche parapluie pour certaines méthodes d'opération de coordonnées qui ne sont pas encore prises en charge dans Apache SIS. Les opérations de coordination comprennent projections cartographiques (ex. Transverse Mercator, Lambert Conic Conformal, etc.), décalages de référence (ex. transformations de NAD27 à NAD83 aux États-Unis), transformation de coordonnées verticales, etc. Nous ne pouvons bien sûr pas lister toutes les formules possibles que nous ne prenons pas en charge, mais cette tâche JIRA liste au moins certaines des opérations répertoriées dans les notes d'orientation de l'EPSG.

Les principaux éléments de ce travail sont les notes d'orientation de l'EPSG, qui peuvent être téléchargées gratuitement sur le site suivant :

Publication IOGP 373-7-2 – Note d'orientation en géomatique numéro 7, partie 2
Coordonner les conversions et les transformations, y compris les formules
http://www.epsg.org/GuidanceNotes

Les étudiants de l'été de code Google intéressés par ce travail devraient être raisonnablement à l'aise avec le langage Java (mais pas nécessairement avec la bibliothèque JDK en général, car ce travail utilise relativement peu de classes JDK en dehors Math) et en mathématiques. En particulier, ce travail nécessite une bonne compréhension des transformations affines: leur représentation sous forme de matrice, et comment mapper un terme dans une formule à un coefficient dans la matrice de transformation affine.

Apache SIS possède une fonctionnalité avancée que l'on ne trouve pas facilement dans les logiciels ou manuels géospatiaux populaires : la capacité de calculer le dérivé (ou plus précisément, le Jacobien) d'une transformation en un point donné. La mise en œuvre de cette fonctionnalité nécessite la capacité de trouver la dérivée analytique d'une formule non linéaire et de la simplifier.


Coordonner les modes opératoires à mettre en œuvre

Il s'agit d'une tâche parapluie pour certaines méthodes d'opération de coordonnées qui ne sont pas encore prises en charge dans Apache SIS. Les opérations de coordination comprennent projections cartographiques (ex. Transverse Mercator, Lambert Conic Conformal, etc.), décalages de référence (par exemple transformations de NAD27 à NAD83 aux États-Unis), transformation de coordonnées verticales, etc. Nous ne pouvons bien sûr pas lister toutes les formules possibles que nous ne prenons pas en charge, mais cette tâche JIRA liste au moins certaines des opérations répertoriées dans les notes d'orientation de l'EPSG.

Les principaux éléments de ce travail sont les notes d'orientation de l'EPSG, qui peuvent être téléchargées gratuitement sur le site suivant :

Publication IOGP 373-7-2 – Note d'orientation en géomatique numéro 7, partie 2
Coordonner les conversions et les transformations, y compris les formules
http://www.epsg.org/GuidanceNotes

Les étudiants de l'été de code Google intéressés par ce travail devraient être raisonnablement à l'aise avec le langage Java (mais pas nécessairement avec la bibliothèque JDK en général, car ce travail utilise relativement peu de classes JDK en dehors Math) et en mathématiques. En particulier, ce travail nécessite une bonne compréhension des transformations affines: leur représentation sous forme de matrice, et comment mapper un terme dans une formule à un coefficient dans la matrice de transformation affine.

Apache SIS possède une fonctionnalité avancée que l'on ne trouve pas facilement dans les logiciels ou manuels géospatiaux populaires : la capacité de calculer le dérivé (ou plus précisément, le Jacobien) d'une transformation en un point donné. La mise en œuvre de cette fonctionnalité nécessite la capacité de trouver la dérivée analytique d'une formule non linéaire et de la simplifier.