Marmap getNOAA.bathy- NOAA 서버가 아래로 이동 했습니까?

Question

MARMAP에서 getNOAA.bathy를 사용하여 해저의 수심지도를 만들려고합니다.

scotia <- getNOAA.bathy(lon1=-42, lon2=-34, lat1=-51, lat2=-55, resolution = 10, keep=TRUE) 

나는 성공적으로 11 월에 높은 해상도로 그것을 실행으로 코딩이 올바른지 알아,하지만 지금 얻을 모든 오류입니다 ... 내가 시도

Querying NOAA database ... 
This may take seconds to minutes, depending on grid size 
Error in getNOAA.bathy(lon1 = -42, lon2 = -34, lat1 = -51, lat2 = -55, : 
The NOAA server cannot be reached 

솔루션은 here 언급 그러나, 그것은 나를 위해 작동하지 않았다. 나는 NOAA가 12 월 (as it says here)에 자신의 웹 사이트를 업그레이드했지만 R 코드와 연결하는 방법을 잘 모르겠다. 어떤 아이디어?

Answer 1

그래서 나는 NOAA에서 매우 도움이되는 사람들을 통해 내 자신의 질문에 대한 답을 찾을 수있었습니다. 웹 주소가 더 이상 내가 그것을 getNOAA.bathy에 다른 무언가를 호출하는 데 필요한 것을 깨닫게 나에게 잠시했다,하지만 업데이트 기능은 이제 here입니다

을 일하지 않도록

NOAA는 페이지를 이동했다.

getNOAA.bathy <- 
function(lon1,lon2,lat1,lat2, resolution = 4, keep=FALSE, antimeridian=FALSE) { 

x1=x2=y1=y2 = NULL 

if (lon1 < lon2) {lon1->x1 ; lon2->x2} else {lon1->x2 ; lon2->x1} 
if (lat1 < lat2) {lat1->y1 ; lat2->y2} else {lat1->y2 ; lat2->y1} 

res = resolution * 0.016666666666666667 

fetch <- function(x1,y1,x2,y2,res) { 
    WEB.REQUEST <- paste("https://gis.ngdc.noaa.gov/cgi-bin/public/wcs/etopo1.xyz?filename=etopo1.xyz&request=getcoverage&version=1.0.0&service=wcs&coverage=etopo1&CRS=EPSG:4326&format=xyz&resx=", res, "&resy=", res, "&bbox=", x1, ",", y1, ",", x2, ",", y2, sep = "") 
    dat <- suppressWarnings(try(read.table(WEB.REQUEST),silent=TRUE)) 
    return(dat) 
} 

# Naming the file 
if (antimeridian) { 
    FILE <- paste("marmap_coord_",x1,";",y1,";",x2,";",y2,"_res_",resolution,"_anti",".csv", sep="") 
} else { 
    FILE <- paste("marmap_coord_",x1,";",y1,";",x2,";",y2,"_res_",resolution,".csv", sep="") 
} 

# If file exists in the working directory, load it, 
if(FILE %in% list.files()) { 
    cat("File already exists ; loading \'", FILE,"\'", sep="") 
    read.bathy(FILE, header=T) -> exisiting.bathy 
    return(exisiting.bathy) 
} else { # otherwise, fetch it on NOAA server 

    if (antimeridian) { 

     l1 <- x2 ; l2 <- 180 ; l3 <- -180 ; l4 <- x1 

     cat("Querying NOAA database ...\n") 
     cat("This may take seconds to minutes, depending on grid size\n") 
     left <- fetch(l1,y1,l2,y2,res) 
     right <- fetch(l3,y1,l4,y2,res) 

     if (is(left,"try-error")|is(right,"try-error")) { 
      stop("The NOAA server cannot be reached\n") 
     } else { 
      cat("Building bathy matrix ...\n") 
      left <- as.bathy(left) ; left <- left[-nrow(left),] 
      right <- as.bathy(right) 
      rownames(right) <- as.numeric(rownames(right)) + 360 
      bath2 <- rbind(left,right) 
      class(bath2) <- "bathy" 
      bath <- as.xyz(bath2) 
     } 

    } else { 

     cat("Querying NOAA database ...\n") 
     cat("This may take seconds to minutes, depending on grid size\n") 
     bath <- fetch(x1,y1,x2,y2,res) 
     if (is(bath,"try-error")) { 
      stop("The NOAA server cannot be reached\n") 
     } else { 
      cat("Building bathy matrix ...\n") 
      bath2 <- as.bathy(bath) 
     } 
    } 

    if (keep) { 
     write.table(bath, file=FILE, sep=",", quote=FALSE, row.names=FALSE) 
    } 

    return(bath2) 
} 
}