Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
HDF to Geotiff используя Python API GDALКод: Выделить всё
from osgeo import gdal
filename = 'MOD07_L2.A2013036.0535.005.2013079062.hdf'
ds = gdal.Open(filename)
#Open output format driver, see gdal_translate --formats for list
format = "GTiff"
driver = gdal.GetDriverByName( format )
#Output to new format
dst_ds = driver.CreateCopy( 'test.tif', ds, 1 )
Код: Выделить всё
# -*- coding: utf-8 -*-
from osgeo import gdal
hdf = gdal.Open("MOD07_L2.A2004026.0110.004.2004026151642.hdf")
sdslist = hdf.GetSubDatasets()
names = []
for s in sdslist:
names.append(s[0])
sds = []
for n in names:
sds.append(gdal.Open(n))
sds[0].GetMetadata()
format = "GTiff"
driver = gdal.GetDriverByName( format )
dst_ds = driver.CreateCopy( "test.tif", sds[0], 1 )
Код: Выделить всё
ERROR 1: CPLMalloc(-482798304): Silly size requested.
ERROR 1: CPLMalloc(-482798304): Silly size requested.
ERROR 1: CPLMalloc(-1116485888): Silly size requested.
Ошибка сегментирования (core dumped)
для Модиса есть специальная тулза от производителя, лучше пользоваться ей (на сайте гислаба есть описание, ссылку не помню). Она знает, как устроены именно продукты MODIS, и позволяет выдавать трансформированные файлы.Rumato писал(а):В чём может быть проблема?)
Да, есть MRT/MRTSwath/HEG и прочие. Но вопрос был про Python+GDAL.gamm писал(а): для Модиса есть специальная тулза от производителя, лучше пользоваться ей (на сайте гислаба есть описание, ссылку не помню). Она знает, как устроены именно продукты MODIS, и позволяет выдавать трансформированные файлы.
Так просто не скажешь. Какая у вас версия GDAL? Какой именно код вы запускаете? Можете выложить ваш HDF-файл?Rumato писал(а):В чём может быть проблема?)
Код: Выделить всё
# -*- coding: utf-8 -*-
from osgeo import gdal
hdf = gdal.Open("MOD07_L2.A2013036.0535.005.2013036133243.hdf")
sdslist = hdf.GetSubDatasets()
names = []
for s in sdslist:
names.append(s[0])
sds = []
for n in names:
sds.append(gdal.Open(n))
sds[0].GetMetadata()
geotrans = gdal.GCPsToGeoTransform(sds[0].GetGCPs())
dst_name = 'test.tif'
format = "GTiff"
driver = gdal.GetDriverByName( format )
dst_ds = driver.CreateCopy( dst_name, sds[0], 1 )
dst_ds.SetGeoTransform( geotrans )
test.tif.zipКод: Выделить всё
for n in names:
... sds.append(gdal.Open(n))
...
ERROR 1: CPLMalloc(-375692544): Silly size requested.
ERROR 1: CPLMalloc(-375692544): Silly size requested.
ERROR 1: CPLMalloc(-1059362816): Silly size requested.
Ошибка сегментирования (сделан дамп памяти)
Код: Выделить всё
# -*- coding: utf-8 -*-
from osgeo import gdal
hdf = gdal.Open('MOD07_L2.A2013036.0535.005.2013079062316.hdf')
sdslist = hdf.GetSubDatasets()
names = []
for s in sdslist:
names.append(s[0])
sds = []
for n in names:
sds.append(gdal.Open(n))
out_names = [s.rsplit(':', 1)[1] for s in names]
for ind, name in enumerate(out_names):
geotrans = gdal.GCPsToGeoTransform(sds[ind].GetGCPs())
if geotrans:
dst_name = name + '.tif'
format = "GTiff"
driver = gdal.GetDriverByName( format )
dst_ds = driver.CreateCopy( dst_name, sds[ind], 1 )
dst_ds.SetGeoTransform( geotrans )
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 16 гостей
© GIS-Lab и авторы, 2002-2017. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).