-
Notifications
You must be signed in to change notification settings - Fork 13
Open
Description
I am writing tiff using the following scripts:
begin
ga = GeoArray(rand(360, 150))
bbox!(ga, (min_x=-180.0, min_y=-60.0, max_x=180.0, max_y=90.0))
epsg!(ga, 4326) # in WGS84
GeoArrays.write!("test.tif", ga)
endWhen reading in R, the crsTransform is error. The latitude range is changed to [-59, 91].
library(raster)
b <- brick("test.tif")Warning message:
In .rasterFromGDAL(x, band = band, objecttype, ...) :
data seems flipped. Consider using: flip(x, direction='y')
b
class : RasterBrick
dimensions : 150, 360, 54000, 1 (nrow, ncol, ncell, nlayers)
resolution : 1, 1 (x, y)
extent : -180, 180, -59, 91 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +no_defs
source : N:/Research/GEE_repos/GEE-latest/test.tif
names : layer Metadata
Metadata
Assignees
Labels
No labels