search

Global-Flood-Assessment / MoMProduction

Production setup for ModelofModels
1 stars 1 forks source link

Compare the results with the development code #14

Closed webgisdeveloper closed 1 year ago

webgisdeveloper commented 2 years ago

especially for GloFAS code . The PDC produced final csv has 70 more records and slightly lower scores for GloFAS but is identifying the same areas. Slight differences in GFMS scores. DFO scores appear identical. HWRF almost identical. PDC VIIRS scores trend slightly lower and some different areas identified. Good agreement with final watershed severity scores and geographic areas of watersheds identified. PDC has 346 individual Warnings vs 393 for IU.

webgisdeveloper commented 2 years ago

first step, check the watershed data, the IU one may use two different version wastersheds shapefile

webgisdeveloper commented 2 years ago

on the server: two shape files first one:
/ModelofModels/data/Watersheds_032020 also in DFO processing 

ogrinfo -so wastershed_prj_latlon.shp wastershed_prj_latlon  

Layer name: wastershed_prj_latlon
Metadata:
  DBF_DATE_LAST_UPDATE=2020-04-23
Geometry: Polygon
Feature Count: 16396
Extent: (-180.000000, -55.987500) - (180.000000, 83.625643)
Layer SRS WKT:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AXIS["Latitude",NORTH],
    AXIS["Longitude",EAST],
    AUTHORITY["EPSG","4326"]]
Data axis to CRS axis mapping: 2,1
pfaf_id: Integer64 (10.0)
areakm2: Real (11.3)

second one: in HWRF and VIIRS 

ogrinfo -so Watershed_pfaf_id.shp Watershed_pfaf_id 

Layer name: Watershed_pfaf_id
Metadata:
  DBF_DATE_LAST_UPDATE=2021-02-18
Geometry: Polygon
Feature Count: 16396
Extent: (-180.000000, -55.987500) - (180.000618, 83.625643)
Layer SRS WKT:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AXIS["Latitude",NORTH],
    AXIS["Longitude",EAST],
    AUTHORITY["EPSG","4326"]]
Data axis to CRS axis mapping: 2,1
pfaf_id: Integer64 (10.0)
area_km2: Real (19.11)
ISO: String (254.0)
Admin0: String (254.0)
Admin1: String (254.0)
rfr_score: Real (19.11)
cfr_score: Real (19.11)
webgisdeveloper commented 2 years ago

the production server use the second one

webgisdeveloper commented 2 years ago

area difference less than 1%

pfaf_id areakm2     pfaf_id area_km2 ISO Admin0 Admin1 rfr_score cfr_score   area_differnce %
111011 1888.469     111011 1894.93587 EGY Egypt As Suways 4.18067408 0   0.34243991
111012 2922.8     111012 2932.50625 EGY Egypt As Suways 0 0   0.33208754
111013 891.315     111013 895.121761 EGY Egypt As Suways 4.0081118 0   0.42709491
111014 4212.89     111014 4226.77571 EGY Egypt Al Bahr al Ahmar 0 0   0.32960054
111015 16606.605     111015 16672.3233 EGY Egypt Al Bahr al Ahmar 4.15950606 1.33295875   0.3957361
111016 1829.616     111016 1836.63245 EGY Egypt Al Bahr al Ahmar 4.1051073 0   0.38349301
111017 8628.285     111017 8665.54055 EGY Egypt Al Bahr al Ahmar 4.00795834 0   0.43178393
111018 1950.134     111018 1960.07436 EGY Egypt Al Bahr al Ahmar 0 0   0.50972702
111019 6550.057     111019 6579.74772 EGY Egypt Al Bahr al Ahmar 4.10117093 0   0.45328945
111020 11651.153     111020 11704.6189 EGY Egypt Al Bahr al Ahmar 4.05183492 0   0.45888927