Velika količina kvantitativnih geoprostornih podataka prikuplja se i spaja u diskretne popisne jedinice (npr. zemlje ili države). Glatka piknofilaktička interpolacija ima za cilj pronaći glatku, nenegativnu funkciju tako da je integral površine nad svakom popisnom jedinicom jednak agregiranim podacima. Konvencionalno, glatka piknofilaktička interpolacija dobiva se algoritmom staničnog automata koji pretvara po dijelovima konstantnu funkciju u približno glatku funkciju definiranu na mreži koordinata na karti u ekvivalentnoj projekciji. Alternativni pristup, kojeg je predložio Tobler (1976), je konstruirati kartografsku projekciju s ujednačavanjem gustoće u kojoj su područja popisnih jedinica proporcionalna agregiranim podacima. Piknofilaktička interpolacija može se dobiti iz Jakobijana te projekcije. U ovom radu opisujemo primjenu softvera za tu metodu. Iako rješenja nisu nužno optimalna u smislu unaprijed definiranih kvantitativnih mjera glatkoće, naša metoda je računski učinkovita i potencijalno se može koristiti zajedno s drugim metodama da bi se ubrzala konvergencija prema optimalnom rješenju.