Запитання з тегом «python»

У мене є растр, з яким я хотів би зробити деякі точкові інтерполяції. Ось де я перебуваю: from osgeo import gdal from numpy import array # Read raster source = gdal.Open('my_raster.tif') nx, ny = source.RasterXSize, source.RasterYSize gt = source.GetGeoTransform() band_array = source.GetRasterBand(1).ReadAsArray() # Close raster source = None # Compute …

12 python  algorithm  interpolation 

Як я можу переключити видимість шару в QGIS за допомогою API Python?

12 python  qgis 

Я хотів би отримати перетин декількох багатокутників. Використовуючи shapelyпакет Python , я можу знайти перетин двох многокутників за допомогою intersectionфункції. Чи існує подібна ефективна функція для отримання перетину декількох багатокутників? Ось фрагмент коду, щоб зрозуміти, що я маю на увазі: from shapely.geometry import Point coord1 = ( 0,0 ) point1 …

12 python  polygon  intersection  polygon-creation  shapely 

Я запускаю Ubuntu 14.4 і встановив віртуальне середовище python. Потім я встановив декілька бібліотек python (gdal, lxml, psycopg2 тощо) Також у моїй системі є модуль для QGIS. Тому коли я ввожу python, я можу імпортувати цей модуль як: import qgis.core Але коли я намагаюся зробити те саме через свій virtualenv, …

12 qgis  python 

Чи існує пакет Python, який забезпечує реалізацію алгоритму Прямий скелет? Я знаю, що проект CGAL з відкритим кодом (C ++) містить реалізацію, але, схоже, cgal-прив'язки не включають цей пакет CGAL . У будь-якому випадку я віддаю перевагу чистій реалізації Python, яку я міг би змінити / розширити відповідно до моїх …

12 python  polygon  geometry  skeleton  computational-geometry 

Я хотів би оновити один атрибут функції. Однак я не знаю оновлювати його за допомогою об’єкта функції. Я маю використовувати провайдера даних для його оновлення. layers = QgsMapLayerRegistry.instance().mapLayersByName('my_line') layer = layers[0] dp = layer.dataProvider() it = dp.getFeatures() for i in range(0, dp.featureCount()): feat = it.next() attrs = { 2 : …

12 qgis  python  pyqgis  features 

У мене виникають проблеми з використанням функції Polygonize в python. Приклад кулінарної книги для цього можна знайти тут . Відповідна частина мого коду: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized" driver = ogr.GetDriverByName("ESRI Shapefile") if os.path.exists(outShapefile+".shp"): driver.DeleteDataSource(outShapefile+".shp") outDatasource = driver.CreateDataSource(outShapefile+ ".shp") outLayer = outDatasource.CreateLayer("polygonized", srs=None) …

12 python  gdal 

Я читаю форму файлів за допомогою Fiona: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] це дає мені дикту, що містить перелік списків координат і введіть "MultiPolygon". Однак я не впевнений, як використовувати дані координат для створення Shapely MultiPolygon, оскільки це вимагає a sequence …

12 python  shapely  fiona 

Чи є спосіб виконати те саме завдання, що й утиліта gdalbuildvrt, використовуючи прив’язки GDAL Python? Поки що я не знайшов іншого способу зробити це, окрім як створити сад одного набору даних та вручну редагувати xml. Я хотів би створити сад з декількох растрових (по суті виконують мозаїку). Це можливо за …

12 python  raster  gdal  vrt  gdalbuildvrt 

[EDIT: рішенням цього було просто використовувати OGR для читання файлів. Дивіться приклад geographika.] У файлі форм ESRI немає різниці між полігонами та мультиполігонами. Крім того, немає чіткого розрізнення між внутрішніми отворами та зовнішніми кільцями (окрім "ручності" даного багатокутника). Отже, прочитавши файл форми, у мене є список послідовностей координат, що описують …

12 python  shapefile  polygon  shapely 

У мене є форма файлів підрозділів країни, і я хотів би отримати масив обмежувальних lat і lon для кожного підрозділу .. чи можна це зробити?

12 python  shapefile  ogr 

Які драйвери має пакет python fiona? Коли я перевіряю посібник, він пише, [...] and the possible formats are enumerated in the fiona.drivers list. щоправда, коли я набираю python from fiona import drivers print drivers > <function drivers at 0x108763050> Як я можу заглянути в це?

12 python  fiona 

Я працюю над інструментом, який перетворює растрові шари в масиви для обробки з NumPy, і в ідеалі я хотів би мати можливість працювати з растрами, які поставляються упакованими в .gdb, не експортуючи їх усіх (особливо, якщо для цього потрібно залучати ArcGIS або ArcPy ). Чи можливо це за допомогою драйвера …

12 python  gdal  file-geodatabase 

використовуючи Gdal / OGR, чи є спосіб перетворити рівень лінії shp, kml або PostGIS в мережу вузлів та сегментів, щоб його можна було використовувати в таких пакетах, як networkx?

12 python  topology  ogr  networkx 

Я намагаюся відкрити форм-файл із ноутбука Юпітера, але не можу його відкрити. Файл форми розміщений у тій самій папці, що і моя робоча книга, і Python не може його визначити, навіть якщо я використовую повний шлях. import shapefile as shp import matplotlib.pyplot as plt sf = shp.Reader(r'C:\Users\Public\Documents\1-11-99n.shp') ShapefileException: Unable to …

11 python  shapefile