Function PIP(byval X as single , byval Y as single) as byte
  ' index is the polygon point array count
  ' fencelon() and fencelatt() global arrays of the x,y coords    -  fencelon is the "Y" axis
  ' x,y are the test point coords - x=latt y=lon

   local f1 as single, f2 as single,f3 as single,f4 as single,f5 as single  ,f6 as single
   local cn as byte
   local i as byte, j as byte    ,q as byte,r as byte


   i=1  :   cn = 0
 index = index-1
   j=index
   for i = 1 to index
      q=0  :    r=0
      if  fencelon(i)< Y And fencelon(J)>= Y then q = 1
      if  fencelon(J)< Y And fencelon(i)>= Y then r=1
      If q=1 or r=1 then
         f1 = Y - fencelon(i) :  f2 = fencelon(j) - fencelon(i)  :  f3 = fencelatt(j) -  fencelatt(i)
         f4 = f1/f2
         f5 = f4 * f3
         f6 = fencelatt(i) + f5
         if f6 < X Then  cn = cn+1
      End If
      j=i
   next i
 index = index+1

   cn = cn mod 2
   if cn > 0 then cn =1  ' even numbers are outside the polygon

   PIP = cn
End Function


'-----------------------------------------------------------------------------------------------------------------------------------------
'-----------------------------------------------------------------------------------------------------------------------------------------

sub PIPtestKML(byval  OffSet as  single, byval resolution as single)
'decimal place / grid resolution / offset
'2        0.01        1.1132 km
'3        0.001        111.32 m
'4        0.0001        11.132 m
'5        0.00001        1.1132 m
'go easy with these - easy to create many millions of points, google earth will have a hissy fit.
' check the terminal program for max number of displayed lines.
'have been using YAT terminal for this test, - need more than 80 columns and can change the buffered displayed lines easy.

'find  max lat and min latt, max lon min lon from stored fencelatt,lon array and add (offset) to them
'this test wont work accross the equator or the prime meridian

   local n as byte    ,m as byte
   local maxlon as single, minlon as single, maxlatt as single, minlatt as single  , original as single

 'set the focus point for initial Earth view and kml header stuff
   print  "<?xml version=";chr(34);"1.0";chr(34);" encoding=";chr(34);"UTF-8";chr(34);"?>"
   print  "<kml xmlns=";chr(34);"http://www.opengis.net/kml/2.2";chr(34);
   print " xmlns:gx=";chr(34);"http://www.google.com/kml/ext/2.2";chr(34);">"
   print  "<Document>"
   print  "<LookAt> <longitude>";
   print fencelon(1);
   print  "</longitude> <latitude>";
   print fencelatt(1);
   print  "</latitude> <range>1300.000000</range> </LookAt>"
   print  "<Folder>"

  ''' gegerate the FENCEPOINTS    /  POLYPOINTS
   for n = 1 to index
      print  "<Placemark>"
      print "<name>";"FP";n;"</name>"
      print"<Point> <coordinates>"
      print fencelon(n);
      print ",";
      print fencelatt(n);
      print  "</coordinates></Point></Placemark>  "
   next n

 'find max and mins then add the offset value to generate the test box
   maxlatt = fencelatt(1) :   minlatt = fencelatt(1):   maxlon = fencelon(1) :  minlon = fencelon(1)
   for n = 1 to index
      if maxlatt <= fencelatt(n) then maxlatt = fencelatt(n)
      if minlatt  >= fencelatt(n) then minlatt  = fencelatt(n)
      if maxlon <= fencelon(n)  then maxlon = fencelon(n)
      if minlon  >=  fencelon(n) then minlon = fencelon(n)
   next n
   minlatt = minlatt - offset
   maxlatt = maxlatt + offset
   minlon = minlon - offset
   maxlon = maxlon + offset

  ''' run the test grid through the point in polygon(pip) function and generate placemark pins
   '' original = minlatt
   while minlatt < maxlatt and minlon < maxlon    'this sets the test boundry box - it's just easier to do a test square
      original = minlatt
      while minlatt < maxlatt   '   print the x axis
         m =  pip(minlatt, minlon)
         if m< 1 then    ' just display the outside points
            print  "<Placemark>"
            print "<name>";m;"</name>"
            print" <Point><coordinates>"
            print minlon;
            print ",";
            print minlatt;
            print  "</coordinates></Point></Placemark>"
         end if
         minlatt = minlatt + resolution
      wend
      minlatt = original
      minlon = minlon + resolution
   wend
   print  "</Folder></Document></kml>"

   do
   loop   '  stop here -  need time to save to file in terminal
end sub

sub   preload()
   fencelatt(1)=-37.702547
   fencelatt(2)=-37.702413
   fencelatt(3)=-37.701884
   fencelatt(4)=-37.702027
   fencelatt(5)=-37.702547    'start repeat


   fencelon(1)=147.840474
   fencelon(2)=147.839278
   fencelon(3)=147.839317
   fencelon(4)=147.840353
   fencelon(5)=147.840474   'start repeat

   index = 5

end sub

 

Function PIP(byval X as single , byval Y as single) as byte

index = index-1
J = index