[Qgis-user] Rotated pole projection: Netcdf

[Etienne Tourigny]

Better post your question on the cdo support forum

https://code.zmaw.de/projects/cdo/boards

On Mon, Feb 18, 2013 at 1:57 PM, Carlos Camino González
<caminoccg at gmail.com> wrote:
> 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