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