Thank you,

Setting the environmental variable worked....

V

On 01/06/2015 09:11 AM, Julien Demaria wrote:
Dear Victoria,

You can try
export GDAL_NETCDF_BOTTOMUP=NO
just before launching gdal_translate
but note that this works only if your dataset is oriented north-south

For GDAL developpers: I have implemented a path on my version to force 
blockYSize to 1 if bottomup dataset :

--------------------- gdal/frmts/netcdf/netcdfdataset.cpp --------------------- 
index f9de39f..5b7a16a 100644 @@ -457,6 +457,15 @@ 
netCDFRasterBand::netCDFRasterBand( netCDFDataset *poNCDFDS,
          }
        }               
  #endif
+
+/* -------------------------------------------------------------------- */
+/*      Force block size to 1 scanline for bottom-up datasets if        */
+/*      nBlockYSize != 1                                                */
+/* -------------------------------------------------------------------- */
+    if( poNCDFDS->bBottomUp && nBlockYSize != 1 ) {
+        nBlockXSize = nRasterXSize;
+        nBlockYSize = 1;
+    }
  }

I have also implemented several other improvements on the netCDF Driver 
(including full support of groups and netCDF4), I will provide them to the 
community ASAP.


Best Regards,

Julien

-----Message d'origine-----
De : [email protected] [mailto:[email protected]] 
De la part de VictoriaH
Envoyé : mardi 6 janvier 2015 16:11
À : [email protected]
Objet : [gdal-dev] Problems with GDAL 1.10.1 converting from NetCDF to Geotiff

Hello,

I am running GDAL 1.10.1 on an Ubuntu 14.04 machine. When I try to create a 
Geotiff from a NetCDF file, I get the following error. The problem does not 
occur when using an earlier version of GDAL (e.g. 1.9.2). I have seen some 
discussion of this issue on the mailing lists, but it is unclear what the 
workaround is. I would appreciate any information about how to get around this 
issue.


ERROR 1: nBlockYSize = 663, only 1 supported when reading bottom-up dataset 
ERROR 1: NETCDF:PM25_2004_06_average.nc:Data, band 1: IReadBlock failed at X 
offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block 
offset 0

I use the following as an example:

gdal_translate -ot Byte -of GTiff NETCDF:"PM25_2004_06_average.nc":Data
pm_2004.tif

Thanks for any help...



--
View this message in context: 
http://osgeo-org.1560.x6.nabble.com/Problems-with-GDAL-1-10-1-converting-from-NetCDF-to-Geotiff-tp5180197.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

Reply via email to