print(filename) ##!#if(composite): ##!# day_filenames = glob(filename[:50]+'*') ##!# cmpst_add = '_composite' ##!#else: day_filenames = glob(filename) cmpst_add = '' # Extract the modis true-color plot limits # ---------------------------------------- ##lat_lims = plot_limits_dict[dt_date_str.strftime('%Y-%m-%d')][dt_date_str.strftime('%H%M')]['modis_Lat'] ##lon_lims = plot_limits_dict[dt_date_str.strftime('%Y-%m-%d')][dt_date_str.strftime('%H%M')]['modis_Lon'] # Use satpy (Scene) to open the file # ---------------------------------- scn = Scene(reader='modis_l1b', filenames=day_filenames) # Load true-color data scn.load(['true_color']) scn.save_dataset('true_color', 'test_image_true2.png') scn.show('true_color') plt.show() ##!#if(save): ##!# outname = 'modis_true_color_' + date_str + zoom_add + cmpst_add + '_satpy.png' ##!# plt.savefig(outname,dpi=300) ##!# print("Saved image",outname) ##!#else: ##!# plt.show()
cot = inputs[8] reff = inputs[9] cwp = inputs[10] lat = inputs[11] lon = inputs[12] msg_con_quick = image.ImageContainerQuick(ir108.squeeze(), area_def) area_con_quick = msg_con_quick.resample(euro_areadef) result_data_quick = area_con_quick.image_data # Create satpy scene testscene = Scene(platform_name="msg", sensor="seviri", start_time=datetime(2013, 11, 12, 8, 30), end_time=datetime(2013, 11, 12, 8, 45), area=area_def) array_kwargs = {'area': area_def} testscene['ir108'] = Dataset(ir108.squeeze(), **array_kwargs) print(testscene['ir108']) testscene.show( 'ir108', overlay={'coast_dir': '/home/mastho/data/', 'color': 'gray'}) resampscene = testscene.resample('germ') print(resampscene.shape) # Define custom fog colormap fogcol = Colormap((0., (250 / 255.0, 200 / 255.0, 40 / 255.0)), (1., (1.0, 1.0, 229 / 255.0))) maskcol = (250 / 255.0, 200 / 255.0, 40 / 255.0)