g.cimis check for bad images before adding to solar.

[qjhart]

The g.cimis program could check for image badness in the goes17 dataset, before deciding whether to project that into the solar location. The way to do that is to project the state boundary into the goes17 mapset, expand it a bit ~10 pixels, and then use that as a mask to do stats on the image. If the stats are bad, then throw an error message, and don't project to solar.

[qjhart]

So, we can create a state boundary in the GOES-17 projection with:

r.proj location=cimis mapset=500m input=state output=state
 r.buffer input=state output=state_plus distance=2000

This leads to this region of interest,

Now, we want use this this, but only if it's needed, first look at the min value of the stats, and only after look at the state boundary

b2=1331PST-B2@20190523
 if [[ $((  $(eval $(r.info -r $b2); echo $min ) >  1000 || 
           $(r.mask --quiet --overwrite  state_plus; 
              r.stats --quiet -c  1331PST-B2@20190523 | head -1 | cut -d' ' -f  1; 
              r.mask --quiet -d ) > 1000)) == 1 ]]; 
then 
 echo YES; 
else 
  echo NO; 
fi;
YES

[qjhart]

And this is the change to the goes.mk file, (see branch check-solar)

# We try and validate the data here, so that we do not propogate bad imagery
$(call f_solar_rast,$1):$(call f_goes_rast,$1)
        @$(call g.mapset,${goes.loc},$(call f_mapset,$1));\
        if [[ $$((  $$(eval $$(r.info -r $(call f_rastname,$1) ); echo $$min ) >  100 ||\
          $$(r.mask --quiet --overwrite  state_plus@500m;\
            r.stats --quiet -c  $(call f_rastname,$1) | head -1 | cut -d' ' -f  1;\
            r.mask --quiet -d ) > 100 )) == 1 ]];  \
        then \
          $(call g.mapset-c,${solar.loc},$(call f_mapset,$1));\
          g.region --quiet ${solar.region};\
          r.proj --quiet location=${goes.loc} mapset=$(call f_mapset,$1) \
            input=$(call f_rastname,$1) method=${solar.proj.method};\
          echo ${solar.loc}/$(call f_mapset,$1)/$(call f_rastname,$1);\
        else \
          echo $(call f_mapset,$1)/$(call f_rastname,$1) BAD DATA;\
        fi;

[qjhart]

@gjscheer-ucd This is ready to test. What was the day with bad data?