ArcGIS + ENVI реализует реализацию оценки точности классификации дистанционного зондирования (метод стратифицированной выборки) 您所在的位置:网站首页 rastervalu ArcGIS + ENVI реализует реализацию оценки точности классификации дистанционного зондирования (метод стратифицированной выборки)

ArcGIS + ENVI реализует реализацию оценки точности классификации дистанционного зондирования (метод стратифицированной выборки)

2023-04-22 05:35| 来源: 网络整理| 查看: 265

Недавно отредактированный документ достиг окончательного порога, и результат ставится под сомнение в отношении метода оценки точности классификации.

Мой первоначальный метод заключался в использовании ArcGIS для создания сетки с центральной точкой сетки в качестве точки визуальной интерпретации, а затем с использованием ArcGIS для сбора результатов классификации изображений этой точки и вычисления матрицы неточностей.

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

Затем точки каждой категории на классифицированном изображении случайным образом разбрасываются в позиции категории, а затем визуально интерпретируются вручную, и матрица неточности вычисляется вместе с результатом классификации изображения.

Я сначала попробовал Erdas, сначала будут проблемы после установки erdas на компьютер с ArcGIS, а потом я напишу обзорную статью. И, наконец, обнаружил, что сложно импортировать точки, выбранные erdas, в ArcGIS или существует таблица Excel для облегчения вычисления матрицы путаницы, поэтому я сдался.

После поиска в Интернете удобнее искать ArcMap 10.4 и выше и выполнять стратифицированную выборку и точки распространения, но в соответствии с более высокой версией это рискованно. В настоящее время я использую 10.2

В ArcGIS 10.2 нет возможности разбросать точки в определенное положение на классифицированном изображении, но я наконец обнаружил, что ENVI может использовать эту функцию.

Текст такой: ---------------------------------------- Вы можете увидеть, если вы не хотите читать мое описание здесь ------------------------------

1. Классификация

Прежде всего, результат моей классификации - это изображение tif. Значения 1, 2 и 3 представляют каждую категорию. Tif не может быть непосредственно выбрана в ENVI. Для классификации в ENVI используйте дерево решений, которое является решением классификации дерево. Прежде всего Сначала создайте дерево решений. Обратите внимание, что это для вас, чтобы выбрать полосу b1. Не выбирайте ее. Я выбрал ее раньше, а затем заменил на другое изображение. После сохранения используйте классификацию-дерево решений-выполнить выход из дерева решений, чтобы выбрать ранее сохраненное дерево решений txt. Выберите изображение результата классификации, а затем вы можете получить изображение классификации.

2. Место отбора проб

Используйте классификацию envi - пост-классификацию - сгенерируйте случайную выборку с использованием достоверного изображения

Затем выберите нужное изображение, нажмите ОК

Затем вы можете выбрать категорию (обратите внимание, что я присвоил категорию изображению, которое не классифицировано, то есть «нет», вы можете не выбирать ее)

После нажатия кнопки ОК выберите стратифицированную выборку и дипропорциональную выборку, а затем нажмите установить размеры выборки классов, вы можете указать, сколько точек будет разбросано в каждой категории.

Остальные варианты остаются неизменными для сохранения рентабельности инвестиций.

Откройте изображение потом

Выберите интересующую область наложения, вы можете увидеть сгенерированные точки (фактически, ROI). Обратите внимание, что ROI здесь может быть сгенерирован раньше, а не ваш на этот раз, удалите его. Выберите в следующем окне, выберите File-output roi to shapefile,

Как показано ниже, выберите категорию, которую нужно сохранить, в основном выберите все

3. Шейп-файл для получения точек классификации изображений + визуальная интерпретация.

Откройте файл shp в ArcGIS, затем выберите инструменты пространственного анализа-извлечение-извлечение значений для точек, выберите файл и изображение shp, и, наконец, вы можете получить векторный файл.Существует столбец атрибута RASTERVALU, который является значением Изображение.

После этого просто откройте таблицу атрибутов, добавьте новый столбец в поле добавления в результате визуальной интерпретации, а затем щелкните правой кнопкой мыши и используйте растровый калькулятор, чтобы сделать его равным столбцу RASTERVALU (этот шаг предназначен для уменьшения рабочей нагрузки, конечно, вы можете присвоить значение самостоятельно), а затем вы можете визуально интерпретировать и изменять значение.

4. Вычислить матрицу неточностей.

Чтобы получить векторный файл точек после визуальной интерпретации, вы можете использовать python, чтобы открыть файл dbf, который поставляется с файлом shp (фактически, это таблица атрибутов), чтобы вычислить матрицу путаницы и коэффициент каппа.

код показать, как показано ниже

from sklearn.metrics import confusion_matrix from sklearn.metrics import accuracy_score from sklearn.metrics import cohen_kappa_score from sklearn.metrics import precision_recall_fscore_support from sklearn.metrics import classification_report from sklearn.metrics import precision_score from sklearn.metrics import recall_score import pandas import os from dbfread import DBF rootdir = "H:\\result_2_type\\points" list = os.listdir(rootdir) colnames = ['OID_', 'RASTERVALU', 'true'] true=[] result=[] for i in range(0,len(list)): if not list[i].endswith(".dbf"): continue folder = os.path.join(rootdir,list[i]) print(folder) table = DBF(folder, encoding='GBK') df = pandas.DataFrame(iter(table)) result_temp = df['RASTERVALU'].tolist() true_temp = df['true'].tolist() score=accuracy_score(true_temp, result_temp) print('score=',score) kappa=cohen_kappa_score(true_temp, result_temp) print('kappa=',kappa) report=classification_report(true_temp, result_temp,digits=4) print('report=',report) true.extend(true_temp) result.extend(result_temp)

 



【本文地址】

公司简介

联系我们

今日新闻

    推荐新闻

    专题文章
      CopyRight 2018-2019 实验室设备网 版权所有