Comparing the shapefile exported from VERDI with a shapefile that defines the border of the United States in QGIS was successful.
VERDI obtains projection information from each netCDF file that it loads and saves it to the shapefile.prj file.
cat 2022hc_nonpt_12US1_PM2_5.prj
PROJCS[“unknown”, GEOGCS[“GCS_unknown”, DATUM[“D_unknown”, SPHEROID[“unknown”, 6370000.0, 0.0]], PRIMEM[“Greenwich”, 0.0], UNIT[“degree”, 0.017453292519943295], AXIS[“Geodetic longitude”, EAST], AXIS[“Geodetic latitude”, NORTH]], PROJECTION[“Lambert_Conformal_Conic_2SP”], PARAMETER[“central_meridian”, -97.0], PARAMETER[“latitude_of_origin”, 40.0], PARAMETER[“standard_parallel_1”, 45.0], PARAMETER[“false_easting”, 0.0], PARAMETER[“false_northing”, 0.0], PARAMETER[“scale_factor”, 1.0], PARAMETER[“standard_parallel_2”, 33.0], UNIT[“m”, 1.0], AXIS[“x”, EAST], AXIS[“y”, NORTH]]
Method:
Download files from emismod - Google Drive
verdi.sh -f $cwd/2022hc_nonpt_12US1.ncf -s "PM2_5[1]" -g tile
Select File> Export as Image/GIS
Save as: 2022hc_nonpt_12US1_PM2_5
Change file format: shapefile
Start QGIS, select Layer: add layer: add vector layer, click on … next to Vector Dataset(s), and select shapefiles from your directory.
Once selected, click add.
I added three shapefiles, two that were output as shapefiles from VERDI, and a third that is provided with VERDI for administrative boundary layers of North America.
Boundary layer file of countries is found under the directory:
VERDI_2.1.5/plugins/bootstrap/data/ne_50m_admin_0_countries/ne_50m_admin_0_countries.shp
They all seem to load into QGIS and are properly geo-referenced.
Perhaps the following layer definition files for each of the shapefiles will be helpful to you.
2022hc_nonpt_12US1_PM2_5_layer_definition_file.qlr.txt (30.0 KB)
2022hc_VMT_12US2_layer_definition_file.qlr.txt (108.4 KB)
layer_definition_file_ne_50m.qlr.txt (59.4 KB)
