意匠の時空間的分析

基礎な手法を用いて日本における登録意匠の時空間的分析などを探検する目的で、このブログによって、手法などについてシェアしたり、意見交換したりします。

福岡市付近における意匠の登録者

f:id:ttcrossroads09:20181011160747p:plain

 

佐賀県と福岡県における意匠の登録者:境界を含むプロット - 意匠の時空間的分析の続きです。

 


## use extents
# Change these parameters
scale.parameter = 0.25  # scaling parameter. less than 1 is zooming in, more than 1 zooming out. 
xshift = -0.05  # Shift to right in map units. 
yshift = 0.1  # Shift to left in map units. 
original.bbox = spSagaFukuoka@bbox  # Pass bbox of your Spatial* Object. 

# Just copy-paste the following
edges = original.bbox

edges[1, ] <- (edges[1, ] - mean(edges[1, ])) * scale.parameter + mean(edges[1, 
    ]) + xshift
edges[2, ] <- (edges[2, ] - mean(edges[2, ])) * scale.parameter + mean(edges[2, 
    ]) + yshift

png(file="DesignRegistrantsFukuokaCityArea.png", width=600, height=480)
plot(spSagaFukuoka, axes=TRUE, pch=1,xlim = edges[1, ], ylim = edges[2, ],
main="Design Registrants in Fukuoka City Area, 2001-2017" )
plot(saga, add=TRUE,xlim = edges[1, ], ylim = edges[2, ], border="red")
plot(fukuoka, add=TRUE,xlim = edges[1, ], ylim = edges[2, ], border="blue")
dev.off()