Suite

Comment convertir un fichier ASCII en TIN où XY est YX dans ArcGIS ?

Comment convertir un fichier ASCII en TIN où XY est YX dans ArcGIS ?


J'ai des fichiers ASCII (vous pouvez les télécharger à partir de : http://analizasan.kvprojekty.pl/underwork/gisse/M3451Ca1.zip) que j'aimerais créer TIN.

J'ai d'abord utilisé l'outil "ASCII3DToFeatureClass" pour créer des fichiers de formes avec des points. Puis l'outil "CreateTin" pour créer un TIN. Tout va bien mais mes coordonnées XY sont inversées et mon TIN est ailleurs par rapport au réel.

Puis-je indiquer X et Y lorsque je crée un TIN avec l'outil "CreateTin" ? Je sais que je peux éditer le fichier .asc dans TXT Editor mais c'est beaucoup de travail.

Ces fichiers sont dans le système de coordonnées PUWG_1992 (Pologne) : Projection : Gauss_Kruger False_Easting : 500000,0 False_Northing : -5300000,0 Central_Meridian : 19,0 Scale_Factor : 0,9993 Latitude_Of_Origin : 0,0 Unité linéaire : Mètre (1,0)

Système de coordonnées géographiques : GCS_ETRS_1989 Unité angulaire : Degré (0,0174532925199433) Méridien principal : Greenwich (0,0) Référence : D_ETRS_1989 Sphéroïde : GRS_1980 Axe semi-majeur : 6378137,0 Axe semi-mineur : 6356752,314140356 Aplatissement inverse : 298 257222101

J'attends de l'aide.


Exportez votre shapefile quelque part. Supprimer le fichier prj. Ouvrez ce fichier dans une nouvelle session ArcGIS. Assurez-vous que le cadre de données n'a pas de projection définie. Exécutez ce script sur le champ Forme :

def SwapXY(shp): fp=shp.firstPoint point=arcpy.Point(fp.Y,fp.X,fp.Z) point de retour

Expression:

SwapXY(!Forme!)

Vous pouvez également exécuter Ajouter des attributs de géométrie (gestion des données) sur le fichier de formes. Exporter le tableau vers Excel et le ramener sous forme de tableau Excel XY


ESRI TIN (fichiers .adf) vers CAD (.dwg ou .dxf)

