Форумы GIS-Lab.info

Natalia Novoselova писал(а): Нда.. понятно, что для открытия такого файла - нужен, как минимум, хороший базис знаний в математике и программировании. Так просто это не сделаешь.

bolotoved писал(а):Когда вам в руки попадает неведомый вам файл геоданных, первое, что вы должны сделать, это использовать утилиту GDAL для ознакомления с этим форматом и что содержит конкретный файл:

gdalinfo filename.nc

Denis Rykov писал(а): Libs -> gdal

trir писал(а):а ведь не надо было ничего устанавливать - через PowerShell это делается на раз!

C:\GIS_Doutorado_QS\Study>gdalinfo 2016121233557Forest_Aboveground_Biomassv3.nc
Driver: HDF5Image/HDF5 Dataset
Files: 2016121233557Forest_Aboveground_Biomassv3.nc
Size is 36000, 14002
Coordinate System is `'
Metadata:
  CDI=Climate Data Interface version 1.6.9 (http://mpimet.mpg.de/cdi)
  CDO=Climate Data Operators version 1.6.9 (http://mpimet.mpg.de/cdo)
  Conventions=CF-1.4
  created_by=R, packages ncdf and raster (version 2.2-31)
  date=2015-08-31 14:39:27
  history=Wed Sep 09 15:23:36 2015: cdo -f nc4c -z zip copy Forest_Aboveground_B
iomass_v3.nc Forest_Aboveground_Biomass_v3_2.nc
  _nc3_strict=1
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,14002.0)
Upper Right (36000.0,    0.0)
Lower Right (36000.0,14002.0)
Center      (18000.0, 7001.0)
Band 1 Block=36000x14002 Type=Float32, ColorInterp=Undefined
  Metadata:
    Forest_Aboveground_Biomass_v3_long_name=Forest_Aboveground_Biomass_v3
    Forest_Aboveground_Biomass_v3_max=611.399963378906
    Forest_Aboveground_Biomass_v3_min=0
    Forest_Aboveground_Biomass_v3_missing_value=-3.4e+038
    Forest_Aboveground_Biomass_v3_projection=+proj=longlat +datum=WGS84 +no_defs
 +ellps=WGS84 +towgs84=0,0,0
    Forest_Aboveground_Biomass_v3_projection_format=PROJ.4
    Forest_Aboveground_Biomass_v3__FillValue=-3.4e+038

NcDump
netcdf Forest_Aboveground_Biomassv3 { dimensions: longitude = 36000 ; latitude = 14002 ; variables: double longitude(longitude) ; longitude:standard_name = "longitude" ; longitude:long_name = "longitude" ; longitude:units = "degrees_east" ; longitude:axis = "X" ; double latitude(latitude) ; latitude:standard_name = "latitude" ; latitude:long_name = "latitude" ; latitude:units = "degrees_north" ; latitude:axis = "Y" ; float Forest_Aboveground_Biomass_v3(latitude, longitude) ; Forest_Aboveground_Biomass_v3:long_name = "Forest_Aboveground_Biomass_v3" ; Forest_Aboveground_Biomass_v3:_FillValue = -3.4e+38f ; Forest_Aboveground_Biomass_v3:missing_value = -3.4e+38f ; Forest_Aboveground_Biomass_v3:projection = "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0" ; Forest_Aboveground_Biomass_v3:projection_format = "PROJ.4" ; Forest_Aboveground_Biomass_v3:min = 0. ; Forest_Aboveground_Biomass_v3:max = 611.399963378906 ; // global attributes: :CDI = "Climate Data Interface version 1.6.9 (http://mpimet.mpg.de/cdi)" ; :history = "Wed Sep 09 15:23:36 2015: cdo -f nc4c -z zip copy Forest_Aboveground_Biomass_v3.nc Forest_Aboveground_Biomass_v3_2.nc" ; :Conventions = "CF-1.4" ; :created_by = "R, packages ncdf and raster (version 2.2-31)" ; :date = "2015-08-31 14:39:27" ; :CDO = "Climate Data Operators version 1.6.9 (http://mpimet.mpg.de/cdo)" ; }

В нашем случае, файл содержит всего один канал (переменная Band 1 в выводе gdalinfo). Далее нас интересует название переменной, которую мы хотим экстрагировать (NETCDF_VARNAME). В случае с температурой на высоте 2-х метров это NETCDF_VARNAME=T2MMEAN

bolotoved писал(а): Итак, экстрагируем конкретный слой c помощью утилиты gdalwarp:

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

gdalwarp -of GTiff -t_srs EPSG:4326 \
netCDF:"MERRA2_100.statD_2d_slv_Nx.19800102.nc4.nc":T2MMEAN output.tif

Где
-of GTiff -- формат выходного файла
-t_srs EPSG:4326 -- целевая система координат (в данном случае +proj=longlat +datum=WGS84 +no_defs)
далее идет конструкция тип_файла:название_файла:название_слоя
может быть и так: тип_файла:название_файла:название_набора:название_слоя (это все смотрим в выводе gdalinfo)
output.tif -- путь и название выходного файла TIF

nickleb писал(а):Здравствуйте. А, где взять этот Ваш, Natalia, конкретный nc-файл? Ну, так, чтобы по этажам-проектам институтов М.Planka не сбить лапки мышке Тоже повозиться, м.б. ... Поучаствовать в online-учениях...

:created_by = "R, packages ncdf and raster (version 2.2-31)" ;

nickleb писал(а):Без регистрации это сделать не удастся? - Так я понял...

nickleb писал(а):Почему-то, кажется, что, если Вы в R навыки имеете, то, задействовав package'и ncdf4 и raster, Вы прорвётесь, в geotiff'ы nc, конвертировав... Там же в R и rgdal есть...