[gdal-dev] NetCDF: how to read coordinates from "lon" and "lat" variables?

andy aborruso at gmail.com
Tue Dec 3 10:33:26 PST 2024


Hi Even,
thank you.

I have this

{
dimensions:
        Time = UNLIMITED ; // (25 currently)
        bottom_top = 1 ;
        south_north = 229 ;
        west_east = 174 ;
        DateStrLen = 19 ;
variables:
        float NO2(Time, bottom_top, south_north, west_east) ;
                NO2:units = "ppb vol" ;
                NO2:long_name = "NO2 Concentration" ;
        float O3(Time, bottom_top, south_north, west_east) ;
                O3:units = "ppb vol" ;
                O3:long_name = "O3 Concentration" ;
        float PM10(Time, bottom_top, south_north, west_east) ;
                PM10:units = "ug/m3" ;
                PM10:long_name = "PM10 Concentration" ;
        float PM25(Time, bottom_top, south_north, west_east) ;
                PM25:units = "ug/m3" ;
                PM25:long_name = "PM25 Concentration" ;
        char Times(Time, DateStrLen) ;
                Times:units = "Date string" ;
                Times:long_name = "WRF-style Times" ;
        double lat(south_north, west_east) ;
                lat:units = "degrees_north" ;
                lat:long_name = "Latitude" ;
        double lon(south_north, west_east) ;
                lon:units = "degrees_east" ;
                lon:long_name = "Longitude" ;
        float pDUST(Time, bottom_top, south_north, west_east) ;
                pDUST:units = "ug/m3" ;
                pDUST:long_name = "pDUST Concentration" ;
        double time(Time) ;
                time:units = "hours since 2000-01-01 00:00:00" ;
        float level(bottom_top) ;
                level:units = "level" ;

// global attributes:
                :Title = "SUITE" ;
                :Sub-title = "Hourly output concentrations file" ;
                :Suite_type = "out" ;
                :Generating_process = "Generated by SUITE" ;
                :Conventions = "None" ;
                :Domain = "BTA7" ;
                :Suite_version = "chimere2017r4v1_maialinux" ;
                :history = "File /fs/out.intput.nc was generated on Sun Jun
16 03:53:11 2024" ;
}


-- 
___________________

Andrea Borruso
website: https://medium.com/tantotanto
38° 7' 48" N, 13° 21' 9" E, EPSG:4326
___________________

"cercare e saper riconoscere chi e cosa,
 in mezzo all’inferno, non è inferno,
e farlo durare, e dargli spazio"

Italo Calvino
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20241203/625880ea/attachment.htm>


More information about the gdal-dev mailing list