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)