rgdal-dev / rasterwise

Hard-won lessons! Don't lose 'em!
4 stars 0 forks source link

GDAL geotransform change to wrapped longitude #16

Open mdsumner opened 4 years ago

mdsumner commented 4 years ago

ncdump below, this from GDAL 2.2.3

sds <- vapour::vapour_sds_names("~/adaptor.cams_regional_fc.retrieve-1587939061.8606012-22402-6-64f0e215-609d-4ca6-b65c-a9ef7dc53370.nc")
vapour::vapour_raster_info(sds$subdataset[1])
$geotransform
[1] 335.25749851  -0.41502144   0.00000000  71.99999695   0.00000000  -0.09999999

$dimXY
[1] 700 420

$minmax
[1] NA NA

$tilesXY
[1] 700   1

GDAL 3.0.4

vapour::vapour_raster_info(sds$subdataset[1])
$geotransform
[1] -25.00001222   0.10000002   0.00000000  71.99999695   0.00000000  -0.09999999

$dimXY
ncdump -h adaptor.cams_regional_fc.retrieve-1587939061.8606012-22402-6-64f0e215-609d-4ca6-b65c-a9ef7dc53370.nc 
netcdf adaptor.cams_regional_fc.retrieve-1587939061.8606012-22402-6-64f0e215-609d-4ca6-b65c-a9ef7dc53370 {
dimensions:
        longitude = 700 ;
        latitude = 420 ;
        level = 1 ;
        time = UNLIMITED ; // (2 currently)
variables:
        float longitude(longitude) ;
                longitude:long_name = "longitude" ;
                longitude:units = "degrees_east" ;
        float latitude(latitude) ;
                latitude:long_name = "latitude" ;
                latitude:units = "degrees_north" ;
        float level(level) ;
                level:long_name = "level" ;
                level:units = "m" ;
        float time(time) ;
                time:long_name = "ANALYSIS time from 20200101" ;
                time:units = "hours" ;
        float pm10_conc(time, level, latitude, longitude) ;
                pm10_conc:_FillValue = -999.f ;
                pm10_conc:species = "PM10 Aerosol" ;
                pm10_conc:units = "µg/m3" ;
                pm10_conc:value = "hourly values" ;
                pm10_conc:standard_name = "mass_concentration_of_pm10_ambient_aerosol_in_air" ;
        float pm2p5_conc(time, level, latitude, longitude) ;
                pm2p5_conc:_FillValue = -999.f ;
                pm2p5_conc:species = "PM2.5 Aerosol" ;
                pm2p5_conc:units = "µg/m3" ;
                pm2p5_conc:value = "hourly values" ;
                pm2p5_conc:standard_name = "mass_concentration_of_pm2p5_ambient_aerosol_in_air" ;

// global attributes:
                :title = "PM10/PM25 Air Pollutant ANALYSIS at the Surface" ;
                :institution = "Data produced by Meteo France" ;
                :source = "Data from ENSEMBLE model" ;
                :history = "Model ENSEMBLE ANALYSIS" ;
                :ANALYSIS = "Europe, 20200101-20200102+[12H_12H]" ;
                :summary = "ENSEMBLE model hourly ANALYSIS of PM10/PM25 concentration at the Surface from 20200101-20200102+[12H_12H] on Europe" ;
                :project = "MACC-RAQ (http://macc-raq.gmes-atmosphere.eu)" ;
mdsumner commented 4 years ago

image