search

pytroll / fogpy

Satellite based fog and low stratus detection and nowcasting
http://fogpy.readthedocs.io/en/latest
GNU General Public License v3.0
20 stars 10 forks source link

Fogpy for ABI produces wrong image #65

Open gerritholl opened 4 years ago

gerritholl commented 4 years ago

Using fogpy 1.2.0+dev52.g86a9cf7.fogpy.for.abi and fogtools 0.1.0.post0.dev43+g12198b9:

show-fog $PLOTDIR/abi-ne.tif -m abi_l1b --sat import/Sat_data/OR_ABI-L1b-RadF-M3C*_G16_s20170731606102_e*.nc --nwcsaf export/CMIC/S_NWC_CMIC_GOES16_NEW-ENGLAND-NR_20170314T160610Z.nc

The result is somewhat artsy but not useful:

abi-ne-failed

gerritholl commented 4 years ago

From the log:

2020-04-07 14:53:35,329 - fogpy.algorithms - DEBUG - Number of spatial coherent fog cloud clusters: --

and

2020-04-07 14:53:47,402 - fogpy.algorithms - INFO - LCTH algorithm results for LowCloudHeightAlgorithm

                    Array size:              5242880
                    Valid cells:             0
                    Assigend cells           0
                      DEM extracted cells    0
                      Lapse rate cells       0
                    Interpolated cells       0
                    Remaining NaN cells      0
                    Excluded negative cells  0
                    Min height:              nan
                    Mean height:             nan
                    Max height:              nan
gerritholl commented 4 years ago

It was trying to resample onto eurol using ABI data, hence everything nan. The black part probably the edge of the ABI full disk (but well out of NWCSAF data). When I correctly use New England:

show-fog $PLOTDIR/abi-ne.tif -m abi_l1b --sat import/Sat_data/OR_ABI-L1b-RadF-M3C*_G16_s20170731606102_e*.nc --nwcsaf export/CMIC/S_NWC_CMIC_GOES16_NEW-ENGLAND-NR_20170314T160610Z.nc -a new-england-500

I get:

abi-ne-1500px