[R-sig-Geo] (no subject)
Roger Bivand
Roger.Bivand at nhh.no
Thu Sep 28 13:09:24 CEST 2006
On Thu, 28 Sep 2006, pir2.jv wrote:
>
> Le 28 sept. 06 à 12:13, Roger Bivand a écrit :
>
> > (Reverting to list)
> >
> >
> > So possibly there is an endian issue, or a contaminated shapefile.
> >
> > In a fresh R session, please do:
> >
> > library(maptools)
> > sessionInfo()
> > fn <- system.file("shapes/sids.shp", package="maptools")[1]
> > file.info(fn)
> > xx <- readShapePoly(fn, proj4string=CRS("+proj=longlat
> > +ellps=clrk66"))
> > class(xx)
> > bbox(xx)
> > summary(xx)
> >
> > and post the output to the list?
OK, the file size at least looks correct, and the output of the summary
looks good too. By the way, why the load()?
Could you now try plot(xx, axes=TRUE), and then traceback() immediately
after the error occurs?
Roger
>
> > load("/Users/jver/Desktop/HMP/.RData")
> > library(maptools)
> Le chargement a nécessité le package : foreign
> Le chargement a nécessité le package : sp
>
> > sessionInfo()
> Version 2.3.1 (2006-06-01)
> powerpc-apple-darwin8.6.0
>
> attached base packages:
> [1] "methods" "stats" "graphics" "grDevices" "utils"
> "datasets"
> [7] "base"
>
> other attached packages:
> maptools sp foreign
> "0.5-15" "0.8-20" "0.8-15"
>
> > fn <- system.file("shapes/sids.shp", package="maptools")[1]
> > str(fn)
> chr "/Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp"
>
>
> > file.info(fn)
>
> size
> /Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp 46196
>
> isdir
> /Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp FALSE
>
> mode
> /Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp 644
>
> mtime
> /Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp 2006-09-23 10:15:57
>
> ctime
> /Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp 2006-09-28 11:39:39
>
> atime
> /Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp 2006-09-28 11:39:39
>
> uid
> /Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp 501
>
> gid
> /Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp 80
>
> uname
> /Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp jver
>
> grname
> /Library/Frameworks/R.framework/Versions/2.3/Resources/library/
> maptools/shapes/sids.shp admin
>
> > xx <- readShapePoly(fn, proj4string=CRS("+proj=longlat
> +ellps=clrk66"))
> > class(xx)
> [1] "SpatialPolygonsDataFrame"
> attr(,"package")
> [1] "sp"
> > bbox(xx)
> min max
> r1 -84.32385 -75.45698
> r2 33.88199 36.58965
> > summary(xx)
> Object of class SpatialPolygonsDataFrame
> Coordinates:
> min max
> r1 -84.32385 -75.45698
> r2 33.88199 36.58965
> Is projected: FALSE
> proj4string : [+proj=longlat +ellps=clrk66]
> Data attributes:
> AREA PERIMETER CNTY_
> CNTY_ID NAME
> Min. :0.0420 Min. :0.999 Min. :1825 Min. :1825
> Alamance : 1
> 1st Qu.:0.0910 1st Qu.:1.324 1st Qu.:1902 1st Qu.:1902
> Alexander: 1
> Median :0.1205 Median :1.609 Median :1982 Median :1982
> Alleghany: 1
> Mean :0.1263 Mean :1.673 Mean :1986 Mean :1986
> Anson : 1
> 3rd Qu.:0.1542 3rd Qu.:1.859 3rd Qu.:2067 3rd Qu.:2067
> Ashe : 1
> Max. :0.2410 Max. :3.640 Max. :2241 Max. :2241
> Avery : 1
>
> (Other) :94
> FIPS FIPSNO CRESS_ID BIR74
> SID74
> 37001 : 1 Min. :37001 Min. : 1.00 Min. : 248
> Min. : 0.00
> 37003 : 1 1st Qu.:37050 1st Qu.: 25.75 1st Qu.: 1077 1st
> Qu.: 2.00
> 37005 : 1 Median :37100 Median : 50.50 Median : 2180
> Median : 4.00
> 37007 : 1 Mean :37100 Mean : 50.50 Mean : 3300
> Mean : 6.67
> 37009 : 1 3rd Qu.:37150 3rd Qu.: 75.25 3rd Qu.: 3936 3rd
> Qu.: 8.25
> 37011 : 1 Max. :37199 Max. :100.00 Max. :21588 Max. :
> 44.00
> (Other):94
> NWBIR74 BIR79 SID79 NWBIR79
> Min. : 1.0 Min. : 319 Min. : 0.00 Min. : 3.0
> 1st Qu.: 190.0 1st Qu.: 1336 1st Qu.: 2.00 1st Qu.: 250.5
> Median : 697.5 Median : 2636 Median : 5.00 Median : 874.5
> Mean :1050.8 Mean : 4224 Mean : 8.36 Mean : 1352.8
> 3rd Qu.:1168.5 3rd Qu.: 4889 3rd Qu.:10.25 3rd Qu.: 1406.8
> Max. :8027.0 Max. :30757 Max. :57.00 Max. :11631.0
>
> >
> >
> >
>
> Jacques Vernin
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand at nhh.no
More information about the R-sig-Geo
mailing list