Re: [netcdf-hdf] NetCDF to ShapeFile

Arnaud,

Take a look at the GDAL OGR library:

http://www.gdal.org/
http://www.gdal.org/ogr/

which has python bindings. This link has installation instructions and links to a python-based tutorial:

http://trac.osgeo.org/gdal/wiki/GdalOgrInPython

 -Hal

Arnaud Vandecasteele wrote:
Hi all,

I would like to know if it's possible to read data from a netcdf file and export it into a EsriShapefile File. The ShapeFile i a common file format in the GIS World. I'm trying to get the latitude, longitude and a determinate value of a netcdf file. But I don't know exactly how to do it. I succeed to open and read a netcdf file in python but i don't know how to export the data.
Here is my simple script :

import Scientific.IO.NetCDF as nc
from Numeric import *
import sys

try :
    ncFile = nc.NetCDFFile("tos_O1_2001-
2002.nc <http://2002.nc>","r")
except :
    print "can't open the file"
    sys.exit(1)

try :
    print "################# Dimensions #################"
    print ncFile.dimensions.keys()

    print "################# Variables #################"
    print ncFile.variables.keys()
#return ['time_bnds', 'lat_bnds', 'lon', 'lon_bnds', 'time', 'lat', 'tos'

    print "################# Var Dim #################"
    tos = ncFile.variables["tos"]
    print tos.dimensions
    #return : ('time', 'lat', 'lon')


    tosValue = tos.getValue()

except :
    ncFile.close()


I also tried to use the *ncdump* and *ncks* utilities and I've been able to get a TXT file of the data but i can't use directly to transform them in shapefile.

Do you know how I could do that?


Best regards


------------------------------------------------------------------------

_______________________________________________
netcdf-hdf mailing list
netcdf-hdf@xxxxxxxxxxxxxxxx
For list information or to unsubscribe, visit: http://www.unidata.ucar.edu/mailing_lists/