Работа с растрами при помощи GDAL и Python

[updates-bot]

Обсуждение статьи "Работа с растрами при помощи GDAL и Python"

http://gis-lab.info/qa/gdal-python.html

[Aleksander]

Есть ли в составе gdal функции для уменьшения разрешения растра (не передискретизации с изменением разрешения, а именно уменьшения разрешения путем усреднения в пределах квадратов пикселей, например, 4*4) ?

[Дмитрий Барышников]

В самом GDAL в чистом виде нет.
Можно посмотреть http://www.gdal.org/gdal__alg_8h.html#4 ... 009078426f но цель его другая - апроксимация для уменьшения кол-ва вычислений.
И еще в Google нашел http://www.jenningsplanet.com/fileadmin ... lter3x3.py

[Aleksander]

Спасибо!

[Aleksander]

Нашел данную функцию в R (пакет raster, функция aggregate)

[ilmir35]

Всем привет. Я импортирую модули:

Код: Выделить всё

import gdal
import _gdal
from gdalconst import *
from _gdal import ptrcreate

Но питон ругается на ptrcreate:
AttributeError: 'module' object has no attribute 'ptrcreate'

Этот компонент нужно отдельно скачивать чтоли?

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

...питон ругается на ptrcreate:
AttributeError: 'module' object has no attribute 'ptrcreate'

Откуда вы взяли "_gdal" и "ptrcreate"? В API такого не находится.

[ilmir35]

Я много где видел использование _gdal, например, http://svn.osgeo.org/gdal/tags/gdal_1_2_0/pymod/gdal.py
И у меня его питон принимает.

[ilmir35]

И в вашей ссылке он есть. Если нажать на osgeo.gdal и прокрутить вниз.

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

Я много где видел использование _gdal, например, http://svn.osgeo.org/gdal/tags/gdal_1_2_0/pymod/gdal.py
И у меня его питон принимает.

А у меня не понимает.Вы видели дату последнего изменения в указанном вами скрипте?

# Revision 1.58 2004/02/25 09:04:33

Что-то мне кажется, что с 2004 года много воды утекло

[ilmir35]

У меня есть источник 2009 года, но там работают через FWTools. Возможно что эти компоненты туда встроены. Но я сомневаюсь что их нет в osgeo4w.

[ilmir35]

Вообще мне нужно реализовать вот эту функцию:
def GetBlockSize(band):
x = _gdal.ptrcreate('int', 0, 2)
_gdal.GDALGetBlockSize(band._o, x, _gdal.ptradd(x, 1))
result = (_gdal.ptrvalue(x, 0), _gdal.ptrvalue(x, 1))
_gdal.ptrfree(x)
return result

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

У меня есть источник 2009 года, но там работают через FWTools. Возможно что эти компоненты туда встроены. Но я сомневаюсь что их нет в osgeo4w.

Надо посмотреть, что запаковано в OSGeo4W. У меня GDAL 1.9.0 на Linux, там никаких "ptrcreate" нет.
Но "_gdal" импортируется, да.

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

Вообще мне нужно реализовать вот эту функцию:
def GetBlockSize(band):

Такая функция уже есть. Вам надо переписать её заново "с нуля"? В чём состоит задача?

[ilmir35]

Переписывать не надо. Она мне нужна для считывание полей из .img файла (get the band and block sizes).