die Berechnung der minimalen Distanz zwischen einem Punkt und der Küste

Ich versuche, die minimale Entfernung zwischen einem gegebenen Punkt und die Küste. Mein Beispiel ist die Entfernung von Madrid an die Küste:

library(rgeos)
library(maptools)
coast <- readShapeLines("Natural_Earth_quick_start/10m_physical/ne_10m_coastline.shp")
MAD = readWKT("POINT(-3.716667 40.383333)")
gDistance(MAD,coast)
[1] 3.021808

Ich habe Schwierigkeiten zu verstehen, was gDistance() gibt. Die docs sagen, es ist die Einheiten der Projektion. Heißt das, es ist in latlong Grad? Wie kann ich konvertieren, dass zu Kilometer?

  • Können Sie uns einen link zu deiner map? Als start geben Sie proj4string(coast) zu sehen, was die Projektion ist. Wenn es ist ein `+units=' string, dann sind das die Einheiten.
  • Die Karte, die ich verwende, ist die Küste von naturalearthdata.com/downloads/10m-physical-vectors. proj4string(coast) zurück NA
  • Haben Sie schaute auf die definition der Hausdorff-Distanz im Zusammenhang von der Hilfe-Seite?. Ich bekomme die Idee, dass diese Antwort wird in "Grad", aber ich bin überhaupt nicht davon überzeugt, dass sinnvoll ist, da die Anzahl der km in einem "Grad" variiert von etwa 100 am äquator, null an den Polen
InformationsquelleAutor user525602 | 2014-01-22
  1. 12

    gDistance(...) gibt die minimale kartesischen (euklidische) Distanz zwischen dem Punkt und feature-set zur Verfügung gestellt, die als Argumente. Da deine Karte ist im long/lat-Koordinaten, die Sie bekommen Abstand in "Grad", z.B.

    d = sqrt { (lange1 lange2)2 + (lat1 - lat2)2 }

    wo long und lat in decimal degrees. Als darauf hingewiesen wurde, das bedeutet nicht viel, weil die Umwandlung auf planare Abstand (z.B. km) hängt davon ab, wo Sie sind. So brauchen wir für die Transformation von Daten in ein CRS, die ungefähr planare in der region von Interesse. Es stellt sich heraus, dass die entsprechenden CRS für Spanien ist EPSG-2062. Die Projektion string für EPSG-2062 ist:

    +proj=lcc +lat_1=40 +lat_0=40 +lon_0=0 +k_0=0.9988085293 +x_0=600000 +y_0=600000 +a=6378298.3 +b=6356657.142669561 +pm=madrid +units=m +no_defs

    die +units=m (Meter). Also müssen wir reproject sowohl der Punkt (MAD) und die Grenzen zu den EPSG-2062.

    library(rgeos)
library(maptools)

epsg.2062 <- "+proj=lcc +lat_1=40 +lat_0=40 +lon_0=0 +k_0=0.9988085293 +x_0=600000 +y_0=600000 +a=6378298.3 +b=6356657.142669561 +pm=madrid +units=m +no_defs"
wgs.84    <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

coast <- readShapeLines("ne_10m_coastline",CRS(wgs.84))
MAD   <- readWKT("POINT(-3.716667 40.383333)",p4s=CRS(wgs.84))
gDistance(MAD,coast)            # WGS-84 (long/lat) projection, units in "degrees"
# [1] 3.021808

coast.proj <- spTransform(coast,CRS(epsg.2062))
MAD.proj   <- spTransform(MAD,CRS(epsg.2062))
gDistance(MAD.proj,coast.proj)  #EPSG-2062 projection, units are in meters.
# [1] 305171.2

    So dass der minimale Abstand ist ~305.2 km.

    Schließlich, beachten Sie, dass Ihre Küste Datei hat, alle Küsten der Welt, so ist dies der minimale Abstand zu einige Küste, nicht unbedingt der spanischen Küste (obwohl in diesem Fall erweisen sich auf der nördlichen Küste von Spanien). Wenn Ihr Bezugspunkt war sehr in der Nähe von der Grenze mit Portugal, das die nächste Küsten-Punkt sein würde, an der westlichen Küste von Portugal.

    • Dank! das ist genau das, was ich wollte. Woher weißt du, dass es die Nordküste?
    • Die kurze Antwort ist, dass ich schaute auf eine Karte, aber wenn Sie wieder in Madrid sein könnte näher an der Mittelmeer-Küste in der Nähe von Castellon. Es ist ein Weg, um die Koordinaten des nächsten Punktes, aber es ist kompliziert. Wenn Sie brauchen, um zu wissen, dass Sie sollten als eine separate Frage.
    • Das sieht toll aus-aber welche Projektion soll ich verwenden bei der Berechnung der Entfernung der Punkte (in km) zum nächsten Ufer auf einer globalen Skala? Ich habe ~250 Punkte!
    • Diese Antwort nicht verallgemeinern.
    InformationsquelleAutor jlhoward
Schreibe einen Kommentar