I currently have shapefile of US PUMAS and another feature shapefile (call it feature X). I am using ArcPy to project both of these files to the following projection coordinate system (using Project_management):
After projecting, I use intersect_analysis and then CalculateAreas_stat to compute the intersection areas. Finally I output a .csv file, which I then use to compute area weighted means.
The issue is that every time I run my simple script, the intersection areas vary by a small amount (somewhere to the 8th decimal point). These small changes are impacting my analysis.
If this is a GIS precision issue or if I making a mistake, how can I go about fixing it?