I need to create a script that determines core habitat areas. Take a vegetation polygon, with a total area of $260$ hectares. From an ecological perspective, not all areas of the polygon are of equal interest since some species prefer the $$core$$ areas away from the outside edge. We can remove the edge buffer from the polygon to obtain just the core areas. Using a $50-$meter buffer, we end up with a core area of only $120$ hectares. The core habitat is therefore $120$ hectares, or $46\mathrm{.}2\%$ of the original area.

The script needs to determine the core area for every polygon in a feature class, both in terms of total area and in terms of the percentage. The specific script requirements are as follows: hardcode the following values at the top of the script: Input features $-$ A polygon feature class, which can be shapefile or geodatabase feature class. Input edge buffer $-$ A numerical value $($positive integer$)$ that defines the buffer to be eliminated from each vegetation polygon to obtain the core area. It is assumed this value will be in the units of the coordinate system of the input features so there is no need for units.

The result of the script should be two new fields added to the existing polygon feature class. $$ The first field should be called COREAREA and should be of type Double. This should contain the calculated core area, in units of the coordinate system of the data.

The second field should be called COREPERC and should be of type Double. This should contain the core area of a polygon as a percentage of the total area of that polygon. Logically, this number will always be less than $100($since the buffer distance should be greater than $0)$ and greater than or equal to $0.$

Other requirements for the script:

The script should work on projected data only. If the coordinate system is geographic, the script should not perform the calculation, but instead report an error message saying that $$the core area could not be calculated since the input data are not projected$.$

The output should consist of the same polygon feature class $($not a copy$),$ with only two fields added. All the calculations should be performed by the script.

The script should ideally not write any intermediate data to disk, for example to a default workspace. All the processing and calculations should be performed by the script in memory.

Finally please no AI use in the example$/$solution provided