Hi,
Our office is working with GIS maps using python and we've
hit a snag. Due
to the way the GIS software behaves, it will take a land
parcel that shares
a border with another parcel and call them just one parcel,
instead of 2.
The land parcels are represented by 2-dimensional polygons.
Our solution is to shrink or scale the parcel down by a unit
of measurement
so that the 2 parcels no longer share a line. Currently we
use a Python
version of the "inpoly" algorithm listed on this
page:
http://www.visibone.
com/inpoly/. It works great until we hit some of our
weird parcels that have 700 vectors.
Since we have to move each vector in four directions and
then test each
vector with the inpoly algorithm, this means that we do 4 X
700 or 2800
calculations. We also have to loop over the number of
vectors, so that means
700 X 2800 which takes too long as this is for a web app.
Does anyone have an alternate algorithm for this kind of
thing? Thank you
for any help!
Mike Driscoll
Applications Specialist
MCIS - Technology Center
1 E. Main St.
Marshalltown, IA 50158
mdriscoll co.marshall.ia.us
_______________________________________________
Image-SIG maillist - Image-SIGpython.org
htt
p://mail.python.org/mailman/listinfo/image-sig
|
