<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Markus,<br><br>I can read a raster from the file into GRASS using:<br><br>r.in.gdal input="NETCDF:test.nc:V10" output="V10_test3"<br><br>and I get the following:<br><br><font face="Courier New"><div>GRASS 6.4.3 (bali_d03):~ > r.info V10_test3</div><div> +----------------------------------------------------------------------------+</div><div> | Layer:    V10_test3                      Date: Sat Nov 23 20:56:28 2013    |</div><div> | Mapset:   PERMANENT                      Login of Creator: Lee             |</div><div> | Location: bali_d03                                                         |</div><div> | DataBase: /Users/Lee/grassdata                                             |</div><div> | Title:     ( V10_test3 )                                                   |</div><div> | Timestamp: none                                                            |</div><div> |----------------------------------------------------------------------------|</div><div> |                                                                            |</div><div> |   Type of Map:  raster               Number of Categories: 255             |</div><div> |   Data Type:    FCELL                                                      |</div><div> |   Rows:         153                                                        |</div><div> |   Columns:      198                                                        |</div><div> |   Total Cells:  30294                                                      |</div><div> |        Projection: x,y                                                     |</div><div> |            N:        153    S:          0   Res:     1                     |</div><div> |            E:        198    W:          0   Res:     1                     |</div><div> |   Range of data:    min = -12.64637  max = 7.255654                        |</div><div> |                                                                            |</div><div> |   Data Description:                                                        |</div><div> |    generated by r.in.gdal                                                  |</div><div> |                                                                            |</div><div> |   Comments:                                                                |</div><div> |    r.in.gdal input="NETCDF:test.nc:V10" output="V10_test3"                 |</div><div> |                                                                            |</div><div> +----------------------------------------------------------------------------+</div><div><br></div><div>gdalinfo on test.nc gives the following (Iíve stripped out a lot of the metadata, but left the info the describes the Mercator projection):</div><div><br></div><div><div>lees-mbp:bali Lee$ gdalinfo test.nc </div><div>Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute</div><div>Driver: netCDF/Network Common Data Format</div><div>Files: test.nc</div><div>Size is 512, 512</div><div>Coordinate System is `'</div><div>Metadata:</div><div>.</div><div>.</div><div>.</div><div>  NC_GLOBAL#CEN_LAT=-8.4095154</div><div>  NC_GLOBAL#CEN_LON=115.02</div><div>.</div></div><div>.</div><div><div>.</div><div>  NC_GLOBAL#MAP_PROJ=3</div><div>.</div><div>.</div><div>.</div><div>  NC_GLOBAL#STAND_LON=115.02</div><div>.</div><div>.</div><div>  NC_GLOBAL#TRUELAT1=-8.4095001</div><div>  NC_GLOBAL#TRUELAT2=0</div><div>.</div><div>.</div><div>.</div><div>Subdatasets:</div><div>  SUBDATASET_1_NAME=NETCDF:"test.nc":U10</div><div>  SUBDATASET_1_DESC=[1x153x198] U10 (32-bit floating-point)</div><div>  SUBDATASET_2_NAME=NETCDF:"test.nc":V10</div><div>  SUBDATASET_2_DESC=[1x153x198] V10 (32-bit floating-point)</div><div>Corner Coordinates:</div><div>Upper Left  (    0.0,    0.0)</div><div>Lower Left  (    0.0,  512.0)</div><div>Upper Right (  512.0,    0.0)</div><div>Lower Right (  512.0,  512.0)</div><div>Center      (  256.0,  256.0)</div></div><div><br></div></font>I can display the data in GRASS (see attachment).<img height="219" width="320" apple-width="yes" apple-height="yes" apple-inline="yes" id="022FF31D-BECF-46D1-89E8-3BB7F7F2CD45" src="cid:008EC921-189D-4E41-A7D7-7AA1280E0EA7@home"><div><br></div><div>I have the latitude and longitude values for every raster point from another file and can read them into collocated rasters.  I was able to interactively georectify the raster, but my goal is produce automated scripts and it doesnít appear that thereís a way to automate the georectify process (please correct me if Iím wrong).  I tried using g.setproj to change the projection as I know what it is, but it wouldnít allow me to because the current projection is x,y.</div><div><br></div><div>Thanks for sticking with me on this,</div><div><br></div><div>Lee</div><div><br><br>On Nov 23, 2013, at 7:43 AM, Markus Neteler <<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>> wrote:<br><br><blockquote type="cite">Lee,<br><br>I am not much familiar with non-map netCDF files. I do not know how to<br>deal with this specific file...<br>Let's hope that an answer from the community arrives. To me it doesn't<br>look like a file containing maps.<br><br>Best<br>Markus<br><br><br>On Sat, Nov 23, 2013 at 9:58 AM, Lee Eddington<br><<a href="mailto:lee.w.eddington@gmail.com">lee.w.eddington@gmail.com</a>> wrote:<br><blockquote type="cite">Markus,<br><br>I tried and got the following:<br><br>lees-mbp:bali Lee$ gdalwarp test.nc test.tif<br><br>Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute<br><br>Input file test.nc has no raster bands.<br><br><br>I've attached the file in case you'd be willing to see what you can do with<br>it.<br><br>Thanks,<br>Lee<br><br><br><br>On Sun, Nov 17, 2013 at 1:29 AM, Markus Neteler <<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>> wrote:<br><blockquote type="cite"><br>On Wed, Nov 13, 2013 at 11:58 PM, Lee Eddington<br><<a href="mailto:lee.w.eddington@gmail.com">lee.w.eddington@gmail.com</a>> wrote:<br><blockquote type="cite">Markus,<br><br>Could you provide me an example of using gdalwarp converting a NetCDF<br>field<br>to a GeoTIFF?<br></blockquote><br>I often due multiple runs here. Initially<br><br>gdalwarp thefile.nc thefile.tif<br><br>then I check if the target resolution is as expected (if not, see the<br>-tr parameter) and so on.<br>Here is an older blog page which may be of interest:<br><br><a href="http://blog.neteler.org/gdal-raster-data-tips-and-tricks/">http://blog.neteler.org/gdal-raster-data-tips-and-tricks/</a><br><br>Markus<br></blockquote><br><br></blockquote></blockquote><br></div></body></html>