Форумы GIS-Lab.info

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

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

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

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

Спасибо!

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

Всем привет. Я импортирую модули: Но питон ругается на ptrcreate:
AttributeError: 'module' object has no attribute 'ptrcreate'

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

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

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

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

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

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

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

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

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

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

Вообще мне нужно реализовать вот эту функцию:
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

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

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

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

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

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