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
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
