[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