nz.r

library(tidyverse)
library(terra)

tas <- vect("E:/geodata/Aus_Coastline/tasmania/csttascd_r.shp")
bclm <- rast(list.files("E:/geodata/bioclim/2_5m/",full.names=TRUE))

bclm <- crop(bclm,tas)
ext <- extract(bclm,tas)

ext <- ext[complete.cases(ext),]
ext <- ext[,-1]
pc <- prcomp(ext)

bclm2 <- rast(list.files("E:/geodata/bioclim/2_5m/",full.names=TRUE))
nz <- vect("nz.shp")
bclm_nz <- crop(bclm2,nz)

ext_nz <- extract(bclm2,nz)
ext_nz <- ext_nz[complete.cases(ext_nz),]
ext_nz <- ext_nz[,-1]
pc_nz <- predict(pc,newdata=ext_nz)
pc_dist <- sqrt(pc_nz[,1]^2 + pc_nz[,2]^2)


ext_nz <- extract(bclm2,nz,xy=TRUE)
ext_nz$prdist = pc_dist

tst <- vect(ext_nz,geom=c("x","y"),crs="EPSG:4326")
plot(tst,"prdist")