Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

GDAL geotransform change to wrapped longitude #16

Open
mdsumner opened this issue Apr 27, 2020 · 1 comment
Open

GDAL geotransform change to wrapped longitude #16

mdsumner opened this issue Apr 27, 2020 · 1 comment

Comments

@mdsumner
Copy link
Collaborator

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
Copy link
Collaborator Author

image

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant