IM_Prec_monthly_anomalies.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import xuleta2, os\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "folder = '/media/denis/GEO/DROUGHT_processed/PRECIPITATION/monthly/'\n",
    "folderout = '/media/denis/GEO/DROUGHT_processed/PRECIPITATION/monthly_mean/'\n",
    "if not os.path.exists(folderout): os.makedirs(folderout)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "lista = xuleta.listfiles(folder)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "meses = ['01','02','03','04','05','06','07','08','09','10','11','12']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "-c:9: RuntimeWarning: invalid value encountered in less\n",
      "/usr/lib/python2.7/dist-packages/numpy/lib/nanfunctions.py:607: RuntimeWarning: Mean of empty slice\n",
      "  warnings.warn(\"Mean of empty slice\", RuntimeWarning)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "calculando media do mes 01\n",
      "mean01 saved.\n",
      " \n",
      "stdv01 saved.\n",
      " \n",
      "calculando media do mes 02\n",
      "mean02 saved.\n",
      " \n",
      "stdv02 saved.\n",
      " \n",
      "calculando media do mes 03\n",
      "mean03 saved.\n",
      " \n",
      "stdv03 saved.\n",
      " \n",
      "calculando media do mes 04\n",
      "mean04 saved.\n",
      " \n",
      "stdv04 saved.\n",
      " \n",
      "calculando media do mes 05\n",
      "mean05 saved.\n",
      " \n",
      "stdv05 saved.\n",
      " \n",
      "calculando media do mes 06\n",
      "mean06 saved.\n",
      " \n",
      "stdv06 saved.\n",
      " \n",
      "calculando media do mes 07\n",
      "mean07 saved.\n",
      " \n",
      "stdv07 saved.\n",
      " \n",
      "calculando media do mes 08\n",
      "mean08 saved.\n",
      " \n",
      "stdv08 saved.\n",
      " \n",
      "calculando media do mes 09\n",
      "mean09 saved.\n",
      " \n",
      "stdv09 saved.\n",
      " \n",
      "calculando media do mes 10\n",
      "mean10 saved.\n",
      " \n",
      "stdv10 saved.\n",
      " \n",
      "calculando media do mes 11\n",
      "mean11 saved.\n",
      " \n",
      "stdv11 saved.\n",
      " \n",
      "calculando media do mes 12\n",
      "mean12 saved.\n",
      " \n",
      "stdv12 saved.\n",
      " \n"
     ]
    }
   ],
   "source": [
    "for mes in meses:\n",
    "    toOpen = []\n",
    "    for mon in lista:\n",
    "        if mon[-6:-4] == mes:\n",
    "            toOpen.append(mon)\n",
    "    arrays = []\n",
    "    for j in toOpen:\n",
    "        imagem, meta = xuleta.TifToArray(folder+j)\n",
    "        imagem = np.where(imagem<0,np.nan,imagem)\n",
    "        arrays.append(imagem)\n",
    "        arrays2 = np.dstack(arrays)\n",
    "    \n",
    "    print 'calculando media do mes %s' %mes\n",
    "    xuleta.ArrayToTif(np.nanmean(arrays2, axis=2), 'mean'+mes, folderout, meta, Type=3)\n",
    "    xuleta.ArrayToTif(np.std(arrays2, axis=2), 'stdv'+mes, folderout, meta, Type=3)\n",
    "    toOpen, arrays, arrays2 = None, None, None"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "folderz = '/media/denis/GEO/DROUGHT_processed/PRECIPITATION/monthly_z/'\n",
    "if not os.path.exists(folderz): os.makedirs(folderz)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "mean = xuleta.listfiles(folderout+'mean') \n",
    "stdv = xuleta.listfiles(folderout+'stdv') "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "writing z-score for sum2000-03.tif\n",
      "sum2000-03.tif saved.\n",
      " \n",
      "writing z-score for sum2000-04.tif\n",
      "sum2000-04.tif saved.\n",
      " \n",
      "writing z-score for sum2000-05.tif\n",
      "sum2000-05.tif saved.\n",
      " \n",
      "writing z-score for sum2000-06.tif\n",
      "sum2000-06.tif saved.\n",
      " \n",
      "writing z-score for sum2000-07.tif\n",
      "sum2000-07.tif saved.\n",
      " \n",
      "writing z-score for sum2000-08.tif\n",
      "sum2000-08.tif saved.\n",
      " \n",
      "writing z-score for sum2000-09.tif\n",
      "sum2000-09.tif saved.\n",
      " \n",
      "writing z-score for sum2000-10.tif\n",
      "sum2000-10.tif saved.\n",
      " \n",
      "writing z-score for sum2000-11.tif\n",
      "sum2000-11.tif saved.\n",
      " \n",
      "writing z-score for sum2000-12.tif\n",
      "sum2000-12.tif saved.\n",
      " \n",
      "writing z-score for sum2001-01.tif\n",
      "sum2001-01.tif saved.\n",
      " \n",
      "writing z-score for sum2001-02.tif\n",
      "sum2001-02.tif saved.\n",
      " \n",
      "writing z-score for sum2001-03.tif\n",
      "sum2001-03.tif saved.\n",
      " \n",
      "writing z-score for sum2001-04.tif\n",
      "sum2001-04.tif saved.\n",
      " \n",
      "writing z-score for sum2001-05.tif\n",
      "sum2001-05.tif saved.\n",
      " \n",
      "writing z-score for sum2001-06.tif\n",
      "sum2001-06.tif saved.\n",
      " \n",
      "writing z-score for sum2001-07.tif\n",
      "sum2001-07.tif saved.\n",
      " \n",
      "writing z-score for sum2001-08.tif\n",
      "sum2001-08.tif saved.\n",
      " \n",
      "writing z-score for sum2001-09.tif\n",
      "sum2001-09.tif saved.\n",
      " \n",
      "writing z-score for sum2001-10.tif\n",
      "sum2001-10.tif saved.\n",
      " \n",
      "writing z-score for sum2001-11.tif\n",
      "sum2001-11.tif saved.\n",
      " \n",
      "writing z-score for sum2001-12.tif\n",
      "sum2001-12.tif saved.\n",
      " \n",
      "writing z-score for sum2002-01.tif\n",
      "sum2002-01.tif saved.\n",
      " \n",
      "writing z-score for sum2002-02.tif\n",
      "sum2002-02.tif saved.\n",
      " \n",
      "writing z-score for sum2002-03.tif\n",
      "sum2002-03.tif saved.\n",
      " \n",
      "writing z-score for sum2002-04.tif\n",
      "sum2002-04.tif saved.\n",
      " \n",
      "writing z-score for sum2002-05.tif\n",
      "sum2002-05.tif saved.\n",
      " \n",
      "writing z-score for sum2002-06.tif\n",
      "sum2002-06.tif saved.\n",
      " \n",
      "writing z-score for sum2002-07.tif\n",
      "sum2002-07.tif saved.\n",
      " \n",
      "writing z-score for sum2002-08.tif\n",
      "sum2002-08.tif saved.\n",
      " \n",
      "writing z-score for sum2002-09.tif\n",
      "sum2002-09.tif saved.\n",
      " \n",
      "writing z-score for sum2002-10.tif\n",
      "sum2002-10.tif saved.\n",
      " \n",
      "writing z-score for sum2002-11.tif\n",
      "sum2002-11.tif saved.\n",
      " \n",
      "writing z-score for sum2002-12.tif\n",
      "sum2002-12.tif saved.\n",
      " \n",
      "writing z-score for sum2003-01.tif\n",
      "sum2003-01.tif saved.\n",
      " \n",
      "writing z-score for sum2003-02.tif\n",
      "sum2003-02.tif saved.\n",
      " \n",
      "writing z-score for sum2003-03.tif\n",
      "sum2003-03.tif saved.\n",
      " \n",
      "writing z-score for sum2003-04.tif\n",
      "sum2003-04.tif saved.\n",
      " \n",
      "writing z-score for sum2003-05.tif\n",
      "sum2003-05.tif saved.\n",
      " \n",
      "writing z-score for sum2003-06.tif\n",
      "sum2003-06.tif saved.\n",
      " \n",
      "writing z-score for sum2003-07.tif\n",
      "sum2003-07.tif saved.\n",
      " \n",
      "writing z-score for sum2003-08.tif\n",
      "sum2003-08.tif saved.\n",
      " \n",
      "writing z-score for sum2003-09.tif\n",
      "sum2003-09.tif saved.\n",
      " \n",
      "writing z-score for sum2003-10.tif\n",
      "sum2003-10.tif saved.\n",
      " \n",
      "writing z-score for sum2003-11.tif\n",
      "sum2003-11.tif saved.\n",
      " \n",
      "writing z-score for sum2003-12.tif\n",
      "sum2003-12.tif saved.\n",
      " \n",
      "writing z-score for sum2004-01.tif\n",
      "sum2004-01.tif saved.\n",
      " \n",
      "writing z-score for sum2004-02.tif\n",
      "sum2004-02.tif saved.\n",
      " \n",
      "writing z-score for sum2004-03.tif\n",
      "sum2004-03.tif saved.\n",
      " \n",
      "writing z-score for sum2004-04.tif\n",
      "sum2004-04.tif saved.\n",
      " \n",
      "writing z-score for sum2004-05.tif\n",
      "sum2004-05.tif saved.\n",
      " \n",
      "writing z-score for sum2004-06.tif\n",
      "sum2004-06.tif saved.\n",
      " \n",
      "writing z-score for sum2004-07.tif\n",
      "sum2004-07.tif saved.\n",
      " \n",
      "writing z-score for sum2004-08.tif\n",
      "sum2004-08.tif saved.\n",
      " \n",
      "writing z-score for sum2004-09.tif\n",
      "sum2004-09.tif saved.\n",
      " \n",
      "writing z-score for sum2004-10.tif\n",
      "sum2004-10.tif saved.\n",
      " \n",
      "writing z-score for sum2004-11.tif\n",
      "sum2004-11.tif saved.\n",
      " \n",
      "writing z-score for sum2004-12.tif\n",
      "sum2004-12.tif saved.\n",
      " \n",
      "writing z-score for sum2005-01.tif\n",
      "sum2005-01.tif saved.\n",
      " \n",
      "writing z-score for sum2005-02.tif\n",
      "sum2005-02.tif saved.\n",
      " \n",
      "writing z-score for sum2005-03.tif\n",
      "sum2005-03.tif saved.\n",
      " \n",
      "writing z-score for sum2005-04.tif\n",
      "sum2005-04.tif saved.\n",
      " \n",
      "writing z-score for sum2005-05.tif\n",
      "sum2005-05.tif saved.\n",
      " \n",
      "writing z-score for sum2005-06.tif\n",
      "sum2005-06.tif saved.\n",
      " \n",
      "writing z-score for sum2005-07.tif\n",
      "sum2005-07.tif saved.\n",
      " \n",
      "writing z-score for sum2005-08.tif\n",
      "sum2005-08.tif saved.\n",
      " \n",
      "writing z-score for sum2005-09.tif\n",
      "sum2005-09.tif saved.\n",
      " \n",
      "writing z-score for sum2005-10.tif\n",
      "sum2005-10.tif saved.\n",
      " \n",
      "writing z-score for sum2005-11.tif\n",
      "sum2005-11.tif saved.\n",
      " \n",
      "writing z-score for sum2005-12.tif\n",
      "sum2005-12.tif saved.\n",
      " \n",
      "writing z-score for sum2006-01.tif\n",
      "sum2006-01.tif saved.\n",
      " \n",
      "writing z-score for sum2006-02.tif\n",
      "sum2006-02.tif saved.\n",
      " \n",
      "writing z-score for sum2006-03.tif\n",
      "sum2006-03.tif saved.\n",
      " \n",
      "writing z-score for sum2006-04.tif\n",
      "sum2006-04.tif saved.\n",
      " \n",
      "writing z-score for sum2006-05.tif\n",
      "sum2006-05.tif saved.\n",
      " \n",
      "writing z-score for sum2006-06.tif\n",
      "sum2006-06.tif saved.\n",
      " \n",
      "writing z-score for sum2006-07.tif\n",
      "sum2006-07.tif saved.\n",
      " \n",
      "writing z-score for sum2006-08.tif\n",
      "sum2006-08.tif saved.\n",
      " \n",
      "writing z-score for sum2006-09.tif\n",
      "sum2006-09.tif saved.\n",
      " \n",
      "writing z-score for sum2006-10.tif\n",
      "sum2006-10.tif saved.\n",
      " \n",
      "writing z-score for sum2006-11.tif\n",
      "sum2006-11.tif saved.\n",
      " \n",
      "writing z-score for sum2006-12.tif\n",
      "sum2006-12.tif saved.\n",
      " \n",
      "writing z-score for sum2007-01.tif\n",
      "sum2007-01.tif saved.\n",
      " \n",
      "writing z-score for sum2007-02.tif\n",
      "sum2007-02.tif saved.\n",
      " \n",
      "writing z-score for sum2007-03.tif\n",
      "sum2007-03.tif saved.\n",
      " \n",
      "writing z-score for sum2007-04.tif\n",
      "sum2007-04.tif saved.\n",
      " \n",
      "writing z-score for sum2007-05.tif\n",
      "sum2007-05.tif saved.\n",
      " \n",
      "writing z-score for sum2007-06.tif\n",
      "sum2007-06.tif saved.\n",
      " \n",
      "writing z-score for sum2007-07.tif\n",
      "sum2007-07.tif saved.\n",
      " \n",
      "writing z-score for sum2007-08.tif\n",
      "sum2007-08.tif saved.\n",
      " \n",
      "writing z-score for sum2007-09.tif\n",
      "sum2007-09.tif saved.\n",
      " \n",
      "writing z-score for sum2007-10.tif\n",
      "sum2007-10.tif saved.\n",
      " \n",
      "writing z-score for sum2007-11.tif\n",
      "sum2007-11.tif saved.\n",
      " \n",
      "writing z-score for sum2007-12.tif\n",
      "sum2007-12.tif saved.\n",
      " \n",
      "writing z-score for sum2008-01.tif\n",
      "sum2008-01.tif saved.\n",
      " \n",
      "writing z-score for sum2008-02.tif\n",
      "sum2008-02.tif saved.\n",
      " \n",
      "writing z-score for sum2008-03.tif\n",
      "sum2008-03.tif saved.\n",
      " \n",
      "writing z-score for sum2008-04.tif\n",
      "sum2008-04.tif saved.\n",
      " \n",
      "writing z-score for sum2008-05.tif\n",
      "sum2008-05.tif saved.\n",
      " \n",
      "writing z-score for sum2008-06.tif\n",
      "sum2008-06.tif saved.\n",
      " \n",
      "writing z-score for sum2008-07.tif\n",
      "sum2008-07.tif saved.\n",
      " \n",
      "writing z-score for sum2008-08.tif\n",
      "sum2008-08.tif saved.\n",
      " \n",
      "writing z-score for sum2008-09.tif\n",
      "sum2008-09.tif saved.\n",
      " \n",
      "writing z-score for sum2008-10.tif\n",
      "sum2008-10.tif saved.\n",
      " \n",
      "writing z-score for sum2008-11.tif\n",
      "sum2008-11.tif saved.\n",
      " \n",
      "writing z-score for sum2008-12.tif\n",
      "sum2008-12.tif saved.\n",
      " \n",
      "writing z-score for sum2009-01.tif\n",
      "sum2009-01.tif saved.\n",
      " \n",
      "writing z-score for sum2009-02.tif\n",
      "sum2009-02.tif saved.\n",
      " \n",
      "writing z-score for sum2009-03.tif\n",
      "sum2009-03.tif saved.\n",
      " \n",
      "writing z-score for sum2009-04.tif\n",
      "sum2009-04.tif saved.\n",
      " \n",
      "writing z-score for sum2009-05.tif\n",
      "sum2009-05.tif saved.\n",
      " \n",
      "writing z-score for sum2009-06.tif\n",
      "sum2009-06.tif saved.\n",
      " \n",
      "writing z-score for sum2009-07.tif\n",
      "sum2009-07.tif saved.\n",
      " \n",
      "writing z-score for sum2009-08.tif\n",
      "sum2009-08.tif saved.\n",
      " \n",
      "writing z-score for sum2009-09.tif\n",
      "sum2009-09.tif saved.\n",
      " \n",
      "writing z-score for sum2009-10.tif\n",
      "sum2009-10.tif saved.\n",
      " \n",
      "writing z-score for sum2009-11.tif\n",
      "sum2009-11.tif saved.\n",
      " \n",
      "writing z-score for sum2009-12.tif\n",
      "sum2009-12.tif saved.\n",
      " \n",
      "writing z-score for sum2010-01.tif\n",
      "sum2010-01.tif saved.\n",
      " \n",
      "writing z-score for sum2010-02.tif\n",
      "sum2010-02.tif saved.\n",
      " \n",
      "writing z-score for sum2010-03.tif\n",
      "sum2010-03.tif saved.\n",
      " \n",
      "writing z-score for sum2010-04.tif\n",
      "sum2010-04.tif saved.\n",
      " \n",
      "writing z-score for sum2010-05.tif\n",
      "sum2010-05.tif saved.\n",
      " \n",
      "writing z-score for sum2010-06.tif\n",
      "sum2010-06.tif saved.\n",
      " \n",
      "writing z-score for sum2010-07.tif\n",
      "sum2010-07.tif saved.\n",
      " \n",
      "writing z-score for sum2010-08.tif\n",
      "sum2010-08.tif saved.\n",
      " \n",
      "writing z-score for sum2010-09.tif\n",
      "sum2010-09.tif saved.\n",
      " \n",
      "writing z-score for sum2010-10.tif\n",
      "sum2010-10.tif saved.\n",
      " \n",
      "writing z-score for sum2010-11.tif\n",
      "sum2010-11.tif saved.\n",
      " \n",
      "writing z-score for sum2010-12.tif\n",
      "sum2010-12.tif saved.\n",
      " \n",
      "writing z-score for sum2011-01.tif\n",
      "sum2011-01.tif saved.\n",
      " \n",
      "writing z-score for sum2011-02.tif\n",
      "sum2011-02.tif saved.\n",
      " \n",
      "writing z-score for sum2011-03.tif\n",
      "sum2011-03.tif saved.\n",
      " \n",
      "writing z-score for sum2011-04.tif\n",
      "sum2011-04.tif saved.\n",
      " \n",
      "writing z-score for sum2011-05.tif\n",
      "sum2011-05.tif saved.\n",
      " \n",
      "writing z-score for sum2011-06.tif\n",
      "sum2011-06.tif saved.\n",
      " \n",
      "writing z-score for sum2011-07.tif\n",
      "sum2011-07.tif saved.\n",
      " \n",
      "writing z-score for sum2011-08.tif\n",
      "sum2011-08.tif saved.\n",
      " \n",
      "writing z-score for sum2011-09.tif\n",
      "sum2011-09.tif saved.\n",
      " \n",
      "writing z-score for sum2011-10.tif\n",
      "sum2011-10.tif saved.\n",
      " \n",
      "writing z-score for sum2011-11.tif\n",
      "sum2011-11.tif saved.\n",
      " \n",
      "writing z-score for sum2011-12.tif\n",
      "sum2011-12.tif saved.\n",
      " \n",
      "writing z-score for sum2012-01.tif\n",
      "sum2012-01.tif saved.\n",
      " \n",
      "writing z-score for sum2012-02.tif\n",
      "sum2012-02.tif saved.\n",
      " \n",
      "writing z-score for sum2012-03.tif\n",
      "sum2012-03.tif saved.\n",
      " \n",
      "writing z-score for sum2012-04.tif\n",
      "sum2012-04.tif saved.\n",
      " \n",
      "writing z-score for sum2012-05.tif\n",
      "sum2012-05.tif saved.\n",
      " \n",
      "writing z-score for sum2012-06.tif\n",
      "sum2012-06.tif saved.\n",
      " \n",
      "writing z-score for sum2012-07.tif\n",
      "sum2012-07.tif saved.\n",
      " \n",
      "writing z-score for sum2012-08.tif\n",
      "sum2012-08.tif saved.\n",
      " \n",
      "writing z-score for sum2012-09.tif\n",
      "sum2012-09.tif saved.\n",
      " \n",
      "writing z-score for sum2012-10.tif\n",
      "sum2012-10.tif saved.\n",
      " \n",
      "writing z-score for sum2012-11.tif\n",
      "sum2012-11.tif saved.\n",
      " \n",
      "writing z-score for sum2012-12.tif\n",
      "sum2012-12.tif saved.\n",
      " \n",
      "writing z-score for sum2013-01.tif\n",
      "sum2013-01.tif saved.\n",
      " \n",
      "writing z-score for sum2013-02.tif\n",
      "sum2013-02.tif saved.\n",
      " \n",
      "writing z-score for sum2013-03.tif\n",
      "sum2013-03.tif saved.\n",
      " \n",
      "writing z-score for sum2013-04.tif\n",
      "sum2013-04.tif saved.\n",
      " \n",
      "writing z-score for sum2013-05.tif\n",
      "sum2013-05.tif saved.\n",
      " \n",
      "writing z-score for sum2013-06.tif\n",
      "sum2013-06.tif saved.\n",
      " \n",
      "writing z-score for sum2013-07.tif\n",
      "sum2013-07.tif saved.\n",
      " \n",
      "writing z-score for sum2013-08.tif\n",
      "sum2013-08.tif saved.\n",
      " \n",
      "writing z-score for sum2013-09.tif\n",
      "sum2013-09.tif saved.\n",
      " \n",
      "writing z-score for sum2013-10.tif\n",
      "sum2013-10.tif saved.\n",
      " \n",
      "writing z-score for sum2013-11.tif\n",
      "sum2013-11.tif saved.\n",
      " \n",
      "writing z-score for sum2013-12.tif\n",
      "sum2013-12.tif saved.\n",
      " \n",
      "writing z-score for sum2014-01.tif\n",
      "sum2014-01.tif saved.\n",
      " \n",
      "writing z-score for sum2014-02.tif\n",
      "sum2014-02.tif saved.\n",
      " \n"
     ]
    }
   ],
   "source": [
    "for mon in lista:\n",
    "    for mean_,stdv_ in zip(mean,stdv):\n",
    "        if mon[-6:-4] == mean_[-6:-4]:\n",
    "            imagem, meta = xuleta.TifToArray(folder+mon)\n",
    "            media, meta = xuleta.TifToArray(folderout+mean_)\n",
    "            desvpad, meta = xuleta.TifToArray(folderout+stdv_)\n",
    "            z = (imagem - media)/desvpad\n",
    "            print 'writing z-score for %s' %mon\n",
    "            xuleta.ArrayToTif(z, mon, folderz, meta, Type=3)\n",
    "        imagem, z, media, desvpad, meta = None, None, None, None, None"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}