Data+ International - Вопросы и Ответы - Вопросы с тегом python

Ошибка cv2.cvtColor в модуле ScannedMapDigitizer

Mon, 18 Apr 2022 08:37:41 +0000

Пытаюсь повторить https://developers.arcgis.com/python/guide/geo-referencing-and-digitization-of-scanned-maps/

from arcgis.learn import ScannedMapDigitizer

smd = ScannedMapDigitizer(r"D:\Arcgis\TEST\RASTER", r"D:\Arcgis\TEST\RASTER")
smd.create_mask(color_list=[[115, 178, 115]], color_delta=5, kernel_size=[2, 2], kernel_type="rect", show_result=False)
smd.create_template_image([100, 100,100], 10, 4, show_result= False)
extent = {
'spatialReference': {'wkid': 4326},
'xmin': -180,
'ymin': -89,
'xmax': 180,
'ymax': 85
}

smd.set_search_region_extent(extent)
ssearch_extent = {
'spatialReference': {'wkid': 4326},
'xmin': -180,
'ymin': -89,
'xmax': 180,
'ymax': 85
}
smd.prepare_search_region(r"D:\Arcgis\TEST\RASTER\M1\cover2_admin_Copy.shp", [255, 255, 255], search_extent, 2068, 3744, show_result=True)
smd.match_template_multiscale(0.2, 2.0, 40, show_result=False)
smd.georeference_image((150, 50), show_result=False)
smd.digitize_image(show_result=False)

Однако на этапе создания маски smd.create_mask(color_list=[[115, 178, 115]], color_delta=5, kernel_size=[2, 2], kernel_type="rect", show_result=False) появляется ошибка :

error: OpenCV(4.5.5) D:\a\opencv-python\opencv-python\opencv\modules\imgproc\src\color.cpp:182: error: (-215:Assertion failed) !_src.empty() in function 'cv::cvtColor

Подскажите что может быть не так по той информации что нашел ошибка означает что нету доступа к изображению хотя он есть и путь прописан правильно

Добавление параметра в скрипт

Wed, 03 Feb 2021 12:51:22 +0000

Имеется скрипт

import arcpy

domains = arcpy.da.ListDomains("Database Connections\sde.sde") #база из которой берутся доменные (эталонные) значения
for dom in domains:
... arcpy.GetParameterAsText(0)

как правильно прописать чтобы передать список доменов которые есть в базе ?

Что нужно выбрать в закладке Parameters чтобы они отобразились ?

Пробовал String но там ничего не показывает

Версия 10.5.1

Кирилица в наименование пути

Wed, 27 May 2020 14:20:49 +0000

Имеется скрипт

mxd=arcpy.mapping.MapDocument(u"R:/Проект/Проверка.mxd")

>>> for lyr in arcpy.mapping.ListLayers(mxd):

... print lyr.dataSource

Выдаёт ошибку

Runtime error

Traceback (most recent call last):

File "<string>", line 1, in <module>

File "c:\program files (x86)\arcgis\desktop10.5\arcpy\arcpy\arcobjects\mixins.py", line 652, in __init__

super(MapDocumentMethods, self).__init__(mxd)

File "c:\program files (x86)\arcgis\desktop10.5\arcpy\arcpy\arcobjects\_base.py", line 47, in __init__

for arg in args))

RuntimeError: Object: CreateObject cannot open map document

С mxd файлом имя которого на английском такую ошибку не выдаёт. В чем может быть проблема

Как узнать какие таблицы присоединены к слою в mxd

Fri, 17 Apr 2020 08:41:53 +0000

Необходимо написать скрипт который бы показывал какие таблицы присоединены (Join) к слою/слоям в mxd файле. Через supports не нашёл ничего подходящего. Поиск привел к этому https://gis.stackexchange.com/questions/7703/detecting-join-programmatically-using-arcpy где пользователь пишет что в arcpy.Layer class нет такой функции. В связи с этим вопрос как организовать код чтобы можно было вытащить наименование таблицы

