[R-sig-Geo] Resoultion in LAT/LON when using WorldHires
Huse, Irene
irene.huse at imr.no
Sun Jun 17 03:25:14 CEST 2007
Hi R-geo users,
I do hope this is a simple task for you, for me, a beginner, it's very difficult. When I plot my positions (of recaptured fish) in a map, the resolution is completely wrong, and the points "clings up in lines" (horizontal and vertical), and looks very funny. The positions I try to plot are not at all at lines.
library(maps)
library(mapdata)
par(bg="white",mar=c(3,1,1,1))
map('worldHires',xlim=c(30.40,30.65),ylim=c(70.60,70.70),type="n")
rist1 <-read.table("C:\\rist1.txt")
colnames(rist1)<-c("lon","lat")
points(rist1$lon,rist1$lat,col="black")
map('worldHires',xlim=c(30.45,30.33),ylim=c(70.60,70.70),col="peach puff",fill=T,add=T) axis(1,c(30.40,30.45,30.50,30.55,30.60),c("30°24'E","30°27'E","30°30'E","30°33'E","30°36'E"))
axis(2,c(70.60,70.65,70.70),c("70°36'N","70°39'N","70°42'N"),las=1)
The rist1 looks like this:
30.500172222 70.666672222
30.500186111 70.666727778
30.450075 70.650263889
30.466930556 70.650161111
30.466930556 70.650161111
30.483411111 70.650144444
30.483505556 70.650125
30.4836 70.650105556
30.500266667 70.650052778
30.516844444 70.650016667
30.566852778 70.633472222
30.583436111 70.633441667
Etc etc
Is there a setting I have misused?
I would be very grateful for som help!
Yours sincerely,
Irene Huse
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Dr. Irene Huse <><
Forsker / Scientist
Havforskningsinstituttet / Institute of Marine Research
Postboks 1870 Nordnes, N-5817 Bergen, Norway
Office: Nordnesgt. 33
Tlf: +47 55 23 68 22 ><>
Faks/Fax +47 55 23 68 30
Mobil/Cell +47 4125 1661 ><>
Internett: www.imr.no
irene.huse at imr.no <><
More information about the R-sig-Geo
mailing list