[Qgis-user] Rotated pole projection: Netcdf

[b.j.kobben at utwente.nl]

As an added comment: trying to see map poles in WGS 84 is going to be
difficult: in WGS84 "projection" where latlons are treated as cartesian
coordinates, a pole (a point) is depicted as a line, making distortions
enormous...

Barend

On 15-02-13 17:06, "Etienne Tourigny" <etourigny.dev at gmail.com> wrote:

>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
>>
>_______________________________________________
>Qgis-user mailing list
>Qgis-user at lists.osgeo.org
>http://lists.osgeo.org/mailman/listinfo/qgis-user