nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Der Herzschlag einer Region

RiccardoKlinger2

<\/HEAD>

Lassen Sie uns mit dem Ergebnis beginnen:<\/P>

Der Herzschlag von Berlin Eine der wichtigsten Fragen an einen Standort ist: Wen oder was kann ich in Minuten erreichen? Daher bietet ArcGIS den Network Analyst mit seinem Werkzeug "Service Areas"<\/A>. Wir können das Gebiet analysieren, das mit dem Auto (andere Verkehrsmittel sind möglich) in einem bestimmten Zeitraum erreichbar ist. Durch die Verwendung mehrerer Eingabepunkte können wir auch Muster in einem gesamten Gebiet analysieren.<\/X><\/P>

Die Verwendung eines normalen Netzwerks würde die Servicebereiche im Tagesverlauf nicht verändern, aber durch die Nutzung des ArcGIS Online-Netzwerkdatensatzes basierend auf HERE-Daten können wir auch die Verkehrsprofile für einen Wochentag und eine bestimmte Uhrzeit verwenden. Verkehrsprofile können in 15-Minuten-Intervallen analysiert werden. Die Analyse der Servicebereiche im Tagesverlauf liefert wertvolle Einblicke in die Verkehrssituation einer Region: Große Unterschiede deuten in erster Linie auf ein höheres Verkehrsaufkommen auf den Straßen hin. Die Gesamtfläche zeigt auch die Qualität des Netzwerks an: Große Flächen deuten auf eine bessere Erreichbarkeit mit dem Auto hin (möglicherweise verursacht durch ein dichtes Netz von Schnellstraßen).<\/P>

