[Qgis-user] Rotated pole projection: Netcdf

[Carlos Camino González]

Hi Etienne and Barend,
Thank you so much for your quickly answer,
I used this command line with cdo:
cdo remapbil,grid_netcdf.txt dust_sconc_00_02.nc test_ncdf.nc
where grid_netcdf.txt correspond with:
gridtype = lonlat
xsize = 382
ysize = 278
xfirst = -57.1926117
xinc = 0.5
yfirst = -10.6783781
yinc = 0.5

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).

Do you know the reason ?

Best Regards

Carlos Camino

2013/2/15 Etienne Tourigny <etourigny.dev at gmail.com>

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



-- 
Carlos Luis Camino González
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/qgis-user/attachments/20130218/9669b4ac/attachment.html>