<div><div>Hi Etienne and Barend,</div><div>Thank you so much for your quickly answer, </div><div>I used this command line with cdo:</div><div>cdo remapbil,grid_netcdf.txt <a href="http://dust_sconc_00_02.nc/" target="_blank">dust_sconc_00_02.nc</a> <a href="http://test_ncdf.nc/" target="_blank">test_ncdf.nc</a></div>
<div>where grid_netcdf.txt correspond with:</div><div>gridtype = lonlat</div><div>xsize = 382 </div><div>ysize = 278</div><div>xfirst = -57.1926117</div><div>xinc = 0.5</div><div>yfirst = -10.6783781</div><div>yinc = 0.5</div>
<div><br></div><div>Now I have a netcdf with a double domain, because in the xsize and ysize, I multiplied the rows and column number (191 & 139, respectively) by 2. I put this xsize and ysize, because with the original size (191 & 139), cdo only processed a part of the netcdf (half exactly). </div>
<div><br></div><div>Do you know the reason ?</div><div><br></div><div>Best Regards</div><div><br></div>Carlos Camino</div><div><br><div class="gmail_quote">2013/2/15 Etienne Tourigny <span dir="ltr"><<a href="mailto:etourigny.dev@gmail.com" target="_blank">etourigny.dev@gmail.com</a>></span><br>


<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Carlos,<br>
<br>
You will get more information on the grid itself if you try<br>
<br>
gdalinfo NETCDF:"<a href="http://dust_sconc_00_02.nc" target="_blank">dust_sconc_00_02.nc</a>":dust_sconc10<br>
<br>
<br>
Unfortunately, the rotated pole projection is not supported by the<br>
gdal netcdf driver.<br>
<br>
Please try to reproject/remap using the cdo command-line utility[1]<br>
<br>
You should try something like<br>
<br>
cdo -f nc remapbil,<grid> <a href="http://in.nc" target="_blank">in.nc</a> <a href="http://out.nc" target="_blank">out.nc</a><br>
<br>
where <grid> should be replaced by a valid grid description such as<br>
"r360x180" for a global 1x1 degree grid.<br>
See the cdo reference for more details (1.3.2. Grid description)<br>
You can also try the other remapping operators, depending on your data.<br>
<br>
[1] <a href="https://code.zmaw.de/projects/cdo" target="_blank">https://code.zmaw.de/projects/cdo</a><br>
<div><div><br>
On Fri, Feb 15, 2013 at 10:40 AM, Carlos Camino González<br>
<<a href="mailto:caminoccg@gmail.com" target="_blank">caminoccg@gmail.com</a>> wrote:<br>
> Hi all,<br>
><br>
> I have many netcdf with rotated pole projection, I would like to proyect<br>
> these netcdf to a WGS84 or similar projection, in order to see correctly<br>
> into QGIS.   I tried with gdal (gdal_translate), r.in.gdal with grass and<br>
> different libraries with R, such as raster, ncdf, etc .. But nothing  ....<br>
> is it possible ???<br>
> Could someone give me  the correct code or the solution for this problem?<br>
><br>
> For your information, this is the header with gdalinfo -h<br>
><br>
> Size is 512, 512<br>
> Coordinate System is `'<br>
> Metadata:<br>
>   NC_GLOBAL#_FillValue=-32767<br>
>   NC_GLOBAL#Conventions=CF-1.2<br>
>   NC_GLOBAL#Grid_type=B-grid: vectors interpolated to scalar positions<br>
>   NC_GLOBAL#Map_Proj=Rotated latitude longitude<br>
>   NC_GLOBAL#missing_value=-32767<br>
>   NC_GLOBAL#NCO=20130201<br>
> Subdatasets:<br>
>   SUBDATASET_1_NAME=NETCDF:"<a href="http://dust_sconc_00_02.nc" target="_blank">dust_sconc_00_02.nc</a>":dust_sconc10<br>
>   SUBDATASET_1_DESC=[29x1x139x191] dust_load (32-bit floating-point)<br>
>  ....<br>
> Corner Coordinates:<br>
> Upper Left  (    0.0,    0.0)<br>
> Lower Left  (    0.0,  512.0)<br>
> Upper Right (  512.0,    0.0)<br>
> Lower Right (  512.0,  512.0)<br>
> Center      (  256.0,  256.0)<br>
><br>
> Thanks in advance for your answers !!!<br>
><br>
> --<br>
> Carlos Camino<br>
><br>
><br>
><br>
> --<br>
> Carlos Luis Camino González<br>
><br>
><br>
</div></div>> _______________________________________________<br>
> Qgis-user mailing list<br>
> <a href="mailto:Qgis-user@lists.osgeo.org" target="_blank">Qgis-user@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/qgis-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/qgis-user</a><br>
><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Carlos Luis Camino González <br>
</div>