Automatisierte Analyse<\/H3>Da ich Geoverarbeitung liebe, interessierte mich die Automatisierung der Analyse und die Visualisierung als Animation. Der Prozess ist daher zweistufig.<\/P>Zuerst habe ich eine einfache Punkt-Feature-Class namens "StartPoints" erstellt, die die Standorte der "facilities" enthält, und diese zur Service Area Analysis Layer hinzugefügt:<\/P>#Erstelle Analyse-Layer<\/SPAN> arcpy.<\/SPAN>na.<\/SPAN>MakeServiceAreaAnalysisLayer(<\/SPAN>"<\/SPAN><A class="jive-link-external-small" href="https:\/\/community.esri.com\/external-link.jspa?url=https%3A%2F%2Fwww.arcgis.com%2F" target="_blank">https:\/\/www.arcgis.com\/<\/A>"<\/SPAN><\/SPAN>, <\/SPAN> "BeatCity Service Layer"<\/SPAN>, <\/SPAN> "Driving Time"<\/SPAN>, <\/SPAN> "FROM_FACILITIES"<\/SPAN>, <\/SPAN> "5;10;15;20;25;30"<\/SPAN>, <\/SPAN> "05.01.1900 00:01:00"<\/SPAN>, <\/SPAN> "LOCAL_TIME_AT_LOCATIONS"<\/SPAN>, <\/SPAN> "POLYGONS"<\/SPAN>, <\/SPAN> "STANDARD"<\/SPAN>, <\/SPAN> "DISSOLVE"<\/SPAN>, <\/SPAN> "RINGS"<\/SPAN>, <\/SPAN> "100 Meters"<\/SPAN>, <\/SPAN> None, <\/SPAN> None)<\/SPAN> #Füge Standorte aus einer Punkt-Feature-Class hinzu<\/Span> arcpy.<\/span>na.<\span>AddLocations(<\span>\"BeatCity Service Layer\"<\span>,<\span> \"Facilities\"<\span>,<\span> \"StartPoints\"<\span>, \"Name Description #;CurbApproach # 0;Attr_Minutes # 0;Attr_TravelTime # 0;Attr_Miles # 0;Attr_Kilometers # 0;Attr_TimeAt1KPH # 0;Attr_WalkTime # 0;Attr_TruckMinutes # 0;Attr_TruckTravelTime # 0;Breaks_Minutes # #;Breaks_TravelTime # #;Breaks_Miles # #;Breaks_Kilometers # #;Breaks_TimeAt1KPH # #;Breaks_WalkTime # #;Breaks_TruckMinutes # #;Breaks_TruckTravelTime # #\", \"5000 Meters\", None, None, \"MATCH_TO_CLOSEST\", \"APPEND\", \"SNAP\", \"5 Meters\", \"EXCLUDE\", None)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<\/span>‍<\\/Span>‍<\\/Span>‍<\\/Span>‍<\\/Span>Jetzt müssen wir die Zeiten automatisch anpassen, um die Tageszeiten widerzuspiegeln. Die Startzeit ist '05.01.1900 00:00:00', was ein Freitag ist und die Zeit Mitternacht.<\/P>Wir können auf die Eigenschaften recht einfach zugreifen. Da ich mit einer 3D-Szene arbeite, ist dies mein einziges Kartenobjekt in meinem ArcGIS Pro-Projekt:<\/P>#Projekt holen<\\/span> doc = arcpy.mp.ArcGISProject('current') #Karte holen<\\/span> map = doc.listMaps()[0] #Servicebereich-Layer holen<\\/span> sa_layer = map.listLayers("BeatCity Service Layer")[0] #Solver-Eigenschaften vom Servicebereich-Layer holen<\\/span> solver_props = arcpy.na.GetSolverProperties(sa_layer)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍Da wir nun Zugriff auf die Eigenschaften haben, können wir einfach iterieren und den Service-Layer lösen. Das Ergebnis wird als 3D-Feature-Class exportiert:<\/P>import os for hour in range(0, 24): #Erstelle das datetime-Objekt<\\/span> date = datetime.datetime(1900, 01, 05, hour, 0, 0) #Wir verwenden später die String-Version<\\/span> datestring = date.strftime("%d.%m.%Y %H:%M") #Setze Tageszeit und Datum:<\\/span> solver_props.timeOfDay = date #Löse den Netzwerk-Layer<\\/span> arcpy.na.Solve("BeatCity Service Layer", "SKIP", "TERMINATE", None, None) arcpy.AddMessage("kopiere Stunde " + str(hour)) #Da wir dies in 3D visualisieren möchten<\\/span> arcpy.ddd.FeatureTo3DByAttribute(r"BeatCity Service Layer\\Polygons") , arcpy.env.workspace + os.sep + "Friday_" + str(hour), "ToBreak", "ToBreak") #und ein Feld mit dem Datum für spätere Verwendung hinzufügen, falls gewünscht: arcpy.AddField_management(arcpy.env.workspace + os.sep + "Friday_" + str(hour), "datetime", "Date", "", "", 8, "datetime", "NULLABLE", "REQUIRED") arcpy.management.CalculateField(arcpy.env.workspace + os.sep + "Friday_" + str(hour), "datetime", '"' + datestring + '"', "PYTHON_9.3") Warnung: Lösung der Ebene verbraucht credits! Die Ergebnisse werden in separaten Feature-Classes gespeichert, was noch nicht optimal ist, aber dies ist momentan meine Lösung. Visualisierung und Export Da wir jetzt separate Layer haben, wende ich denselben Stil auf alle an und steuere die Sichtbarkeit jedes Layers und exportiere das Layout, das vor dem Export erstellt werden muss! Mein Layout sieht ungefähr so aus: Ich habe ein Textsymbol mit dem Standardtext "Friday" im Layout hinzugefügt, um den Zeitstempel auch im PNG-Ergebnis zu exportieren. Außerdem habe ich ein Markersymbol mit dem Städtenamen zum Layout hinzugefügt und die Ansicht etwas geneigt, um den Tiefeneffekt für die extrudierte Ebene zu verstärken. Aber zuerst müssen wir alle Layer für unseren automatisierten Export mit einem ziemlich einfachen Ansatz deaktivieren: #Alle Layer abwählen: for layer in map_obj.listLayers(): if layer.name[0] == "F": layer.visible = False Jetzt können wir durch die Layer iterieren, den Unique Value Renderer anwenden, die Symbole ändern, um Polygone ohne Ränder zu erhalten, die Extrusion einstellen, das Textelement ändern, um den richtigen Zeitstempel zu erhalten und das Layout als PNG exportieren: for hour in range(0, 24): poly_layer = map_obj.listLayers("Friday_" + str(hour))[0] poly_layer.visible = True symbology = poly_layer.symbology symbology.updateRenderer('UniqueValueRenderer') symbology.renderer.fields = ['ToBreak'] symbology.renderer.colorRamp = doc.listColorRamps('Green-Blue (6 Classes)')[0] #Alle Symbole ändern: for sym in symbology.renderer.groups[0].items: sym.symbol.outlineColor = {'RGB': [255, 255, 255, 0]} sym.symbol.size =<\/SPAN> 0.0<\/SPAN>\n     poly_layer.<\/SPAN>symbology =<\/SPAN> symbology\n     #extrude:<\/SPAN>\n     poly_layer.<\/SPAN>extrusion(<\/SPAN>'BASE_HEIGHT'<\/SPAN>,<\/SPAN> "3000-100*[ToBreak]"<\/SPAN>)<\/SPAN>\n     #export:<\/SPAN>\n     datestring =<\/SPAN> "2017\/12\/08 "<\/SPAN> +<\/SPAN> str(<\/SPAN>hour)<\/SPAN> +<\/SPAN> ":00"<\/SPAN>\n     for<\/SPAN> lyt in<\/SPAN> doc.<\/SPAN>listLayouts(<\/SPAN>)<\/SPAN>:<\/SPAN>\n          for<\/SPAN> elm in<\/SPAN> lyt.<\/SPAN>listElements(<\/SPAN>"TEXT_ELEMENT"<\/SPAN>)<\/SPAN>:<\/SPAN>\n               if<\/SPAN> elm.<\/Span>text[0]== "2"\n                   elmtext = datestring\n     lytexportToPNGarcpyenvworkspace + ossep + friday_seattle_final" + strhour+ .png"      poly_layervisible = False<<EMOJI_0>><EMOJI_1><<EMOJI_2><<EMOJI_3><<EMOJI_4><<EMOJI_5><<EMOJI_6><<EMOJI_7><<EMOJI_8><<EMOJI_9><<EMOJI_10><<EMOJI_11><<EMOJI_12><<EMOJI_13><<EMOJI_14><<EMOJI_15><<EMOJI_16><<EMOJI_17><<EMOJI_18><<EMOJI_19><<EMOJI_20><<EMOJI_21><<EMOJI_22><</spa n>Leider dauert der gesamte Rendering- und Exportvorgang auf meinem Lenovo X270 ohne dedizierte GPU etwa 2,5 Minuten pro Szene.<\/P>Das Ergebnis sind eine Reihe von PNGs, die mit verschiedenen GIF-Erstellungsseiten im Web oder dem Open-Source-Bildeditor <\/> <A href=\"https:\/\/www.gimp.org\/tutorials\/Simple_Animations\/\

Finde mehr Beiträge mit dem Tag

Knowhow

Kommentare

There are no comments yet

Nehmen Sie an der Diskussion teil

Angemeldete Mitglieder können Beiträge verfassen, Updates folgen und mehr. Neu hier? Registrieren Sie ein kostenloses Konto.

Schnellzugriffe

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.