GIS-LAB

Географические информационные системы и дистанционное зондирование


Форумы GIS-Lab.info

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


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

Обсуждение материалов сайта: вопросы, замечания, предложения

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

Сообщение updates-bot » 01 апр 2010, 08:56

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

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

updates-bot
Bot
 
Зарегистрирован: 04 фев 2008
 
Сообщения: 276
Репутация: 2

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

Сообщение Aleksander » 14 июн 2010, 17:06

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

Aleksander
Участник
 
Зарегистрирован: 12 май 2009
Откуда: Минск (на карте)
 
Сообщения: 60
Репутация: 0

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

Сообщение Дмитрий Барышников » 14 июн 2010, 17:35

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

Аватара пользователя
Дмитрий Барышников
Гуру
 
Зарегистрирован: 17 ноя 2009
Откуда: Москва (на карте)
 
Сообщения: 2571
Репутация: 258
Статьи: 7
Проекты: 2

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

Сообщение Aleksander » 14 июн 2010, 17:51

Спасибо!
  • 0

Aleksander
Участник
 
Зарегистрирован: 12 май 2009
Откуда: Минск (на карте)
 
Сообщения: 60
Репутация: 0

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

Сообщение Aleksander » 14 июн 2010, 18:42

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

Aleksander
Участник
 
Зарегистрирован: 12 май 2009
Откуда: Минск (на карте)
 
Сообщения: 60
Репутация: 0

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

Сообщение ilmir35 » 06 июн 2013, 15:05

Всем привет. Я импортирую модули:
Код: Выделить всё
import gdal
import _gdal
from gdalconst import *
from _gdal import ptrcreate

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

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

ilmir35
Интересующийся
 
Зарегистрирован: 18 апр 2013
 
Сообщения: 26
Репутация: 0

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

Сообщение Александр Мурый » 06 июн 2013, 15:17

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

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

Редактор материалов, модератор форума
Александр Мурый
званий не имею
 
Зарегистрирован: 26 сен 2009
Откуда: Москва (на карте)
 
Сообщения: 4854
Репутация: 661
Статьи: 3
Проекты: 5/1

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

Сообщение ilmir35 » 06 июн 2013, 15:19

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

ilmir35
Интересующийся
 
Зарегистрирован: 18 апр 2013
 
Сообщения: 26
Репутация: 0

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

Сообщение ilmir35 » 06 июн 2013, 15:21

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

ilmir35
Интересующийся
 
Зарегистрирован: 18 апр 2013
 
Сообщения: 26
Репутация: 0

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

Сообщение Александр Мурый » 06 июн 2013, 15:24

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 года много воды утекло :)
  • 0

Редактор материалов, модератор форума
Александр Мурый
званий не имею
 
Зарегистрирован: 26 сен 2009
Откуда: Москва (на карте)
 
Сообщения: 4854
Репутация: 661
Статьи: 3
Проекты: 5/1

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

Сообщение ilmir35 » 06 июн 2013, 15:27

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

ilmir35
Интересующийся
 
Зарегистрирован: 18 апр 2013
 
Сообщения: 26
Репутация: 0

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

Сообщение ilmir35 » 06 июн 2013, 15:30

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

ilmir35
Интересующийся
 
Зарегистрирован: 18 апр 2013
 
Сообщения: 26
Репутация: 0

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

Сообщение Александр Мурый » 06 июн 2013, 15:32

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

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

Редактор материалов, модератор форума
Александр Мурый
званий не имею
 
Зарегистрирован: 26 сен 2009
Откуда: Москва (на карте)
 
Сообщения: 4854
Репутация: 661
Статьи: 3
Проекты: 5/1

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

Сообщение Александр Мурый » 06 июн 2013, 16:14

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

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

Редактор материалов, модератор форума
Александр Мурый
званий не имею
 
Зарегистрирован: 26 сен 2009
Откуда: Москва (на карте)
 
Сообщения: 4854
Репутация: 661
Статьи: 3
Проекты: 5/1

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

Сообщение ilmir35 » 06 июн 2013, 16:22

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

ilmir35
Интересующийся
 
Зарегистрирован: 18 апр 2013
 
Сообщения: 26
Репутация: 0

След.

Вернуться в Материалы сайта

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 2


(Геокруг)

© GIS-Lab и авторы, 2002-2013. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).