if len(sys.argv) < 6:
print("*** ")
print("*** Warning, please specify date and time completely, e.g.")
print("*** python plot_hsaf.py 2014 07 23 16 10 ")
print("*** ")
quit() # quit at this point
else:
year = int(sys.argv[1])
month = int(sys.argv[2])
day = int(sys.argv[3])
hour = int(sys.argv[4])
minute = int(sys.argv[5])
else:
if True: # automatic choise of last 5min
from my_msg_module import get_last_SEVIRI_date
datetime1 = get_last_SEVIRI_date(False)
if delay != 0:
datetime1 -= timedelta(minutes=delay)
year = datetime1.year
month = datetime1.month
day = datetime1.day
hour = datetime1.hour
minute = datetime1.minute
else: # fixed date for text reasons
year = 2015
month = 12
day = 16
hour = 13
minute = 30
#prop_str='DBZH'
product_ids["cloud_top_pressure"] = ['ctth_pres', 'ctth_pres_pal']
product_ids['cloud_top_height'] = ['ctth_alti', 'ctth_alti_pal']
product_ids['cloudtype'] = ['ct', 'ct_pal']
product_ids['cloudmask'] = ['cma', 'cma_pal']
for p_ in product_list:
fixed_date = False
if fixed_date:
start_time = datetime(2017, 7, 7, 12, 0)
end_time = datetime(2017, 7, 7, 12, 0)
#base_dir="/data/COALITION2/database/meteosat/SAFNWC_v2013/2015/07/07/",
base_dir = "/data/COALITION2/database/meteosat/SAFNWC_v2016/2017/07/07/" + p_ + "/",
else:
from my_msg_module import get_last_SEVIRI_date
start_time = get_last_SEVIRI_date(False, delay=3)
end_time = start_time
base_dir = "/data/cinesat/in/eumetcast1/"
#file_pattern='S_NWC_' + product + '_' + sat_id + '_*_' + timestamp_NWCSAF + extension
#if glob.glob(file_pattern):
# print ("NWCSAF data product " + product + " exists and is readable")
# # creates list of available products
# product_list_to_process.append(nwcsaf.product.get(product)[:])
#else:
# print ("NWCSAF data product " + product + " is missing or is not readable")
print("... search files in " + base_dir)
files_nwc = find_files_and_readers(sensor='seviri',
start_time=start_time,
end_time=end_time,
print("*** Error, please specify date and time completely, e.g.")
print("*** python " + inspect.getfile(inspect.currentframe()) +
" 2014 07 23 16 10 ")
print("*** ")
quit() # quit at this point
else:
print("*** get date from command line")
year = int(sys.argv[1])
month = int(sys.argv[2])
day = int(sys.argv[3])
hour = int(sys.argv[4])
minute = int(sys.argv[5])
else:
if True: # automatic choise of last 5min
from my_msg_module import get_last_SEVIRI_date
datetime1 = get_last_SEVIRI_date(False, delay=10)
year = datetime1.year
month = datetime1.month
day = datetime1.day
hour = datetime1.hour
minute = datetime1.minute
else: # fixed date for text reasons
year = 2015
month = 12
day = 16
hour = 13
minute = 30
datetime1 = datetime.datetime(year, month, day, hour, minute)
#prop_str='DBZH'
print("*** Warning, please specify date and time completely, e.g.")
print("*** python plot_radar.py 2014 07 23 16 10 ")
print("*** ")
quit() # quit at this point
else:
year = int(sys.argv[1])
month = int(sys.argv[2])
day = int(sys.argv[3])
hour = int(sys.argv[4])
minute = int(sys.argv[5])
if len(sys.argv) > 6:
cell = int(sys.argv[6])
else:
if True: # automatic choise of last 5min
from my_msg_module import get_last_SEVIRI_date
datetime1 = get_last_SEVIRI_date(True)
year = datetime1.year
month = datetime1.month
day = datetime1.day
hour = datetime1.hour
minute = datetime1.minute
else: # fixed date for text reasons
year = 2014 # 2014 09 15 21 35
month = 7 # 2014 07 23 18 30
day = 23
hour = 18
minute = 00
time_slot = datetime(year, month, day, hour, minute)
global_data = GeostationaryFactory.create_scene("swisstrt", "04", "radar",
time_slot)
#import os
#os.chdir(base_dir)
#filenames = glob("*201507071200*__")
#print base_dir
#print filenames
##global_scene = Scene(reader="hrit_msg", filenames=filenames, base_dir=base_dir, ppp_config_dir="/opt/users/hau/PyTroll//cfg_offline/")
#global_scene = Scene(reader="hrit_msg", filenames=filenames, base_dir=base_dir, ppp_config_dir="/opt/users/hau/PyTroll/packages/satpy/satpy/etc")
# new version of satpy after 0.8
#################################
from satpy import Scene, find_files_and_readers
from datetime import datetime
from my_msg_module import get_last_SEVIRI_date
RSS = True
starttime = get_last_SEVIRI_date(RSS, delay=3)
#files = find_files_and_readers(sensor='olci',
# start_time=datetime(2017, 10, 11, 12, 0),
# end_time=datetime(2017, 10, 11, 12, 59),
# base_dir="/home/a001673/data/satellite/Sentinel-3",
# reader='nc_olci_l1b')
#scn = Scene(filenames=files)"
#def find_files_and_readers(start_time=None, end_time=None, base_dir=None,
# reader=None, sensor=None, ppp_config_dir=get_environ_config_dir(),
# filter_parameters=None, reader_kwargs=None):
make_rgb = False
if make_rgb:
# along with this program. If not, see <http://www.gnu.org/licenses/>.
"""Test deriving the 3.9 relfectance of real data and make some imagery
"""
from mpop.satellites import GeostationaryFactory
from datetime import datetime
from mpop.projector import get_area_def
import sys
europe = get_area_def("EuropeCanary")
if len(sys.argv) == 1:
from my_msg_module import get_last_SEVIRI_date
tslot = get_last_SEVIRI_date(True, delay=5)
#tslot = datetime(2013, 10, 11, 11, 30)
#tslot = datetime(2015, 7, 7, 15, 10)
#tslot = datetime(2015, 10, 15, 11, 00)
elif len(sys.argv) == 6:
year = int(sys.argv[1])
month = int(sys.argv[2])
day = int(sys.argv[3])
hour = int(sys.argv[4])
minute = int(sys.argv[5])
tslot = datetime(year, month, day, hour, minute)
else:
print("\n*** Error, wrong number of input arguments")
print(" usage:")
print(" python demo_refl039.py")
print(" or")
from __future__ import division
from __future__ import print_function
from mpop.satellites import GeostationaryFactory
from mpop.projector import get_area_def
import datetime
## uncomment these two lines for more debugging information
from mpop.utils import debug_on
debug_on()
if False:
from my_msg_module import get_last_SEVIRI_date
time_slot = get_last_SEVIRI_date(True, delay=10)
else:
out_path = None
import sys
if len(sys.argv) <= 2:
time_slot = datetime.datetime(2015, 7, 7, 12, 00)
elif len(sys.argv) <= 7:
# python
year = int(sys.argv[1])
month = int(sys.argv[2])
day = int(sys.argv[3])
hour = int(sys.argv[4])
minute = int(sys.argv[5])
time_slot = datetime.datetime(year, month, day, hour, minute)
if len(sys.argv) == 7:
out_path = sys.argv[6]
elif len(sys.argv) >= 7:
print("*** Error: Incorrect number of arguments ***")
#area="cosmo7"
#area="EuropeCanaryS"
elif sat_nr == "10":
Rapid_Scan = False
delay = 10
area = "EuropeCanaryS"
elif sat_nr == "11":
Rapid_Scan = False
delay = 10
area = "EuropeCanaryS"
else:
print("Unknown satellite")
quit()
print("... Meteosat ", sat_nr, ", Rapid Scan mode: ", Rapid_Scan)
time_slot = get_last_SEVIRI_date(Rapid_Scan, delay=delay)
#time_slot = datetime.datetime(2015, 11, 26, 15, 30)
print(str(time_slot))
global_data = GeostationaryFactory.create_scene("meteosat", sat_nr, "seviri",
time_slot)
#prop="CMa"
#prop='CT'
prop = "CTH"
# e.g. return "CTTH" for prop "CTH"
from my_msg_module import get_NWC_pge_name
pge = get_NWC_pge_name(prop)
print("read pge: ", pge)
#chn='IR_108'
chn = 'HRV'
channels = [chn, 'IR_108']
if not verbose:
channels = [
'VIS006', 'VIS008', 'IR_016', 'IR_039', 'WV_062', 'WV_073',
'IR_087', 'IR_097', 'IR_108', 'IR_120', 'IR_134', 'HRV'
]
colors = "rainbow"
#colors="greys"
if near_real_time:
from my_msg_module import get_last_SEVIRI_date
time_slot = get_last_SEVIRI_date(False, delay=3)
else:
if len(sys.argv) == 6:
yyyy = int(sys.argv[1])
month = int(sys.argv[2])
dd = int(sys.argv[3])
hh = int(sys.argv[4])
mm = int(sys.argv[5])
else:
yyyy = int(sys.argv[1])
month = int(sys.argv[2])
dd = int(sys.argv[3])
hh = int(sys.argv[4])
mm = int(sys.argv[5])
time_slot = datetime.datetime(yyyy, month, dd, hh, mm)
minute = int(sys.argv[5])
t1 = datetime(year, month, day, hour, minute)
if len(sys.argv) > 6:
yearStop = int(sys.argv[6])
monthStop = int(sys.argv[7])
dayStop = int(sys.argv[8])
hourStop = int(sys.argv[9])
minuteStop = int(sys.argv[10])
tStop = datetime(yearStop, monthStop, dayStop, hourStop,
minuteStop)
else:
tStop = deepcopy(t1)
data_new, first_time_step = properties_cells(t1, tStop)
else:
if True: # automatic choise of last 5min
from my_msg_module import get_last_SEVIRI_date
time_slot = get_last_SEVIRI_date(True)
if delay != 0:
time_slot -= timedelta(minutes=delay)
nrt = True
tStop = time_slot
print("... chose time (automatically): ", str(tStop))
else:
print("*** ")
print("*** Warning, please specify date and time completely, e.g.")
print("*** python plot_radar.py 2014 07 23 16 10 ")
print("*** ")
quit() # quit at this point
data_new, first_time_step = properties_cells(t1, tStop, None)
import datetime
from PIL import Image
#from mpop.utils import debug_on
#debug_on()
from my_msg_module import get_last_SEVIRI_date
#sat="Meteosat-9"
sat = "Meteosat-10"
if sat == "Meteosat-9":
rapid_scan = True
else:
rapid_scan = False
time_slot = get_last_SEVIRI_date(rapid_scan, delay=5)
#time_slot = datetime.datetime(2015, 11, 26, 19, 30)
print(str(time_slot))
chn = "IR_108"
global_data = GeostationaryFactory.create_scene(sat, "", "seviri", time_slot)
global_data.load([10.8]) #, reader_level="seviri-level2"
area = "EuropeCanaryS95" # only 'PLAT', 'MERC', stere (in Ninjo called 'NPOL', 'SPOL') possible
area = "ccs4" # only 'PLAT', 'MERC', stere (in Ninjo called 'NPOL', 'SPOL') possible
#area="EuroMercator"
area = "germ"
chan_id = 9600015
#chan_id=8800015
local_data = global_data.project(area)
