Ошибка загрузки растра в NetCDF

[Максим Дубинин]

Попробуйте загрузить, пожалуйста в QGIS в такой растр:

http://eca.knmi.nl/download/ensembles/d ... v4.0.nc.gz

У меня Bad allocation, другие nc загружаются более-менее нормально.
Gdalinfo на растре тоже работает.

Win7, c9bb721

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

Подтверждаю ошибку:

2011-07-29-190617_1366x768_scrot.png (19.72 КБ) 5169 просмотров

Linux, QGIS 1.8.0+git20110725+33c766c

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

В общем, оно и не удивительно: gdalinfo в этом netCDF находит аж 22280 каналов. Т.е., видимо, надо либо вытаскивать данные поканально, либо как-то ещё.

Если вытащить первый канал в геотифф,

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

gdal_translate -of GTiff -a_srs "EPSG:4326" tn_0.50deg_reg_v4.0.nc tn_0.50deg_reg_v4.0_band1.tif -b 1

, то на выходе имеем такой растр:

Спойлер

Driver: GTiff/GeoTIFF
Files: tn_0.50deg_reg_v4.0_band1.tif
Size is 232, 101
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (-40.500000000000000,75.500000000000000)
Pixel Size = (0.500000000000000,-0.500000000000000)
Metadata:
NC_GLOBAL#Ensembles_ECAD=4.0
NC_GLOBAL#Conventions=CF-1.4
NC_GLOBAL#References=http://eca.knmi.nl\nhttp://eca.knmi.nl/download/ensembles/ensembles.php\nhttp://eca.knmi.nl/download/ensembles/Haylock_et_al_2007.pdf
tn#long_name=minimum temperature
tn#units=Celsius
tn#standard_name=air_temperature
tn#_FillValue=-9999
tn#scale_factor=1.000000e-02
longitude#long_name=Longitude values
longitude#units=degrees_E
longitude#standard_name=longitude
latitude#long_name=Latitude values
latitude#units=degrees_N
latitude#standard_name=latitude
time#long_name=Time in days
time#units=days since 1950-01-01 00:00
time#standard_name=time
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -40.5000000, 75.5000000) ( 40d30' 0.00"W, 75d30' 0.00"N)
Lower Left ( -40.5000000, 25.0000000) ( 40d30' 0.00"W, 25d 0' 0.00"N)
Upper Right ( 75.5000000, 75.5000000) ( 75d30' 0.00"E, 75d30' 0.00"N)
Lower Right ( 75.5000000, 25.0000000) ( 75d30' 0.00"E, 25d 0' 0.00"N)
Center ( 17.5000000, 50.2500000) ( 17d30' 0.00"E, 50d15' 0.00"N)
Band 1 Block=232x17 Type=Int16, ColorInterp=Gray
NoData Value=-9999
Offset: 0, Scale:0.00999999977648258
Metadata:
NETCDF_VARNAME=tn
NETCDF_DIMENSION_time=0
NETCDF_time_units=days since 1950-01-01 00:00

2011-07-29-194736_1366x768_scrot.png (57.23 КБ) 5165 просмотров

[Максим Дубинин]

amuriy, спасибо за проверку, извлекать то я умею тремя простыми и четырьмя сложными способами

теперь надо кугисик научить с такими растрами работать нормально.