Command line works fine. I didn't have time to test well though. I will try to get on a windows machine and download tamas package and give it another shot.
kss On Fri, Jun 8, 2012 at 9:31 AM, Etienne Tourigny <[email protected]>wrote: > with commandline and/or python? > > I can do both fine using gdal-1.9.1 from osgeo4w (32bit), haven't > tried Tamas' build yet > > On Fri, Jun 8, 2012 at 12:08 PM, Kyle Shannon <[email protected]> > wrote: > > Etienne, Rutger, > > I can reproduce on windows with 32 bit 1.9 branch and 1.9.1 from Tamas' > site. In fact, I couldn't open any datasets. > > > > kss > > > > -----Original Message----- > > From: [email protected] [mailto: > [email protected]] On Behalf Of Etienne Tourigny > > Sent: Friday, June 08, 2012 9:10 AM > > To: Rutger > > Cc: [email protected] > > Subject: Re: [gdal-dev] Reading a NetCDF file with python (not > recognised as a supported dataset name > > > > Rutger, > > > > There should not be any difference when using the python bindings. If > gdalinfo can read it the python bindings should also. > > > > I cannot reproduce your problem using gdal-1.9.1 under linux. Which > exact version are you using? > > > > By the way, the vrtfile in unnecessary, because GDAL's Open accepts the > NETCDF:file.nc:variable syntax. > > > > Try adding this to your script and send the result back here > gdal.SetConfigOption( 'CPL_DEBUG', 'ON' ) > > > > also try gdal_info.exe -stats > > > NETCDF:LPRM-TMI_L3_DY_SOILM3_V001-20120411T144114Z_20120102.nc:soil_moisture_x > > this actually reads the file contents and not just the metadata. If it > fails then something is wrong. > > > > you should see > > STATISTICS_MAXIMUM=99 > > STATISTICS_MEAN=21.746993313128 > > STATISTICS_MINIMUM=0 > > STATISTICS_STDDEV=20.100893016522 > > > > > > Etienne > > > > On Fri, Jun 8, 2012 at 9:00 AM, Rutger <[email protected]> wrote: > >> Dear list, > >> > >> I am trying to open a NetCDF file with the Python bindings. My GDAL > >> 1.9 build (Win32 from Tamas) shows the NetCDF driver when i run > >> "gdalinfo --formats", so it should be possible. > >> > >> When using the pre-build command-line utilities all works well. But > >> when i move to Python i get the message: > >> /" ... does not exist in the file system, and is not recognised as a > >> supported dataset name."/ I dont have the slightest clue about what > >> i'm doing wrong. > >> > >> I have simplified the problem to 3 steps: > >> 1. Use gdal_translate to select a subdataset and output to VRT 2. Read > >> the VRT in Python with gdal.Open() 3. Use gdal_translate to convert > >> the VRT to a Geotiff > >> > >> Step 1 & 3 work well, but 2 fails. Is there some limitation of the > >> Python bindings regarding NetCDF? The script runs fine if i replace > >> the NetCDF with a Geotiff file. > >> > >> Any help would be appreciated. > >> > >> An example NetCDF i used can be found at: > >> ftp://hydro1.sci.gsfc.nasa.gov/data/s4pa/WAOB/LPRM_TMI_DY_SOILM3.001/2 > >> 012/05/ > >> ftp://hydro1.sci.gsfc.nasa.gov/data/s4pa/WAOB/LPRM_TMI_DY_SOILM3.001/2 > >> 012/05/ > >> > >> Below is the Python script i have used for exploring the issue: > >> > >> import os > >> try: > >> from osgeo import gdal > >> except ImportError: > >> import gdal > >> > >> netcdffile = > >> > "NETCDF:LPRM-TMI_L3_DY_SOILM3_V001-20120411T144114Z_20120102.nc:soil_moisture_x" > >> vrtfile = "LPRM-TMI_L3_DY_SOILM3_V001-20120411T144114Z_20120102.vrt" > >> tiffile = "LPRM-TMI_L3_DY_SOILM3_V001-20120411T144114Z_20120102.tif" > >> > >> # 1: convert the netcdf to vrt externally (works well) > >> os.system('gdal_translate -of "VRT" ' + netcdffile + ' ' + vrtfile) > >> > >> # 2: try to open the vrt in Python (fails at .ReadAsArray) dsin = > >> gdal.Open(vrtfile) > >> > >> if not dsin is None: > >> indata = dsin.ReadAsArray(0, 0, dsin.RasterXSize, dsin.RasterYSize) > >> > >> # 3: convert the same vrt to a tif externally (works well) > >> os.system('gdal_translate ' + vrtfile + ' ' + tiffile) > >> > >> > >> Regards, > >> Rutger > >> > >> > >> -- > >> View this message in context: > >> http://osgeo-org.1560.n6.nabble.com/Reading-a-NetCDF-file-with-python- > >> not-recognised-as-a-supported-dataset-name-tp4979943.html > >> Sent from the GDAL - Dev mailing list archive at Nabble.com. > >> _______________________________________________ > >> gdal-dev mailing list > >> [email protected] > >> http://lists.osgeo.org/mailman/listinfo/gdal-dev > > _______________________________________________ > > gdal-dev mailing list > > [email protected] > > http://lists.osgeo.org/mailman/listinfo/gdal-dev > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev >
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
