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/
