/usr/share/doc/mpop/polar_aapp1b.py is in python-mpop 1.0.0-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 | #!/usr/bin/env python
# -*- coding: utf-8 -*-
# Copyright (c) 2011 SMHI
# Author(s):
# Martin Raspaud <martin.raspaud@smhi.se>
# This program 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.
# This program 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 this program. If not, see <http://www.gnu.org/licenses/>.
"""Minimal script for polar data (from aapp level 1b).
We take the case of level 1b data (calibrated and geolocalized) from noaa 19,
as output from AAPP.
- Install mpop and pyresample
- Don't forget to set up the PPP_CONFIG_DIR variable to point to your
configuration files.
- Edit the noaa19.cfg configuration file (a template is provided in case
you don't have one) with your data directory:
.. code-block:: ini
[avhrr-level2]
filename = hrpt_%(satellite)s_%Y%m%d_%H%M_%(orbit)s.l1b
dir = /data/polar/
format = aapp1b
- Here is an example of a minimal script that has to be called as soon as a new
swath has arrived
"""
from mpop.satellites import PolarFactory
import sys
from datetime import datetime
if sys.version_info < (2, 5):
import time
def strptime(string, fmt=None):
"""This function is available in the datetime module only
from Python >= 2.5.
"""
return datetime(*time.strptime(string, fmt)[:6])
else:
strptime = datetime.strptime
if __name__ == '__main__':
if len(sys.argv) < 3:
print "Usage: " + sys.argv[0] + " time_string orbit"
sys.exit()
time_string = sys.argv[1]
orbit = sys.argv[2]
time_slot = strptime(time_string, "%Y%m%d%H%M")
global_data = PolarFactory.create_scene("noaa", "19",
"avhrr", time_slot, orbit)
global_data.load()
areas = ["euro4", "scan2"]
for area in areas:
local_data = global_data.project(area)
img = local_data.image.overview()
img.save("overview_" + area + "_" + time_string + ".png")
img = local_data.image.cloudtop()
img.save("cloudtop_" + area + "_" + time_string + ".png")
|