Координаты углов клеток в очень большом растре

[Odd-Bird]

Доброй ночи всем! Я пытаюсь решить очень сложную задачу, а именно - извлечь координаты углов клеток в растре GeoTiff. Проблема, собственно, состоит только в том, что растр этот очень большой; конечно, его можно разбить на части, преобразовать их в полигоны функцией rasterToPolygon(){raster} и извлечь из них координаты, но это займет чрезвычайно много времени.

Если кто-нибудь сталкивался с такой задачей - посоветуйте, пожалуйста: стоит ли пытаться решить ее в R или лучше воспользоваться другими программами, и какими именно?..

[Максим Дубинин]

Ваша очень сложная задача не выглядит как очень сложная

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

Ну и непонятно в каком виде нужен результат.

[Odd-Bird]

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

Ну и непонятно в каком виде нужен результат.

Я построил простейший цикл, он обсчитывает координаты со скоростью пикселя в секунду; 8 миллиардов пикселей (300-метровая карта мира) будут обсчитываться примерно 130 лет .

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

[Odd-Bird]

Ваша очень сложная задача не выглядит как очень сложная

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

Ну и непонятно в каком виде нужен результат.

Да, проблема еще в том, что пиксели на сфере не идеально квадратные, и их размер варьирует в разных широтах.

[Александр Мурый]

Распишите подробнее, что за растр, какая система координат, какие именно "углы клеток" надо найти?

[sergsh]

8 миллиардов пикселей будут иметь 32 миллиарда координат углов ...

а может проще все же вычислять каждый раз 4 координаты для нужного пикселя, и не хранить их все ?

[Odd-Bird]

Растр Land Cover Maps отсюда: http://maps.elie.ucl.ac.be/CCI/viewer/download.php (GeoTiff ~4 Гб)

Система координат longlat, эллипсоид wgs84.

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

[Odd-Bird]

а может проще все же вычислять каждый раз 4 координаты для нужного пикселя, и не хранить их все ?

Пиксели неравного размера, вычислять их координаты не получится - только считывать.

[Odd-Bird]

В общем, я решил ограничиться координатами центров пикселей.

[sergsh]

Растр Land Cover Maps отсюда: http://maps.elie.ucl.ac.be/CCI/viewer/download.php (GeoTiff ~4 Гб)

Система координат longlat, эллипсоид wgs84.

Спасибо за ссылку, интересные данные.

[gamm]

Проблема, собственно, состоит только в том, что растр этот очень большой;

проблема собственно в том, что никакие координаты ячеек вам хранить не нужно, поскольку растр и придуман для того, чтобы хранить только два 2-3 угла, и шаг сетки, а все остальное вычислять 4 действиями арифметики вместо SQL запросов. Если растр у вас в lon/lat, то достаточно просто перевести координаты в lon/lat, и вычислить номер ячейки/ячеек, куда они попадают. Т.е. от всей сетки вам нужно хранить 8 чисел, независимо от размера растра ... ну, может еще строку проекции. Пересечения ячеек с полигонами тоже быстро вычисляются без запросов (заметанием), с линиями - та же картина.