Ich habe eine Postgis-Tabelle mit Polygonen. Ich möchte meine Polygone in 4-Punkte-Polygone umwandeln: Den kleinsten MBR abrufen. Hier kommt eine Illustration von dem, was ich will: Mein ursprüngliches Polygon ist in rosa Farbe. Ich denke, es wird ziemlich einfach sein, den MBR zu bekommen (mit...
Ich habe gepostet Gibt es ein Projektionsproblem mit den Hydrologiedaten von Natural Resources Canada (NRC) Canada Atlas? Beschreibung einer Verschiebung einiger NRC-Daten (nicht aller). Anscheinend hat diese Verschiebung mit der Umrechnung zwischen geografischen Koordinaten, Breiten- und...
Laut Dokumentation stehen folgende Optionen zur Verfügung: Syntax ExportToTIFF (map_document, out_tiff, {data_frame}, {df_export_width}, {df_export_height}, {Auflösung}, {world_file}, {color_mode}, {tiff_compression}, {geoTIFF_tags}) Dies funktioniert jedoch: import arcpy.mapping mxd =...
Was ist der einfachste Weg, um Openstreetmap-Daten zu drehen und zu rendern und trotzdem ein schönes Rendering zu erhalten? Zum Beispiel würde Toronto im Querformat am besten leicht gedreht auf eine Seite passen, sodass die nördliche Grenze der Stadt parallel zum oberen Rand der gedruckten Seite...
Ich versuche, die Bertin 1953-Kartenprojektion zu verwenden, wie im Buch "Semiology of Graphics" von Jacques Bertin beschrieben. Ich kann die Einstellungen jedoch in keiner meiner GIS-Software finden. Hat jemand eine Lösung? Vielen Dank,
Wenn ich nach Bildern der Erde suche, die sich in Längen- und Breitengrade teilen, scheinen das zwei Haupttypen zu sein: ist ein perfektes Quadrat und die Breiten darin wachsen unterschiedlich schnell vom Äquator nach oben (und unten) [Mittelwert scheint, dass die Länge von 0 bis 66 kürzer ist als...
Ich versuche, den Begrenzungsrahmen des Kartenbereichs von OpenLayers anzuzeigen. Ich habe versucht, map.getExtent () zu verwenden, da es genau das sein sollte, was ich brauche, aber die zurückgegebenen Werte sind "seltsam": Wenn ich zum Beispiel die Karte mit Zoomstufe 14 um Grenoble (Frankreich,...
Ich suche Rat bei der Arbeit mit Weltdatensätzen, bei denen das GCS einen positiven Bereich aufweist (0 bis 360 Grad). Die Daten, mit denen ich arbeite, sind ozeanografische NetCDF- Daten und haben wie erwähnt positive Koordinatenwerte. Wenn es in einem normalen GCS WGS84 in ArcGIS 10 von ESRI...
Ich sehe hier, dass es Fortran-Quellcode gibt, um von USNG nach lon / lat zu wechseln: http://www.ngs.noaa.gov/PC_PROD/USNG/ Gibt es Python- oder C / C ++ - Code, der dasselbe tut? Mein Verständnis, dass dies keine Funktion von Proj4 / pyproj / GDAL ist,
Soweit ich die Literatur verstehe, gibt es verschiedene Möglichkeiten, einen Satz geozentrischer kartesischer Koordinaten in geodätische Koordinaten umzuwandeln Bowring, 1976 Borkowski, 1989 Lin & Wang, 1995 Toms, 1996 ... Was wäre der effizienteste Weg, um kartesische Koordinaten in...
Ich habe derzeit ein Projekt, in dem ich Shapefiles in einfache Raster für Anwesenheit / Abwesenheit (1 & 0) rastere und sie dann zu einer laufenden Summe hinzufüge, um am Ende ein Raster für die Gesamtfülle zu erstellen. Die Shapefiles sind in WGS84-Projektion, aber ich brauche das endgültige...
Ich habe eine Datenbank mit vielen Punkten in WGS84. Jetzt erstelle ich einen Cache, der NN ausführt und mithilfe eines KDtree in Bereichsabfragen zeigt. Der Punkt ist, dass der Suchradius in Metern angegeben wird und dass lat / lon kein netter SRS für diese geometrischen Abfragen ist. Ich suche...
Beim Zeichnen der gadmKarte von Russland in habe Rich ein Problem am 180-Grad-Meridian, wo ein Teil der russischen Region separat dargestellt wird. Sie können es auf dem Plot mit dem folgenden Code sehen (es ist die Kante links und der Teil ist rechts). Kann jemand vorschlagen, wie dieses Problem...
Ich habe ungefähr 100-200 Bodenaufnahmen, die ich in ArcGIS 10 überlagern möchte und die keine räumlichen Informationen enthalten, die ich mithilfe einer photogrammetrischen Formel irgendwie entwickeln muss. Die Kamera wurde 2 Meter über dem Boden an einem RTK-GPS-Mast befestigt....
Ich möchte Karten einer "geneigten Erde" erstellen, wie zum Beispiel das folgende Beispiel, in dem Australien am Nordpol liegt: http://i1048.photobucket.com/albums/s379/laskaris_mods/tilted_earth.jpg Mit der von mir verwendeten Software Geocart 3 kann ich Karten in Bildformaten wie .jpeg wie diesen...
Ich habe ein Shapefile der Weltländer, das als Azimuthal Equidistant mit Chicago, USA als Zentrum projiziert wird. Ich möchte die Projektion so ändern, dass Edinburgh, Großbritannien, das Zentrum ist. Man würde denken, dass dies so einfach sein sollte wie das Ändern von 2 Zahlen in den...
Ich bin ein bisschen neu in der gesamten GIS-Welt, und obwohl ich große Fortschritte beim Zeichnen von Polygonen und Ähnlichem in Bing Maps mache, stecke ich derzeit fest. In unseren SQL-Geodaten bin ich auf MULTIPOLYGON-Daten gestoßen, bei denen alle Koordinatenwerte sehr große Zahlen sind (6...
Ich habe im ganzen Internet Dinge gelesen, die behaupten, dass NAD83 und WGS84 effektiv gleichwertig sind (mit geringfügigen Ellipsoidunterschieden). Wenn ich jedoch geodätische (WGS84) -> ECEF -> geodätische (NAD83) Konvertierungen durchführe, bei denen nur die Ellipsoidparameter WGS84 und...
Ich verwende TileStache, um meine eigenen Vektorkacheln mithilfe einer Postgis-Datenquelle bereitzustellen. Zuerst habe ich die Shapefiles, die ich mit ArcMap hatte (neu projiziert in Web Mercator (Auxilary Sphere)), neu projiziert und sie mit shp2pgsql (SRID: 900913) in Postgres geladen. Ich kann...
1. Frage Ich habe eine Warnung mit der Funktion projectRaster () im Rasterpaket in R festgestellt. Ein vollständig reproduzierbares Beispiel wird unten eingefügt. Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 33940 projected point(s) not finite Meine Frage ist:...