佐賀県と福岡県における意匠の登録者 - Rを用いてプロット

R Code

setwd("{.....}/R Projects/Saga-Fukuoka")
### read locations from SQL Server database
library(RODBC)
dc <- odbcConnect("{database}",uid = "{userid}", pwd = "{password}")
## [Location] is a geography column in SQL Server 2008 that contains the coordinates of 1193 design registrants
dfSagaFukuoka <- sqlQuery(dc, "select [Location].STAsText() as [ShapeWKT] FROM [DesignV3Holders].[dbo].[PointDataSet]")
###
library(sp)
## adapted from https://gis.stackexchange.com/questions/57935/reading-ogc-geometry-data-from-sql-server-using-r
library(rgeos)
points <- vector("list",length(dfSagaFukuoka$ShapeWKT))
z = 0
for(line in dfSagaFukuoka$ShapeWKT)
{
   {
     points[[z+1]]<-readWKT(line, id = NULL, p4s = llCRS)
   }
   z = z + 1
}
library(raster)
spSagaFukuoka <- do.call(bind,points)
plot(spSagaFukuoka)
png("SagaFukuokaDesignRegistrants.png")
plot(spSagaFukuoka)
dev.off()

        

Bivand, Pebesma, Gómez-Rubio (2008)