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非常乐于助人的人找到了我自己问题的答案。

NOAA确实移动了网页，因此网址不再有效

更新的功能现在是here但是，我花了一段时间才意识到我需要将它称为与getNOAA.bathy不同的东西。

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)
}
}

Marmap getNOAA.bathy- NOAA服务器移动/关闭？

1 个答案: