Tag: gdal

Как изменить разрешение растра с помощью GDAL?

Я ищу лучший способ изменить разрешение растрового набора данных GDAL. Например, у меня есть растровый размер пикселя (30, -30), и я хотел бы изменить размер пикселя на (5, -5), интерполируя все значения для данного пикселя в выходной растровый. Поэтому для каждого пикселя входного растра, я хотел бы иметь 36 пикселей в выходном растре, все они […]

Ошибка при установке GDAL

Я пытаюсь установить GDAL через pip. Но я получаю эту ошибку: extensions/gdal_wrap.cpp:3089:27: fatal error: cpl_vsi_error.h: No such file or directory #include "cpl_vsi_error.h" ^ compilation terminated. error: command 'x86_64-linux-gnu-gcc' failed with exit status 1 Я использовал эти команды: sudo apt-get install libgdal-dev export CPLUS_INCLUDE_PATH=/usr/include/gdal export C_INCLUDE_PATH=/usr/include/gdal pip install GDAL Может ли кто-нибудь сказать мне, как его […]

Установка GDAL на Mac OS X Lion

Я пытаюсь установить GDAL 1.7.1 на Mac OS X Lion, используя: python setup.py build python setup.py install и получить ошибку: running build running build_py running build_ext building 'osgeo._gdal' extension llvm-gcc-4.2 -fno-strict-aliasing -fno-common -dynamic -g -Os -pipe -fno-common -fno-strict-aliasing -fwrapv -mno-fused-madd -DENABLE_DTRACE -DMACOSX -DNDEBUG -Wall -Wstrict-prototypes -Wshorten-64-to-32 -DNDEBUG -g -fwrapv -Os -Wall -Wstrict-prototypes -DENABLE_DTRACE -pipe -arch […]

Эффективный способ ввода больших растровых данных в PyTables

Я ищу эффективный способ загрузить файл растровых данных (GeoTiff) размером 20 ГБ в PyTables для дальнейшего вычисления. В настоящее время я читаю его как массив numpy с использованием Gdal и записываю массив numpy в pytables, используя следующий код: import gdal, numpy as np, tables as tb inraster = gdal.Open('infile.tif').ReadAsArray().astype(np.float32) f = tb.openFile('myhdf.h5','w') dataset = f.createCArray(f.root, […]

Использование pip для установки модулей в python с ошибкой

У меня возникли проблемы с установкой модулей python с помощью pip. Ниже выведено окно команд: Обратите внимание, что я установил программу непосредственно перед установкой модуля GDAL. Я нахожусь на 64-битной машине w7 с python 2.7 Microsoft Windows [Version 6.1.7601] Copyright (c) 2009 Microsoft Corporation. All rights reserved. C:\Users\jnunn\Desktop>python get-pip.py Downloading/unpacking pip Downloading pip-1.2.1.tar.gz (102Kb): 102Kb […]

Как проецировать и перемасштабировать сетку для соответствия другой сетке с помощью python GDAL?

Уточнение: я как-то не учитывал ключевой аспект: не использовать os.system или subprocess – просто API-интерфейс python. Я пытаюсь преобразовать раздел сетки смещения NOAA GTX для вертикальных преобразований нулевой точки и не полностью следовать тому, как это сделать в GDAL с помощью python. Я хотел бы взять сетку (в этом случае атрибут Attributed Grid, но это […]

Массивный массив в Python

У меня есть базовые 2-D массивы numpy, и я бы хотел «уменьшить их» до более грубого разрешения. Есть ли простой модуль numpy или scipy, который может легко сделать это? Следует также отметить, что этот массив отображается географически через модули Basemap. ОБРАЗЕЦ:

Растеризация слоя GDAL

редактировать Вот правильный способ сделать это и документацию : import random from osgeo import gdal, ogr RASTERIZE_COLOR_FIELD = "__color__" def rasterize(pixel_size=25) # Open the data source orig_data_source = ogr.Open("test.shp") # Make a copy of the layer's data source because we'll need to # modify its attributes table source_ds = ogr.GetDriverByName("Memory").CopyDataSource( orig_data_source, "") source_layer = source_ds.GetLayer(0) […]

Работа с растрами в базе геоданных файлов (.gdb) с помощью GDAL

Я работаю над инструментом, который преобразует растровые слои в массивы для обработки с помощью NumPy, и в идеале я хотел бы работать с растрами, которые поставляются в формате .gdb, не экспортируя их все (особенно, если для этого требуется привлечение ArcGIS или ArcPy ). Возможно ли это с драйвером OpenFileGDB? Из того, что я могу сказать, […]

Сохранение большого цветного изображения как `GTiff` с помощью` gdal`

Я пытаюсь сохранить большое изображение размера (15000, 80000, 3) . Этот массив представляет собой массив numpy, который я инициализировал как im_final = np.zeros((15000,80000,,3)) . Чтобы сделать сохранение, я использую gdal следующим образом: dst_ds = gdal.GetDriverByName('GTiff').Create('val.tif', 80000, 15000, 3, gdal.GDT_Byte) dst_ds.GetRasterBand(1).WriteArray(im_final[:,:,0]) # write r-band to the raster dst_ds.GetRasterBand(2).WriteArray(im_final[:,:,1]) # write g-band to the raster dst_ds.GetRasterBand(3).WriteArray(im_final[:,:,2]) # […]

Python - лучший язык программирования в мире.