WHEN I CHECK MY IMAGE AND SHAPE FILE IN ERDAS, IT IS OVERALYING.
BUT WHEN I WANTED TO OVERLAY IT IN R.
I GET THE FOLLOWING ERROR.
Error in .local(x, y, ...) : extents do not overlap In addition: Warning
message: In intersect(extent(x), extent(y)) : Objects do not overlap
BUT RASTER EXTENTS SHOW UP IN METERS AND POLYGON IS IN LAT/LON
THE FOLLOWING IS MY BOUNDING BOX RANGE
> bbox(raster) min max
s1 571045 1006045
s2 3405030 3788030> bbox(map) min max
x 3663120 3875470
y 4845817 5013605
THE FOLLOWING ARE MY PROJECTIONS
> proj4string(map)[1] "+proj=lcc +lat_1=35.17280444444445 +lat_2=12.472955
> +lat_0=24 +lon_0=80 +x_0=4000000 +y_0=4000000 +datum=WGS84 +units=m +no_defs
> +ellps=WGS84 +towgs84=0,0,0"> proj4string(raster)[1] "+proj=utm +zone=43
> +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
