Рассчитать минимальное расстояние от полигона до пространственной точки - PullRequest
Новая
       79

Рассчитать минимальное расстояние от полигона до пространственной точки

0 голосов
/ 27 февраля 2019

У меня проблемы с поиском расстояния между полигонами и пространственными точками.Я пытаюсь измерить расстояние между районами (полигоны) и бывшими оловянными шахтами (пространственные точки).Я подозреваю, что, возможно, что-то не так с проекцией.Кажется, я могу рассчитать расстояние между районами, но не для районов и шахт.

#Open file
mapping<- readOGR(dsn="C:/Users/noble/OneDrive - London School of Economics/LSE/EC465/Extended Essay/Stack Exchange Question/gadm2.Peninsula.dbf")


#Define the districts
Kinta <- mapping[mapping@data$NAME_1 == "Perak" &
             mapping@data$NAME_2 == "Kinta", ]
Gombak <- mapping[mapping@data$NAME_1 == "Selangor" &
             mapping@data$NAME_2 == "Gombak", ]

#Set a projection
EPSG.3375<-"+proj=omerc +lat_0=4 +lonc=102.25 +alpha=323.0257964666666 +k=0.99984 +x_0=804671 +y_0=0 +ellps=GRS80 +units=m +no_defs"
Gombak.km<-spTransform(Gombak,CRS(EPSG.3375))
Kinta.km<-spTransform(Kinta,CRS(EPSG.3375))

#Calculate the distance
gDistance(Kinta.km, Gombak.km, byid=TRUE)

#Open data on tin mines and define as spatial points
tin.01<-read.csv("C:/Users/noble/OneDrive - London School of Economics/LSE/EC465/Extended Essay/Stack Exchange Question/Tin Mine Location_01.csv")
coordinates(tin.01)<-8:9

proj4string(tin.01) <- CRS("+proj=omerc +lat_0=4 +lonc=102.25 +alpha=323.0257964666666 +k=0.99984 +x_0=804671 +y_0=0 +ellps=GRS80 +units=m +no_defs")

#Find distance between district and mines
gDistance(Kinta.km,tin.01,byid=TRUE)

Мой вывод для расстояния между полигонами выглядит правильно: 

> gDistance(Kinta.km, Gombak.km, byid=TRUE)
   59
71 100676

Но мой вывод для расстояния между районами и минами определенно неверен: 

> gDistance(Kinta.km,tin.01,byid=TRUE)
    59
1 661153.5
2 661152.6
3 661153.0
4 661152.7
5 661151.8
6 661152.9
7 661153.1
8 661153.3
9 661153.2

Что я собираюсь сделать, это: 1) рассчитать расстояние между всеми районами Малайзии со всеми бывшими оловянными копями;и 2) извлеките расстояние от ближайшей жестяной шахты до каждого района.

Вот ссылка на данные, которые я использую.Я использую данные GADM для районных полигонов, и я сам вручную кодировал местоположение исторических жестяных шахт.Буду признателен за оказанную помощь.Спасибо.

1 Ответ

0 голосов
/ 27 февраля 2019

Есть два вопроса с вашим текущим подходом.1.) Назначение coordinates не совсем корректно.Это должно быть long / lat, но вы назначаете его как lat / long.2.) Непосредственная настройка CRS текущим способом фактически не изменит точки нужным образом.Сначала необходимо назначить соответствующий CRS long / lat, а затем выполнить действие spTransform.

#Open data on tin mines and define as spatial points
tin.01 <- read.csv("random/Tin Mine Location_01.csv")
coordinates(tin.01) <- c("Longitude","Latitude")  ## Should be `9:8` if you wanted to stick with indexes, but using the names here is generally lower risk.
proj4string(tin.01) <- CRS(proj4string(Kinta)) ## Setting the initial projection to the same one the polygons are using. You should change this if your original data source uses some other known long/lat projection.

tin.km <- spTransform(tin.01,CRS(EPSG.3375)) ## Creating a transformed set of points for the distance calculation.

#Find distance between district and mines
gDistance(Kinta.km,tin.km,byid=TRUE) ## '0' distance means the mine is inside the district.

          59
1 194384.372
2 223773.999
3      0.000
4  36649.914
5 102944.361
6      0.000
7      0.000
8   6246.066
9      0.000
...