<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>