将测深分配给光栅 map

Question

I have a fishing zone between Latitudes 20°N and 26°N and I want to assign bathymetry to each square of the resolution 5 nautical miles (Or 1 nautical mile). I tried the next R code but I get positive bathymetry inside the ocean which is wrong.

marmap::getNOAA.bathy(xmin ,
                      xmax ,
                      ymin ,
                      ymax ,
                      resolution = 5)

Where xmin, xmax, ymin, ymax represent the limits of the zone.

Is there any other solution that you suggest? Thanks

Answer 1

I suppose you reversed longitudes and latitudes. the syntax of getNOAA.bathy is as follows:

getNOAA.bathy(lon1, lon2, lat1, lat2, resolution)

If your longitudes are located to the West of the greenwich meridian, they should also be negative.

Finally, you might get a "could not reach NOAA servers" error message when using getNOAA.bathy() since the NOAA removed the server on which the function used to make its queries up until very recently. We are working on a fix and should be releasing it on github and on CRAN servers in the next few days.