This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#!/usr/bin/env python | |
""" | |
Code extracts MODIS derived products (1km) from the HDF, rescales them, checks values against | |
acceptable QA and outputs to a binary file. Code seems generic enough that one can just | |
add another method for a new MODIS product. Currently works on LAI, NDVI/EVI and LST. Code | |
also writes out an ENVI .hdr file...nice. | |
*** NOTE...I don't vouch for the reliability of the default QA, please supply your own on | |
the cmd line. |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
""" | |
Zonal Statistics | |
Vector-Raster Analysis | |
Copyright 2013 Matthew Perry | |
Usage: | |
zonal_stats.py VECTOR RASTER | |
zonal_stats.py -h | --help | |
zonal_stats.py --version |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
def fleiss_kappa(M): | |
""" | |
See `Fleiss' Kappa <https://en.wikipedia.org/wiki/Fleiss%27_kappa>`_. | |
:param M: a matrix of shape (:attr:`N`, :attr:`k`) where `N` is the number of subjects and `k` is the number of categories into which assignments are made. `M[i, j]` represent the number of raters who assigned the `i`th subject to the `j`th category. | |
:type M: numpy matrix | |
""" | |
N, k = M.shape # N is # of items, k is # of categories | |
n_annotators = float(np.sum(M[0, :])) # # of annotators |