Habitat model for marten Martes americana




Дата канвертавання25.04.2016
Памер15.89 Kb.

HABIT@ APL Scripts



Marten

; Habitat model for marten (Martes americana)

; In RMLANDS workspace, use 10 800 RUNLOOP 'marten.txt'

; Full version includes roads, but for HRV there aren’t any roads, so we’ve simplified this model


gridpath = 'z:\RMLANDS\SJNF\hrv_a1\run1\elk' ; Location of zipped input grids (uses same grid as elk)

BLOCK 140 140 36 ; 1100x1100 blocks, with 69 cell edge buffer


covcond = READ 'covcond.asc' ; read grids
lri1 = (GET 'mcover.txt') CATEGORY covcond ; lri1 = cover type

SETWINDOW 50 10 ; logistic function determines window size

lri2 = (GET 'medge.txt') 'LOGISTIC' (50 10 .5) (50 10 .6) (50 10 .8) BADEDGE covcond ; lri2 = adverse edge effects

lra = lri1 * lri2 ; lra = product of lris


hrc = 892 ('MEAN' 'L') COMBINE lra ; hrc1 = lra over homerange

SAVE 'hrc'

END


Three-toed Woodpecker

; Habitat model for Three-toed Woodpecker (Picoides tridactylus)

; 13 May 2003, 2 Mar 2004

; In RMLANDS workspace, use 10 800 RUNLOOP 'TTWO.TXT'

; 8 March 2004, Modifications to include insect disturbances
gridpath = 'd:\habit@\input\ttwo'
BLOCK 1100 1100 14 ; 1100x1100 blocks, with 14 cell edge buffer

covercond = READ 'covcond.asc' ; read grids

covercondt1 = -10 READ 'covcond.asc' ; covercond for step t-1

agewfire = READ 'wfashm.asc'

agesbeetleh = READ 'sbashm.asc'

agesbeetlel = READ 'sbaslm.asc'

agefbeetleh = READ 'fbashm.asc'

agefbeetlel = READ 'fbaslm.asc'


nodist = (GET 'ttwocover.txt') CATEGORY covercond ; undisturbed based on cover type

withdist = (GET 'ttwocover.txt') CATEGORY covercondt1 ; burned based on cover type in step t-1


dist1 = 1.33 * withdist ; it's always better after a disturbance

sbh = (-9999 0, 0 1, 10 .75, 20 0) ARBFN agesbeetleh

sbl = (-9999 0, 0 .75, 10 .25, 20 0) ARBFN agesbeetlel

fbh = (-9999 0, 0 1, 10 .75, 20 0) ARBFN agefbeetleh

fbl = (-9999 0, 0 .75, 10 .25, 20 0) ARBFN agefbeetlel

fire = (-9999 0, 0 1, 10 1, 20 0) ARBFN agewfire ; <10 y since burn or bugs is optimal;

; >20 y has no additional value
dist = MAX (dist1 * fire) (dist1 * sbh) (dist1 * sbl) (dist1 * fbh) (dist1 * fbl)
hrc1 = 125 COMBINE dist ; burned/bug-killed hrc

hrc2 = 350 COMBINE nodist ; undisturbed hrc (larger home range)

hrc = MAX hrc1 hrc2 ; hrc is best of disturbed/undisturbed
SAVE 'hrc' ; save files with year attached to names

END



Olive-sided Flycatcher

; Habitat model for Olive-Sided Flycatcher (Contopus cooperi)

; "Quick, three beers!"

; In RMLANDS workspace, use 10 800 RUNLOOP 'osfl.txt'


;test = 1 ; uncomment this line to do a one-block time test

gridpath = 'z:\RMLANDS\SJNF\hrv_a1\run1\osfl' ; Location of zipped input grids

BLOCK 1100 1100 6 ; 1100x1100 blocks, with 6 cell edge buffer
covcond = READ 'covcond.asc' ; read grids

cond = READ 'cond.asc'


f1 = (GET 'osflforage.txt') CATEGORY covcond

edgemax = THRESHOLD 100 -20

f2 = (GET 'osfledge.txt') (1 4 1, 1 3 .75, 1 2 .05, 2 4 .75, 2 3 .5, 3 4 .1) 'LOGISTIC' (100 -20) GOODEDGE cond

forage = MIN f1 f2


nest = (GET 'osflnest.txt') CATEGORY covcond
hrcf = 150 ('MEAN' 'S') COMBINE forage ; food based on squared distance-weighted mean

hrcn = 150 ('AREA' .5 .1 .6 .75 .1 1) COMBINE nest ; at least 10% of nest > .75 for prime nesting;

; or 10% better than .5 for value of 0.6

hrc = GEOMEAN hrcf hrcn


SAVE 'hrc'

END



Elk

; Habitat model for Elk (Cervus elaphus)

; In RMLANDS workspace, use 10 800 RUNLOOP 'elk.txt'

gridpath = 'z:\RMLANDS\SJNF\hrv_a1\run1\elk' ; Location of zipped input grids

BLOCK 300 300 69 ; 1100x1100 blocks, with 69 cell edge buffer
covcond = READ 'covcond.asc' ; read grids

slope = 1 READ 'slope.asc'


forage1 = (GET 'elkforage.txt') CATEGORY covcond ; cover types for foraging

steep = 25 -3 .8 LOGISTIC slope ; like shallow slopes

forage = forage1 * steep
cov1 = (GET 'elkcover.txt') CATEGORY covcond ; cover type for cov

cov = cov1 * steep


SETWINDOW 300 -40

forage_r = forage * 300 -40 LOGISTIC JUXTAPOSE cov >= .75 ; realized forage must be near cover

SETWINDOW 500 -70

cov_r = cov * 500 -70 LOGISTIC JUXTAPOSE forage >= .5 ; realized cover must be near forage


hrc1 = 1725 ('AREA' .5 .66 1) COMBINE forage_r ; at least 2/3 of forage must be > 0.5

hrc2 = 1725 ('AREA' .75 .5 1 .85 .25 1) COMBINE cov_r ; at least half of cov must be > 0.75, or

; a quarter better than 0.85.

hrc = GEOMEAN hrc1 hrc2 ; geometric mean


SAVE 'hrc'

END




База данных защищена авторским правом ©shkola.of.by 2016
звярнуцца да адміністрацыі

    Галоўная старонка