Systemy Informacji Geograficznej

Czy w QGIS można policzyć liczbę wierzchołków wielokątów i linii w danej warstwie i utworzyć osobne pole zawierające te wartości? Podejrzewam, że jest to funkcja gdzieś w kalkulatorze polowym, ale nie mogę jej znaleźć.

14 qgis  vertices  count 

Chcę zaimportować dużą liczbę polilinii i / lub wielokątów do jednej lub więcej warstw QGIS w celu wykreślenia danych. Każdy będzie miał tylko proste atrybuty (np. Nazwę). Muszę wiedzieć: Wymagany format pliku CSV, aby umożliwić krok 2. Polecenie lub wtyczka są wymagane do utworzenia polilinii lub warstwy wielokątów w QGIS. …

14 qgis  csv  import 

Próbuję dokończyć wybór według atrybutu w Pythonie, ale na podstawie zapytania, czy atrybut jest obecny na liście. Takie zapytanie w najprostszej formie powinno być takie: qry = " \"OBJECTID\" in oid_list" arcpy.SelectLayersByAttribute_management(inft, "NEW_SELECTION", qry) ale to podejście zwraca nieprawidłowy błąd wyrażenia. W przeszłości musiałem używać bardziej skomplikowanej sytax dla tego …

14 arcpy  arcgis-10.0  sql  select-by-attribute  where-clause 

Regularnie tworzę własne rastry GeoTIFF z GDAL w Pythonie, np .: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds = driver.Create('MyRaster.tif', 200, 300) band = ds.GetRasterBand(1) band.WriteArray(data) ds = band = None # save, close Jednak gdy wynik jest wyświetlany …

14 python  gdal  geotiff-tiff  statistics 

ArcGIS 10 oferuje nową możliwość dodawania załączników do poszczególnych obiektów do klas obiektów przechowywanych w geobazie. Jednym z zastosowań tego byłoby dodanie wielu zdjęć do określonej lokalizacji, hydranty przeciwpożarowe wydają się być popularnym przykładem (ważne, aby wiedzieć, czy jesteś psem w Internecie, który potrzebuje postoju). W bazie danych struktura klasy …

14 arcgis-desktop  python  convert  esri-geodatabase  attachments 

Mam proste pytanie dotyczące obliczania odległości w PostGIS. Chciałbym uzyskać odległość między dwiema geometriami. Używam tego sid: 4269 w metrach, co teraz robię, to: ST_Distance((a.geom,b.geom)) FROM ...ale otrzymuję wynik w stopniach. Myślę, że powinienem pracować z geografią, ale jak mogę rzucić geom na geografię? Próbowałem z (a.geom :: geography), ale …

14 postgis  coordinate-system  postgresql  distance 

Niektóre pliki shapefile są powiązane z plikiem .prj, a plik .prj zawiera informacje o projekcji pliku shapefile w formacie WKT. Czasami muszę przekonwertować WKT na ciąg proj4, a czasem muszę go przekonwertować. Czy jest na to gotowa biblioteka?

14 shapefile  coordinate-system  wkt  proj 

Używam plotrixw Rcelu wykreślenia mapy stanu USA. Ma doskonałą funkcję floating.pieumieszczania wykresów kołowych w każdym stanie. Zastanawiam się, czy w plotrixpakiecie znajdują się podobne funkcje do wyświetlania wykresów słupkowych w każdym stanie? (Rzuciłem okiem na dokumentację i funkcje, które radzą sobie z wykresami słupkowymi, nie wydają się mieć takiej możliwości, …

14 r  statistics  chart 

Zbudowałem prostą aplikację do map, w której mam warstwę na mapie. Jak mogę sprawdzić, czy warstwa jest już obecna w Ulotce?

14 leaflet 

Chcę otworzyć wyskakujące okienko u dołu ikony znacznika w Ulotce. mój kod: var mymap = L.map('mapid').setView([51.505, -0.09], 13); L.tileLayer('https://api.tiles.mapbox.com/v4/{id}/{z}/{x}/{y}.png?access_token=pk.eyJ1IjoibWFwYm94IiwiYSI6ImNpandmbXliNDBjZWd2M2x6bDk3c2ZtOTkifQ._QA7i5Mpkd_m30IGElHziw', { maxZoom: 18, attribution: 'Map data © <a href="http://openstreetmap.org">OpenStreetMap</a> contributors, ' + '<a href="http://creativecommons.org/licenses/by-sa/2.0/">CC-BY-SA</a>, ' + 'Imagery © <a href="http://mapbox.com">Mapbox</a>', id: 'mapbox.streets' }).addTo(mymap); var lati = 51.51; var longi = -0.09; …

14 leaflet  popup 

Korzystam z poniższego kodu, aby znaleźć kraj (a czasem stan) dla milionów punktów GPS. Kod zajmuje obecnie około jednej sekundy na punkt, co jest niezwykle powolne. Plik shapefile ma 6 MB. Czytałem, że geopandy wykorzystują rtree do połączeń przestrzennych, co czyni je niesamowicie wydajnymi, ale wydaje się, że to tutaj …

14 python  spatial-join  geopandas  rtree 

Mam duży plik rastrowy (129600 na 64800 pikseli) z globalnymi częściami wód (1 bit wartości 0 i 1) i próbuję wydobyć linie brzegowe oceanów i wód śródlądowych. Próbowałem z ArcGIS i QGIS przekonwertować z rastra na polilinię, ale to trwa wieki. Czy ktoś zna lepszy / szybszy sposób (Python lub …

14 python  r  performance 

Chcę zaimportować plik shp do postgreSQL. Najpierw tworzę plik sql, a następnie uruchamiam PostgreSQL. Aby utworzyć plik sql, uruchamiam polecenie w Windows cmd: shp2pgsql -s 4326 worldCountries.shp worldcountries postgres > worldcountries.sql a następnie uruchom: psql -d postgres -U postgres -p 4321 -f worldcountries.sql ale wynikiem jest: psql: worldcountries.sql: 21: ERROR: …

14 postgis  shapefile  postgresql  import 

Szukam sposobu na utworzenie diagramu voronoi na podstawie segmentów linii zamiast punktów. Zobacz poniższy przykład (zaczerpnięty z tej prezentacji ). Idealnie chciałbym coś, co mogę napisać w skrypcie z Pythona, ale akceptowalne byłoby również rozwiązanie korzystające z ArcGIS lub podobnego. Jedyną biblioteką, która to zrobiła, jaką do tej pory znalazłem …

14 line  voronoi-thiessen 

Stworzyłem obiekt typu SpatialPointsDataFrameza pomocą sppakietu w R. Jednak jestem zdezorientowany co do @, $, . and []operatorów i kiedy ich używać, aby uzyskać dostęp do różnych właściwości mojego obiektu. Oto mój przykładowy kod: library(sp) library(rgdal) #creating a SpatialPointsDataFrame with sample points in UTM x <- c(15.2, 15.3, 15.4, 15.5, …

14 r  syntax