Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
https://gis-lab.info/forum/
1. А как их установить, и нет ли GUI-аналога?Александр Мурый писал(а):Можно попробовать через питоновский модуль shapely, как написано здесь. Или через GRASS (там есть пару дополн. модулей для получения средних точек).
Код: Выделить всё
"xypoint"=(midPoint.to_wkt())У вас точно GRASS 7 стоит? Потому как <v.line.center> работает только под GRASS 6.4.* (и только под Linux, т.к. это шелл-скрипт), а <v.centerpoint> - только под GRASS 7.vsugig писал(а):Спасибо. Но что-то я его не нашёл в списке, есть v.line.center. А в QGIS в панель инструментов GRASS нельзя его установить?
Код: Выделить всё
center_points.py lines.shp cent_pts.shpКод: Выделить всё
#!/usr/bin/env python
# -*- coding: utf-8 -*-
import sys, os
from osgeo import ogr
from shapely import wkt
from shapely.geometry import LineString
in_lns = sys.argv[1]
out_pts = sys.argv[2]
# open shapefile with OGR
input = ogr.Open(in_lns)
in_lyr = input.GetLayer(0)
## create output point shapefile
ogrDriver = ogr.GetDriverByName("ESRI Shapefile")
outDataSource = ogrDriver.CreateDataSource(out_pts)
# copy SRS from the input shapefile
srs = in_lyr.GetSpatialRef()
# create layer
out_lyr = outDataSource.CreateLayer(out_pts, srs, ogr.wkbPoint)
# export to WKT, compute center points and write them to shp
for f in in_lyr:
g = f.geometry().ExportToWkt()
shapelyLine = LineString(wkt.loads(g))
midPoint = shapelyLine.interpolate(shapelyLine.length/2).to_wkt()
# write point features to layer
feature = ogr.Feature(out_lyr.GetLayerDefn())
point = ogr.CreateGeometryFromWkt(midPoint)
feature.SetGeometry(point)
out_lyr.CreateFeature(feature)
feature.Destroy()
# Destroy the data source to free resources
outDataSource.Destroy()
center_points.pyТут вот 7, а десктопная 6,4:Александр Мурый писал(а):У вас точно GRASS 7 стоит? Потому как <v.line.center> работает только под GRASS 6.4.* (и только под Linux, т.к. это шелл-скрипт), а <v.centerpoint> - только под GRASS 7.vsugig писал(а):Спасибо. Но что-то я его не нашёл в списке, есть v.line.center. А в QGIS в панель инструментов GRASS нельзя его установить?
Результат:Александр Мурый писал(а):Вот простой Питон-скрипт с использованием библиотек <OGR> и <shapely>, на входе - линейный шейп-файл, на выходе - точечный с "серединными" точками.
Код: Выделить всё
center_points.py lines.shp cent_pts.shp
С взаимоотношениями QGIS и GRASS 7 без пол-литра не разберешься. Как-то добавлять модули в Processing можно, но как именно, быстро не подскажу. Сам такого пока не делал.vsugig писал(а): меня интересует-то главным образом эта, панелька. Она зависит от десктопной или нет?
vsugig писал(а):Можно пример с конкретным синтаксисом путей в Windows? И в какой консоли его запускать? Консоль питон в кугис/коммандная строка в осгеов/cmd?
Код: Выделить всё
C:\Program Files\GDAL;C:\Python27;C:\Python27\Scripts Код: Выделить всё
No module named osgeoКод: Выделить всё
C:\Program Files\GDAL;C:\OSGeo4W64\apps\Python27\;C:\OSGeo4W64\apps\Python27\Scripts\; ***Код: Выделить всё
from qgis.core import *
import qgis.utils
layer = qgis.utils.iface.mapCanvas().currentLayer()
iter = layer.getFeatures()
crs = layer.crs().toWkt()
out_name = layer.name() + '_centerpoints'
vl = QgsVectorLayer("Point?crs=" + crs, out_name, "memory")
pr = vl.dataProvider()
fet = QgsFeature()
for feature in iter:
geom = feature.geometry()
length = geom.length()
mid = geom.interpolate(length/2)
fet.setGeometry(mid)
pr.addFeatures([fet])
QgsMapLayerRegistry.instance().addMapLayer(vl)
Код: Выделить всё
Traceback (most recent call last):
File "<input>", line 1, in <module>
File "G:/center_points.py", line 6, in <module>
in_lns = sys.argv[1]
IndexError: list index out of rangeКак-то можно, естественно. Надо разбираться. Просто это не входило в изначальную задачу.vsugig писал(а):Однако, в этом случае в слое не сохраняются атрибуты исходного линейного, что желательно получить, и связи с исходным слоем тоже не работают. Можно ли как-то таким методом унаследовать родительские атрибуты?
Это нормально, т.к. первый скрипт не предназначен для запуска в Питон-консоли, только в системном терминале.vsugig писал(а): При попытке запуска первого скрипта в редакторе консоли питон выходит ошибка: