[gdal-dev] WRF NetCDF import

[Lee Eddington]

I'm trying to use r.in.gdal in GRASS GIS to import Weather Research &
Forecasting (WRF) model NetCDF output.  There have been some tips provided
by the GRASS users list, but none have worked for me.  I can import the
data as a simple x,y array with no georeferencing, but information about
georeferencing is in the file as other programs read, georeference and
display the file correctly.

gdalinfo produces the following:

lees-mbp:full Lee$ gdalinfo test_full.nc
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Driver: netCDF/Network Common Data Format
Files: test_full.nc
Size is 512, 512
Coordinate System is `'
Metadata:
  NC_GLOBAL#BL_PBL_PHYSICS=1
  NC_GLOBAL#BOTTOM-TOP_GRID_DIMENSION=28
  NC_GLOBAL#BOTTOM-TOP_PATCH_END_STAG=28
  NC_GLOBAL#BOTTOM-TOP_PATCH_END_UNSTAG=27
  NC_GLOBAL#BOTTOM-TOP_PATCH_START_STAG=1
  NC_GLOBAL#BOTTOM-TOP_PATCH_START_UNSTAG=1
  NC_GLOBAL#BUCKET_J=-1
  NC_GLOBAL#BUCKET_MM=-1
  NC_GLOBAL#CEN_LAT=-8.4095154
  NC_GLOBAL#CEN_LON=115.02
  NC_GLOBAL#CU_PHYSICS=0
  NC_GLOBAL#DAMP_OPT=0
  NC_GLOBAL#DAMPCOEF=0.2
  NC_GLOBAL#DFI_OPT=0
  NC_GLOBAL#DIFF_6TH_FACTOR=0.12
  NC_GLOBAL#DIFF_6TH_OPT=0
  NC_GLOBAL#DIFF_OPT=1
  NC_GLOBAL#DT=16.666666
  NC_GLOBAL#DX=3000
  NC_GLOBAL#DY=3000
  NC_GLOBAL#FEEDBACK=0
  NC_GLOBAL#GFDDA_END_H=0
  NC_GLOBAL#GFDDA_INTERVAL_M=0
  NC_GLOBAL#GMT=12
  NC_GLOBAL#GRID_FDDA=0
  NC_GLOBAL#GRID_ID=3
  NC_GLOBAL#GRID_SFDDA=0
  NC_GLOBAL#GRIDTYPE=C
  NC_GLOBAL#HYPSOMETRIC_OPT=2
  NC_GLOBAL#I_PARENT_START=34
  NC_GLOBAL#ISFTCFLX=0
  NC_GLOBAL#ISHALLOW=0
  NC_GLOBAL#ISICE=24
  NC_GLOBAL#ISLAKE=-1
  NC_GLOBAL#ISOILWATER=14
  NC_GLOBAL#ISURBAN=1
  NC_GLOBAL#ISWATER=16
  NC_GLOBAL#J_PARENT_START=34
  NC_GLOBAL#JULDAY=324
  NC_GLOBAL#JULYR=2013
  NC_GLOBAL#KHDIF=0
  NC_GLOBAL#KM_OPT=4
  NC_GLOBAL#KVDIF=0
  NC_GLOBAL#MAP_PROJ=3
  NC_GLOBAL#MFSHCONV=0
  NC_GLOBAL#MMINLU=USGS
  NC_GLOBAL#MOAD_CEN_LAT=-8.4095078
  NC_GLOBAL#MOIST_ADV_OPT=1
  NC_GLOBAL#MP_PHYSICS=3
  NC_GLOBAL#NUM_LAND_CAT=24
  NC_GLOBAL#OBS_NUDGE_OPT=0
  NC_GLOBAL#OMLCALL=0
  NC_GLOBAL#PARENT_GRID_RATIO=3
  NC_GLOBAL#PARENT_ID=2
  NC_GLOBAL#POLE_LAT=90
  NC_GLOBAL#POLE_LON=0
  NC_GLOBAL#PREC_ACC_DT=0
  NC_GLOBAL#RA_LW_PHYSICS=1
  NC_GLOBAL#RA_SW_PHYSICS=1
  NC_GLOBAL#SCALAR_ADV_OPT=1
  NC_GLOBAL#SF_SFCLAY_PHYSICS=1
  NC_GLOBAL#SF_SURFACE_PHYSICS=2
  NC_GLOBAL#SF_URBAN_PHYSICS=0
  NC_GLOBAL#SGFDDA_END_H=0
  NC_GLOBAL#SGFDDA_INTERVAL_M=0
  NC_GLOBAL#SHCU_PHYSICS=0
  NC_GLOBAL#SIMULATION_START_DATE=2013-11-20_12:00:00
  NC_GLOBAL#SMOOTH_OPTION=0
  NC_GLOBAL#SOUTH-NORTH_GRID_DIMENSION=154
  NC_GLOBAL#SOUTH-NORTH_PATCH_END_STAG=154
  NC_GLOBAL#SOUTH-NORTH_PATCH_END_UNSTAG=153
  NC_GLOBAL#SOUTH-NORTH_PATCH_START_STAG=1
  NC_GLOBAL#SOUTH-NORTH_PATCH_START_UNSTAG=1
  NC_GLOBAL#SST_UPDATE=0
  NC_GLOBAL#STAND_LON=115.02
  NC_GLOBAL#START_DATE=2013-11-20_12:00:00
  NC_GLOBAL#SURFACE_INPUT_SOURCE=1
  NC_GLOBAL#SWRAD_SCAT=1
  NC_GLOBAL#TITLE= OUTPUT FROM WRF V3.4 MODEL
  NC_GLOBAL#TKE_ADV_OPT=1
  NC_GLOBAL#TRUELAT1=-8.4095001
  NC_GLOBAL#TRUELAT2=0
  NC_GLOBAL#W_DAMPING=0
  NC_GLOBAL#WEST-EAST_GRID_DIMENSION=199
  NC_GLOBAL#WEST-EAST_PATCH_END_STAG=199
  NC_GLOBAL#WEST-EAST_PATCH_END_UNSTAG=198
  NC_GLOBAL#WEST-EAST_PATCH_START_STAG=1
  NC_GLOBAL#WEST-EAST_PATCH_START_UNSTAG=1
Subdatasets:
  SUBDATASET_1_NAME=NETCDF:"test_full.nc":Times
  SUBDATASET_1_DESC=[1x19] Times (8-bit character)
  SUBDATASET_2_NAME=NETCDF:"test_full.nc":LU_INDEX
  SUBDATASET_2_DESC=[1x153x198] LU_INDEX (32-bit floating-point)

(more SUBDATASETs)

  SUBDATASET_106_NAME=NETCDF:"test_full.nc":LANDMASK
  SUBDATASET_106_DESC=[1x153x198] LANDMASK (32-bit floating-point)
  SUBDATASET_107_NAME=NETCDF:"test_full.nc":SST
  SUBDATASET_107_DESC=[1x153x198] SST (32-bit floating-point)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)

My first question is why am I getting the following warning?:

Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute

Next, I see what appears is georeferencing info in the following lines:

  NC_GLOBAL#CEN_LAT=-8.4095154
  NC_GLOBAL#CEN_LON=115.02
  NC_GLOBAL#MAP_PROJ=3 (Mercator)
  NC_GLOBAL#MOAD_CEN_LAT=-8.4095078
  NC_GLOBAL#STAND_LON=115.02
  NC_GLOBAL#TRUELAT1=-8.4095001

but obviously this is not being interpreted or used.

Trying to use gdalwarp:

lees-mbp:full Lee$ gdalwarp NETCDF:test_full.nc:HGT test_full_HGT.tif
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
ERROR 1: Unable to compute a transformation between pixel/line
and georeferenced coordinates for NETCDF:test_full.nc:HGT.
There is no affine transformation and no GCPs.

Can anyone explain what I need to do?

Thanks,
Lee
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20131123/a7153c47/attachment.html>