is a Python script for calculating the Multi-order Adjacent order index MAI of urban construction land using raster data The script uses the GDAL library to read and write raster data and numpy and ma

e period. The output is a plot of the MAI values over time.

Here is an example of how to use the script:

1. Install the required libraries: GDAL, numpy, and matplotlib. You can do this using pip:

pip install gdal numpy matplotlib

2. Save the script to a file named 'mai.py'.

3. Prepare your raster data. You will need three raster files: one for the original urban construction land, and one for each time period of new urban construction land. The raster files should be in the same projection and have the same pixel size.

4. Open a Python shell or create a new Python script. Import the mai module:

import mai

5. Call the mai.calculate_mai() function, passing in the paths to the raster data files and the list of new urban construction land file names:

mai.calculate_mai('original.tif', ['new_2000.tif', 'new_2010.tif', 'new_2020.tif'])

This will calculate the MAI values for the three time periods and plot them using matplotlib.

Note: The script assumes that the pixel values in the raster data represent land use categories, with the value 1 indicating urban construction land. If your data uses different values, you may need to modify the script.