Ma rage a lentement augmenté toute la journée jusqu'à ce point. Quelqu'un sait-il comment obtenir un TIN (en l'état) à partir du format propriétaire ESRI et au format .dwg ou .dxf ?

J'ai fait la chose où vous le chargez dans ArcScene et l'exportez en tant que VRML, mais cela ne préserve pas les coordonnées du monde réel et le place à une origine X/Y arbitraire.

Utilisez l'outil Nœud TIN dans 3D Analyst pour générer un fichier de formes de points pour le TIN. Utilisez la classe d'entités Z à ASCII. Importez l'ASCII en tant que fichier de données XYZ et vous pouvez enregistrer la surface dans n'importe quel type de dessin pris en charge. Sachez simplement qu'un grand TIN générera une quantité extrême de points de données. J'ai trouvé que 26 millions est la limite d'importation pour Civil3D et est facilement réalisable si votre TIN est basé sur LIDAR.

Si vous disposez de Civil3D, vous pouvez lire en mode natif la sortie de l'outil Nœud TIN et créer une surface directement à partir des données.

C'est exactement ce que j'ai fait, mais j'ai utilisé Microstation. A très bien fonctionné !

Génial, je ne connais pas Civil3D mais j'y ai accès.

Voudriez-vous simplement utiliser des points d'insertion à partir du fichier et spécifier xyz ?

Puisque vous ne l'avez pas mentionné, j'exporterais votre VRML .. mais lorsque vous choisissez le nom du fichier. Cliquez sur "Options. " sous le champ et assurez-vous que 'Utiliser l'extension GeoVRML' est coché. Cela peut aider avec ce problème.

Malheureusement, je n'ai pas trouvé un excellent support avec les VRML et cela pourrait valoir la peine d'essayer de prendre votre TIN et d'exécuter Interpolate Polygon to Multipatch (qui accepte les TIN), puis de convertir à l'aide de Multipatch en Collada.

Je ne peux pas confirmer comment les informations spatiales sont transférées, mais j'avais l'habitude de l'utiliser. DAE peut être plus pris en charge que VRML. Je n'ai pas trouvé de support idéal pour les formes très avancées/complexes d'Arc, soz. J'espère que vous trouverez une bonne solution de toute façon!


1 réponse 1

Les éléments de cette notation de présentation sont

Cela n'a rien à voir avec les tuples.

Soit $G$ le groupe donné par cette présentation. Pour obtenir un isomorphisme $f : G o mathbb imes mathbb$ vous pourriez :

  • Choisissez deux éléments $u, v in mathbb imes mathbb$ pour être l'image de $x$ et $y$ respectivement, c'est-à-dire pour définir $f(x) = u$ et $f(y) = v$.
  • Vérifiez que toutes les relations sont vérifiées après avoir substitué $u$ pour $x$ et $v$ pour $y$ c'est-à-dire vérifier $f(xy) = f(yx)$, ou de manière équivalente, que $u + v = v + u$.
  • Trouver un homomorphisme $mathbb imes mathbb à G$ qui est l'inverse de $f$

La première puce est suffisante pour définir un homomorphisme du groupe libre sur deux éléments à $mathbb imes mathbb$. Le deuxième point montre que cet homomorphisme est bien défini sur le groupe $G$. Le troisième point montre qu'il s'agit d'un isomorphisme.

En fin de compte, il n'est pas difficile de deviner le choix de $u$ et $v$ la toute première chose qui vient à l'esprit en réponse à la question "Donnez-moi une paire d'éléments de $mathbb imes mathbb$ qui génèrent le groupe !" vont probablement fonctionner.


ArcGMT :Une suite d'outils de conversion entre Arc/INFO® et les outils de cartographie génériques (GMT)

[1] Département de géosciences, Oregon State University, Corvallis, OR 97331 États-Unis
[2] Institute of Geological & Nuclear Sciences Ltd., P.O. Box 30-368, Lower Hutt, NOUVELLE-ZÉLANDE
[3] Collège des sciences océaniques et atmosphériques, Oregon State University, Corvallis, OR 97331 États-Unis

Téléchargez le code source qui accompagne cet article (la version 2 a une conversion de fichier de contour !)

Utiliser ArcGIS 9.x/10.x ?
Assurez-vous que vous utilisez la dernière version de GMT et dans cette version, faites : grdreformat -fg oldformat.grd newformat.nc

ArcGMT est une suite d'outils de conversion Outils de cartographie génériques (GMT) grilles dans un format lisible par les plus anciens Arc/INFO® et ArcView® des progiciels de système d'information géographique (SIG), et inversement, pour la création de grilles et de fichiers GMT à partir des formats Arc/INFO. La GMT est devenue une norme non officielle parmi les géoscientifiques pour la production de cartes, en particulier à partir de données bathymétriques marines, d'imagerie sonar à balayage latéral, de gravité et de données magnétiques. De nombreux géoscientifiques se sont également tournés vers les capacités de cartographie, de gestion de bases de données et d'analyse spatiale d'Arc/INFO et d'ArcView. Par la suite, un besoin s'est fait jour, notamment dans la communauté de la géologie et de la géophysique marines, pour une interface de conversion entre ces deux milieux. Cette version initiale d'ArcGMT est conçue pour les postes de travail UNIX exécutant Arc/INFO version 7.0 ou supérieure, GMT version 3.0 ou supérieure et ArcView version 3.0 avec l'extension Spatial Analyst. ArcGMT se compose d'une série de scripts shell UNIX, de fichiers Arc Macro Language (AML) et de scripts Avenue, regroupés dans une interface de menu facile à utiliser pour une utilisation dans Arc/INFO ou comme extension d'ArcView.

Mots clés: conversion de données système d'information géographique (SIG) géologie et géophysique marines géographie physique

Note ajoutée en preuve, 10/9/00: La version actuelle de MBSystem inclut désormais "mbm_grd2arc" pour les utilisateurs d'ArcView, qui convertit les grilles GMT en grilles ASCII qu'Arc/INFO peut lire.

Ces dernières années, de nombreux géoscientifiques se sont tournés vers le système d'information géographique (SIG) pour la cartographie et la gestion avancées de leurs données, ainsi que pour l'intégration avec d'autres types de données environnementales (Wadge, 1992a Wadge, 1992b Walker et autres, 1996). Cela est particulièrement vrai dans la communauté de la géologie et de la géophysique marines (Bobbitt et autres, 1997 Hatcher, Maher et Orange, 1997 Fox, Bobbitt et Wright, 1996 Wright, 1996 Wright et Goodchild, sous presse) où les scientifiques s'appuient également fortement sur la cartographie générique. Paquet d'outils (GMT). GMT est devenu une norme non officielle pour la production de cartes de bathymétrie marine, d'imagerie sonar à balayage latéral, de gravité et de magnétisme. Sur le front des SIG, beaucoup se sont tournés vers Arc/INFO® et son petit cousin de bureau, ArcView®, tous deux de l'Environmental Systems Research Institute. Les deux peuvent maintenant être considérés comme les principaux progiciels SIG vectoriels au monde (GIS World, 1996) et ont été intégralement adoptés par des agences telles que l'USGS et la NOAA, ainsi que par des laboratoires universitaires et des salles de classe dans tout l'hémisphère occidental. Un besoin s'est fait, encore une fois, en particulier dans la communauté de la géologie et de la géophysique marines, d'une interface entre ces progiciels très utilisés. GMT est souvent utilisé pour traiter et créer des cartes de base, mais fournit peu de moyens pour l'incorporation d'attributs détaillés aux emplacements qui sont cartographiés. Arc/INFO et ArcView, en tant que packages SIG, offrent une combinaison de capacités de cartographie et de bases de données relationnelles qui dépassent l'heure GMT. Cependant, les cartes de base qui sont plus facilement créées en GMT (de concert avec le package MB-System de l'observatoire terrestre de Lamont-Doherty, qui traite la bathymétrie multifaisceaux brute et le balayage latéral haute résolution) fournissent une couche de base cruciale pour de nombreuses applications qui sont souhaitables pour le SIG. Un exemple est la cartographie à petite échelle des emplacements des cheminées hydrothermales, des communautés animales et des coulées de lave dans les centres d'épandage des fonds marins, une tâche qui n'est

5% achevé à l'échelle mondiale (par exemple, Embley et al., 1995 Haymon et al., 1991). Récemment, avec la popularité croissante des SIG, des environnements d'acquisition de données à bord des navires ont été mis en place pour importer des données de submersibles ou de véhicules télécommandés directement dans un SIG (Bobbitt, 1996 Bobbitt, Fox et Wright, 1996 Wright, 1994 Wright, 1996), certains en temps réel (Hatcher, Maher et Orange, 1997). Cela a grandement facilité l'intégration des observations géologiques avec la chimie et la biologie des évents. Cependant, les images bathymétriques ou sonar à balayage latéral, qui permettent l'interprétation de structures géologiques plus importantes qui placent ces observations à petite échelle dans un contexte régional, sont le plus souvent collectées, traitées et archivées dans des grilles GMT. Cela se produit généralement lors d'expéditions distinctes précédant les plongées du submersible ou du véhicule. Il est crucial de pouvoir ajouter ces données en tant que thèmes ou couches de base aux bases de données SIG existantes. Inversement, il peut être nécessaire d'exporter des entités cartographiées dans le SIG pour les inclure dans des cartes et des tracés GMT élaborés.

Nous présentons ici ArcGMT, une suite d'outils pour convertir les grilles GMT dans un format lisible par Arc/INFO et ArcView. ArcGMT créera également des grilles et des fichiers GMT au format Arc/INFO. Nous commençons par fournir de brefs aperçus de GMT, Arc/INFO et ArcView, suivis d'une description fonctionnelle d'ArcGMT.

GMT est un progiciel du domaine public pour la manipulation de données et la génération de cartes et d'illustrations scientifiques de haute qualité. La version initiale est sortie en octobre 1991 (Wessel et Smith, 1991) et rééditée dans une version 3 très améliorée en 1995 (Wessel et Smith, 1995). Une enquête électronique (Wessel et Smith, données non publiées, 1994), ainsi qu'une évaluation du trafic ftp et des formulaires d'inscription soumis, révèlent que plus de 5 000 utilisateurs dans le monde utilisent fréquemment le GMT dans leur travail. Au sein de la communauté de la géologie et de la géophysique marines, le GMT est actuellement utilisé par des scientifiques de tous les continents, à bord de navires de recherche de la flotte du University National Oceanographic Laboratory System, et même à bord de sous-marins et d'avions de la marine américaine en mission scientifique.

GMT fonctionne sur presque toutes les plates-formes informatiques capables d'exécuter UNIX&tm, des superordinateurs Cray aux stations de travail en passant par les PC et les Macintosh. Bien qu'il soit capable de filtrer et de quadriller des ensembles de données à 1 et 2 dimensions, GMT est surtout connu et utilisé pour ses capacités avancées de traçage et de cartographie (Fig. 1)

Figure 1. Exemple des capacités cartographiques avancées des Generic Mapping Tools (GMT). Figure composite montrant un maillage 3-D de la topographie hawaïenne à partir de la bathymétrie ETOPO 5 au-dessus de la carte des contours du géoïde hawaïen. Après Wessel et Smith (1991).

Cliquez sur l'image pour une vue agrandie

Avec des ventes estimées à plus de 200 millions de dollars (Daratech Inc., 1994), l'Environmental Systems Research Institute est la plus grande entreprise de SIG au monde. Il fournit des logiciels autonomes dans des catégories de produits allant du SIG de bureau au SIG analytique avancé haut de gamme aux kits d'outils de développement de logiciels. Le produit phare de la société, ARC/INFO, est le SIG le plus couramment utilisé dans une variété de segments de marché et de recherche scientifique, notamment le pétrole et l'exploitation minière, la cartographie des eaux souterraines, la cartographie géologique et la géologie marine. De plus, ArcView, le nouveau SIG de bureau de l'entreprise, a connu un énorme succès et fournit un complément plus facile à utiliser et rentable à ARC/INFO. En plus des agences susmentionnées, Arc/INFO et ArcView ont été adoptés par 40 des 50 plus grandes compagnies pétrolières du monde (Petroleum Intelligence Weekly, 1995), sont des sites autorisés sur plus de 100 établissements d'enseignement supérieur aux États-Unis (M. Phoenix, ESRI, communication personnelle, 1997), et est utilisé par des sociétés telles que Sea Beam Instruments, le leader mondial des systèmes de sonar pour la cartographie des fonds marins.

  1. fourniture de liens au sein d'un relationnel environnement de gestion de base de données entre les emplacements des points, des lignes, des zones et/ou des grilles et des informations descriptives détaillées les concernant (attributs), ainsi que les métadonnées associées.
  2. fourniture d'algorithmes pour l'analyse spatiale des données (par exemple, autocorrélation spatiale, modélisation d'advection et de dispersion, reconnaissance et analyse de formes, mise en mémoire tampon, interpolation spatiale, superposition, etc.).
  3. « intelligence spatiale » ou la capacité non seulement de modifier les projections cartographiques, mais également de calculer les distances, d'identifier les voisins les plus proches, de calculer les régions de proximité autour des entités existantes et de mettre à jour leurs attributs en conséquence.

DESCRIPTION FONCTIONNELLE DE L'ARCGMT

ArcGMT permet aux utilisateurs de tirer parti des capacités de GMT et d'Arc/INFO ou d'ArcView en fournissant un moyen de conversion de données entre les deux environnements. Cette version initiale convient uniquement aux postes de travail UNIX qui exécutent également Arc/INFO version 7.0 ou supérieure, GMT version 3.0 ou supérieure et ArcView version 3.0 avec l'extension Spatial Analyst. ArcGMT se compose d'une série de scripts shell UNIX, de fichiers Arc Macro Language (AML) et de scripts Avenue. AML et Avenue sont les langages utilisés dans Arc/INFO et ArcView, respectivement, pour programmer et personnaliser ces environnements. Il existe deux options pour l'utilisation du code, selon le progiciel SIG qui sera principalement utilisé, Arc/INFO ou ArcView.

ArcGMT avec Arc/INFO

Dans l'environnement Arc/INFO, ArcGMT est exploité via un simple menu déroulant à partir duquel l'utilisateur peut sélectionner la documentation à lire ou choisir parmi trois types de conversion :

ArcGMT avec ArcView

ArcView a été introduit en 1994 en tant que SIG de bureau à échelle réduite, avec une interface utilisateur graphique non disponible dans Arc/INFO, permettant ainsi à l'utilisateur de visualiser, d'explorer et d'interroger des données beaucoup plus facilement. Avec une courbe d'apprentissage de 1 à 5 jours (contre 6 à 12 mois pour Arc/INFO), de nombreux géoscientifiques se sont tournés vers ArcView comme principale introduction au monde des SIG. Les récentes révisions d'ArcView ont inclus des capacités d'analyse spatiale plus puissantes, bien qu'Arc/INFO reste la boîte à outils SIG standard et à service complet. ArcGMT permet aux utilisateurs d'effectuer les mêmes conversions décrites ci-dessus dans l'environnement ArcView. Ceci est facilité par la mise en œuvre d'appels de procédure à distance (RPC) entre Arc/INFO et ArcView. Il existe plusieurs formes de RPC, mais ArcView et ARC/INFO utilisent la norme Open Network Computing. Les RPC permettent aux applications clientes d'initier des procédures dans une application serveur sur un réseau. Les serveurs UNIX RPC sont identifiés de manière unique par un nom d'hôte, un identifiant de serveur et un numéro de version. Une application cliente RPC utilise ces informations pour établir la connexion RPC. Dans ArcGMT, ArcView envoie des requêtes client à un serveur Arc/INFO sous la forme d'instructions Avenue. Ces instructions appellent les AML de conversion ArcGMT et les scripts shell, puis permettent à l'utilisateur d'importer les résultats dans l'environnement ArcView. Avant d'utiliser ArcGMT dans ArcView, l'utilisateur doit démarrer Arc/INFO en tant que serveur avec un simple "&type [iacopen]" directive, puis démarrez ArcView.

    Des boutons pour effectuer les trois conversions décrites ci-dessus, ainsi que pour découper des sous-ensembles de grilles pour créer de nouvelles grilles et créer des thèmes ombrés à partir de grilles (Fig. 2A).

ArcGMT est également fourni avec un fichier de projet ArcGMT en option ("arcgmt.apr") qui ouvre automatiquement ArcView avec l'extension déjà chargée.

RÉSULTATS ET CONCLUSION

La figure 3 montre le résultat d'une conversion de grille GMT vers Arc/INFO dans ArcGMT. La grille originale était de la bathymétrie Sea Beam 2000 révélant un mont sous-marin dans le sud-ouest du Pacifique découvert en mai 1996 (Wright, Bloomer et Boomerang Leg 8 Shipboard Scientific Party, 1996). Les cellules de la grille n'étaient pas carrées en raison de l'algorithme de grille utilisé lors de la création de la grille avec MB-System, mais ArcGMT était toujours capable d'effectuer une conversion parfaite. La figure 4 montre une carte GMT de la bathymétrie Sea Beam recouverte de linéations de failles qui ont été exportées à partir d'une couverture de ligne Arc/INFO avec ArcGMT. La figure 5 montre les résultats d'une conversion de grille Arc/INFO en GMT dans ArcGMT. La grille d'origine était un modèle d'élévation numérique (DEM) USGS de 1 degré (espacement des données de 3 x 3 secondes d'arc), 1:250 000 de la grande île d'Hawaï téléchargé depuis le centre de données USGS Eros.

Figure 3. Résultats d'une conversion de grille GMT en Arc/INFO dans ArcGMT. (Haut) Carte de contour GMT d'une grille GMT de bathymétrie Sea Beam 2000 du sud-ouest du Pacifique, entrée dans ArcGMT pour conversion. Les données proviennent de Wright, Bloomer et Boomerang Leg 8 Shipboard Scientific Party (1996). Les étoiles sur la carte montrent où les roches ont été draguées au sommet de ce mont sous-marin. (Bas) Capture d'écran d'une session ArcGMT, montrant la même grille après conversion au format Arc/INFO. La grille est affichée en niveaux de gris par défaut dans le module GRID d'Arc/INFO.

Cliquez sur l'image pour une vue agrandie

Figure 4. Carte de contour GMT d'une grille GMT de bathymétrie Sea Beam du Pacifique NE, recouverte de linéations de failles qui ont été exportées à partir d'une couverture de ligne Arc/INFO par ArcGMT. Les données proviennent du site Web de RIDGE Multibeam Synthesis.

Cliquez sur l'image pour une vue agrandie

Figure 5. Résultats d'une conversion de grille Arc/INFO en GMT dans ArcGMT. (Haut) Un DEM USGS de 1 degré, 1:250 000 au format Arc/INFO, affiché en niveaux de gris par défaut dans le module GRID d'Arc/INFO. Les données proviennent du site Web de l'USGS Eros Data Center.
(Milieu) Construction TIN de la grille à l'aide de la commande LATTICETIN avec ombrage analytique dans Arc/INFO avant conversion finale au format GMT. Cette option peut être utilisée dans ArcGMT si vous le souhaitez. Ce n'est pas la valeur par défaut.
(Bas) Carte de contour créée en GMT de la grille GMT résultante (contours de 100 m, annotés tous les 500 m).

Cliquez sur l'image pour une vue agrandie

Quelques mises en garde s'appliquent aux conversions de grille Arc/INFO vers GMT dans ArcGMT. La conversion est conçue principalement pour les données d'élévation ou de bathymétrie en unités de m positif. Si l'utilisateur essaie de convertir un autre type de grille (par exemple, la gravité en mgals), la valeur z peut être décalée de 1 à 2 ordres de grandeur, ce qui peut facilement être résolu avec "grdmath" en GMT. Les DEM USGS peuvent être exprimés en unités x-y de secondes décimales, qui sont des unités inacceptables pour GMT. Avant de convertir, l'utilisateur doit toujours revérifier les caractéristiques de la grille avec la commande DESCRIBE dans Arc/INFO. Les limites de la grille produites par DESCRIBE doivent être notées par l'utilisateur. Il peut être nécessaire de projeter les unités de la grille en degrés décimaux et/ou de faire une transformation de datum. Dans le cas de l'USGS DEM d'Hawaï, le datum d'origine WGS72 a été modifié pour correspondre à une base de données stockée dans NAD83 avec les sous-commandes suivantes dans Arc/INFO PROJECT : sortie, projection géographique, unités dd, datum nar_d trois, paramètres, fin. Une telle transformation à trois paramètres peut entraîner des erreurs ou des incertitudes de l'ordre de

Les améliorations futures d'ArcGMT incluent des scripts Avenue supplémentaires pour lire automatiquement le nom d'hôte du serveur Arc/INFO pour l'extension ArcView "arcgmt.avx" et permettre à l'utilisateur de tester de manière interactive les projections, les transformations de datum et divers intervalles de maillage avant de convertir les grilles d'Arc en GMT.

Le code source et un fichier README pour ArcGMT sont disponibles ici (clic droit et "Enregistrer sous").

Remerciements -- Ce travail est soutenu par la National Science Foundation sous la subvention OCE-9521039. Le document a été amélioré par les critiques utiles de deux arbitres anonymes.

Bobbitt, A. M. (1996) Exploration de l'observatoire RIDGE et intégration de données à l'aide d'ALVIN et de systèmes d'information géographique. Événements RIDGE 7(1), 10-12.

Bobbitt, A. M., Dziak, R. P., Stafford, K. M. et Fox, C. G. (1997) Analyse SIG des données océanographiques de télédétection et d'observation sur le terrain. Géodésie marine 20(2-3), 153-161.

Bobbitt, A. M., Fox, C. G. et Wright, D. J. (1996) ArcView prend en charge la recherche en haute mer avec des sous-marins célèbres. Actes de la 16e conférence annuelle des utilisateurs d'ESRI, Palm Springs, Californie, document 082.

Daratech Inc. (1994) Nouvelle étude de marché SIG. GIM : International Journal for Surveying, Mapping, and Applied GIS 8(5), 36-37.

Embley, R. W., Chadwick, W. W., Jr., Jonasson, I. R., Butterfield, D. A. et Baker, E. T. (1995) Résultats initiaux de la réponse rapide à l'événement CoAxial de 1993 : relations entre les processus hydrothermaux et volcaniques. Lettres de recherche géophysique 22(2), 143-146.

GIS World (1996) 1996 International GIS Sourcebook, GIS World, Inc., Fort Collins, Colorado, 771 pages.

Hatcher, G. A., Maher, N. M. et Orange, D. L. (1997) La personnalisation d'ArcView comme outil en temps réel pour la recherche océanographique. Actes de la 17e conférence annuelle des utilisateurs d'ESRI, San Diego, Californie, article 676.

Haymon, RM, Fornari, DJ, Edwards, MH, Carbotte, S., Wright, D. et Macdonald, KC (1991) Distribution des évents hydrothermaux le long de la crête est du Pacifique (9&# 176 09'-54'N) et ses relation avec les processus magmatiques et tectoniques sur les dorsales médio-océaniques à propagation rapide. Earth and Planetary Science Letters 104, 513-534.

Fox, C. G., Bobbitt, A. M. et Wright, D. J. (1996) Integration and distribution of deepsea oceanographic data from the NE Pacific using Arc/Info et ArcView. Actes de la 16e conférence annuelle des utilisateurs d'ESRI, Palm Springs, Californie, article 293.

Petroleum Intelligence Weekly (1995) Top 50 des compagnies pétrolières. Petroleum Intelligence Weekly, décembre 1995.

Wadge, G. (1992a) Applications géologiques des SIG. Journal de la Société géologique, Londres 149, 672.

Wadge, G. (1992b) SIG pour la géologie. Terra Nova 3(1), 93-98.

Walker, J. D., Black, R. A., Linn, J. K., Thomas, A. J., Wiseman, R. et D'Attilio, M. G. (1996) Développement de systèmes d'information géographique : bases de données orientées pour des applications géologiques et géophysiques intégrées. GSA Aujourd'hui 6(3), 1-7.

Wessel, P. et Smith, W. H. F. (1991) Un logiciel gratuit permet de cartographier et d'afficher les données. Eos, Transactions, Union géophysique américaine 72, 441.

Wessel, P. et Smith, W. H. F. (1995) Publication d'une nouvelle version des outils de cartographie génériques. Eos, Transactions, American Geophysical Union 76(33), 329.

Wright, D. J. (1994) From Pattern to Process on the Deep Ocean Floor: A Geographic Information System Approach. doctorat thèse, Université de Californie, Santa Barbara, Californie, 341 pp.

Wright, D. J. (1996) Rumblings on the ocean floor : le SIG soutient la recherche en haute mer. Systèmes d'informations géographiques 6 (1), 22-29.

Wright, D. J., Bloomer, S. H. et Boomerang Leg 8 Shipboard Scientific Party (1996) New Sea Beam 2000 bathymetry and side scan data from the Tonga Forearc and Trench: Results from Boomerang Leg 8, an ODP site survey. Eos, Transactions de l'American Geophysical Union 77(46), F326.

Wright, D. J. et Goodchild, M. F. (1997) Data from the deep: Implications for the GIS community. Journal international des systèmes d'information géographique 11(6), 523-528.


Travail (résolu)

On travaille dans un cadrage orthonormé $< E_i >$ sur le fibré tangent, à base duale $< omega^i >$ . On définit la forme unique $omega^j_i$ par $omega^j_i (X) := langle abla_X E_i, E_j angle$ . Ensuite, en travaillant sur la première équation, et en laissant $langle X,Y angle := g(X,Y),$ on obtient

$commencer langle abla_XY, Z angle &= ldots = Z^j(X^iE_i(Y^j) + Y^iomega^_(X)), langle Y, abla_X Z angle &= ldots = Y^j(X^iE_i(Z^j) + Z^i omega^j_i (X)), X langle Y,Z angle &= ldots = X^i Y^j E_i (Z^j) + X^i Z^j E_i(Y^j), implique 0 &= (Z^jY^i omega^j_i + Y^jZ^iomega^j_i)(X), forall X,Y,Z in Gamma(TM) implies 0 &= (Z^jY^i)( omega^j_i + omega^i_j)(X), forall X,Y,Z in Gamma(TM) implies 0 &= omega^j_i + omega^i_j. finir$

Cette partie est bien. Maintenant, en définissant $[E_i,E_j]=c_^k E_k$ et en insérant les coordonnées dans la condition de torsion, j'obtiens $ Gamma^k_ - Gamma^k_ - c^k_ =0, $ où $Gamma^i_E_i:= abla_E_k$ . Maintenant clairement $Gamma^_ = omega^k_j(E_i)$ et $-c^_=domega^k(E_i,E_j)$ , nous avons donc $ egin omega^k_j(E_i) - omega^k_i (E_j) + domega^k(E_i,E_j) &=0, forall i,j. finir$ Cependant, $ (omega^p wedge omega^_p)(E_i,E_j) = egin omega^p(E_i) & omega^p(E_j) omega^k_p(E_i) & omega^k_p(E_j) end = ldots = omega^k_i(E_j) - omega^k_j(E_i).$ En mettant tout cela ensemble, nous obtenons $ egin domega^k(E_i,E_j) &= omega^p wedge omega^k_p (E_i,E_j), forall i,j implies domega^k &= omega^p coin omega^k_p. finir $

Notre système différentiel correspondant est donc

$ début omega^j_i + omega^i_j &= 0, domega^i - omega^j wedge omega^i_j &= 0.end $


Maintenant, $cos hetasin heta(cos heta-2sin heta)$ est limité entre, disons, $A$ et $B$.

Presser avec $r o0$ donne $f(r, heta) o0$ comme requis.

$0leqslant frac<|xy(x-2y)|>leqslant frac<(x^2+y^2)|(x-2y)|><2(x^2+y^2)>=frac<|x-2y|> <2> o 0 mbox < as >(x, y) o (0,0).$ Si $f(0,0)=0$, alors $f$ est continu à $(0,0)$. Sinon, $f$ n'est pas continu à $(0,0)$.

La fonction $f(x,y)$ est indéfinie à l'origine, elle n'est donc pas continue à $(0,0)$.

Notez qu'il n'a pas de limite comme $(x,y) o (0,0)$ car la limite de $f$ est différente selon le chemin que $(x,y) o (0,0) $. Par exemple,


Utilisez la boîte de dialogue Options d'importation de fichier texte ASCII générique pour spécifier comment formater le(s) fichier(s) texte sélectionné(s) pour l'importation. La boîte de dialogue apparaît une fois que vous avez sélectionné un fichier ASCII à charger.

  1. Sélectionnez le Type d'importation, pour spécifier comment traiter les données importées :
    • Point uniquement (toutes les entités sont des points): les lignes avec des données de coordonnées donneront lieu à une entité ponctuelle unique.
    • Entités ponctuelles, linéaires et surfaciques (par défaut) : des chaînes de deux lignes consécutives ou plus avec des données de coordonnées se traduiront par une entité linéaire ou surfacique. Toutes les lignes de coordonnées isolées se traduiront par une entité ponctuelle.
    • Grille d'élévation à partir de données de point 3D: Des lignes avec des données de coordonnées 3D seront utilisées pour générer un terrain triangulé, qui est quadrillé pour créer une grille d'élévation. La grille possède toutes les capacités d'un MNT importé, y compris la génération de contours, l'analyse de la ligne de visée et du champ de vision, et le drapage de trame. Lorsque les paramètres d'importation de fichiers ASCII sont terminés, la boîte de dialogue Options de création de grille d'élévation s'affiche pour configurer le processus de quadrillage.
  2. Sélectionner un Ordre/format des colonnes de coordonnées pour spécifier l'ordre des lignes de coordonnées dans le fichier d'import.
    • Les coordonnées peuvent être x suivi de y (longitude puis latitude), ou l'inverse. La valeur par défaut est Coordonnée X/Est/Longitude en premier. Les valeurs d'altitude, si elles sont présentes, sont toujours supposées venir après les valeurs x et y. Vous pouvez également utiliser des coordonnées au format WKT ou MGRS/USNG.
    • Utilisation Champs à ignorer en début de ligne pour contrôler dans quelle colonne les coordonnées commencent. Par exemple, si les coordonnées x et y se trouvent dans les colonnes 3 et 4, définissez cette valeur sur 2 (ignorez les deux premières colonnes) pour vous assurer que les coordonnées sont lues dans la bonne colonne.
    • Sélectionnez le Format des coordonnées à utiliser dans la liste déroulante.
  3. Sélectionnez le Préfixe de ligne de coordonnées pour spécifier si les coordonnées commencent au début de la ligne (Rien, la valeur par défaut) ou si les lignes de coordonnées commencent par une autre séquence de caractères. Par exemple, certains formats peuvent commencer les lignes de coordonnées par la séquence « XY ».
  4. Spécifie le Lignes à ignorer au début du fichier avant d'extraire les données. Ceci est utile s'il y a des lignes d'en-tête au début du fichier que vous souhaitez ignorer.
  5. Sélectionner un Délimiteur de coordonnées pour spécifier ce qui sépare les coordonnées sur les lignes de coordonnées. Détection automatique est la méthode par défaut et recommandée.
  6. Sélectionnez le Classification des caractéristiques dans les listes déroulantes pour affecter l'entité aux entités de surface, de ligne et de point importées à partir du fichier.
  7. Activez ou désactivez les cases à cocher suivantes selon vos besoins. Les paramètres par défaut dépendent du Type d'importation vous avez sélectionné :
    • Inclure les attributs des lignes avec l'option de données de coordonnées: lorsqu'elle est activée, les lignes contenant du texte après les données de coordonnées seront incluses en tant qu'attributs pour l'entité de cette coordonnée. Si cette option est désactivée, seules les lignes qui ne contiennent pas de données de coordonnées seront utilisées comme attributs.
    • En-têtes de colonne dans la première ligne du fichier: Si activé et le Type d'importation est Point seulement, les valeurs de la première ligne du fichier seront utilisées aux noms d'attributs pour les attributs trouvés dans les lignes de données de coordonnées. Ceci est utile pour les fichiers CSV.
    • Traiter la 3e valeur de coordonnée comme une élévation: lorsqu'elle est activée, et si une valeur numérique est trouvée immédiatement après les valeurs de coordonnées x et y (ou lat et lon), cette valeur sera traitée comme une élévation. Sinon, si Inclure les attributs des lignes avec l'option de données de coordonnées est activé, la valeur sera incluse en tant qu'attribut. Laissez cette option activée, sauf si vous importez des données de point où la 3e colonne est un attribut qui peut contenir toutes les valeurs numériques ainsi que les noms.
    • Si le fichier contient des données de ligne, des données de zone ou les deux qui sont délimitées par un changement dans un champ ou une colonne de données particulier (au lieu d'avoir des lignes sans coordonnées les séparant), spécifiez le Fonctionnalités de ligne de rupture/surface lors d'un changement de champ pour définir le champ (utiliser un index basé sur 1) pour vérifier la division des données en entités de ligne/surface distinctes.
  8. Clique le Sélectionnez le décalage/l'échelle de coordonnées pour ouvrir la boîte de dialogue Sélectionner l'échelle/le décalage de coordonnées. Utilisez cette boîte de dialogue pour définir un décalage et un facteur d'échelle à appliquer à chaque coordonnée. Le décalage que vous entrez dans la boîte de dialogue sera d'abord ajouté à chaque coordonnée, puis chaque coordonnée sera multipliée par le facteur d'échelle.

Abstrait

La capacité d'extraire des propriétés des modèles de cadre géologique 3D pour une utilisation dans la construction de modèles conceptuels et mathématiques est considérée comme de plus en plus importante, cependant, des outils et des techniques sont nécessaires pour prendre en charge ces flux d'informations. Le développement de telles méthodologies maximisera les opportunités d'utilisation et de réutilisation des informations, ce qui est particulièrement important car la valeur réelle de ces actifs n'est pas toujours connue lors de leur première acquisition. Cet article décrit brièvement les défis culturels et techniques associés à l'application des informations dérivées des modèles de cadre géologique 3D par des modèles de processus hydrogéologiques. Nous examinons comment ces problèmes sont traités et présentons un outil, SurfGrid, qui permet à un utilisateur de générer des grilles 3D (voxels) de données paramétrées à partir d'une série de surfaces géologiques.

Les procédures et outils décrits offrent la possibilité de réutiliser des actifs créés de manière coûteuse en fournissant des techniques conviviales qui permettent aux scientifiques multidisciplinaires d'extrapoler les distributions de propriétés à partir de modèles géologiques.


Modèle spatial et analyse dans ArcGIS 10.0+

Voici un didacticiel simple dans ArcGIS sur l'utilisation des outils de statistiques spatiales. Pas grand chose de nouveau ici, mais je l'ai trouvé utile et intéressant car il utilise les outils pour aider à éclairer l'analyse avec d'autres outils de l'ensemble d'outils.

  • Aller sur le site du blog :
    • Suivez le lien à http://www.arcgis.com/home/item.html?id=020e1d9e9d2d44f9ae7da8663cf3bac5
    • Blog de statistiques spatiales de Google arcgis->Tutorial/Training->Analyse des modèles spatiaux dans ArcGIS 10

    Notre zone d'analyse est la réserve naturelle de Morasko Meteorite située dans la voïvodie de Grande-Pologne. It was established in 1976 in order to protect the area of impact craters, which, according to researchers, were formed in the fall of the Morasko meteorite about 5,000 years ago. In addition, the oak-hornbeam forest with rare species of plants (lilium martagon, ceratophyllum submersum) and birds (european nightjar, black woodpecker) is protected.

    Vector data

    The centroid (geometric center) of the Morasko Meteorite nature reserve has X = 16.895 and Y = 52.487 coordinates in World Geodetic System 1984 (EPSG 4326). Let’s start by creating this point with the sf package.

    Now our point is embedded in space (has a spatial reference). In the next step, let’s create an approximate zone that will include the area of the reserve. The function st_buffer() will be used for this. Before this operation, we need to transform the coordinate system to a system with metric units, e.g. Poland CS92 (EPSG 2180), using st_transform() function.

    We have created a buffer with a radius of 400 meters. Let’s visualize it.

    Of course, the area shown above is not exactly the reserve area. The exact area can be determined from the polygon layer as in orthophotomap example using the General Geographic Database.

    Raster data

    Now we can search for available elevation data for this area using the DEM_request() function (it is analogous to the ortho_request() function). The only argument of the function is our reserve buffer.

    Let’s check the obtained results.

    We have received metadata with many types of data of different formats, timeliness, resolution, and accuracy. For our analysis, we need digital terrain model (DTM) and digital surface model (DSM) in the “ARC/INFO ASCII GRID” format. Let’s make data selection by creating two tables and combining them together.

    Now we can download the data using the tile_download() function with our filtered data frame as input.

    If you run into any problem with the download, remember that you can pass another download method from download.file() as a function argument.

    Traitement

    Let’s load the downloaded numerical models using the read_stars() function from the stars package, which allows working on spatiotemporal arrays. We have two files, one represents DTM and second represents DSM.

    You probably noticed the four-fold difference in their sizes. It is due to the difference between their cells resolutions. We need to unify them to a common resolution to be able to combine them into one stack. It is much better to use a lower resolution than to increase it, because we cannot get more information and the processing will be faster. Let’s use the st_warp() function to do this.

    Now, both models have the same dimensions (the number of rows and columns) and spatial resolution. Thus, we can combine them into one object ( DEM ).

    Now we have a DEM object that consists of two attributes (DTM and DSM). In fact, both attributes contains same type of data as they are representing elevation. Therefore, we can collapse the attributes into a new dimension. Let’s do that using st_redimension() .

    After this operation, our elevation attribute consists of the DTM and DSM layers (dimensions). Then let’s crop the rasters to our buffer.

    Let’s check what the result looks like.

    In the first quadrant of the circle, we can see five smaller circles. These are the craters formed after the impact of the Morasko meteorite. The largest fragment found weighs 272 kg and it is the largest meteorite found in Poland. The collection of found meteorites can be seen at the Earth Museum in Poznań.

    Let’s calculate the crater width using the terrain transverse profile. We can use our centroid and add a second example point 30 degrees towards N. Next, we connect these points into a line ( st_linestring() ) and then sample this line every 1 m ( st_line_sample() ), because our DEM has this resolution. As a result, we get one complex geometry (MULTIPOINT), which we have to convert into a simple geometry (INDIQUER) consisting of many points. The function st_cast() is used for this.

    In the last step, we extract the elevation values for these points using st_extract() .

    Now we can see how our transverse profile looks like.

    The largest width of the impact crater is about 90 m.

    Okay, we checked the terrain. In the last step, let’s examine the height of the objects on it. For this purpose, we calculate the height of the trees by subtracting the DTM from the DSM. The product of this difference is called normalized DSM, because it takes the terrain elevation as a reference.


    Voir la vidéo: Reclassify u0026 Convert Raster to Shapefile Polygon u0026 Calculate Area Using ArcGIS