Геоинформационные системы

Я создал коллекцию изображений в GEE, и с помощью функции я вычислил индекс NDVI и сопоставил его, чтобы создать другую коллекцию с NDVI в качестве группы. Теперь я хочу создать сложенное изображение с группами NDVI всей коллекции изображений в одно изображение. Так должно быть, как NDVI_1, NDVI_2 и так далее …

10 google-earth-engine  stack 

Это таблица очков. ~ 1 млн. Записей SELECT COUNT(*) as value FROM alasarr_social_mv s; Output: 976270 Похоже, что st_intersects заставляет использовать пространственные индексы, но && этого не делает. Пример использования ST_Intersects(282 мс) SELECT COUNT(*) as value FROM alasarr_social_mv WHERE ST_Intersects( the_geom_webmercator, ST_MakeEnvelope(-410961,4920492,-402305,4926887,3857) ) Aggregate (cost=34370.18..34370.19 rows=1 width=0) (actual time=282.715..282.715 rows=1 …

10 postgis  carto  postgis-2.2 

Я хотел бы знать, есть ли способ добавить поле динамической системы координат в компоновщик компоновщика QGIS. Я могу использовать выражение для добавления многих элементов, таких как имя документа, путь к файлу и т. Д., Но я не могу понять, как добавить систему координат?

10 qgis  coordinate-system  print-composer  qgis-3 

Закрыто . Этот вопрос должен быть более сфокусированным . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он был сосредоточен только на одной проблеме, отредактировав этот пост . Закрыто 2 года назад . Я ищу некоторые ресурсы для изучения PyQGIS. Было бы интересно иметь …

10 pyqgis  references 

У меня есть геоданные с точками и некоторыми связанными данными. Я хочу нанести его на карту, используя геопанды, и чтобы размер точек соответствовал одному из столбцов в геоданных. Пока у меня есть следующий код: base = world.plot(color='white', figsize=(20,10)) geo_df.plot(ax=base, marker='.', color='red', markersize = geo_df['Pop_2005']) plt.xlim([-85, -60]) plt.ylim([-5, 12.5]); Но я …

10 python  geopandas  matplotlib  pandas 

Используя Geometry generatorтип слоя символов, я рисую прямоугольники размеров @nv_bg_w(ширина) и @nv_bg_h(высота) (переменные проекта) для линейных объектов, либо в координатах text_x, text_y(атрибуты, если нет NULL), либо, в качестве альтернативы, в центре линии следующим выражением: geom_from_wkt( 'POLYGON((' || COALESCE("text_x", x(point_on_surface($geometry))) ||' '|| COALESCE("text_y", y(point_on_surface($geometry))) || ','|| (to_real(COALESCE("text_x", x(point_on_surface($geometry)))+ @nv_bg_w )) ||' …

10 qgis  symbology  expression  geometry-generator 

Я новичок в python и ломаю голову над тем, как использовать gdal_merge в (spyder) скрипте python. Я использую Windows 10, Python 3.6 и у меня установлен инструмент gdal из osgeo4w. Я понимаю, что многие другие сообщения описывают эту проблему, но ни одна не может помочь мне решить эту проблему. Когда …

10 python  dem  windows  gdal-merge  spyder 

Я создал несколько линий в EPSG: 32632 (WGS 84 / UTM зона 32N) длиной ровно 10 000 м. QGIS рассчитывает все длины линий (горизонтальные, вертикальные, диагональные) как 10 001,9 м. Инструмент измерения объясняет: «Преобразование Project CRS включено и выбран эллипсоидальный расчет. Координаты преобразуются в выбранный эллипсоид (WGS84), а расстояние рассчитывается …

10 qgis  measurements 

Есть ли способ получить список установленных форматов gdal в системе с python? В командной строке это легко сделать с помощью: gdalinfo --formats ... но можно ли это сделать на python?

10 python  raster  gdal 

Когда я хочу выполнить кригинг, он срабатывает только иногда, в зависимости от того, какие значения я использую в моей таблице данных. В результате функции krige, которую я получаю для var1.pred: NA NA NA ...и var1.var: NA NA NA ...(но только когда я использую «неправильные» значения в моей таблице данных). Например: …

10 r  interpolation  kriging  gstat 

У меня есть очень большой шейп-файл ESRI (> 8000 полигонов, файл .shp> 32 МБ), который я хочу прочитать в R в Windows 7. Я знаком с rgdal и readOGR (). Однако при таком размере шейп-файла считывание всего шейп-файла в память является очень медленным процессом. Кроме того, меня интересуют только несколько …

10 shapefile  r  rgdal 

У меня нет большого опыта работы с QGIS, и я работаю над картой с линиями, соединяющими местоположения, созданные с помощью плагина mmqgis. Есть ли простой способ сделать эти прямые линии несколько пышными для эстетических целей? Я пытаюсь достичь чего-то вроде карты ниже, но с несколькими центрами. Я прочитал предложения о …

10 qgis  qgis-plugins 

Я хочу использовать PyQGIS, чтобы добавить новое поле и рассчитать значения для каждой функции. Аналогично опции Калькулятор поля. Мое выражение «Калькулятор поля», например: y(start_point($geometry)) from PyQt4.QtCore import QVariant from qgis.core import QgsField, QgsExpression, QgsFeature vl = iface.activeLayer() vl.startEditing() #step 1 myField = QgsField( 'myNewColumn', QVariant.Float ) vl.addAttribute( myField ) idx …

10 pyqgis  fields-attributes 

Я немного изо всех сил пытался сгенерировать плитки для изображения с высоким разрешением, которое у нас есть. В настоящее время у нас есть очень большое (+ 20 ГБ) изображение, сохраненное в виде файла GeoTiff. Я хотел бы сгенерировать плитки с помощью утилиты командной строки gdal2tiles, а затем открыть и просмотреть …

10 gdal  gdalwarp  cesium  gdal2tiles 

Я делаю карты для книги, где карты бегут к краю страницы. Принтер посоветовал удалить пространство вокруг изображений на 3 мм, поэтому, если моя страница имеет размер 127x177 мм, изображения должны иметь размер 133x183 мм. Это пространство для отбраковки должно допускать ошибки при вырезании страниц книги, а в случае карт не …

10 qgis  labeling  print-composer