migurski

var LambertProjection = function(lat0, lon0, lat1, lat2, xmin, ymin, xmax, ymax)
{
    var pi = Math.PI, ln = Math.log, pow = Math.pow,
        sin = Math.sin, cos = Math.cos, tan = Math.tan,
        atan = Math.atan, sqrt = Math.sqrt;
    
    function sec(t) { return 1 / cos(t); }
    function cot(t) { return 1 / tan(t); }

    function deg2rad(deg) { return pi * deg / 180; }

migurski

from sys import argv
from time import time
from glob import glob
from os import stat, kill, getuid
from os.path import basename, dirname, join
from datetime import datetime
from random import choice
from signal import SIGTERM

if __name__ == '__main__':

migurski

"""

>>> usa = USA()

>>> loc = Location(48.38544, -124.72916)
>>> pt = usa.locationProj(loc)
>>> abs(-2109470 - pt.x) < 1
True
>>> abs( 1532790 - pt.y) < 1
True

migurski

from time import time
from urlparse import urlparse
from httplib import HTTPConnection
from os.path import basename
from cStringIO import StringIO
from datetime import timedelta
from os import SEEK_SET, SEEK_CUR, SEEK_END

class RemoteFileObject:
    """ Implement enough of this to be useful:

migurski

from sys import argv, stdout, stderr
from subprocess import Popen
from os import stat, unlink
from math import log, ceil
from stat import ST_SIZE
from time import time

from ModestMaps.OpenStreetMap import Provider
from ModestMaps.Core import Coordinate

migurski

''' Extract route relations from osm2pgsql-generated database to CSV.
'''
from sys import argv
from csv import writer
from psycopg2 import connect

if __name__ == '__main__':
    
    db = connect(database='raster_osm', user='osm').cursor()
    

migurski

'''

ogr2ogr -overwrite -t_srs '+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs' PG:"user=nhd dbname=nhd" NHDFlowline.shp

create table nhdflowline_snapped ( ogc_fid integer primary key );

select addgeometrycolumn('', 'nhdflowline_snapped', 'geometry', 900914, 'LINESTRING', 2);
create index nhdflowline_snapped_geom_idx on nhdflowline_snapped using gist (geometry);

insert into nhdflowline_snapped

migurski

'''

#!/bin/sh -x

apt-get -y update
apt-get -y install htop osmosis git python-pip python-shapely python-pyproj python-networkx qhull-bin

git clone git://github.com/migurski/Skeletron.git ~/Skeletron
cd ~/Skeletron
python setup.py install

migurski

from itertools import product
from multiprocessing import Pool

from numpy import zeros, ubyte
from psycopg2 import connect
from osgeo import gdal, osr
from PIL import Image

def measure_quad(x, y, cellsize):
    '''

migurski

#!/usr/bin/env python
''' Convert binary data to UTF8.

If you want binary data in Javascript over XMLHttpRequest but don't
have ArrayBuffer available, one option is to read strings of text and
use xhr.responseText.charCodeAt() to read them. Because this expects to
work with text, it will return the unicode of the character at that
position. This script converts raw binary bytes to unicode text that
evaluates to the same 8-bit unsigned int at each position. The overall
length will increase slightly.