From 4568f894c517b7299da6e095db263939cd5bec9b Mon Sep 17 00:00:00 2001 From: "Adam.Dybbroe" Date: Thu, 3 Dec 2015 16:40:26 +0100 Subject: [PATCH] Remove old FY3 mersi reader Signed-off-by: Adam.Dybbroe --- mpop/satin/fy3_mersi_aggr1km.py | 162 -------------------------------- 1 file changed, 162 deletions(-) delete mode 100644 mpop/satin/fy3_mersi_aggr1km.py diff --git a/mpop/satin/fy3_mersi_aggr1km.py b/mpop/satin/fy3_mersi_aggr1km.py deleted file mode 100644 index cadb222c..00000000 --- a/mpop/satin/fy3_mersi_aggr1km.py +++ /dev/null @@ -1,162 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# Copyright (c) 2010. - -# SMHI, -# Folkborgsvägen 1, -# Norrköping, -# Sweden - -# Author(s): - -# Adam Dybbroe - -# This file is part of mpop. - -# mpop is free software: you can redistribute it and/or modify it under the -# terms of the GNU General Public License as published by the Free Software -# Foundation, either version 3 of the License, or (at your option) any later -# version. - -# mpop is distributed in the hope that it will be useful, but WITHOUT ANY -# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR -# A PARTICULAR PURPOSE. See the GNU General Public License for more details. - -# You should have received a copy of the GNU General Public License along with -# mpop. If not, see . - -"""Interface to 1km agregated MERSI hdf5 level 1 format. -""" -import os.path -from ConfigParser import ConfigParser - -import _pyhl -import numpy as np - -from mpop import BASE_PATH - - -def load(satscene): - """Read data from file and load it into *satscene*. - """ - conf = ConfigParser() - conf.read(os.path.join(BASE_PATH, "etc", satscene.fullname + ".cfg")) - options = {} - for option, value in conf.items(satscene.instrument_name+'-level2', raw = True): - options[option] = value - CASES[satscene.instrument_name](satscene, options) - - -def load_1km_aggregated_mersi(satscene, options): - """Read 1km agregated mersi data from file and load it into *satscene*. - """ - # Example: FY3A_MERSI_GBAL_L1_20100308_0915_1000M_MS.HDF - filename = satscene.time_slot.strftime("FY3A_MERSI_GBAL_L1_%Y%m%d_%H%M_1000M_MS.HDF") - filename = os.path.join(options["dir"], filename) - - a = _pyhl.read_nodelist(filename) - b = a.getNodeNames() - # Should only select/fetch the datasets needed. FIXME! - a.selectAll() - a.fetch() - - # MERSI Channel 1-4: EV_250_Aggr.1KM_RefSB - # MERSI Channel 5: EV_250_Aggr.1KM_Emissive - # MERSI Channel 6-20: EV_1KM_RefSB - - datasets = ['/EV_250_Aggr.1KM_RefSB', - '/EV_250_Aggr.1KM_Emissive', - '/EV_1KM_RefSB'] - - for nodename in datasets: - band_data = a.getNode(nodename).data() - valid_range = a.getNode('%s/valid_range' % (nodename)).data() - band_names = a.getNode('%s/band_name' % (nodename)).data().split(",") - # Special treatment of the situation where the bandnames are stored - # as '6~20' (quite inconvenient): - if '6~20' in band_names: - band_names = ['6','7','8','9','10','11','12','13', - '14','15','16','17','18','19','20'] - - for (i, band) in enumerate(band_names): - if band not in satscene.channels_to_load: - continue - - satscene[band] = np.ma.masked_outside(band_data[i], - valid_range[0], - valid_range[1], - copy = False) - satscene[band].info = { - 'var_name' : 'ch'+str(band), - 'var_data' : satscene[band].data, - 'var_dim_names': ('x','y'), - '_FillValue' : 32767, - 'standard_name' : '', - 'short_name' : band, - 'scale_factor' : 1.0, - 'add_offset' : 0.0, - } - - - satscene.info = { - 'var_children' : [ #{'var_name' : 'lat', 'var_callback': Functor(satscene.get_lat, low_res), 'var_dim_names': ('x','y') }, - #{'var_name' : 'lon', 'var_callback' : Functor(satscene.get_lon, low_res) , 'var_dim_names': ('x','y') }, - ## {'var_name' : 'lat_hrvis', 'var_data' : satscene.lat[high_res]}, - ## {'var_name' : 'lon_hrvis', 'var_data' : satscene.lon[high_res]}, - ], - 'Satellite' : satscene.fullname, - 'Antenna' : 'None', - 'Receiver' : 'SMHI' , - 'Time' : satscene.time_slot.strftime("%Y-%m-%d %H:%M:%S UTC"), - #'Area_Name' : satscene.area_id, - 'Area_Name' : "swath", - 'Projection' : 'satproj', - 'Platform' : 'fengyun', - 'Number' : '3a', - 'Service' : '', - #'Columns' : satscene.channels[0].shape[1], - #'Lines' : satscene.channels[0].shape[0], - 'SampleX' : 1.0, - 'SampleY' : 1.0, - 'title' : 'MERSI Level 1', - } - - -def get_lat_lon(satscene, resolution): - """Read lat and lon. - """ - del resolution - - conf = ConfigParser() - conf.read(os.path.join(BASE_PATH, "etc", satscene.fullname + ".cfg")) - options = {} - for option, value in conf.items(satscene.instrument_name+'-level2', raw = True): - options[option] = value - - return LAT_LON_CASES[satscene.instrument_name](satscene, options) - -def get_lat_lon_1km_aggregated_mersi(satscene, options): - """Read latitude and longitude for each (aggregated) pixel. - """ - # Example: FY3A_MERSI_GBAL_L1_20100308_0915_1000M_MS.HDF - filename = satscene.time_slot.strftime("FY3A_MERSI_GBAL_L1_%Y%M%D_%H%M_1000M_MS.HDF") - filename = os.path.join(options["dir"], filename) - - a = _pyhl.read_nodelist(filename) - b = a.getNodeNames() - # Should only select/fetch the datasets needed. FIXME! - a.selectAll() - a.fetch() - - lat = a.getNode("/Latitude").data() - lon = a.getNode("/Longitude").data() - - return lat, lon - -CASES = { - "mersi": load_1km_aggregated_mersi - } - -LAT_LON_CASES = { - "mersi": get_lat_lon_1km_aggregated_mersi - }