refactoring everything
something is f*!k@d with the polygon index. see line 42, and the plots that come out. abpolyID!=ab index?
new rstudio added something here