address nan's in pslg

[krober10nd]

• Fixes NaN-delimited polygons failure  #8
• Simple switch from np.amin and np.amax to np.nanmin and np.nanmax.
• Need to still add a unit test.

[krober10nd]

Okay here's a complex NaN-delimited polygon to test with.

import numpy as np
import matplotlib.pyplot as plt
 
from inpoly import inpoly2
 
nan = np.nan
 
points = np.loadtxt("qp.txt")
node = np.loadtxt("poly.txt")
edge = np.loadtxt("e.txt", dtype=int)
 
IN, ON = inpoly2(points, node, edge)
 
fig, ax = plt.subplots()
plt.plot(points[IN == 1, 0], points[IN == 1, 1], "b.")
plt.plot(points[IN == 0, 0], points[IN == 0, 1], "r.")
plt.plot(points[ON == 1, 0], points[ON == 1, 1], "ms")
 
ax.set_aspect("equal", adjustable="box")
plt.show()

This should produce (if the changes in this PR are not used the answer is incorrect).

qp.txt

poly.txt

e.txt

[krober10nd]

@dengwirda do you have comments on this?