[R-sig-Geo] Even simpler TM_WORLD_BORDERS_SIMPLE (AND SP???)

Jim Burke j.burke at earthlink.net
Thu Jan 29 23:15:06 CET 2009


OK I am on decaf coffee now. I reinstalled R and
all is well now.

Now Hans, below is what I get when I run your code.
Is this expected? Quick and certainly DIRTY.

 > library("maptools")
Loading required package: foreign
 > load(url("http://spatial.nhh.no/R/etc/TM_WORLD_BORDERS_SIMPL-0.2.RData"))
 > wrld_simpl <- wrld_simpl[wrld_simpl$NAME != "Antarctica",]
Warning message:
closing unused connection 3 
(gzcon(http://spatial.nhh.no/R/etc/TM_WORLD_BORDERS_SIMPL-0.2.RData))
 > for (i in (1:length(wrld_simpl at polygons))) {
+ polyCnt <- length(wrld_simpl at polygons[[i]]@Polygons)
+ if(polyCnt > 2) {
+ pO <- sapply(wrld_simpl at polygons[[i]]@Polygons, function(x) 
{slot(x,'area')>0.1})
+ new_pO <- pO[ wrld_simpl at polygons[[i]]@plotOrder ]
+ wrld_simpl at polygons[[i]]@plotOrder <- 
wrld_simpl at polygons[[i]]@plotOrder[ new_pO ]
+ }
+ }
 > save(wrld_simpl,'SIMPLEMAP.Rdata')
Error in save(wrld_simpl, "SIMPLEMAP.Rdata") :
object ‘SIMPLEMAP.Rdata’ not found




Thanks,
Jim Burke

Jim Burke wrote:
> Thanks Hans,
>
> Its been a day or so since I have used R graphics.
> Now sp seems to lack grid. I cannot run my tried and
> true graphics workspace at all.
>
> Are you folks having any trouble? What am I doing
> wrong?
>
> Below
> 1. try to load library sp. No go.
> 2. download sp from Berkeley (US1). Different errors
>
>  > library(sp)
> Error : package 'grid' does not have a name space
> Error: package/namespace load failed for 'sp'
>  > chooseCRANmirror()
>  > utils:::menuInstallPkgs()
> trying URL 
> 'http://cran.cnr.Berkeley.edu/bin/windows/contrib/2.8/sp_0.9-29.zip'
> Content type 'application/zip' length 845890 bytes (826 Kb)
> opened URL
> downloaded 826 Kb
>
> package 'sp' successfully unpacked and MD5 sums checked
>
> The downloaded packages are in
>         C:\Documents and Settings\Jim Burke\Local 
> Settings\Temp\RtmpqCMAZS\downloaded_packages
> updating HTML package descriptions
> Warning message:
> dependency 'grid' is not available
>  >
>
> Thanks,
> Jim Burke
>
>
>
> Hans-Jörg Bibiko wrote:
>   
>> On 29.01.2009, at 16:49, Jonathan Callahan wrote:
>>     
>>> Is anyone aware of an even more simplified version of the data that went
>>> into creating TM_WORLD_BORDERS_SIMPL-0.2.RData so that we might have
>>> TM_WORLD_BORDERS_LOWRES.RData for use cases like my own.
>>>       
>> A very quick solution would be to filter out all irrelevant small 
>> polygons (islands).
>> [to reduce the number of polygon points would lead to very unshaped 
>> boundaries, I mean]
>>
>> What is an irrelevant polygon? My first approach is if a country 
>> consists of more than 2 polygons then draw only those which are 
>> greater than a given area let's say 0.1 .
>>
>> To achieve that without recreating the entire SpatialPolygonsDataFrame 
>> one could simply change the plotOrder of each PolygonsList and save it.
>>
>> Like that quick&dirty script:
>>
>> library("maptools")
>> load(url("http://spatial.nhh.no/R/etc/TM_WORLD_BORDERS_SIMPL-0.2.RData"))
>> wrld_simpl <- wrld_simpl[wrld_simpl$NAME != "Antarctica",]
>> for (i in (1:length(wrld_simpl at polygons))) {
>>   polyCnt <- length(wrld_simpl at polygons[[i]]@Polygons)
>>   if(polyCnt > 2) {
>>     pO <- sapply(wrld_simpl at polygons[[i]]@Polygons, function(x) 
>> {slot(x,'area')>0.1})
>>     new_pO <- pO[ wrld_simpl at polygons[[i]]@plotOrder ]
>>     wrld_simpl at polygons[[i]]@plotOrder <- 
>> wrld_simpl at polygons[[i]]@plotOrder[ new_pO ]
>>   }
>> }
>> save(wrld_simpl,'SIMPLEMAP.Rdata')
>>
>>
>>
>> --Hans
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>>
>>     
>
>
> 	[[alternative HTML version deleted]]
>
>   
> ------------------------------------------------------------------------
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>




More information about the R-sig-Geo mailing list