Как отцентрировать надпись

Fri, 10 Nov 2017 11:30:36 +0000

Задача заключается в следующем. В ArcMap 10.2.2 необходимо надписать точку. При этом надпись должна быть представлена в виде дроби. Дробь должна быть выровнена по центру. Значения для дроби берутся из разных полей.

Выражение на Phyton я смог написать: u'<UND>{0:^5}</UND>\n{1:^5}'.format([Field1], [Field2]), однако после написания одних объектов все нормально, после написания других gis просто вылетает и даже не выдает ошибки.

Вопрос, есть ли аналог выражения на VBScript? Или может кто то знает другой способ выранивания дроби по центру?

P.s. В аннотации переводить не хочу, дабы не захломлять базу.

arcpy.da.UpdateCursor повторно походит по записям в таблице

Wed, 18 Oct 2017 01:56:15 +0000

Здравствуйте. Столкнулся со странным поведением курсора, который обновляет записи в таблице. Возможно, кто-то сталкивался с похожей проблемой и подскажет решение.

На сервере установлен ArcGIS for Server. В СУБД MS SQL создана база пространственных данных. В базе есть полигональный слой, который необходимо периодически обновлять на основании внешних шейп-файлов. Был подготовлен скрипт на python, который читает шейп-файлы и обновляет записи в базе данных, если выполняются определенные условия. Запускается скрипт на том же сервере.

Код метода для обновления:

class Dataset(object):
    def update(self, new_features):
        with arcpy.da.Editor(self.workspace):
            with arcpy.da.UpdateCursor(self.path, self.fields) as cursor:
                for row in cursor:
                    feature = Feature.create(self.fields, row)
                    if feature.id in new_features:
                        new_feature = new_features[feature.id]
                        if new_feature.date > feature.date:
                            cursor.updateRow(new_feature.as_row(self.fields))
                        else:
                            pass
                    else:
                        cursor.deleteRow()

Проблема заключается в том, что курсор проходит по всем записям в таблице (обновляя, удаляя или пропуская каждую запись), а потом проходит по некоторым из обновленных записей повторно. Поясню на примере: в таблице всего 10000 записей. Курсор прошел по всем 10000, и еще раз по 200 обновленным записям. Причем иногда это может быть 10 лишних записей, иногда 800, а иногда и 0. Закономерность уследить не удалось. Почему курсор проходит по обновленным записям повторно?

Версия ArcGIS for Server 10.3.1, СУБД MS SQL Server 2012, скрипт подключается к базе под учетной записью владельца данных. Версионность для обновляемого слоя не включена, т.к. никто кроме владельца данных (скрипта) изменять данные не должен. Спасибо

Экспорт вложений из ГБД с кириллическими символами

Tue, 04 Jul 2017 18:45:33 +0000

Нашел скрипт .py экспортирующий вложения из ГБД.

С латинскими названиями вложений работает, но при попадании на кириллические символы стопорится и вылетает.

Как указать правильную кодировку для работы скрипта с лат и кириллическими символами?

    import arcpy
    from arcpy import da
    import os

    inTable = arcpy.GetParameterAsText(0)
    fileLocation = arcpy.GetParameterAsText(1)

    with da.SearchCursor(inTable, ['DATA', 'ATT_NAME', 'ATTACHMENTID']) as cursor:
    for item in cursor:
    attachment = item[0]
    filenum = "ATT" + str(item[2]) + "_"
    filename = filenum + str(item[1])
    open(fileLocation + os.sep + filename, 'wb').write(attachment.tobytes())
    del item
    del filenum
    del filename
    del attachment

Изменение части строк в атрибутивной таблице

Mon, 29 May 2017 12:57:49 +0000

Для добавление русских символов в строку, например, названия рек, я использовала следующую команду u'р. ' + !Name!, а как убрать лишние символы из строки чтобы осталось только название, т.е.из поля Name в строке "город Болгар", убрать слово "город").