ODW_example

import sys
sys.path.append('/Users/toshan/dev/pysal/pysal/weights')
from spintW import ODW
import pysal as ps

With an equal number of origins and destinations (n=16)

origins = ps.weights.lat2W(4,4)
dests = ps.weights.lat2W(4,4)

origins.n

dests.n

ODw = ODW(origins, dests)

print ODw.n, 16*16

256 256

ODw.full()[0].shape

(256, 256)

With non-equal number of origins (n=9) and destinations (m=25)

origins = ps.weights.lat2W(3,3)

dests = ps.weights.lat2W(5,5)

origins.n

dests.n

ODw = ODW(origins, dests)

print ODw.n, 9*25

225 225

ODw.full()[0].shape

(225, 225)