adehabitat бумасы үшін орын деректерінің пішімі

Менде осы форматта файл бар:

ASCII пішімі

Бірінші жолдар келесідей:

ncols 1440
nrows 720
xllcorner -180.0
yllcorner -90
cellsize 0.25
NODATA_value -9999

Мен негізінен 1440 'плиткамен' x бағытта (бойлық) және 720 'тақта' бағыты бойынша (ендік) әлем бар. Әрбір «плитка» - ұзындығы 0,25 градус болатын шаршы. Менің ойымша, xllcorner және yllcorner дұрыс. Мен бұл картаны R:

library("adehabitat")
bio1 <- import.asc("D:/ENFA/data.asc")
maps <- as.kasc(list(data = bio1))
image(maps, col = cm.colors(256), clfac = list(Aspect = cl))

Карта жақсы көрінеді.

Adehabitat бумасын пайдаланып, кейбір экологиялық экологиялық факторларды талдауды (ENFA) орындағым келеді және орын туралы деректерге қатысты сенімді емеспін. Негізінде, мен оларды бойлықтар мен ендік ретінде ұстап тұрамын, бірақ сол кезде «қаңылтыр индексі» ретінде пайда бола аламын (мысалы, төменгі сол жақ бұрышы ендік -90 және бойлық -180 болғандықтан, «қаңылтыр индексі» 0, 0 болса керек)? ). Орналасқан жері туралы деректер пішімі қандай? Мен ENFA кодын келесідей қолданамын:

locs <- read.table("D:/ENFA/Locs.txt", header = TRUE, sep="\t")
dataenfa1 <- data2enfa(maps, locs)
pc <- dudi.pca(dataenfa1$tab, scannf = FALSE)
enfa1 <- enfa(pc, dataenfa1$pr,scannf = FALSE)
hist(enfa1)

Мен кез келген ескертулерді бағалаймын. Алдын-ала рақмет.

2

2 жауаптар

Сіздің координаталарыңызды ұзақ уақытқа қалдырып кету проблемасы, жер бетіндегі көп жерлерде ұзындықтың ендік ендік дәрежесінен басқа ұзындыққа ие. Бұл басқаларға қатысты кейбір бағыттардағы қашықтықты асқақтатып, ENFA-ны бұрмалауы мүмкін.

Әсіресе сіздің деректеріңіз салыстырмалы түрде кішігірім аймақтан болса, мен W/E x-осі бойымен метрлердегі координаталарды және У/У-у осінің қайтадан білдірілуін ұсынамын. Егер барлық ұпайларыңыз бір UTM аймағы ішінде болса, онда rgdal бумасында project() арқылы R

Бір мысал: алынды мұнда:

library(rgdal)

# Make a two-column matrix, col1 = long, col2 = lat
xy <- cbind(c(118, 119), c(10, 50))
# Convert it to UTM coordinates (in units of meters)
project(xy, "+proj=utm +zone=51 ellps=WGS84")
          [,1]    [,2]
[1,] -48636.65 1109577
[2,] 213372.05 5546301 

Кеңістіктік деректерді қалай басқаруға болатынын білу үшін Bivand, Pebesma және Gomez-Rubio қолданбаларымен «Қолданбалы кеңістіктік деректерді талдау» бөлімінде қол жетімді. Егер сізге нақты көмек қажет болса, R-sig-Geo тарату тізімі .

Бұл көмектеседі деп үміттенемін.

7
қосылды
Джошқа рахмет. Бұл мәселе әлемдегі кез-келген жерде болуы мүмкін (яғни, 60 УТМ бойлық аймағынан артық) болуы мүмкін. Бұл туралы не істеу керектігін білмеймін ...
қосылды автор cs0815, көзі
Рахмет, мен бұл туралы ойлайтынмын.
қосылды автор cs0815, көзі
@Josh - бұны түзете аласыз ба: «бойлықтың ұзақтығы бойлықтың ұзақтығынан өзгеше ұзындыққа ие» Мен бұл бұрмаланулардың мәнін әлдеқайда жақсырақ түсінемін және менің талдауыма әсер етуі мүмкін. Рахмет!
қосылды автор cs0815, көзі
Балама түрде, сіз өзіңізді конверсиялауды жүзеге асыра аласыз. Ендік дәрежесі - 111 км. lat радда көрсетілген болса, бойлықтың дәрежесі cos (lat) * 111 km-ге жақын. Градустарда градустарда cos (((2 * pi)/360) * lat) орнына пайдаланыңыз * 111 . Оны бірнеше ендік үшін қолданып көріңіз және берілген мәндермен салыстырыңыз [ en.wikipedia.org/ wiki/Longitude # Degree_length]] {{Бұл Wikipedia мақаласында] бұл жақсы жұмыс істейтінін көру үшін. Бұл деректер жиынтығы көптеген градустарға дейін созылмағанша (немесе полюстердің жанында үлкен аумақтарды қосқанда) жұмыс істеуі керек. Егер бұл пайдалы болса, мен оны негізгі жауапқа қосамын.
@csetzkorn - Аулау үшін рахмет. Негізгі мәселе - экватордан алыстағанда, дәреженің ұзақтығы қысқа және қысқа болады. Полюсте сіз бірнеше қадаммен 360 градусқа өтуіңізге болады; экваторға 40 000 км жол жүруге тура келеді.
қосылды автор Josh O'Brien, көзі
Кітапқа сілтеме үшін рахмет. Мен осы салаға көбірек кірісуді ұсындым.
қосылды автор John Colby, көзі

Мүмкін сіз координаталарды түрлендіргіңіз келуі мүмкін

GHAM (жаһандық, иерархиялық, алфавиттік-цифрлық және мортон-кодталған)

кез келген lat/lon-де альфа-сандық мекен-жайы бар, ол сұрыпталып қалады.

GHAM: сұрыптауға жарамды ықшам жаһандық геокоде авторефераты. Дункан Агнеу :

GHAM коды - географиялық орналасуды белгілеу әдісі   олардың лауазымдарына. Ол шектелген тең аймақ ұяшықтары үшін мекенжайларды анықтайды   тұрақты ендік пен бойлық арқылы, еркін дәлдікпен.   Ұяшық кодтары Мортонды рекурсивтеуге тапсырыс беру арқылы анықталады   16-16 торға бөлу, нәтижесінде алынған кодтар кодталады   әріп-санды жұптар. Белгіленген пункттердің лексикалық түрлерінің тізбесі болады   көршілерге жақындастырыңыз (әдетте) бірге; әртүрлі сынақтар   жаһандық деректер жиынтығы көптеген жағдайларда нақты ең жақын нүкте екенін көрсетеді   тізімде 50%, ал 5 жазба ішінде 80%   уақыт.

Бастапқы код - IAMG репозиторийі, бірақ егер мүмкін болмаса оған қол жеткізетініне сенімдімін.

1
қосылды
Қызықты көрінеді, бірақ ол осы тұрғыда қолданылмайды деп ойлаймын
қосылды автор cs0815, көзі