Вопросы с тегом «python»

Я репроектирую растры в Python, используя GDAL. Мне нужно спроецировать несколько tiff из географических координат WGS 84 в WGS 1984 Web Mercator (Вспомогательная сфера), чтобы потом использовать их в Openlayers вместе с OpenStreetMap и, возможно, картами Google. Я использую Python 2.7.5 и GDAL 1.10.1 отсюда , и преобразовываю координаты, используя …

17 python  gdal  coordinate-system  epsg  web-mercator 

Я часто вижу в коде Python GDAL, что люди закрывают наборы данных в конце своего скрипта. Почему имеет смысл закрывать набор данных в Python GDAL? Есть ли последствия, если я не сделаю это? import gdal # open dataset ds = gdal.Open('test.tif') # close dataset ds = None

17 python  gdal 

В настоящее время я работаю над тем, чтобы импортировать климатические данные CANGRID (предоставленные как Surfer Grid ascii, файлы .grd) в ArcGIS. Сетка размером 95 строк на 125 столбцов. Метаданные обеспечивают широту / долготу происхождения (нижний левый угол), размер ячейки (50 км), а также проекцию нот в виде полярной стереографии с …

16 arcgis-desktop  python  gdal  coordinate-system 

Закрыто . Этот вопрос должен быть более сфокусированным . В настоящее время не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он был сосредоточен только на одной проблеме, отредактировав этот пост . Закрыто 9 месяцев назад . За исключением ArcPy, есть ли какая-либо библиотека python, которая может выполнять геообработку, …

16 python  buffer  software-recommendations 

Мне интересно узнать, как использовать массивы NumPy для оптимизации геообработки. Большая часть моей работы связана с «большими данными», где для геообработки часто требуются дни для выполнения определенных задач. Излишне говорить, что я очень заинтересован в оптимизации этих процедур. ArcGIS 10.1 имеет ряд функций NumPy, к которым можно получить доступ через …

16 arcgis-desktop  python  arcpy  numpy  big-data 

Я хочу программно создать файл GeoJSON, используя полигоны из шейп-файла, но добавляя атрибуты из моего собственного приложения. Это легко сделать для шейп-файла: def create_data_dayer(self,varlist, data): """ Creates a new shape to contain data about nodes. varlist is the list of fields names associated with the nodes. data is a list …

16 python  geojson  ogr  fiona 

У меня есть часть программного обеспечения (не Arc), которая работает ночью на непроизводственной машине, которая обновляет функции шейп-файла из внешней базы данных. Периодически я хотел бы удалить все функции в шейп-файле (а не сам файл, который должен остаться) и позволить программному обеспечению "перестроить" шейп-файл с нуля. Я хотел бы автоматизировать …

16 python  shapefile  open-source-gis 

Я работаю с надстройкой ArcMap в C #. Из кода C # я выполнил несколько скриптов Python. Теперь, чтобы запустить этот скрипт, у меня есть жестко запрограммированный путь к Python. Но это не портативно. Итак, я хочу получить путь к исполняемому файлу Python из кода и использовать его. Вопрос: Как …

16 arcgis-10.0  python  c# 

Я действительно новичок в Python и хотел бы знать, существует ли быстрый метод получения значений ячеек растрового пикселя за пикселем и координат (координата XY центра каждого пикселя) с использованием Python в ArcGIS 10? Чтобы описать это далее, мне нужно получить карту X, карту Y и значение ячейки первого пикселя, назначить …

16 arcgis-10.0  python  raster 

Извиняюсь, если следующий вопрос несколько глуп, но я только ОЧЕНЬ новичок во всей этой ГИС. Я пытаюсь преобразовать некоторые проецируемые изображения GeoTiff в WGS84, используя gdal в python. Я нашел сообщение, в котором описывается процесс преобразования точек в проецируемых геотиффах с использованием чего-то похожего на следующее: from osgeo import osr, …

16 python  gdal  coordinate-system  wgs84 

Пакет geopy для Python содержит две методики измерения расстояний: формулы Большого круга и Винсенти . >>> from geopy.distance import great_circle >>> from geopy.distance import vincenty >>> p1 = (31.8300167,35.0662833) # (lat, lon) - https://goo.gl/maps/TQwDd >>> p2 = (31.8300000,35.0708167) # (lat, lon) - https://goo.gl/maps/lHrrg >>> vincenty(p1, p2).meters 429.16765838976664 >>> great_circle(p3, p4).meters …

16 python  distance  great-circle  geopy  vincenty-formulae 

Учитывая наличие geopandas GeoDataFrameряда полигонов, я хотел бы получить площадь в км кв каждой функции в моем списке. Это довольно распространенная проблема, и обычно предлагаемое решение в прошлом было использовать shapelyи pyprojнапрямую (например, здесь и здесь ). Есть ли способ сделать это в чистом виде geopandas?

16 python  area  shapely  geopandas  pyproj 

я использую данные ГИС моего округа, чтобы обновить дороги в OpenStreetMap. У меня есть копия, которую они опубликовали в прошлом году, и копия, которая была только что опубликована, и я хотел бы найти все строки строк, у которых либо были изменены свои атрибуты, либо изменена их геометрия. Это поможет мне …

15 qgis  python  linestring  change-detection 

Коротко: я хочу изменить разрешение растра и сгладить серую рампу, как показано на изображениях ниже. Предпочтение отдается использованию GDAL, PIL или Numpy. Описание: Я делаю точки в растры с выходным разрешением 20 метров с помощью высокопроизводительной геостатистической библиотеки . Я не хочу менять выходное разрешение, потому что время интерполяции увеличивается …

15 python  raster  gdal 

Я ищу другое, более элегантное решение проблемы пространственной статистики. Необработанные данные состоят из координаты xy для каждого отдельного дерева (т.е. преобразуются в точечный файл .shp). Хотя это и не используется в этом примере, каждое дерево также имеет соответствующий многоугольник (то есть как .shp), который представляет диаметр кроны. На двух изображениях …

15 arcgis-desktop  python  r  remote-sensing  spatial-statistics