danclewley

#!/usr/bin/env python
#
# A script topographically correct SAR data using RIOS applier.
#
# Dan Clewley (clewley@usc.edu) - 05/02/2013
#
# The correction uses:
#
# sigma0_norm = sigma0 * (A_flat / A_slope) * (cos(theta_ref) / cos(theta_loc))^n
#

danclewley

import sys
from rios import applier
from rios import cuiprogress
import numpy

def calcstats(info, inputs, outputs):
    # Calc standard deviation of bands for each pixel
    stdev_pixels = numpy.std(inputs.inimage, axis=0)
    # Reshape to make 3-dimensional array (required for output)
    outputs.outimage = stdev_pixels.reshape((1,stdev_pixels.shape[0],stdev_pixels.shape[1]))

danclewley

import sys
from rios import applier
from rios import cuiprogress
import numpy

def calcstats(info, inputs, outputs):
    # Calc standard deviation of bands for each pixel
    stdev_pixels = numpy.std(inputs.inimage, axis=0)
    # Reshape to make 3-dimensional array (required for output)
    outputs.outimage = stdev_pixels.reshape((1,stdev_pixels.shape[0],stdev_pixels.shape[1]))

danclewley

#!/usr/bin/env python
"""
A script to create a GIF from a sequence of images using
ImageMagick.
Applies natural sorting to file names.
Dan Clewley (http://spectraldifferences.wordpress.com)
"""

import argparse
import re

danclewley

#!/usr/bin/env python
# -*- coding: utf-8 -*-
"""
Mosaic files within a zip archive using RSGISLib

Author: Dan Clewley
Creation Date: 02/09/2015

Designed to mosaic Environemnt Agency LiDAR-derived DTM and DSM files
downloaded in ASCII format from http://environment.data.gov.uk/ds/survey

danclewley

#!/usr/bin/env python

#############################################
# ZonalStats.py
#
# A script to calculate zonal statistics
# using RSGISLib python bindings
# Requires RSGISLib > Mercurial version 588
#
# Dan Clewley (daniel.clewley@gmail.com) 18/08/2013

danclewley

DAY=`date +%Y%m%d`
curl http://soilscape.usc.edu/webmapping/tonzi/tonzi_data_last7days.zip > $1/tonzi_$DAY.zip
curl http://soilscape.usc.edu/webmapping/newhogan/newhogan_data_last7days.zip > $1/newhogan_$DAY.zip
curl http://soilscape.usc.edu/webmapping/canton/canton_data_last7days.zip > $1/canton_$DAY.zip
curl http://soilscape.usc.edu/webmapping/matthaeigardens/matthaeigardens_data_last7days.zip > $1/matthaeigardens_$DAY.zip
curl http://soilscape.usc.edu/webmapping/terradoro/terradoro_data_last7days.zip > $1/terradoro_$DAY.zip
curl http://soilscape.usc.edu/webmapping/blm/blm_data_last7days.zip > $1/blm_$DAY.zip

danclewley

<?xml version="1.0" encoding="UTF-8" ?>
<!--
    Co-register two images
-->

<rsgis:commands xmlns:rsgis="http://www.rsgislib.org/xml/">

    <!-- Generate GCPs -->
    <rsgis:command algor="registration" option="basic" 
        reference="ref_image.kea" 

danclewley

# A script to convert PALSAR data to dB using the RIOS applier.
# Dan Clewley (daniel.clewley@gmail.com) - 05/02/2013

import sys
from rios import applier
from rios import cuiprogress
import numpy as np

def converttodB(info, inputs, outputs):
    """

danclewley

"""
Example of reading the spectra from a single pixel using
arsf_envi_reader (https://github.com/pmlrsg/arsf_tools)

Dan Clewley, NERC-ARF-DAN
2018-05-01
"""
from arsf_envi_reader import numpy_bin_reader
from arsf_envi_reader import envi_header
from matplotlib import pyplot as plt