Перевести растр в точки с их действительными координатами

[Variant]

А теперь делаем таблицам обычный sql join по столбцу с координатами. И получаем шиш с маслом (ну то есть несколько случайных точек, центры которых совпали) потому что исходные растры - в разных проекциях. Чтобы join получился, нужно перепроецировать оба растра в одну проекцию, сделать ресемплинг для уравнивания разрешений, выровнять центры пикселей этих растров и только тогда конвертить их оба в таблицы, которые можно будет склеить по координатам.

Про sql join и речи не было. Только честный ресэмплинг. Либо по растровому пути, как предложен выше, либо перейдя сначала в облако точек и сформировав итоговый "растр" уже на конечной стадии. И так, и так, немного морочный путь. Вот мне и интересно, нет ли готовых инструментов, позволяющих автоматизировать такие манипуляции. Или вся возможная автоматизация здесь - это только через unix shell?

[gamm]

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

у вас странные представления о ГИС. Я бы даже сказал, никакие.

Многофакторные модели - стандартная вещь в ГИС, как и подключение современных статистических систем для работы с ними (S-PLUS в ArcInfo или ее полный открытый клон R в любых системах).

То, что вы описываете, является не структурированным описанием проблемы, а самой общей типовой постановкой задачи, типичной для геопроцессинга в ГИС, причем описание дается на уровне школьного реферата.

И уж совсем непонятно, что такое процедуру многофакторного распознавания, особенно при отсутствии исходников. Потому что все известные процедуры многофакторного распознавания в R есть, и какое отношение имеет к этому отсутствие исходников, вообще непонятно, reverse engineering - это о другом.

Есть ощущение, что все, что вам нужно - это обычный сэмплер, снимающий данные со слоев в заданных точках. Он есть в любой ГИС.

[Variant]

У GDAL есть формат XYZ, не что иное как две координаты (центр точки) и значение в ней

В любом случае, спасибо за идею с XYZ в GDAL, о такой возможности я не знал. Начну, видимо именно с неё, прогоню всё хотя бы руками, а уж как автоматизировать - там видно будет

[Variant]

у вас странные представления о ГИС. Я бы даже сказал, никакие.

Здесь я позволю себе не согласиться, какие-то представления уж точно есть, я гарантирую это. А странность - понятие относительное. То, что я видел что-то кроме ГИС, считаю растр абстракцией, и пытаюсь излагать свои мысли "языком школьного реферата" - чтобы они были общедоступны, ещё не означает, что я не знаком с рассматриваемым предметом

То, что вы описываете, является не структурированным описанием проблемы, а самой общей типовой постановкой задачи, типичной для геопроцессинга в ГИС

Я постарался описать всё как можно более структурированно и доступно. Но если даже так не нравится, то я уж даже и не знаю как ещё написать

S-PLUS в ArcInfo или ее полный открытый клон R в любых системах

Процедура, реализованная на некотором специализированном языке и среда для анализа это всё же немного разные вещи. Я же говорю не про какую-то одну уникальную процедуру, а именно про среду. Её и не всякий опытный программист возьмётся портировать, даже имея исходники, уж не говоря про стороннего пользователя

Потому что все известные процедуры многофакторного распознавания в R есть

Довольно спорное утверждение. Что-то мне подсказывает, что, скорее всего, не все

Есть ощущение, что все, что вам нужно - это обычный сэмплер, снимающий данные со слоев в заданных точках. Он есть в любой ГИС.

Близко, но не совсем. Нужны не данные со слоёв в точках, а точки, попадающие в определённые пространственные ячейки. Причём, как раз скорее это можно было бы назвать сэмплингом. Только вот действительно ли они есть именно во всех ГИС?

[gamm]

Вы так и не смогли сформулировать задачу. Скорее всего, вам нужны элементарные манипуляции с пространственными данными. Но об этом можно только догадываться, нам это не удалось, Вы победили.

Успехов!

[ericsson]

Ну вот охота человеку через по-своему делать, пусть делает... может пользу какую народному хозяйству принесет?
Автоматизация через шелл, надо просто атомарные операции описать один раз, а потом их комбинировать.

[gamm]

да я только за, пусть шаманит

Может, велосипеды еще не все изобретены.

[Variant]

Я знаю, что способ обмена данными, да и вообще построения всей системы, которым пока приходится пользоваться - неправильный. И некрасивый, и неудобный. Как костыль для лошади. Но пока что, я не могу ничего изменить.

Спасибо за помощь!