Netcdf clip using lat and longitude import xarray as xr ds =xr.open_dataset("") lat_bnds,lon_bnds = [20.2,25.2] , [81.2,85.5] ds1=ds.sel()
Posts
Showing posts with the label netcdf
- Get link
- X
- Other Apps
NetCDF to PostGIS Script Script # Import modules import psycopg2, time, datetime from scipy.io import netcdf # Establish connection db1 = psycopg2.connect("host=xxxx dbname=priogrid user=xxxx password=xxxx") cur = db1.cursor() # Create Table in postgis print str(time.ctime())+ " Creating precipud19002008 table." cur.execute("DROP TABLE IF EXISTS precipud19002008;") cur.execute("CREATE TABLE precipud19002008 (gid serial PRIMARY KEY not null, year int, month int, lon decimal, lat decimal, prec decimal);") # Read netcdf file using the netcdf_file class. 'r' for read. f = netcdf.netcdf_file('/mnt/pc258/prio_grid/source/ClimateData/precipud19002008.nc', 'r') # Create lathash. Read the f.variables['lat'].data.tolist() creates a list of lat coodinates. print str(time.ctime())+ " Looping through lat coords." temp = f.variables['lat'].data.tolist() # Create list of lat coords lathash = {} f...