Qgis не работает tin интерполяция

Интерполяция точечных данных¶

This tutorial is now obsolete. A new and updated version is available at Interpolating Point Data (QGIS3)

Интерполяция — это широко используемая методика ГИС для создания непрерывной поверхности из дискретных точек. Многие реальные явления в мире являются непрерывным: высоты, почвы, температуры и т. д. Если нам нужно смоделировать эти поверхности для анализа, невозможно провести измерения по всей поверхности. Следовательно, полевые измерения проводятся в различных точках поверхности, а промежуточные значения рассчитываются в процессе под названием „интерполяция“. В QGIS, интерполяция осуществляется с помощью встроенного модуля Интерполяция .

Обзор задачи¶

Мы возьмем полевые измерения глубин для озера Арлингтон в Техасе и создадим карту рельефа дна и изолинии из этих измерений.

Вы научитесь также¶

создавать изолинии из точечных данных;

скрывать пустые значения из растрового слоя;

добавлять подписи к векторному слою.

Получение данных¶

Texas Water Development Board provides the shapefiles for completed lake surveys.

For convenience, you can directly download the sample data used in this tutorial from link below.

Методика¶

Откройте QGIS. Перейдите к Layer ‣ Add Layer ‣ Add Vector Layer..

Перейдите к загруженному файлу Shapefiles.zip и выберите его. Нажмите Open .

In the Select layers to add… dialog, hold the Shift key and select Arlington_Soundings_2007_stpl83.shp and Boundary2004_550_stpl83.shp layers. Click OK .

You will see the 2 layers loaded in QGIS. The Boundary2004_550_stpl83 layer represents the boundary of the lake. Un-check the box next to it in the Table of Contents.

This will reveal the data from the second layer Arlington_Soundings_2007_stpl83 . Though the data looks like lines, it is a series of points that are very close.

Click the Zoom icon and select a small area on the screen. As you zoom closer, you will see the points. Each point represents a reading taken by a Depth Sounder at the location recorded by a DGPS equipment.

Select the Identify tool and click on a point. You will see the Identify Results panel show up on the left with the attribute value of the point. In this case, the ELEVATION attribute contains the depth of the lake at the location. As our task is to create a depth profile and elevation contours, we will use this values as input for the interpolation.

Make sure you have the Interpolation plugin enabled. See Использование модулей расширения for how to enable plugins. Once enabled, go to Raster ‣ Interpolation ‣ Interpolation .

In the Interpolation dialog, select Arlington_Soundings_2007_stpl83 as the Vector layers in the Input panel. Select ELEVATION as the Interpolation attribute . Click Add . Change the Cellsize X and Cellsize Y values to 5 . This value is the size of each pixel in the output grid. Since our source data is in a projected CRS with Feet-US as units, based on our selection, the grid size will be 5 feet. Click on the … button next to Output file and name the output file as elevation_tin.tif . CLick OK .

Interpolation results can vary significantly based on the method and parameters you choose. QGIS interpolation supports Triagulated Irregular Network (TIN) and Inverse Distance Weighting (IDW) methods for interpolation. TIN method is commonly used for elevation data whereas IDW method is used for interpolating other types of data such as mineral concentrations, populations etc. See the Spatial Analysis module of the QGIS documentation for more details.



You will see the new later elevation_tin loaded in QGIS. Right-click the layer and select Zoom to layer .

Now you will see the full extent of the created surface. Interpolation does not give accurate results outside the collection area. Let’s clip the resulting surface with the lake boundary. Go to Raster ‣ Extraction ‣ Clipper .

