{
 "metadata": {
  "name": "TRMM Month to 8 Day interval_additional code"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import numpy as np\n",
      "import numpy.ma as ma\n",
      "import gdal\n",
      "from gdalconst import *"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 1
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "def hr2_8day(Month,outFile):\n",
      "    ## -- For single Rasters -- ##\n",
      "    # Load Raster 1\n",
      "\n",
      "    dataset_1 = gdal.Open(Month, GA_ReadOnly)\n",
      "    band_1 = dataset_1.GetRasterBand(1)\n",
      "    data_1 = band_1.ReadAsArray()\n",
      "    NA_1 = band_1.GetNoDataValue()\n",
      "\n",
      "    # Do Computation on NumPy Array 1\n",
      "    data_mask_1 = ma.masked_where(data_1 == NA_1, data_1)\n",
      "    data_mask_1 = data_mask_1 * 192\n",
      "    #data_mask_1 = 0.62 * data_mask_1\n",
      "    data_1 = data_mask_1.data\n",
      "\n",
      "    # Set Driver\n",
      "    format = \"GTiff\"\n",
      "    driver = gdal.GetDriverByName( format )\n",
      "    driver.Register()\n",
      "\n",
      "    # Set Metadata for Raster output\n",
      "    cols = dataset_1.RasterXSize\n",
      "    rows = dataset_1.RasterYSize\n",
      "    bands = dataset_1.RasterCount\n",
      "    datatype = band_1.DataType\n",
      "\n",
      "    # Set Projection for Raster\n",
      "    outDataset = driver.Create(outFile, cols, rows, bands, datatype)\n",
      "    geoTransform = dataset_1.GetGeoTransform()\n",
      "    outDataset.SetGeoTransform(geoTransform)\n",
      "    proj = dataset_1.GetProjection()\n",
      "    outDataset.SetProjection(proj)\n",
      "\n",
      "    # Write output to band 1 of new Raster\n",
      "    outBand = outDataset.GetRasterBand(1)\n",
      "    outBand.WriteArray(data_1,0,0)\n",
      "\n",
      "    # Close and finalise newly created Raster\n",
      "    F_M01 = None\n",
      "    outBand = None\n",
      "    outDataset = None\n",
      "    driver = None\n",
      "    geTransform = None\n",
      "    proj = None\n",
      "    cols = None\n",
      "    rows = None\n",
      "    bands = None\n",
      "    datatype = None\n",
      "\n",
      "    dataset_1 = None\n",
      "    band_1 = None\n",
      "    data_1 = None\n",
      "    NA_1 = None\n",
      "\n",
      "    data_mask_1 = None"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 4
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "M01 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20090101.7A.tif'\n",
      "hr2_8day(M01,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009001.tif')\n",
      "hr2_8day(M01,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009009.tif')\n",
      "hr2_8day(M01,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009017.tif')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 6
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "M02 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20090201.7A.tif'\n",
      "hr2_8day(M02,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009033.tif')\n",
      "hr2_8day(M02,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009041.tif')\n",
      "hr2_8day(M02,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009049.tif')\n",
      "M03 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20090301.7A.tif'\n",
      "hr2_8day(M03,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009065.tif')\n",
      "hr2_8day(M03,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009073.tif')\n",
      "hr2_8day(M03,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009081.tif')\n",
      "M04 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20090401.7A.tif'\n",
      "hr2_8day(M04,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009097.tif')\n",
      "hr2_8day(M04,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009105.tif')\n",
      "hr2_8day(M04,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009113.tif')\n",
      "M05 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20090501.7A.tif'\n",
      "hr2_8day(M05,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009121.tif')\n",
      "hr2_8day(M05,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009129.tif')\n",
      "hr2_8day(M05,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009137.tif')\n",
      "M06 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20090601.7A.tif'\n",
      "hr2_8day(M06,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009153.tif')\n",
      "hr2_8day(M06,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009161.tif')\n",
      "hr2_8day(M06,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009169.tif')\n",
      "M07 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20090701.7A.tif'\n",
      "hr2_8day(M07,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009185.tif')\n",
      "hr2_8day(M07,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009193.tif')\n",
      "hr2_8day(M07,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009201.tif')\n",
      "M08 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20090801.7A.tif'\n",
      "hr2_8day(M08,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009217.tif')\n",
      "hr2_8day(M08,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009225.tif')\n",
      "hr2_8day(M08,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009233.tif')\n",
      "M09 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20090901.7A.tif'\n",
      "hr2_8day(M09,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009249.tif')\n",
      "hr2_8day(M09,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009257.tif')\n",
      "hr2_8day(M09,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009265.tif')\n",
      "M10 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20091001.7A.tif'\n",
      "hr2_8day(M10,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009281.tif')\n",
      "hr2_8day(M10,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009289.tif')\n",
      "hr2_8day(M10,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009297.tif')\n",
      "M11 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20091101.7A.tif'\n",
      "hr2_8day(M11,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009305.tif')\n",
      "hr2_8day(M11,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009313.tif')\n",
      "hr2_8day(M11,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009321.tif')\n",
      "M12 = r'I:\\Data\\TRMM 3B43\\2009_tif\\resampled_extract\\rspr3B43.20091201.7A.tif'\n",
      "hr2_8day(M12,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009337.tif')\n",
      "hr2_8day(M12,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009345.tif')\n",
      "hr2_8day(M12,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009353.tif')\n",
      "#hr2_8day(M01,r'I:\\Data\\TRMM 3B43\\2009_tif\\convert_juliandate\\TRMM2009361.tif')\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 8
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}