[GRASS-user] trouble with r.in.gdal to georeference a NetCDF file

Lee Eddington lee.w.eddington at gmail.com
Tue Nov 12 06:43:58 PST 2013 

OK, here it is:

Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Driver: netCDF/Network Common Data Format
Files: none associated
Size is 198, 153
Coordinate System is `'
Metadata:
  HGT_M#description=Topography height
  HGT_M#FieldType=104
  HGT_M#MemoryOrder=XY
  HGT_M#sr_x=1
  HGT_M#sr_y=1
  HGT_M#stagger=M
  HGT_M#units=meters MSL
  NC_GLOBAL#BOTTOM-TOP_GRID_DIMENSION=0
  NC_GLOBAL#CEN_LAT=-8.4095154
  NC_GLOBAL#CEN_LON=115.02
  NC_GLOBAL#corner_lats={ -10.454437, -6.3537445, -6.3537445, -10.454437,
-10.454437, -6.3537445, -6.3537445, -10.454437, -10.46785, -6.3401871,
-6.3401871, -10.46785, -10.46785, -6.3401871, -6.3401871, -10.46785 }
  NC_GLOBAL#corner_lons={ 112.3332, 112.3332, 117.70679, 117.70679,
112.31956, 112.31956, 117.72044, 117.72044, 112.3332, 112.3332, 117.70679,
117.70679, 112.31956, 112.31956, 117.72044, 117.72044 }
  NC_GLOBAL#DX=3000
  NC_GLOBAL#DY=3000
  NC_GLOBAL#DYN_OPT=2
  NC_GLOBAL#FLAG_MF_XY=1
  NC_GLOBAL#grid_id=3
  NC_GLOBAL#GRIDTYPE=C
  NC_GLOBAL#i_parent_end=99
  NC_GLOBAL#i_parent_start=34
  NC_GLOBAL#ISICE=24
  NC_GLOBAL#ISLAKE=-1
  NC_GLOBAL#ISOILWATER=14
  NC_GLOBAL#ISURBAN=1
  NC_GLOBAL#ISWATER=16
  NC_GLOBAL#j_parent_end=84
  NC_GLOBAL#j_parent_start=34
  NC_GLOBAL#MAP_PROJ=3
  NC_GLOBAL#MMINLU=USGS
  NC_GLOBAL#MOAD_CEN_LAT=-8.4095078
  NC_GLOBAL#NUM_LAND_CAT=24
  NC_GLOBAL#parent_grid_ratio=3
  NC_GLOBAL#parent_id=2
  NC_GLOBAL#POLE_LAT=90
  NC_GLOBAL#POLE_LON=0
  NC_GLOBAL#SIMULATION_START_DATE=0000-00-00_00:00:00
  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#sr_x=1
  NC_GLOBAL#sr_y=1
  NC_GLOBAL#STAND_LON=115.02
  NC_GLOBAL#TITLE=OUTPUT FROM GEOGRID V3.4
  NC_GLOBAL#TRUELAT1=-8.4095001
  NC_GLOBAL#TRUELAT2=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
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  153.0)
Upper Right (  198.0,    0.0)
Lower Right (  198.0,  153.0)
Center      (   99.0,   76.5)
Band 1 Block=198x1 Type=Float32, ColorInterp=Undefined
  NoData Value=9.96920996838686905e+36
  Metadata:
    description=Topography height
    FieldType=104
    MemoryOrder=XY
    NETCDF_DIMENSION_Time=1
    NETCDF_VARNAME=HGT_M
    sr_x=1
    sr_y=1
    stagger=M
    units=meters MSL



On Tue, Nov 12, 2013 at 5:44 AM, Anna Petrášová <kratochanna at gmail.com>wrote:

>
>
>
> On Tue, Nov 12, 2013 at 8:38 AM, Lee Eddington <lee.w.eddington at gmail.com>wrote:
>
>>  SUBDATASET_26_NAME=NETCDF:"geo_em.d03.nc":HGT_M
>>  SUBDATASET_26_DESC=[1x153x198] HGT_M (32-bit floating-point)
>>
>>
> no, I mean when you run this command (in the directory with the netcdf
> file):
> gdalinfo NETCDF:geo_em.d03.nc:HGT_M
>
>
>
>
>>
>> On Tue, Nov 12, 2013 at 5:28 AM, Anna Petrášová <kratochanna at gmail.com>wrote:
>>
>>>
>>>
>>>
>>> On Tue, Nov 12, 2013 at 12:55 AM, Lee Eddington <
>>> lee.w.eddington at gmail.com> wrote:
>>>
>>>> Anna,
>>>>
>>>> I tried your suggestion, but now I'm unable to import the data.  Here's
>>>> some output showing the region and the error messages:
>>>>
>>>
>>>
>>> What does the gdalinfo says about the subdataset?
>>> gdalinfo NETCDF:geo_em.d03.nc:HGT_M
>>>
>>>
>>>
>>>> GRASS 6.4.3 (bali_d03_ll):~ > g.proj -p
>>>>
>>>> -PROJ_INFO-------------------------------------------------
>>>>
>>>> name       : Lat/Lon
>>>>
>>>> proj       : ll
>>>>
>>>> datum      : wgs84
>>>>
>>>> ellps      : wgs84
>>>>
>>>> no_defs    : defined
>>>>
>>>> towgs84    : 0.000,0.000,0.000
>>>>
>>>> -PROJ_UNITS------------------------------------------------
>>>>
>>>> unit       : degree
>>>>
>>>> units      : degrees
>>>>
>>>> meters     : 1.0
>>>>
>>>> GRASS 6.4.3 (bali_d03_ll):~ > g.region -p
>>>>
>>>> projection: 3 (Latitude-Longitude)
>>>>
>>>> zone:       0
>>>>
>>>> datum:      wgs84
>>>>
>>>> ellipsoid:  wgs84
>>>>
>>>> north:      6:19:40.44S
>>>>
>>>> south:      10:27:59.04S
>>>>
>>>> west:       112:19:46.2E
>>>>
>>>> east:       117:43:20.64E
>>>>
>>>> nsres:      0:01:36.12
>>>>
>>>> ewres:      0:01:37.56
>>>>
>>>> rows:       155
>>>>
>>>> cols:       199
>>>>
>>>> cells:      30845
>>>>
>>>> GRASS 6.4.3 (bali_d03_ll):~ > r.in.gdal -o
>>>> input="NETCDF:grassdata/bali/geo_em.d03.nc:HGT_M" output="d03_HGT_M"
>>>>
>>>> Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
>>>>
>>>> WARNING: Over-riding projection check
>>>>
>>>> WARNING: G_set_window(): Illegal latitude for North
>>>>
>>>>
>>>> On Mon, Nov 11, 2013 at 9:21 PM, Anna Petrášová <kratochanna at gmail.com>wrote:
>>>>
>>>>> Hi,
>>>>>
>>>>> On Mon, Nov 11, 2013 at 11:37 PM, Lee Eddington <
>>>>> lee.w.eddington at gmail.com> wrote:
>>>>>
>>>>>> I'm trying to import a NetCDF file from the WRF weather forecast
>>>>>> model with r.in.gdal, but I can't get it to georeference.  Instead I end up
>>>>>> with a x,y coordinate system of rows and columns.  The file georeferences
>>>>>> fine in a number of other weather graphics programs.  I tried to follow the
>>>>>> directions at:
>>>>>>
>>>>>> http://www.gdal.org/frmt_netcdf.html
>>>>>>
>>>>>> but when I run gdalinfo I get the following:
>>>>>>
>>>>>> $ gdalinfo geo_em.d03.nc
>>>>>>
>>>>>> Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
>>>>>>
>>>>>> Driver: netCDF/Network Common Data Format
>>>>>>
>>>>>> Files: geo_em.d03.nc
>>>>>>
>>>>>> Size is 512, 512
>>>>>>
>>>>>> Coordinate System is `'
>>>>>>
>>>>>> Metadata:
>>>>>>
>>>>>>   NC_GLOBAL#BOTTOM-TOP_GRID_DIMENSION=0
>>>>>>
>>>>>>   NC_GLOBAL#CEN_LAT=-8.4095154
>>>>>>
>>>>>>   NC_GLOBAL#CEN_LON=115.02
>>>>>>
>>>>>>   NC_GLOBAL#corner_lats={ -10.454437, -6.3537445, -6.3537445,
>>>>>> -10.454437, -10.454437, -6.3537445, -6.3537445, -10.454437, -10.46785,
>>>>>> -6.3401871, -6.3401871, -10.46785, -10.46785, -6.3401871, -6.3401871,
>>>>>> -10.46785 }
>>>>>>
>>>>>>   NC_GLOBAL#corner_lons={ 112.3332, 112.3332, 117.70679, 117.70679,
>>>>>> 112.31956, 112.31956, 117.72044, 117.72044, 112.3332, 112.3332, 117.70679,
>>>>>> 117.70679, 112.31956, 112.31956, 117.72044, 117.72044 }
>>>>>>
>>>>>>   NC_GLOBAL#DX=3000
>>>>>>
>>>>>>   NC_GLOBAL#DY=3000
>>>>>>
>>>>>>   NC_GLOBAL#DYN_OPT=2
>>>>>>
>>>>>>   NC_GLOBAL#FLAG_MF_XY=1
>>>>>>
>>>>>>   NC_GLOBAL#grid_id=3
>>>>>>
>>>>>>   NC_GLOBAL#GRIDTYPE=C
>>>>>>
>>>>>>   NC_GLOBAL#i_parent_end=99
>>>>>>
>>>>>>   NC_GLOBAL#i_parent_start=34
>>>>>>
>>>>>>   NC_GLOBAL#ISICE=24
>>>>>>
>>>>>>   NC_GLOBAL#ISLAKE=-1
>>>>>>
>>>>>>   NC_GLOBAL#ISOILWATER=14
>>>>>>
>>>>>>   NC_GLOBAL#ISURBAN=1
>>>>>>
>>>>>>   NC_GLOBAL#ISWATER=16
>>>>>>
>>>>>>   NC_GLOBAL#j_parent_end=84
>>>>>>
>>>>>>   NC_GLOBAL#j_parent_start=34
>>>>>>
>>>>>>   NC_GLOBAL#MAP_PROJ=3
>>>>>>
>>>>>>   NC_GLOBAL#MMINLU=USGS
>>>>>>
>>>>>>   NC_GLOBAL#MOAD_CEN_LAT=-8.4095078
>>>>>>
>>>>>>   NC_GLOBAL#NUM_LAND_CAT=24
>>>>>>
>>>>>>   NC_GLOBAL#parent_grid_ratio=3
>>>>>>
>>>>>>   NC_GLOBAL#parent_id=2
>>>>>>
>>>>>>   NC_GLOBAL#POLE_LAT=90
>>>>>>
>>>>>>   NC_GLOBAL#POLE_LON=0
>>>>>>
>>>>>>   NC_GLOBAL#SIMULATION_START_DATE=0000-00-00_00:00:00
>>>>>>
>>>>>>   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#sr_x=1
>>>>>>
>>>>>>   NC_GLOBAL#sr_y=1
>>>>>>
>>>>>>   NC_GLOBAL#STAND_LON=115.02
>>>>>>
>>>>>>   NC_GLOBAL#TITLE=OUTPUT FROM GEOGRID V3.4
>>>>>>
>>>>>>   NC_GLOBAL#TRUELAT1=-8.4095001
>>>>>>
>>>>>>   NC_GLOBAL#TRUELAT2=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:"geo_em.d03.nc":Times
>>>>>>
>>>>>>   SUBDATASET_1_DESC=[1x19] Times (8-bit character)
>>>>>>
>>>>>>   SUBDATASET_2_NAME=NETCDF:"geo_em.d03.nc":XLAT_M
>>>>>>
>>>>>>   SUBDATASET_2_DESC=[1x153x198] XLAT_M (32-bit floating-point)
>>>>>>
>>>>>>   SUBDATASET_3_NAME=NETCDF:"geo_em.d03.nc":XLONG_M
>>>>>>
>>>>>>   SUBDATASET_3_DESC=[1x153x198] XLONG_M (32-bit floating-point)
>>>>>>
>>>>>>   SUBDATASET_4_NAME=NETCDF:"geo_em.d03.nc":XLAT_V
>>>>>> .
>>>>>> .
>>>>>> .
>>>>>> .
>>>>>>
>>>>>> 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)
>>>>>>
>>>>>> The warning -
>>>>>>
>>>>>> Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
>>>>>>
>>>>>> leads me to believe that all the projection info is not being
>>>>>> recognized, but I'm not sure why.  Is it possible that GDAL wasn't compiled
>>>>>> with libnetcdf?  But if that's the case I would think that gdalinfo
>>>>>> wouldn't work at all.
>>>>>>
>>>>>> I've imported the data into GRASS using:
>>>>>>
>>>>>> r.in.gdal input="NETCDF:geo_em.d03.nc:HGT_M" output="d03_HGT_M"
>>>>>> location="bali_d03"
>>>>>>
>>>>>> to create a new location.
>>>>>>
>>>>>> I also created a location with the Mercator projection of the data
>>>>>> and tried importing using:
>>>>>>
>>>>>> r.in.gdal -o input="NETCDF:grassdata/bali/geo_em.d03.nc:HGT_M"
>>>>>> output="d03_HGT_M"
>>>>>>
>>>>>> but as far as I can tell still get an x,y row,column coordinate
>>>>>> system.
>>>>>>
>>>>>> I also have the lat and lon of each cell.  Is there anyway to
>>>>>> georeference using that data?
>>>>>>
>>>>>>
>>>>>
>>>>> try to create a WGS84 latitude-longitude projection (EPSG 4326) and
>>>>> then run r.in.gdal with -o flag. At least that worked for me and I had
>>>>> similar gdalinfo output.
>>>>>
>>>>> Anna
>>>>>
>>>>>
>>>>>
>>>>>>  Thanks,
>>>>>> Lee
>>>>>>
>>>>>> _______________________________________________
>>>>>> grass-user mailing list
>>>>>> grass-user at lists.osgeo.org
>>>>>> http://lists.osgeo.org/mailman/listinfo/grass-user
>>>>>>
>>>>>
>>>>>
>>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20131112/1c92e42a/attachment-0001.html>

More information about the grass-user mailing list