到这一天结束时，这是解决问题的最佳方法：

extensions [gis]

globals

[

land-use-map

lotid-patch-map

def-risk-map

market-dist-map

]

patches-own

[

land-use

lotid-patch

def-risk

market-dist

]

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

;; load the maps ;;

;;;;;;;;;;;;;;;;;;;;;

to load-gis

clear-all

set land-use-map gis:load-dataset "area2_lu_black.asc" ;loads the land use map

set lotid-patch-map gis:load-dataset "area2_lot.asc" ;loads the lots map

set def-risk-map gis:load-dataset "area2_risk.asc" ;loads the deforestation risk map

set market-dist-map gis:load-dataset "area2_mkt.asc" ;loads the distance from markets map

gis:set-world-envelope-ds gis:envelope-of land-use-map ;sets the envelope of the world to match that of the GIS dataset

gis:apply-raster land-use-map land-use ;patches in the land-use-map have a specific land-use now

gis:apply-raster lotid-patch-map lotid-patch ;patches in the lot-id-map have a specific lot-id now

gis:apply-raster def-risk-map def-risk ;patches in the def-risk-map have a specific def-risk now

gis:apply-raster market-dist-map market-dist ;patches in the market-dist-map have a specific market-dist now

ask patches [

if land-use = 1 [ set pcolor 64 ] ; Green = Forest

if land-use = 2 [ set pcolor 14 ] ; Dark red = Agriculture

if land-use = 3 [ set pcolor 45 ] ; Yellow = Reforestation

]

let view gis:load-dataset "area2.shp" ;load a shapefile of the properties

gis:set-world-envelope-ds gis:envelope-of view

foreach gis:feature-list-of view

[

gis:set-drawing-color white ;draws the line of the shapefile

gis:draw ? 1

]

end