netcdf

My aim is to access data from a netcdf file and write to a CSV file in the following format. Latitude Longitude Date1 Date2 Date3 100 200 <-- MIN_SFC values --> So far I have accessed the variables, written the header to the file and populated the lat/lons. How can I access the MIN_SFC values for specified lon,lat coordinates and dates and then write to a CSV file. I'm a python newbie if there is a better way to go about this please let me know. NetCDF file info: Dimensions: time = 7 latitude = 292 longitude =341 Variables: float MIN_SFC (time=7, latitude = 292, longitude = 341) Here's what I

I'd like to be able to find the lon/lat coordinate indices of the closest location to a lon/lat tuple. This is already available in the Java API as GridCoordSystem.findXYindexFromLatLon(), but I haven't found anything comparable in the Python API. What I'm hoping to find in the API, or write myself and contribute to the API (if useful), is something like this: def get_indices(netcdf_dataset, lon_value, lat_value): ''' :param netcdf_dataset an open NetCDF data set object :param lon_value a longitude value, in degrees (-180...180) :param lat_value a latitude value, in degrees (-90...90) :return

I have downloaded a nc file from f=open.ncdf("file.nc") [1] "file Lfile.nc has 2 dimensions:" [1] "Longitude Size: 1440" [1] "Latitude Size: 720" [1] "------------------------" [1] "file filr.nc has 8 variables:" [1] "short ts[Latitude,Longitude] Longname:Skin Temperature (2mm) Missval:NA" I then wanted to work with the variable soil_moisture_c A = get.var.ncdf(nc=f,varid="soil_moisture_c",verbose=TRUE) I then plot A with image(A) . I got the map shown below,I even transposed it image(t(a)) but that was changed to other direction,and not how it should be. Anyway,in order to know what is wrong

I have been trying plot the following gridded netcdf file: "air.1999.nc" found at the following website: http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.html I have tried the code below based on answers I have found here and elsewhere, but no luck. library(ncdf); temp.nc <- open.ncdf("air.1999.nc"); temp <- get.var.ncdf(temp.nc,"air"); temp.nc$dim$lon$vals -> lon temp.nc$dim$lat$vals -> lat lat <- rev(lat) temp <- temp[nrow(temp):1,] temp[temp==-32767] <- NA temp <- t(temp) image(lon,lat,temp) library(maptools) data(wrld_simpl) plot(wrld_simpl, add = TRUE) This code was from modified from

I would like to extract a spatial subset of a rather large netcdf file. From Loop through netcdf files and run calculations - Python or R from pylab import * import netCDF4 f = netCDF4.MFDataset('/usgs/data2/rsignell/models/ncep/narr/air.2m.1989.nc') # print variables f.variables.keys() atemp = f.variables['air'] # TODO: extract spatial subset How do I extract just the subset of netcdf file corresponding to a state (say Iowa). Iowa has following boundary lat lon: Longitude: 89° 5' W to 96° 31' W Latitude: 40° 36' N to 43° 30' N Well this is pretty easy, you have to find the index for the upper

题目： 编译C#环境下GDAL（支持HDF4、NetCDF） 需要文件： gdal-1.8.0 //gdal1.8版本 swigwin-1.3.4 //gdal1.8版本支持1.3.39系列，1.3.4亦可 HDF4.2.6_win_x86 //HDF4驱动，如要支持hdf4，必须编译 netcdf-4.1.1-win32-src.zip //netCDF驱动，需libsrc文件夹 netcdf-4.1.1-win32-bin.zip //netCDF驱动，需netcdf.lib文件 过程： 1、全部解压 2、用记事本打开gdal-1.8.0文件夹中，make.opt 3、分别修改 MSVC_VER=1400 //1400为visualstudio版本 GDAL_HOME = "C:\warmerda\bld6" //gdal编译生成目录 SWIG = C:\gdal\swigwin-1.3.40\swig.exe //swig目录（全路径，否则加环境变量（不建议）） *hdf4： # Uncomment the following and update to enable NCSA HDF Release 4 support. HDF4_PLUGIN = NO HDF4_DIR = E:\hdf4 HDF4_LIB = $(HDF4_DIR)\dll\hd426m.lib $