Name the Output file as elevation_tin_clipped.tif . Select the Cliiped mode as Mask layer . Select Boundary2004_550_stpl83 as the Mask layer` . Click OK .

A new raster elevation_tin_clipped will be loaded in QGIS. We will now style this layer to show the difference in elevations. Note the min and max elevation values from the elevation_tin layer. Right-click the elevation_tin_clipped layer and select Properties .

Go to the Style tab. Select Render type as Singleband pseudocolor . In the Generate new color map panel, select Spectral color ramp. As we want to create a depth-map as opposed to a height-map, check the Invert box. This will assign blues to deep areas and reds to shallow areas. Click Classify .

Switch to the Tranparency tab. We want to remove the black-pixels from our output. Enter 0 as the Additional no data value . Click OK .

Now you have a elevation relief map for the lake generated from the individual depth readings. Let’s generate contours now. Go to Raster ‣ Extraction ‣ Contours .

In the Contour dialog, enter contours as the Output file for contour lines . We will generate contour lines at 5ft intervals, so enter 5.00 as the Interval between contour lines . Check the Attribute name box. Click OK .

The contour lines will be loaded as contours layer once the processing is finished. Right-click the layer and select Properties .

Go to the Labels tab. Check the Label this layer with box and select ELEV as the field. Select Curved as the Placement type and click OK .

You will see that each contour line will be appropriately labeled with the elevation along the line.

© Copyright 2019, Ujaval Gandhi.
Обновлено: сент. 04, 2021.
Created using Sphinx 4.0.1.

Источник

Форумы GIS-Lab.info

Геоинформационные системы (ГИС) и Дистанционное зондирование Земли

Улучшение плагина интерполяции. TIN

Улучшение плагина интерполяции. TIN

Сообщение wanderer » 13 сен 2010, 16:49

что для хранения используется формат «ASCII grid (text format)» —
то есть по сути информация о рельефе храниться в текстовом виде, что очень неэффективно.
В чем, собственно и признался Marco. В решениях ESRI для этого используется бинарный растровый формат
«ASCII binary grid». Загвоздка заключается в том, что формат является закрытым и в библиотеке GDAL
доступен только для чтения (я так понимаю, что то как его читать получено эмперическим путем ).
Еще один важный момент, вытекающий из ответа Marco, — он расчитывает на помощь в «расширении»
возможностей плагина интерполяции.

Наконец, я перешел к предложениям и вопросам:
1) Предлагаю на уровне сообщества взяться за написание «патча» для сохранения TIN в каком-то
открытом формате аналоге «ASCII binary grid». Необходимо проанализировать существующие,
поддерживаемые GDAL и при необходимости подумать над созданием нового.
2) Эта тема также возможно относится к улучшению/модификации OGR/GDAL. Вопрос к завсегдатаям
форума (я пока новичок): как в данной ситуации лучше поступить, чтобы тема была «видна» и из форума OGR/GDAL?
3) Вопрос к разработчикам QGIS — членам сообщества: изменилась ли ситуация с плагином интерполяции
с января 2010 (дата моей переписки с Marco)?

Источник

Форумы GIS-Lab.info

Геоинформационные системы (ГИС) и Дистанционное зондирование Земли

Интерполяция в Q gis

Интерполяция в Q gis

Сообщение ogorod » 09 авг 2018, 12:05

Re: Интерполяция в Q gis

Сообщение gamm » 09 авг 2018, 12:58

Re: Интерполяция в Q gis

Сообщение ogorod » 10 авг 2018, 08:47

Re: Интерполяция в Q gis

Сообщение Александр Мурый » 10 авг 2018, 15:49

Re: Интерполяция в Q gis

Сообщение ogorod » 13 авг 2018, 08:24

Из GPS координат сделал файл с точками. Открыл в Q Gis.
Потом этот файл соединил в exel таблицей в которой данные хим. анализа.

Пробовал и в ручную набивать данные, но результат тот же.
делать новый слой и ставить точки в ручную не пробовал.

Re: Интерполяция в Q gis

Сообщение konst555 » 13 авг 2018, 16:04

Re: Интерполяция в Q gis

Сообщение gamm » 13 авг 2018, 18:51

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

Сделайте все в одной проекции (в метрах), там же и интерполируйте, там же и смотрите (судя по смещению, у вас просто где-то глобусы разные).

Re: Интерполяция в Q gis

Сообщение chashin-an » 13 сен 2019, 11:48

Re: Интерполяция в Q gis

Сообщение gamm » 13 сен 2019, 14:32

Re: Интерполяция в Q gis

Сообщение chashin-an » 13 сен 2019, 17:00

Re: Интерполяция в Q gis

Сообщение gamm » 13 сен 2019, 18:52

Re: Интерполяция в Q gis

Сообщение chashin-an » 16 сен 2019, 15:11

Re: Интерполяция в Q gis

Сообщение gamm » 16 сен 2019, 18:17

Использовать лучше SRTM, но нужно его причесывать под гидрологию. И нужны не локальные параметры, а те, что связаны с переносом и потоками. Потом нужно оценивать потерю почвы и ее аккумуляцию. И рельеф нужен в пределах водосбора. Но если вы этим не занимались, то может лучше и не начинать, нужен опыт. Если делать, можно использовать WhiteBox (https://jblindsay.github.io/ghrg/Whitebox/), прожечь рельеф (breaching), и посчитать произведение L*S (http://www.iwr.msu.edu/rusle/factors.htm). Если будет корреляция с гумусом, то можно использовать для регрессии, еще нужно как-то оценить аккумуляцию.

Вместо рельефа можно использовать ДДЗЗ (http://www.50northspatial.org/soil-humu . sing-data/), например Sentinel, который бесплатный. Это проще и понятнее.

В обоих случаях нужно интерполировать остатки от регрессии

Напрямую, таком количестве точек, лучше использовать что-то консервативное, типа tension spline, поиграв напряжением, чтобы результат не задирался при экстраполяции.

Источник

TIN Interpolation plugin in QGIS 3.0

I just downloaded QGIS 3.0 on Mac OS X but I could not find TIN interpolation plugin in the Search & manage plugin. Is there something wrong with the Mac version of QGIS 3?

3 Answers 3

According to the changelog, the TIN interpolation plugin has been removed, as it is now part or QGIS core.

To find it in Qgis 3.0, open the processing toolbox and search for Interpolation -> TIN interpolation.

Some of the plugins are «deprecated» for this version of QGIS, they need to be upgraded/ moved within the QGIS 3.0 environment. Here is some info for Plugin Devs that you can read to maybe help you understand the concept «Ported Over»

Considering QGIS 3.0 is in early release, if you desparately need the functionality of that Plugin I recommend using QGIS 2 considering it is «Long Term Deployment» and the plugins are up to date.

You can try download a zipped version of the plugin from Github or the QGIS Plugin Repo from a web browser (Google the plugin name) and install it via the «Install From Zip» option in the plugin window, if it works, great, but it should then tell you if the plugin is compatible

Источник

Модуль интерполяции¶

Модуль интерполяции может использоваться для интерполяции точечного векторного слоя методом триангуляции (TIN — Triangular Irregular Network) или обратного взвешивания расстояний (IDW — Inverse Distance Weighted). Данная операция довольно несложная и основывается на интуитивно понятном графическом интерфейсе для создания интерполированных растровых слоев (cм. рисунок Figure_interpolation_1). Модуль требует наличия следующих параметров для выполнения:

Исходный векторный слой: Выберите исходный точечный векторый слой из списка загруженых точечных слоев. Если выбраны несколько слоев, для интерполяции используются данные всех слоев. Примечание: существует возможность вставки линий или полигонов в качестве ограничений для триангуляции; для этого необходимо выбрать «Линии структуры» или «Линии разбивки» в выпадающем списке Тип .

Атрибут интерполяции: Выберите необходимый атрибут для интерполяции или установите флаг Использовать для интерполяции Z-координату для того, чтобы задействовать значения Z, хранимые в слоях.

Метод интерполяции: Выберите метод интерполяции. Это может быть либо «Триангуляция (TIN)» или же «Обратное взвешивание расстояний (IDW)».

Количество столбцов/строк: Выберите количество строк и столбцов в результирующем растровом файле.

Файл вывода: Выберите название для выходного растрового файла.

Добавить результат в проект для загрузки результата на в проект.

Figure Interpolation 1:

Модуль интерполяции

Использование модуля¶

Запустите QGIS и загрузите точечный векторый слой (к примеру, elevp.csv ).

Активируйте модуль интерполяции в Менеджере модулей (см. раздел Загрузка основных модулей QGIS), а затем нажмите кнопку Интерполяция , которая появится на панели инструментов QGIS. Откроется диалоговое окно модуля интерполяции, как показано на рисунке Figure_interpolation_1.

Выберите исходный слой (к примеру, elevp ) и колонку (к примеру, ELEV ) для интерполяции.

Укажите метод интеполяции (например, «Триангуляция») и задайте Разрешение по Х и Разрешение по Y равными 5000, а также задайте название растрового файла вывода (например, elevation_tin ).

Источник