Output raster to netcdf with decimal degrees - ArcGIS 10

MattiasCape — Wed, 05 Dec 2012 16:02:52 GMT

I'm working with shapefiles downloaded from

ftp://sidads.colorado.edu/pub/DATASETS/MOA/coastlines/

which contain polygons and lines outlining certain Antarctic features (coastline, groundingline, and islands). To be able to use these in another software program (WIM) that does not handle shapefiles I converted these features to raster (e.g. using the polyline to raster tool) and then to netcdf (using raster to netcdf). While I have been able to open the final netcdf file, the coordinates listed in the file are given as x and y coordinates of projection in meters, as opposed to the decimal degrees I want them in. Unable to find latitude and longitude fields for the data, WIM forces me to define a geolocation file for the netcdf file I created, but I do not have one.

What changes should I make to my processing to change the units of my final netcdf file to decimal degrees? Alternatively, how can I create another netcdf file containing the latitude and longitude coordinates of the points in my raster (lon lat arrays if you wish)?

The header for the netcdf file I create appears below. Thanks in advance for your help.

netcdf file:/Volumes/lemaire/mattias/shapefiles/moa_coastline_PolylineToRast.nc {
dimensions:
   longitude = 2627;
   latitude = 2193;
variables:
   short moa_coastline_PolylineToRast(latitude=2193, longitude=2627);
     :long_name = "moa_coastline_PolylineToRast";
     :esri_pe_string = "PROJCS[\"Polar_Stereographic_South\",GEOGCS[\"GCS_WGS_1984\",DATUM[\"D_WGS_1984\",SPHEROID[\"WGS_1984\",6378137.0,298.257223563]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Stereographic_South_Pole\"],PARAMETER[\"False_Easting\",0.0],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",0.0],PARAMETER[\"Standard_Parallel_1\",-71.0],UNIT[\"Meter\",1.0]]";
     :coordinates = "longitude latitude";
     :grid_mapping = "polar_stereographic";
     :units = "Meter";
     :missing_value = 255S; // short
   int polar_stereographic;
     :grid_mapping_name = "polar_stereographic";
     :straight_vertical_longitude_from_pole = 0.0; // double
     :false_easting = 0.0; // double
     :false_northing = 0.0; // double
     :standard_parallel = -71.0; // double
     :latitude_of_projection_origin = -90.0; // double
     :_CoordinateTransformType = "Projection";
     :_CoordinateAxisTypes = "GeoX GeoY";
   double longitude(longitude=2627);
     :long_name = "x coordinate of projection";
     :standard_name = "projection_x_coordinate";
     :units = "Meter";
     :_CoordinateAxisType = "GeoX";
   double latitude(latitude=2193);
     :long_name = "y coordinate of projection";
     :standard_name = "projection_y_coordinate";
     :units = "Meter";
     :_CoordinateAxisType = "GeoY";

:Conventions = "CF-1.0";
:Source_Software = "Esri ArcGIS";
}

topic Output raster to netcdf with decimal degrees - ArcGIS 10 in Data Management Questions

Output raster to netcdf with decimal degrees - ArcGIS 10