{"id":262,"date":"2012-07-31T04:21:08","date_gmt":"2012-07-31T01:21:08","guid":{"rendered":"http:\/\/ekqvist.goeuropeinfo.com\/rbloggerqvist\/?p=262"},"modified":"2012-07-31T04:21:08","modified_gmt":"2012-07-31T01:21:08","slug":"finding-centroid-from-shape-file","status":"publish","type":"post","link":"https:\/\/science.ekqvist.fi\/blogi\/r_maps_spatial\/finding-centroid-from-shape-file\/","title":{"rendered":"Finding Centroid point from shape file and labelling map"},"content":{"rendered":"<p><html><br \/>\n<head><br \/>\n<meta http-equiv=\"Content-Type\" content=\"text\/html; charset=utf-8\"><\/p>\n<p><title>Finding Centroid point from shape file and labelling map<\/title><\/p>\n<p><base target=\"_blank\"\/><\/p>\n<style type=\"text\/css\">\nbody, td {\n   font-family: sans-serif;\n   background-color: white;\n   font-size: 12px;\n   margin: 8px;\n}<\/p>\n<p>tt, code, pre {\n   font-family: 'DejaVu Sans Mono', 'Droid Sans Mono', 'Lucida Console', Consolas, Monaco, monospace;\n}<\/p>\n<p>h1 { \n   font-size:2.2em;\n}<\/p>\n<p>h2 { \n   font-size:1.8em; \n}<\/p>\n<p>h3 { \n   font-size:1.4em; \n}<\/p>\n<p>h4 { \n   font-size:1.0em; \n}<\/p>\n<p>h5 { \n   font-size:0.9em; \n}<\/p>\n<p>h6 { \n   font-size:0.8em; \n}<\/p>\n<p>a:visited {\n   color: rgb(50%, 0%, 50%);\n}<\/p>\n<p>pre {\t\n   margin-top: 0;\n   max-width: 60%;\n   border: 1px solid #ccc;\n}<\/p>\n<p>pre code {\n   display: block; padding: 0.5em;\nmax-width: 60%;\n}<\/p>\n<p>code.r {\n   background-color: #F8F8F8;\nmax-width: 60%;\n}<\/p>\n<p>table, td, th {\n  border: none;\n}<\/p>\n<p>blockquote {\n   color:#666666;\n   margin:0;\n   padding-left: 1em;\n   border-left: 0.5em #EEE solid;\n}<\/p>\n<p>hr {\n   height: 0px;\n   border-bottom: none;\n   border-top-width: thin;\n   border-top-style: dotted;\n   border-top-color: #999999;\n}<\/p>\n<p>@media print {\n   * { \n      background: transparent !important; \n      color: black !important; \n      filter:none !important; \n      -ms-filter: none !important; \n   }<\/p>\n<p>   body { \n      font-size:12pt; \n      max-width:80%; \n   }<\/p>\n<p>   a, a:visited { \n      text-decoration: underline; \n   }<\/p>\n<p>   hr { \n      visibility: hidden;\n      page-break-before: always;\n   }<\/p>\n<p>   pre, blockquote { \n      padding-right: 1em; \n      page-break-inside: avoid; \n   }<\/p>\n<p>   tr, img { \n      page-break-inside: avoid; \n   }<\/p>\n<p>   img { \n      max-width: 100% !important; \n   }<\/p>\n<p>   @page :left { \n      margin: 15mm 20mm 15mm 10mm; \n   }<\/p>\n<p>   @page :right { \n      margin: 15mm 10mm 15mm 20mm; \n   }<\/p>\n<p>   p, h2, h3 { \n      orphans: 3; widows: 3; \n   }<\/p>\n<p>   h2, h3 { \n      page-break-after: avoid; \n   }\n}<\/p>\n<\/style>\n<p><!-- Styles for R syntax highlighter --><\/p>\n<style type=\"text\/css\">\n   pre .operator,\n   pre .paren {\n     color: rgb(104, 118, 135)\n   }<\/p>\n<p>   pre .literal {\n     color: rgb(88, 72, 246)\n   }<\/p>\n<p>   pre .number {\n     color: rgb(0, 0, 205);\n   }<\/p>\n<p>   pre .comment {\n     color: rgb(76, 136, 107);\n   }<\/p>\n<p>   pre .keyword {\n     color: rgb(0, 0, 255);\n   }<\/p>\n<p>   pre .identifier {\n     color: rgb(0, 0, 0);\n   }<\/p>\n<p>   pre .string {\n     color: rgb(3, 106, 7);\n   }\n<\/style>\n<p><!-- R syntax highlighter --><br \/>\n<script type=\"text\/javascript\">\nvar hljs=new function(){function m(p){return p.replace(\/&\/gm,\"&amp;\").replace(\/<\/gm,\"&lt;\")}function f(r,q,p){return RegExp(q,\"m\"+(r.cI?\"i\":\"\")+(p?\"g\":\"\"))}function b(r){for(var p=0;p<r.childNodes.length;p++){var q=r.childNodes[p];if(q.nodeName==\"CODE\"){return q}if(!(q.nodeType==3&#038;&#038;q.nodeValue.match(\/\\s+\/))){break}}}function h(t,s){var p=\"\";for(var r=0;r<t.childNodes.length;r++){if(t.childNodes[r].nodeType==3){var q=t.childNodes[r].nodeValue;if(s){q=q.replace(\/\\n\/g,\"\")}p+=q}else{if(t.childNodes[r].nodeName==\"BR\"){p+=\"\\n\"}else{p+=h(t.childNodes[r])}}}if(\/MSIE [678]\/.test(navigator.userAgent)){p=p.replace(\/\\r\/g,\"\\n\")}return p}function a(s){var r=s.className.split(\/\\s+\/);r=r.concat(s.parentNode.className.split(\/\\s+\/));for(var q=0;q<r.length;q++){var p=r[q].replace(\/^language-\/,\"\");if(e[p]){return p}}}function c(q){var p=[];(function(s,t){for(var r=0;r<s.childNodes.length;r++){if(s.childNodes[r].nodeType==3){t+=s.childNodes[r].nodeValue.length}else{if(s.childNodes[r].nodeName==\"BR\"){t+=1}else{if(s.childNodes[r].nodeType==1){p.push({event:\"start\",offset:t,node:s.childNodes[r]});t=arguments.callee(s.childNodes[r],t);p.push({event:\"stop\",offset:t,node:s.childNodes[r]})}}}}return t})(q,0);return p}function k(y,w,x){var q=0;var z=\"\";var s=[];function u(){if(y.length&#038;&#038;w.length){if(y[0].offset!=w[0].offset){return(y[0].offset<w[0].offset)?y:w}else{return w[0].event==\"start\"?y:w}}else{return y.length?y:w}}function t(D){var A=\"<\"+D.nodeName.toLowerCase();for(var B=0;B<D.attributes.length;B++){var C=D.attributes[B];A+=\" \"+C.nodeName.toLowerCase();if(C.value!==undefined&#038;&#038;C.value!==false&#038;&#038;C.value!==null){A+='=\"'+m(C.value)+'\"'}}return A+\">\"}while(y.length||w.length){var v=u().splice(0,1)[0];z+=m(x.substr(q,v.offset-q));q=v.offset;if(v.event==\"start\"){z+=t(v.node);s.push(v.node)}else{if(v.event==\"stop\"){var p,r=s.length;do{r--;p=s[r];z+=(\"<\/\"+p.nodeName.toLowerCase()+\">\")}while(p!=v.node);s.splice(r,1);while(r<s.length){z+=t(s[r]);r++}}}}return z+m(x.substr(q))}function j(){function q(x,y,v){if(x.compiled){return}var u;var s=[];if(x.k){x.lR=f(y,x.l||hljs.IR,true);for(var w in x.k){if(!x.k.hasOwnProperty(w)){continue}if(x.k[w] instanceof Object){u=x.k[w]}else{u=x.k;w=\"keyword\"}for(var r in u){if(!u.hasOwnProperty(r)){continue}x.k[r]=[w,u[r]];s.push(r)}}}if(!v){if(x.bWK){x.b=\"\\\\b(\"+s.join(\"|\")+\")\\\\s\"}x.bR=f(y,x.b?x.b:\"\\\\B|\\\\b\");if(!x.e&#038;&#038;!x.eW){x.e=\"\\\\B|\\\\b\"}if(x.e){x.eR=f(y,x.e)}}if(x.i){x.iR=f(y,x.i)}if(x.r===undefined){x.r=1}if(!x.c){x.c=[]}x.compiled=true;for(var t=0;t<x.c.length;t++){if(x.c[t]==\"self\"){x.c[t]=x}q(x.c[t],y,false)}if(x.starts){q(x.starts,y,false)}}for(var p in e){if(!e.hasOwnProperty(p)){continue}q(e[p].dM,e[p],true)}}function d(B,C){if(!j.called){j();j.called=true}function q(r,M){for(var L=0;L<M.c.length;L++){if((M.c[L].bR.exec(r)||[null])[0]==r){return M.c[L]}}}function v(L,r){if(D[L].e&#038;&#038;D[L].eR.test(r)){return 1}if(D[L].eW){var M=v(L-1,r);return M?M+1:0}return 0}function w(r,L){return L.i&#038;&#038;L.iR.test(r)}function K(N,O){var M=[];for(var L=0;L<N.c.length;L++){M.push(N.c[L].b)}var r=D.length-1;do{if(D[r].e){M.push(D[r].e)}r--}while(D[r+1].eW);if(N.i){M.push(N.i)}return f(O,M.join(\"|\"),true)}function p(M,L){var N=D[D.length-1];if(!N.t){N.t=K(N,E)}N.t.lastIndex=L;var r=N.t.exec(M);return r?[M.substr(L,r.index-L),r[0],false]:[M.substr(L),\"\",true]}function z(N,r){var L=E.cI?r[0].toLowerCase():r[0];var M=N.k[L];if(M&#038;&#038;M instanceof Array){return M}return false}function F(L,P){L=m(L);if(!P.k){return L}var r=\"\";var O=0;P.lR.lastIndex=0;var M=P.lR.exec(L);while(M){r+=L.substr(O,M.index-O);var N=z(P,M);if(N){x+=N[1];r+='<span class=\"'+N[0]+'\">'+M[0]+\"<\/span>\"}else{r+=M[0]}O=P.lR.lastIndex;M=P.lR.exec(L)}return r+L.substr(O,L.length-O)}function J(L,M){if(M.sL&&e[M.sL]){var r=d(M.sL,L);x+=r.keyword_count;return r.value}else{return F(L,M)}}function I(M,r){var L=M.cN?'<span class=\"'+M.cN+'\">':\"\";if(M.rB){y+=L;M.buffer=\"\"}else{if(M.eB){y+=m(r)+L;M.buffer=\"\"}else{y+=L;M.buffer=r}}D.push(M);A+=M.r}function G(N,M,Q){var R=D[D.length-1];if(Q){y+=J(R.buffer+N,R);return false}var P=q(M,R);if(P){y+=J(R.buffer+N,R);I(P,M);return P.rB}var L=v(D.length-1,M);if(L){var O=R.cN?\"<\/span>\":\"\";if(R.rE){y+=J(R.buffer+N,R)+O}else{if(R.eE){y+=J(R.buffer+N,R)+O+m(M)}else{y+=J(R.buffer+N+M,R)+O}}while(L>1){O=D[D.length-2].cN?\"<\/span>\":\"\";y+=O;L--;D.length--}var r=D[D.length-1];D.length--;D[D.length-1].buffer=\"\";if(r.starts){I(r.starts,\"\")}return R.rE}if(w(M,R)){throw\"Illegal\"}}var E=e[B];var D=[E.dM];var A=0;var x=0;var y=\"\";try{var s,u=0;E.dM.buffer=\"\";do{s=p(C,u);var t=G(s[0],s[1],s[2]);u+=s[0].length;if(!t){u+=s[1].length}}while(!s[2]);if(D.length>1){throw\"Illegal\"}return{r:A,keyword_count:x,value:y}}catch(H){if(H==\"Illegal\"){return{r:0,keyword_count:0,value:m(C)}}else{throw H}}}function g(t){var p={keyword_count:0,r:0,value:m(t)};var r=p;for(var q in e){if(!e.hasOwnProperty(q)){continue}var s=d(q,t);s.language=q;if(s.keyword_count+s.r>r.keyword_count+r.r){r=s}if(s.keyword_count+s.r>p.keyword_count+p.r){r=p;p=s}}if(r.language){p.second_best=r}return p}function i(r,q,p){if(q){r=r.replace(\/^((<[^>]+>|\\t)+)\/gm,function(t,w,v,u){return w.replace(\/\\t\/g,q)})}if(p){r=r.replace(\/\\n\/g,\"<br \/>\")}return r}function n(t,w,r){var x=h(t,r);var v=a(t);var y,s;if(v){y=d(v,x)}else{return}var q=c(t);if(q.length){s=document.createElement(\"pre\");s.innerHTML=y.value;y.value=k(q,c(s),x)}y.value=i(y.value,w,r);var u=t.className;if(!u.match(\"(\\\\s|^)(language-)?\"+v+\"(\\\\s|$)\")){u=u?(u+\" \"+v):v}if(\/MSIE [678]\/.test(navigator.userAgent)&&t.tagName==\"CODE\"&&t.parentNode.tagName==\"PRE\"){s=t.parentNode;var p=document.createElement(\"div\");p.innerHTML=\"<\/p>\n<pre><code>\"+y.value+\"<\/code><\/pre>\n<p>\";t=p.firstChild.firstChild;p.firstChild.cN=s.cN;s.parentNode.replaceChild(p.firstChild,s)}else{t.innerHTML=y.value}t.className=u;t.result={language:v,kw:y.keyword_count,re:y.r};if(y.second_best){t.second_best={language:y.second_best.language,kw:y.second_best.keyword_count,re:y.second_best.r}}}function o(){if(o.called){return}o.called=true;var r=document.getElementsByTagName(\"pre\");for(var p=0;p<r.length;p++){var q=b(r[p]);if(q){n(q,hljs.tabReplace)}}}function l(){if(window.addEventListener){window.addEventListener(\"DOMContentLoaded\",o,false);window.addEventListener(\"load\",o,false)}else{if(window.attachEvent){window.attachEvent(\"onload\",o)}else{window.onload=o}}}var e={};this.LANGUAGES=e;this.highlight=d;this.highlightAuto=g;this.fixMarkup=i;this.highlightBlock=n;this.initHighlighting=o;this.initHighlightingOnLoad=l;this.IR=\"[a-zA-Z][a-zA-Z0-9_]*\";this.UIR=\"[a-zA-Z_][a-zA-Z0-9_]*\";this.NR=\"\\\\b\\\\d+(\\\\.\\\\d+)?\";this.CNR=\"\\\\b(0[xX][a-fA-F0-9]+|(\\\\d+(\\\\.\\\\d*)?|\\\\.\\\\d+)([eE][-+]?\\\\d+)?)\";this.BNR=\"\\\\b(0b[01]+)\";this.RSR=\"!|!=|!==|%|%=|&#038;|&#038;&#038;|&#038;=|\\\\*|\\\\*=|\\\\+|\\\\+=|,|\\\\.|-|-=|\/|\/=|:|;|<|<<|<<=|<=|=|==|===|>|>=|>>|>>=|>>>|>>>=|\\\\?|\\\\[|\\\\{|\\\\(|\\\\^|\\\\^=|\\\\||\\\\|=|\\\\|\\\\||~\";this.ER=\"(?![\\\\s\\\\S])\";this.BE={b:\"\\\\\\\\.\",r:0};this.ASM={cN:\"string\",b:\"'\",e:\"'\",i:\"\\\\n\",c:[this.BE],r:0};this.QSM={cN:\"string\",b:'\"',e:'\"',i:\"\\\\n\",c:[this.BE],r:0};this.CLCM={cN:\"comment\",b:\"\/\/\",e:\"$\"};this.CBLCLM={cN:\"comment\",b:\"\/\\\\*\",e:\"\\\\*\/\"};this.HCM={cN:\"comment\",b:\"#\",e:\"$\"};this.NM={cN:\"number\",b:this.NR,r:0};this.CNM={cN:\"number\",b:this.CNR,r:0};this.BNM={cN:\"number\",b:this.BNR,r:0};this.inherit=function(r,s){var p={};for(var q in r){p[q]=r[q]}if(s){for(var q in s){p[q]=s[q]}}return p}}();hljs.LANGUAGES.cpp=function(){var a={keyword:{\"false\":1,\"int\":1,\"float\":1,\"while\":1,\"private\":1,\"char\":1,\"catch\":1,\"export\":1,virtual:1,operator:2,sizeof:2,dynamic_cast:2,typedef:2,const_cast:2,\"const\":1,struct:1,\"for\":1,static_cast:2,union:1,namespace:1,unsigned:1,\"long\":1,\"throw\":1,\"volatile\":2,\"static\":1,\"protected\":1,bool:1,template:1,mutable:1,\"if\":1,\"public\":1,friend:2,\"do\":1,\"return\":1,\"goto\":1,auto:1,\"void\":2,\"enum\":1,\"else\":1,\"break\":1,\"new\":1,extern:1,using:1,\"true\":1,\"class\":1,asm:1,\"case\":1,typeid:1,\"short\":1,reinterpret_cast:2,\"default\":1,\"double\":1,register:1,explicit:1,signed:1,typename:1,\"try\":1,\"this\":1,\"switch\":1,\"continue\":1,wchar_t:1,inline:1,\"delete\":1,alignof:1,char16_t:1,char32_t:1,constexpr:1,decltype:1,noexcept:1,nullptr:1,static_assert:1,thread_local:1,restrict:1,_Bool:1,complex:1},built_in:{std:1,string:1,cin:1,cout:1,cerr:1,clog:1,stringstream:1,istringstream:1,ostringstream:1,auto_ptr:1,deque:1,list:1,queue:1,stack:1,vector:1,map:1,set:1,bitset:1,multiset:1,multimap:1,unordered_set:1,unordered_map:1,unordered_multiset:1,unordered_multimap:1,array:1,shared_ptr:1}};return{dM:{k:a,i:\"<\/\",c:[hljs.CLCM,hljs.CBLCLM,hljs.QSM,{cN:\"string\",b:\"'\\\\\\\\?.\",e:\"'\",i:\".\"},{cN:\"number\",b:\"\\\\b(\\\\d+(\\\\.\\\\d*)?|\\\\.\\\\d+)(u|U|l|L|ul|UL|f|F)\"},hljs.CNM,{cN:\"preprocessor\",b:\"#\",e:\"$\"},{cN:\"stl_container\",b:\"\\\\b(deque|list|queue|stack|vector|map|set|bitset|multiset|multimap|unordered_map|unordered_set|unordered_multiset|unordered_multimap|array)\\\\s*<\",e:\">\",k:a,r:10,c:[\"self\"]}]}}}();hljs.LANGUAGES.r={dM:{c:[hljs.HCM,{cN:\"number\",b:\"\\\\b0[xX][0-9a-fA-F]+[Li]?\\\\b\",e:hljs.IMMEDIATE_RE,r:0},{cN:\"number\",b:\"\\\\b\\\\d+(?:[eE][+\\\\-]?\\\\d*)?L\\\\b\",e:hljs.IMMEDIATE_RE,r:0},{cN:\"number\",b:\"\\\\b\\\\d+\\\\.(?!\\\\d)(?:i\\\\b)?\",e:hljs.IMMEDIATE_RE,r:1},{cN:\"number\",b:\"\\\\b\\\\d+(?:\\\\.\\\\d*)?(?:[eE][+\\\\-]?\\\\d*)?i?\\\\b\",e:hljs.IMMEDIATE_RE,r:0},{cN:\"number\",b:\"\\\\.\\\\d+(?:[eE][+\\\\-]?\\\\d*)?i?\\\\b\",e:hljs.IMMEDIATE_RE,r:1},{cN:\"keyword\",b:\"(?:tryCatch|library|setGeneric|setGroupGeneric)\\\\b\",e:hljs.IMMEDIATE_RE,r:10},{cN:\"keyword\",b:\"\\\\.\\\\.\\\\.\",e:hljs.IMMEDIATE_RE,r:10},{cN:\"keyword\",b:\"\\\\.\\\\.\\\\d+(?![\\\\w.])\",e:hljs.IMMEDIATE_RE,r:10},{cN:\"keyword\",b:\"\\\\b(?:function)\",e:hljs.IMMEDIATE_RE,r:2},{cN:\"keyword\",b:\"(?:if|in|break|next|repeat|else|for|return|switch|while|try|stop|warning|require|attach|detach|source|setMethod|setClass)\\\\b\",e:hljs.IMMEDIATE_RE,r:1},{cN:\"literal\",b:\"(?:NA|NA_integer_|NA_real_|NA_character_|NA_complex_)\\\\b\",e:hljs.IMMEDIATE_RE,r:10},{cN:\"literal\",b:\"(?:NULL|TRUE|FALSE|T|F|Inf|NaN)\\\\b\",e:hljs.IMMEDIATE_RE,r:1},{cN:\"identifier\",b:\"[a-zA-Z.][a-zA-Z0-9._]*\\\\b\",e:hljs.IMMEDIATE_RE,r:0},{cN:\"operator\",b:\"<\\\\-(?!\\\\s*\\\\d)\",e:hljs.IMMEDIATE_RE,r:2},{cN:\"operator\",b:\"\\\\->|<\\\\-\",e:hljs.IMMEDIATE_RE,r:1},{cN:\"operator\",b:\"%%|~\",e:hljs.IMMEDIATE_RE},{cN:\"operator\",b:\">=|<=|==|!=|\\\\|\\\\||&#038;&#038;|=|\\\\+|\\\\-|\\\\*|\/|\\\\^|>|<|!|&#038;|\\\\||\\\\$|:\",e:hljs.IMMEDIATE_RE,r:0},{cN:\"operator\",b:\"%\",e:\"%\",i:\"\\\\n\",r:1},{cN:\"identifier\",b:\"`\",e:\"`\",r:0},{cN:\"string\",b:'\"',e:'\"',c:[hljs.BE],r:0},{cN:\"string\",b:\"'\",e:\"'\",c:[hljs.BE],r:0},{cN:\"paren\",b:\"[[({\\\\])}]\",e:hljs.IMMEDIATE_RE,r:0}]}};\nhljs.initHighlightingOnLoad();\n<\/script><\/p>\n<p><\/head><\/p>\n<p><body><\/p>\n<h1>Finding Centroid point from shape file and map data labelling<\/h1>\n<p>In this test I try to find centroid point in shp file. In additionally I will show how to use this centroid point information when to labelling area name into map.<br \/>\nR Library <em>PBSmapping<\/em> offers wide range of functions to do that job. <\/p>\n<p><a href=\"http:\/\/science.ekqvist.fi\/blogi\/wp-content\/uploads\/2012\/07\/unnamed-chunk-7.png\"><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter size-full wp-image-265\" title=\"unnamed-chunk-7\" src=\"http:\/\/science.ekqvist.fi\/blogi\/wp-content\/uploads\/2012\/07\/unnamed-chunk-7.png\" alt=\"\" width=\"504\" height=\"504\" srcset=\"https:\/\/science.ekqvist.fi\/blogi\/wp-content\/uploads\/2012\/07\/unnamed-chunk-7.png 504w, https:\/\/science.ekqvist.fi\/blogi\/wp-content\/uploads\/2012\/07\/unnamed-chunk-7-150x150.png 150w, https:\/\/science.ekqvist.fi\/blogi\/wp-content\/uploads\/2012\/07\/unnamed-chunk-7-300x300.png 300w\" sizes=\"(max-width: 504px) 100vw, 504px\" \/><\/a><\/p>\n<h2>Needed Library<\/h2>\n<pre><code class=\"r\">require(rgdal)\n<\/code><\/pre>\n<pre><code>## Loading required package: rgdal\n<\/code><\/pre>\n<pre><code>## Loading required package: sp\n<\/code><\/pre>\n<pre><code>## Warning: package &#39;sp&#39; was built under R version 2.15.1\n<\/code><\/pre>\n<pre><code>## Geospatial Data Abstraction Library extensions to R successfully loaded\n## Loaded GDAL runtime: GDAL 1.9.0, released 2011\/12\/29 Path to GDAL shared\n## files: C:\/Program Files\/R\/R-2.14.2\/library\/rgdal\/gdal Loaded PROJ.4\n## runtime: Rel. 4.7.1, 23 September 2009, [PJ_VERSION: 470] Path to PROJ.4\n## shared files: C:\/Program Files\/R\/R-2.14.2\/library\/rgdal\/proj\n<\/code><\/pre>\n<pre><code class=\"r\">require(PBSmapping)\n<\/code><\/pre>\n<pre><code>## Loading required package: PBSmapping\n<\/code><\/pre>\n<pre><code>## ----------------------------------------------------------- PBS Mapping\n## 2.62.34 -- Copyright (C) 2003-2012 Fisheries and Oceans Canada\n## \n## PBS Mapping comes with ABSOLUTELY NO WARRANTY; for details see the file\n## COPYING. This is free software, and you are welcome to redistribute it\n## under certain conditions, as outlined in the above file.\n## \n## A complete user guide &#39;PBSmapping-UG.pdf&#39; is located at C:\/Program\n## Files\/R\/R-2.14.2\/library\/PBSmapping\/doc\/PBSmapping-UG.pdf\n## \n## To see demos, type &#39;.PBSfigs()&#39;.\n## \n## Packaged on 2012-03-01 Pacific Biological Station, Nanaimo\n## -----------------------------------------------------------\n<\/code><\/pre>\n<pre><code class=\"r\">require(maptools)\n<\/code><\/pre>\n<pre><code>## Loading required package: maptools\n<\/code><\/pre>\n<pre><code>## Loading required package: foreign\n<\/code><\/pre>\n<pre><code>## Loading required package: lattice\n<\/code><\/pre>\n<pre><code>## Checking rgeos availability: TRUE\n<\/code><\/pre>\n<h2>Importing Shapefile<\/h2>\n<pre><code class=\"r\">FIN_adm2.pb &lt;- importShapefile(&quot;FIN_adm2&quot;)\n<\/code><\/pre>\n<p>ImportShapefile reads the .prj file if it exists, but it<br \/>\ndoes not adopt the proj4 format <\/p>\n<pre><code class=\"r\">proj.abbr &lt;- attr(FIN_adm2.pb, &quot;projection&quot;)\n# abbreviated projection info\n<\/code><\/pre>\n<p>print(proj.abbr)<\/p>\n<h2>Generate map using PBSmapping plotting functions<\/h2>\n<pre><code class=\"r\">plotPolys(FIN_adm2.pb, projection = proj.abbr, \n    border = &quot;gray&quot;, xlab = &quot;Longitude&quot;, ylab = &quot;Latitude&quot;)\n<\/code><\/pre>\n<p><img decoding=\"async\" src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAMAAACR9g9NAAAAkFBMVEX9\/v0AAAAAADkAAGUAOY8AZrU5AAA5AGU5OTk5OWU5OY85ZrU5j9plAABlADllAGVlOQBlOWVlOY9lZjllZmVlZrVltf2POQCPOTmPOWWPZgCP29qP2\/21ZgC1Zjm1ZmW1j4+124+1\/rW1\/tq1\/v2+vr7ajznatWXa\/rXa\/tra\/v39tWX924\/9\/rX9\/tr9\/v2ZXfTVAAAAMHRSTlP\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/wBipdB4AAAACXBIWXMAAAsSAAALEgHS3X78AAAZZUlEQVR4nO2daXvbOJZGm4ldVZM4lZn0jJN025WeKVNdsSTz\/\/+7EXcsF8AFeEFSwns+yDKxiA+OAGIj9bcGFMnftj4BsA0QXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXygQXyhXJP5YtTxufRo3wvWIf\/v2dHk9\/fqy9YncBtcj\/vzl9fL69l0RX10Vm5UcyRWJ\/92q8auWZd1y+ZOaHuJTsa\/xa5ZlPRiH+B2wYvtZT+JTzUN8In1Tr7G+9ibdPMQncv6v\/7w3DmUty7oeruq69+TGHuITufTqj5WuPmdZDq5n+0pIknqIT6Qbzh2r9+v06me5tmbrq8AC4hPpx\/Eaq4ing+NzhHg58pVlqEpD\/KZkK8tgU57Q2EO8HLnKkmM1Wj3Ey5FNPC9WnHuIlyNPWUb4jFEP8XLkKMvIJpxf7SFejizic6WAeDkylGXC1AwzCcTLkaepj0\/CigXxcuxFPCsJxMuxj2s8Mw3Ey7GPazwzEcTLkUN8UiKIX5cMO6\/SdllA\/LrspcZD\/MpkqPBo6q+B3YjnRIJ4OXYiHuP4tdlH3w5TtqsjXpYpFR6rc+szlyWxCzqFjN4hXpCpLKn7HlLItjTXQLwkY1lOpb\/UfHRy7MDZhKEs1RvbFuUX\/72JSADxcozi5yPLqnx06pgEEC+HLX6p+dj4EL8JdlO\/sLGPFh8RF+LlqIaOvFiVh\/jroLpYtp5Ls8Q8mvrroOqqu33\/enKGGM5dB9XwICqj9BeIj0yKGr8N7Q4coujTq3xkStw7txFtWVKFn2w+9u6pmOgQL0dXlqT5tPyiTOI26e3oxRMBiVU+p3eIF2So8SJVPnJdN34ZGOLl6MuSMhDfPa8j01DDCS8QL4e1LDvjckIfT+0NQvw2KBsx7IfRkSnoir2Gd4gXRClLrviG+pasMvqDeDlU8bzpu5oKW2fYD\/FyaGVp2ST63fMDaomj0WACZyt08cRSjb5Sr\/jWvhSrVHiIF0QvS3pYN1Z86iHUw4FVKjzEC2KWpXMM55hXHx9Hn\/bpmLLdDKssHebdk2x1+m58zNVvh12W8Vtoku\/CwFz9dhDiYz2mD+UgfjvIsoyfck8iOh3Ey0GXpcDNkwwgfkMcDz+CeAbXLZ4+XJOzdsJA\/Ib4yrJTn1E+OncbEijLvLUe4rcjWJY5W3yI3w5OWUI8ze2Lz1XnMWW7HbyyzGY+LjrEy8Esy0zmIX4z2OKzmIf4zWCXZR7x2IixFduKx567zdhYfNw1BOJTOVb3h+rdk3KELz7bPA4\/Z4hP5O37y+G+OX95nQ9F1PiMM3h4enVeLsqPD63++VBMWWZdr+FkDvGJLKnxTeZFeo55iE9lwTW+2cFSHcTLUUX9\/tjG5iE+lePdz6\/pNb7J29pDfDbevj0dHtKv8U3e1j6cNcQnclF+eEzv1bfkbOwhPhs\/Xk4fm9OHJTVe9oTiMof4VM6fLp2590qFjyzLzPtuQ9lDvByR4jOdBTN\/iJdjZ+L9HwDxcuxLfMA8xMuxM\/F+8xAvx97Ee81DvBy7E+8zD\/Fy7E+851MgXg6IT6ck8dveNw\/xcuxodS78IRAvxw7Fu5sViJdjh009evVrEFWWWz8TCeLliNt5le001A+B+DXYzY47zodAvByRTX32R2Fhrn4l4styw931EC\/H3sR7QyFejn2Jx0aM1diVeOy5W499iQ+EQ7wc8fvqN7yjAuLl2NG+etxJsyY7WqTBvXNrsp8az8ga4uXYTY3HgxHWZTc1Ho9CWZfoXn2Ws2B+oyBejp3cJs3LFuLliC7LLZ9mDPFyxIvPYZ6ZJ8TLEbUDpyODeO6XCeLlSJirlzePR5quT0JZypuH+PVJKUuI7ylPvLR5iF+fpLKE+I4CxQubh\/j1SStLafEYzq3OLmo8N0OIlyOtc7dRvx7i5djFcA6LNOuzE\/G8PCFejt2Ix3r8uuxFPLZercxuxOM3adZlR+Kxy3ZN9iMe++pXZRfLstx8IV6OXWzEGHMORYB4OaK3XslO22m5QfyKbLy9utbU427Z9dj4hoouu3oiEBvi5di2xk+51TUnc4iXY9Pd1ZZoiF+NbcVHZg7xckB8OkWJFx7MxeUO8XJs3bkzsoP4tdj41jlLPB5puhJb3y0b1dhDvBxbz9VHjeggXo6k1TlB9XbHHs+rz8Dx7ufX6t2TcmTbW6io7xB+oUKet29Ph4fm\/OV1PrSx+LjcIT6Ri\/LDY\/P2\/WU+tHWNj8od4lP58XL62Jw+7KfGU+YhXp7zp6qq3isVfvN75yjxztwhXo4q\/ulHjYx59xI8xIsj1KuXEu\/IB+KlkerVp4k3Krg7D4iXRqpXn2S+q+DLzEN8KkK9+ialY19b6ZzmIV4aoV59kyB+sMyq8hCfn7XEz9sq5yNo6jckWbx3qcYOVH0H91K7wiBejiXiXfK6gZp5yIzjz50+DPFypIt3Ndb90agtVXYe9GGIl2NZWdZ6J11pBfROXOzoD+Kzk16WXb1Wb3zTbdU152ruyZsA4uVYLH4UK7Zw4\/mqQLwcS5v6UZOUd2\/7APFyLCzLuhkmYsXquy8Q4uWQKMu067grL08gxMshVJbrmId4OSA+HYjvWGS+v1o4t2YMQLwcSTuvSBbV+WEiMNBwQLwcgmW5QPzsG+JXQqQs7T0WEWn1eg7xIhyr++dHbwwZ8clTOHE3yEM8k\/Pn\/\/14\/vziiyJQlrX2Jymt54AKxDNZR\/xASo2300C8CMeqqlZo6luSFuGibpKGeEm2HcdHPQ8D4iXZqsbTK\/WBLCCeRbeX2thNbbGJeOdMDcTL0N4zo90wZbNF586zr9o\/dwfxTNoe\/Q579Uk7bFsgnsuqvXqZmKjx67DBsmygwkP8KkiVJW8XDmfTLcQvx75H0mJpWSq7bANGOdsyA5lAfATn\/\/EGL95sqdT15Gu3mp0PiI8gb69+3Fzf\/+PzxuoDQLwMXVN\/742ycAfO3HaHrt6cGo\/fpFmPpTVeu8YvbeshXoh21k574o3N8rKcL\/Ch3h0jL38wxPN47ubq7\/JP2dZ1WDvEr0hgnr5FYgfOoH3BAJ0ZBeJZnL\/8lX8cP+kOz98s791BvBzLO3fEWzKewEAe4plk79wxp+j7u2RY4rH1SoAVOndM7xEZ+toGiGeSv3O3vBbH5Arxcgg+\/MgZSSxXiOfynLVXz73Cx2fsOA7xTM6f\/+\/LX\/lW57g9O7GMIZ7J+fOf3\/\/M16tnDczTbqmjD0M8l8P9IbQ6tyD35asukckgXo7M4tPyxTVegHwTOOE52tQK7wqA+AiyiQ\/vmkx9RAbES5BLfD7vEC9CSPyS55b7IyRmjGv8UrLeNBmuzwueggXx+cn482PybT3Ey5HzBwfFe3cQL0eieF4siN8vOX9NekFbTyeFeDmy\/n78gseZkwkhXo6EsoyyuaDS28cgXo4k8bkih1JCvBwp4jPG1lJaSSFejtiylPoJOdZHmUcgXo7IshT77UDWZ5lHIH4J+uOsd1zj7aQQnwgxe59Q4yMfYRn3Ad6kEJ\/K4X5hjY9\/Mv0S8WZaiE\/meL9AvPoQBH6iiA8IJYX4dI53\/0gXPxuPEC+4SAfxCzj9skA89TaUaMG0rXkA4uWoIrbgqOKZMmve74qFPm4A4lM53v38Wr17Uo7EXeOV97Epot1jylaOt29Phwf9Hlp2WRq1NuFu90j1WKSR46L88KhvvwyX5TCCM7UlPb0sQj2WZSX58XL62Jw+RNX42vE8i8TbKbTW3xdIpIX4VOynGvPEUxKDzywOBqgZ18bHYAdOZkJlaQrRg3wJPVmqf\/UPgPi1YIl3Vd808cqzMF2psOcuO4GyrGu1Rpo6vA848+bqSuyfEYZ4OXjiXU1+ovjGLd77XDSIl4N3jW96T1abnyze6be7sLiSQbwcbPHT\/3qoOyEjW8dYATdNrgCrqZ\/\/NStquO\/uDXbEwTU+P\/6yNHphdiV16130CDxHYoiXw1eWhGb3JKwVtES8IwDi5VgqXp2HIQ4HiOwiQLwcAfHWISra8KLFXnojNfkdg3g5PGVJ9bpdVX5cwpsPsj49cgYI4uXwbMBhL83MQ++598cU7wvDRoyc8L17Zu2Jkf7yGo87abJil2U9Ndv2LJ1nuk2LyN5w4RdvhkK8HGNZqtsjOsXW\/Kx7\/cTu+4uItz8N4uVQxA8X57omF1HGazfzss\/8+NAXRFcP8XLMZaksxoy7os2w7q2ZA6MN8BCMqGYF8XJMZald1Pt\/piZA2xPFmaETFK\/mBfFyqOK1Gq\/\/r6RgTc0Kim8gPgfTNb4eh2sKDVV3rblZKlupa7z+gRAvhy6+MbQ7h+3quJ3Klt2rj5rZhXg5KnV\/s1XfnWKCk7Miwznz4yBeDuUaP7lX6rSzRiq9PTqc9ensTiDES6OLt0frbvPeCEyj\/ImeLh7Ey6GP44ldFb623jNFJ3qNh3h52rK0Vlgaz4E5YLwo0MHMz4+a24V4OVTx5oi9O+YWE7m05ooYEQ3i5VDLklhdSfUuXuW7aBAvhy3eXpKjCS6vMM8goh8I8XJYwzkN79ao0JIq9xS4XXuIl6RSpu5sBQuu4nzx\/EE\/xMvRde6oldimCckLtfXsc2D37CFeDqOp1wOTu28Re3BiegMQL4cq3lIQasyd4RHNPMRvg7UDRyVkxC0+zjw3HsTLYW3EUKwFu2\/Jg\/y0yBAviLoeP70bj\/mTiszc+aZ9ragQL4d2je\/feDbU6vjN8z4\/Yi82xvGSzOKVN9TQjsJf5TlCa+cTtaioEC\/HUJZ1ozfwvGroH+Z7HmajxGFTY65eklF8ze3SqYSihrOC+K2YxU+HBOdaBat8DfGizL36xnoXJljlgxnwP6qBeEks8aIjcEaNj9mwAfFyqItz\/d+Y1Esv8vxOPTZiCKNM2SqvXJbVeNe9Os6cIF4Oq8ZHpV4kPqqRgXhhFooPBNeuf\/rEtSPElRPEy1FN47nhLz9peJJHGyK6t3Wx5wghXg5zrp4zXzfN7YUn8+txNt7exqnOHDDnAyBeDkM8q9UddmqxviOjcWsVLnzHbWOFQ7wc2hMxQssuwytvBUfNVs9Cfzfl6cnEPNl9cBvih3bZHXO8jZq3\/DKlcv1rtvveT7ZOdhfcjnjv74KovyWdvpOSzsL7VZpjQrwc5p479+2vmih2\/g7xMfv75hCIl8No6p1ts9ElZ+dv+XX1Id0tDcTnQClLsnunPyBjPsrM3o7o2mvFaWkgXg79psnGqtrUHXVNxGiOPEgmpfPT4kK8HJZ44mJOKWEP49lnQkbVD0K8HFpZWlWevOyPYc6QKXHEsI\/OzPgyQLwc3qbeO43rH\/xF9AOmBMFDEC+HT3xg2dRfo6O9E3mZeUB8Kseqemzevis\/IG\/36s3JNa9B33xPLOEhHsQn8vbtqXl+dIvvXqf\/WPskcoq3DkB8Iucvr03zxz954ufjviw58y5M3Av2IxCfSFvjLy\/\/QYvXxm7KnHqo6+6IsLjK2xlAfCqnX1+GlxH9Gj9Nq9WGfzfOHl5KlaeXcCcgXo6qMlbk9UrM3whJHF5U56nUEJ\/K8e7n1+rdk3LEnrnr3hHH3AiK9\/csID6RyzX+8NB38UaMBxzaaVgzs\/yjgZz8qSE+kYvyg3sc33BmTUnkxdOJIT6VHy+nj83pg6PGu5ZjgtnKibcHkyoQn8r506Uz916p8MZwLk282HhuniyE+NwY4knzoTxcnbuE0+mHkq7vGsTLYTT1SVXe3J7DTUdl5V36gXg5OOKDmRjip9mf6LOZb74ggXg5jLIk98gFM9FiqE\/MizyZcdrQFQ7xcjDEB\/PQH6CjVNl48f5UEC+HvvWKEh\/OQ9u7oU71x1b5OpAI4uUw9twl1Xh1DX9OVPOfXaim9qWBeDmWi1eXdLSDMVts5w\/ypYF4OSzxZsGHxVNv++oeI55zdYB4OczhnC0+lIOSJLCc7s+G8XEQL4ct3hyUB7MQEM\/c4AfxctjDOXoazge5tMIVr\/UGAokgXo4q0Lvj3SlFxOUO\/2O+LBAvh7EDJ2U4Nz4qI7hlzkxnXlmCSSBeDnM93i57jnji3nre0o4+2Rv6JIiXw+rcmRFipupjxKd8EsTLEVid490Nbb2x\/7GHC1Yu4U+CeDnMXr0RHLe9upuptbOaJvCnTgQ5RRh9sltzQ+KTOvXmTK2VdN5bMfUDiUwgflUCy7KRd9Jo\/bVaUx8x8+sE4uUIde68PqabrhwDeHU\/TcA8a7oH4uUIbcTw+eq11sbzkfR6rszK+U6DN88H8XIwhnOu3dP2\/ovxn\/HNeFEf\/3PCnN+FeDmC6\/GNS4vZG7DeDg1C0wTbeu68PsTLYdw756jc5EHXMp51UR9HcK5zYC\/kQbwcwZsmuwDikKf\/XxsTNv7tOBEbdSBejvC9cw1vQk+r8cQqn+sEYtbtIV4ORo0nNmeQPYH5DRHuFh86QwWIl4NxjbfWWwMdgaide1EbdSBeDrPGE1skeeL1o+x5oBjvEC8IJd6+PDP2R\/nX3kQqPMQLomy9midgqVWXeZsNQzy7qY\/yDvGCqM+8Iuq0PvE+rqlSGYXEO3oPaSe7D25FvOJhWkTVtsTNMW2MeVrWKh9jzc5xsvvgNsQba+m2PL8kJSW3Wx9pvYF4SbSmvpn\/sWdkGatr41IdrZnTCDBPdh\/cgni7mnJn3pTgQMc\/LkcKiJdj+jFp33R6w6jxtaqSjO0f6XOAeDn0Xr0nol9Urey1cUbXexH8c5yAeDnI4RxFUHwdFD9E7P\/wT3EG4uUwHlvujqgN9xyhyl4b54x+3YemnCvECzJ17oLDakdvfQzrK\/LUWfCtv6dd4BuIl4S1ESMYZZrWHaUTaz1qHhC\/OdavUIUwl+qUzoFuW2CrlQnEy0HN3IX7ccqAXZ\/20SNyhnUxQLwcRK8+KKae5+O9i7GOgUKyd4gXxBDv+PFonWl+1hJvz8rah9IrPMQLYooPDecbsydHhFj\/69cGgZPdB7chXlfvS8LfP6ct9g1vEs+zA+LlqOgqHz1rT43OtaaB15oETnZB2gxct3jlfd0o26vY1Xo8TF7P9doO8ftB3YGjefdM0pGHzZ5crVX1ZvwGoKnfCco4ft49M9nybrLSj2pfidqiWdihN052F9yGeM2S+s5M4RFv59WMyp0Jk052H9yUeKuyWikc13jlr\/EdWtqXJ052H9ySeLuRHgKbOR6VjdKBbxTn6Qtx\/pPdB7cknvSuVfyQeH0yX1I7xEtC3Ehj9Mb1jry7sz\/moU3pKVcAge8AxMtB3kFldO8U8Z7hvZ4888nugpsTb1\/v5\/hM8blPdhfcnHjToTZOs7NQI\/KW95af7C64DfEu681Uy12\/LzRGzWjcPNl9cBPi7Rlbrca7vZItQ+aT3QlXJP5YtTwqRybx3Qtd5YcWnMyRuWlHBohP5O3b0+X19OvLfMj+3bl6atBV8Y4pvLXqeg\/EJ3L+8np5fftOi6\/nibbxVZuANTLztQOZgPhEzr\/7a3yPNmi3enpDFNfXISsQn4r7Gq8w9+YV6UO\/3QzIfsYaEC9HNe69Iodzej+vru0u4Lonu+qnBbki8ce7n1+rd0\/KEW0Cx29fmY+3125WAeITufTqDw9DF2\/AscvWUm5+K7bwDvGpXJQfHj29esf0nbfhXxOIT+XHy+ljc\/pA1\/gOlvhVz1kB4lM5f7p05t4rFd78bVljcyxlf+1zVoB4OcynVwfZ7EwbiJeEnLLV6rdyUHorVSwQL4dd49u\/1WRfnaifq3vAwLJgdzjEy1Gpf1rjw7\/9SL3qB+\/1EDyVfEWlZv9Njg7xclRXxdalpXPV4kE6EF8oEF8oEF8oEF8oEF8oEF8oEF8oEF8oEF8oEF8oNyX+7WtV3b32r+7wy5vDgyf4UFX3nuBnfcOnwnO3+dsdvituSvzxIuz5sX91h7c79EnxgeRD8N3rkf5aXcLPn1\/c4fvipsS3HB7nV0f46cO\/SfF98OG\/HTW+Dz59eD06gy\/i253An19cEfbDrYnvpbjVDNXSJf4S\/Hzv\/t60+R4d15GmbeUfGnMn8F65MfHP9\/OrK\/zQLo7T5kfpji9GF3yv37+nc3hAjd+At6+P06s7vHGJ7YOPrhrfB7vFtpkeHnCN34Dnri4\/O2v0FEKLH4IP3tRt397xxXruOgfo1YM9A\/GFAvGFAvGFAvGFAvGFAvGFAvGFAvGFAvGFAvGFAvGFUrB47TlKCu3i27+UoD+uYc0lGoin0J6mB\/E3xiC+W0w9\/fZLv03z\/d+fLjX+ubo\/feyUP7dH2jjOzVpXSvHi+y2Ul\/fDkzN\/b8Vfanwvvt2p9fuTc2\/GFVO8+H6nVa+53Vz7hyb+0B3xbNa6WooXP9T4XrNa4z+89nX90hI492ZeMSWL\/6VqL+\/9Nb4TP13jz5\/aJ2X3V\/d3v3XXeNd+q2ulYPE217JDVgKIV2g3Ut5gq04C8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YUC8YXy\/5vZR66u7QZKAAAAAElFTkSuQmCC\" alt=\"plot of chunk unnamed-chunk-4\"\/> <\/p>\n<h2>Calculate the centroids of adm2 area<\/h2>\n<pre><code class=\"r\">FIN.centroids2 &lt;- calcCentroid(FIN_adm2.pb, rollup = 1)\n<\/code><\/pre>\n<p>Check results<\/p>\n<pre><code class=\"r\">str(FIN.centroids2)\n<\/code><\/pre>\n<pre><code>## Classes &#39;PolyData&#39; and &#39;data.frame&#39;: 21 obs. of  3 variables:\n##  $ PID: int  1 2 3 4 5 6 7 8 9 10 ...\n##  $ X  : num  29.9 27.5 26 27.9 26.3 ...\n##  $ Y  : num  62.9 63.2 61.4 61.9 67.7 ...\n##  - attr(*, &quot;projection&quot;)= chr &quot;LL&quot;\n<\/code><\/pre>\n<h2>Plot the centroids against the source polygons<\/h2>\n<p>In this purpose we use the PBSMapping package graphics routines. If you want you can plot these centroid point in to the map<\/p>\n<pre><code class=\"r\">plotPolys(FIN_adm2.pb, projection = proj.abbr, \n    border = &quot;gray&quot;, xlab = &quot;Longitude&quot;, ylab = &quot;Latitude&quot;)\naddPoints(FIN.centroids2, col = &quot;blue&quot;, pch = 19)\n<\/code><\/pre>\n<p><img decoding=\"async\" src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAMAAACR9g9NAAAAk1BMVEX9\/v0AAAAAADkAAGUAAP8AOY8AZrU5AAA5AGU5OTk5OWU5OY85ZrU5j9plAABlADllAGVlOQBlOWVlOY9lZjllZmVlZrVltf2POQCPOTmPOWWPZgCP29qP2\/21ZgC1Zjm1ZmW1j4+124+1\/rW1\/tq1\/v2+vr7ajznatWXa\/rXa\/tra\/v39tWX924\/9\/rX9\/tr9\/v3LuyHZAAAAMXRSTlP\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/8AH5pOIQAAAAlwSFlzAAALEgAACxIB0t1+\/AAAGj5JREFUeJztnWt72zp2RoeJ43OmiXMyk7ZO0ton046pTizJ\/P+\/ruIdlw1gA9wgKeFdH2yJJEA+WMIdJP\/SgCL5y9YXALYB4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gsF4gvlisQfq5bHrS\/jRrge8W\/fny5\/T7+\/bH0ht8H1iD9\/fb38ffuhiK+uis1SjuSKxP9h5fhV07JuufxLDQ\/xqdh1\/JppWQ\/GIX4HrFh+1pP4VPMQn0hf1Gusr71JNw\/xiZz\/\/rd7Y1PWtKzroVbXvScX9hCfyKVVf6x09TnTcnA921f2JKmH+ES67tyxer9Oq36Wa2u2fgosID6Rvh+vsYp4end8jBAvR760DGVpiN+UbGkZLMoTCnuIlyNXWnKsRquHeDmyiecdFece4uXIk5YRPmPUQ7wcOdIysgjnZ3uIlyOL+FwhIF6ODGmZMDTDDALxcuQp6uODsI6CeDn2Ip4VBOLl2EcdzwwD8XLso45nBoJ4OXKITwoE8euSYeVV2ioLiF+XveR4iF+ZDBkeRf01sBvxnIMgXo6diEc\/fm320bbDkO3qiKdlSobH7Nz6zGlJrIJOIaN3iBdkSkvqvocUsk3NNRAvyZiWU+ovNR8dHCtwNmFIS\/XGtkXxxf9uIgJAvByj+HnLsiwfHTomAMTLYYtfaj72eIjfBLuoX1jYR4uPOBbi5aiGhrxYlof466C6WLaeS7PEPIr666Dqsrt9\/3pyhOjOXQfV8CAqI\/UXiI8Mihy\/De0KHCLp07N8ZEjcO7cRbVpSiZ9sPvbuqZjDIV6OLi1J82nxRZnEbdLb0YsndiRm+ZzeIV6QIceLZPnIed34aWCIl6NPS8pAfPO8jgxDdSe8QLwc1rTsjMsJvT21NQjx26AsxLAfRkeGoDP2Gt4hXhAlLbniG+pXskrvD+LlUMXzhu9qat863X6Il0NLS8sm0e6eH1BLbI0GAzhboYsnpmr0mXrFt\/ajWCXDQ7wgelrS3box41MPoR42rJLhIV4QMy2dfTjHuPr4OPq0s2PIdjOstHSYdw+y1emr8TFWvx12WsYvoUm+CwNj9dtBiI\/1mN6Vg\/jtINMyfsg9iehwEC8HnZYCN08ygPgNcTz8COIZXLd4enNNjtoJA\/Eb4kvLTn1G+WjcbUggLfPmeojfjmBa5izxIX47OGkJ8TS3Lz5XnseQ7Xbw0jKb+bjDIV4OZlpmMg\/xm8EWn8U8xG8GOy3ziMdCjK3YVjzW3G3GxuLj6hCIT+VY3R+qd0\/KFr74bOM4\/JghPpG3Hy+H++b89XXeFJHjM47g4enVebkoPz60+udNMWmZdb6GEznEJ7IkxzeZJ+k55iE+lQV1fLODqTqIl6OKev\/YxuYhPpXjh1\/f0nN8k7e0h\/hsvH1\/Ojyk1\/FN3tI+HDXEJ3JRfnhMb9W35CzsIT4bP19On5rTxyU5XvaC4iKH+FTOny+NufdKho9My8zrbkPRQ7wckeIzXQUzfoiXY2fi\/SeAeDmyib+7u4u8lO4E3jNAvBy5xN\/dZTAP8XLsTbzXPMTLsTvxPvMQL8fO6nj\/WSBejn216gNngXg5djWAM57FtQPi5djR7Fz4JBAvxw7Fu4sViJdjh0U9WvVrEJWWWz8TCeLliFt5le0y1JNA\/BrsZsUd5yQQL0dkUZ\/9UVgYq1+J+LTccHU9xMuxN\/HevRAvx77EYyHGauxKPNbcrce+xAf2Q7wc8evqN7yjAuLl2NG6etxJsyY7mqTBvXNrsp8cz4ga4uXYTY7HgxHWZTc5Ho9CWZfoVn2Wq2D+oiBejp3cJs2LFuLliE7LLZ9mDPFyxIvPYZ4ZJ8TLEbUCpyODeO6PCeLlSBirlzePR5quT0JaypuH+PVJSUuI7ylPvLR5iF+fpLSE+I4CxQubh\/j1SUtLafHozq3OLnI8N0KIlyOtcbdRux7i5dhFdw6TNOuzE\/G8OCFejt2Ix3z8uuxFPJZercxuxOOdNOuyI\/FYZbsm+xGPdfWrsotpWW68EC\/HLhZijDGHDoB4OaKXXskO22mxQfyKbLy8utbU427Z9dj4hoouunoicDTEy7Ftjp9iq2tO5BAvx6arqy3REL8a24qPjBzi5YD4dIoSL9yZi4sd4uXYunFnRAfxa7HxrXOWeDzSdCW2vls2qrCHeDm2HquP6tFBvBxJs3OC6u2GPZ5Xn4Hjh1\/fqndPypZtb6GifkN4Q4U8b9+fDg\/N+evrvGlj8XGxQ3wiF+WHx+btx8u8aescHxU7xKfy8+X0qTl93E+Op8xDvDznz1VVvVcy\/Ob3zlHinbFDvBxV\/NOPGhnz7il4iBdHqFUvJd4RD8RLI9WqTxNvZHB3HBAvjVSrPsl8l8GXmYf4VIRa9U1Kw762wjnNQ7w0Qq36JkH8YJmV5SE+P2uJn5dVzltQ1G9IsnjvVI29U\/UdXEvt2gfxciwR75LXddTMTeYx\/tjpzRAvR7p4V2Hdb41aUmXHQW+GeDmWpWWtN9KVUkBvxMX2\/iA+O+lp2eVr9cY33VZdc2pzT9wEEC\/HYvGjWLGJG89PBeLlWFrUj5qkvHvLB4iXY2Fa1s0wECuW3307IV4OibRMq8ddcXl2QrwcQmm5jnmIlwPi04H4jkXm+9rCuTRjAOLlSFp5RbIozw8DgYGCA+LlEEzLBeJn3xC\/EiJpaa+xiAir53OIF+FY3T8\/eo+QEZ88hBN3gzzEMzl\/+Z9P5y8vvkME0rLW\/iWF9WxQgXgm64gfSMnxdhiIF+FYVdUKRX1L0iRc1E3SEC\/Jtv34qOdhQLwkW+V4eqY+EAXEs+jWUhurqS02Ee8cqYF4Gdp7ZrQbpmy2aNx51lX7x+4gnknbot9hqz5phW0LxHNZtVUvcyRy\/DpsMC0byPAQvwpSaclbhcNZdAvxy7HvkbRYmpbKKtuAUc6yzEAkEB\/B+T+9uxcvtlTyenLdrUbnA+IjyNuqHxfX918sb3d3d+qRrOg8QDyTrqi\/9x6ycAXOXHZTtffd3Wyek+PxTpr1WJrjtTreJ55jHuKFaEfttCfe2CxPy7mC94vnFPYQL8JzN1b\/If+QbV1T2hutjof4FQmM07dIrMAZtC\/ooDMPgXgW56\/\/yt+Pn3SHx2+Wt+4gXo7ljTviI3mcQEce4plkb9wxh+j7u2RY4rH0SoAVGndM7xER+soGiGeSv3G3PBfHxArxcgg+\/Mh5kFisEM\/lOWurnlvDx0fs2A7xTM5f\/vfrv\/LNznFbdmIRQzyT85d\/\/vhnvlY9q2OedksdvRniuRzuD6HZuQWxL591iQwG8XJkFp8WL+p4AfIN4ITHaENHaHM4SjjX8RAfQTbx4VWToR+GPmurBHQFgPgIcolf7h3isxISv+S55f4DglE4xaOOX0TWmybD+ZnRonfU8RjAWYGMrx9Lf\/4ZxOcn5wsHk81DfH4SxfOOgvj9kvNt0gvKejooxMuR9f3xCx5nTgaEeDkS0jLK5oJMb2+DeDmSxOc6OBQS4uVIEZ\/xaC2kFRTi5YhNS6lXyLFOZW6BeDki01Ls3YGsc5lbIH4J+uOsd5zj7aAQnwgxep+Q4yMfYRl3Am9QiE\/lcL8wx8c\/mX6JeDMsxCdzvF8gXn0IAj9QxAlCQSE+neOH\/0oXPxuPEC84SQfxCzj9tkA89TEUaMGwrbkB4uWoIpbgqOKZMmvee8VCpxuA+FSOH359q949KVvi6njlc2yIaPcYspXj7fvT4UG\/h5adlkauTbjbPVI9JmnkuCg\/POrLL8NpOfTgTG1JTy+LUI9pWUl+vpw+NaePUTm+djzPIvF2Cq309+0kwkJ8KvZTjXniKYnBZxYHd6gR18ZpsAInM6G0NIXou3wBPVGq\/\/UTQPxasMS7sm+aeOVZmK5QWHOXnUBa1rWaI00d3geceWN1BfaPCEO8HDzxriI\/UXzjFu99LhrEy8Gr45vek1XmJ4t3+u0qFlcwiJeDLX76ru91B2RE6+gr4KbJFWAV9fPXIaNOdzuG2+7e3Y5jUMfnx5+WRits\/BXM9ze79S56BJ4jMMTL4UtLqywevys3tisjMGZoxtndPxtyK8TLsVS8Og5jHMs4e2QTAeLlCIi3NvX\/7DdPmG2ypTdSk41+iJfDk5ZUq5vSMTXE9WYgh8gRIIiXw7MAhz01M3e95z4AU7xvHxZi5ITv3TNqT\/T0l+d43EmTFTst66nYtkfpPMNt2oHsBRd+8eZeiJdjTEt1eUSn2Bqfdc+f2G1\/EfH22SBeDkX8UDnXNTmJMtbdzGqfefrQD0RXD\/FyzGmpTMaMq6LNfd1HMwZGGeAheKAaFcTLMaWlVqn3X6YiQFsTxRmhExSvxgXxcqjitRyvf1dCsIZmBcU3EJ+DqY6vx+6aQkPlXWtslopWqo7XTwjxcujiG0O7s9uu9tupaNmt+qiRXYiXo1LXN1v53SkmODgr0p0zTwfxcih1\/OReydPOHKm09uj99ibiWdXsRiDES6OLt3vrbvPeA4it1NPp+QM9tXax++A2xGsDbuFptrnx79hvbSFfSwDxW9GmpTXD0szlsqesHysFere1hX4fRdTYLsTLoYpXeuyTJY+YyKk1x\/soIH4b1LRUCvtRfKp32Z78eBjEy2GL70RM4t0hg9MrzCtgm68hXhCrOzcQzPDBHM8fteWqh3hBKmXozlYQWYtH7GafRo0R4uXoGnfUTGzThOSFynr2Neincb2T7HIcxMthFPX6zuTmW8QaHMK7yzzEC6KKt3J4qDB37o8o5iF+G6wVOCohf27xcebVLxC\/DtZCDMVasPmW3Mn3HYw6fhXU+fjp07jNHzR25M55HLM9gFa9JFod33\/wLKjV8ZvnnT9iLTb68ZLM4pUPVtfOUfz6szxHaO18ohZ1KMTLMaRl3egFvG7N+aJ3X8y152E2yjFsaozVSzKKr31NOqf4xeP1EL8Vs\/hpE198en9vOiBqmAfi5Zhb9Y31acbVxQpm+dDpI8Z1G4iXxBK\/oAdO7A9HwJ69bSBeEnVyrv8fE3ppJc9v1GMhhjDKkK3yl4tmjVpLGQhdm3G4j23\/QLwcVo6PCh3q8wX6exGFDMQLs1C88pkQryklpoBqxx7iRBAvTDX154b\/\/KDG0BwlXv1kzfl6+xHmufqL5V\/cGly3+PEDP8fX09iefrDjRhnjbjxtn35uzymNi90HtySeVeoOK7VYv5HRuDULF77jtrH2Q7wc2hMxQtMuw1+i3PYEout5c7O\/GWhe7D64DfFDuew+cryNmjf9MoVyfTXLfe+ZrYvdBbcj3vteEPVd0qkrKZVGJD\/Lz0dCvBzmmjv37a+aKHb8DvEx6\/vmPRAvh1HUO8tmo0nOjt\/y62pDuksaiM+BkpZk805\/QMa8lRm9faBrrRWnpIF4OfSbJhsra9e29X47szdHbiSD0vFpx0K8HJZ4ojKnlLC78ewrIQ\/VN0K8HFpaWlmerPbHfc49U+CIbh8dmfFjgHg5vEW9dxjX3\/mLaAdMAYKbIF4On\/jAtKk\/R0d7J+Iy44D4VI5V9di8\/VBeIG+36s3BNa9B33hPLOEuHsQn8vb9qXl+dIvv\/k7fWOskcoq3NkB8Iuevr03z53\/zxM\/bfVFyxl2YuCfsRyA+kTbHX\/78Gy1e67spY+qhprvjgMVZ3o4A4lM5\/f4y\/BnR6\/hpWK02\/LtxtvBSsjw9hTsB8XJUlTEjr2di\/kJIYvOiPE+FhvhUjh9+favePSlb7JG77hOxzU2UePdzD\/QgVGCIT+RSxx8e+ibeiPGAQzsMa2SWv9XzpJMm+IOD+EQuyg\/ufnzDGTUlkRdPB4b4VH6+nD41p4+OHO+ajglGKyfe7kyqQHwq58+Xxtx7JcMb3bk08VH9OW8dPw0WQnxuDPGk+VAcrsZdwuX0XUnXbw3i5TCK+qQsby7P4YajovJO\/UC8HBzxwUgM8dPoT\/TVzDdfkEC8HEZakmvktO\/k+0WMAOxhP\/tc3uUdEC8HQ7z2NfQuKe0+m3jx\/lAQL4e+9IoSr3+lxSsf1aH+2CxfBwJBvBzGmrukHK\/O4c+Bav6zC9XQvjAQL0eseLuOV6d0tI0xS2znE\/nCQLwclngz4cOdOepjn91jxHNqB4iXw+zOjQk\/ZWzGuB0xqZdYvwdOB\/Fy2OK7pOe8cHAONH0itzJgLvCDeDns7pwhnhEHObXCFa+1BgKBIF6Oim7dxeR4c51W8FWF88mauB8LxMthrMCZkn6u41mVvJJvh98MI5zZmAwGgXg5zPl4O+054tV767ni7dHZ4JkgXg6rcWceELMOI0Z8ypkgXo7A7Bzvbmj9A9UstEcIUn5hEC+H2ao3dsctr+5Gau2opgH8qRFBDhFGX+zW3JB4y0fkunp6Sf68tmJqBxKRQPyqBKZlI++k0dprtaY+YuTXCcTLEWrceX1MN11pnXEjOt4ULWu4B+LlCC3E8PnqtdbG85H0fK6MyumB9Vk+3jgfxMvB6M65Vk\/b6y\/GL+OHsVIfvyno8\/rM8V2IlyM4H9+4tJitAevjUCA0DVXWa+K54\/oQL4dx75wjc5Mba\/OrPrxvTNrVbvHsiTyIlyN402S3g9hkt\/\/VGT2txUcux1EmA9gTuBAvR\/jeuYY3oDeLN1t7vpgjsnsD8ZIwcrw92kq3BCbxlEy3eNZl9kC8HIw63uynOxsCUx0fsXIvaqEOxMth5nhiiSRPvL6VPQ4U4x3iBaHE29UzY32Uf+5NJMNDvCDK0qt5AJaadZmX2TDEs4v6KO8QL4j6zCsiT+sD7+OcavfduLMiJN7Reki72H1wK+IVD9MkqrYkbj6yse6lMsZpWbN8jDk7x8Xug9sQb8yl2\/KMr6b4OSS3WR9pvYF4SbSivpm\/2COy9sibNcsyTtVpgedOnh5BvHeIF2R8p7CVTYMjb3aG72IwV1rOPxDuWF7wYvfCTYh3N9enA3zRTLNwvWflaEW8t6cfcbF74RbE0+NtKn5R9bjWZh63HVCqBL0VseRi98ENifdleW809VizW+KtpTZqlzD1YvfBbYgP9a607p5j72C+bjwj+nW\/d+nF7oFbED+48pb1xGiusq\/PyFNjwR3XvCYn9WL3wi2I72CMp\/hnaGbpxFyPGgfEb471FqoQ5lSd0jjQbXvrhTQgXg5q5C7cjlNG9fVhH\/1A+reU7B3iBSFa9froCxGmnsfjvZOxjo5CsneIF8QQry+XczxcfhqftcTbo7L2pvQMD\/GCmOK1XEqL11tyxB7ru143CFzsPrgN8br6fotXvM00XzttUGMePghd7B64bvF0lp\/NWyFoeVO5rw3REhEvMQ\/xcpivJgkK8ozI1fMszbRNy+0Qvx\/UFTiad88gHbl5bMkpg\/V6XEPbUeRi98FtiB\/fCqLZchXqxNbpJzFMy9oF\/MLs3kC8JGQV36ifzBAe8XZczajcGTDpYvfBTYm3MqsVQunrEVunQl0xvzinWxe7D25JvF1IDzub+bj+P\/lkAy0KaohH6GL3wS2JJ71rGT8kXh\/Ml9QO8ZIQN9IYrXG9Ie8a1VN+ImoVodQAAr8BiJeDvIPKaN4p4p0TOEY5IZvT7YvdBTcn3q7v5+Pd4zpWwyDnxe6CmxNvOtT6aXYU6oHG9F62i90FtyHeZb2Zcrnr\/ULjoRmNmxe7D25CvD1iq+V4t1eyZMh8sTvhisQfq5ZHZcskvvtDZ\/mhBCdjtBbt5ATiE3n7\/nT5e\/r9Zd5kv3eungp0VbxjCG+tvN4D8Ymcv75e\/r79oMXX80Db+FcbgDUi85UDmYD4RM5\/+HN8j9Zpt1p6wyGun0NWID4Vdx2vMLfmFelDu93ckf2KNSBejmpce0V25\/R2Xl3bTcB1L3bVswW5IvHHD7++Ve+elC3aAI7fvjIeb8\/drALEJ3Jp1R8ehibegGOVraXc\/FVs4R3iU7koPzx6WvWO4Ttvwb8mEJ\/Kz5fTp+b0kc7xHSzxq16zAsSncv58acy9VzK8+W5ZY3EsZX\/ta1aAeDnMp1cH2exKG4iXhByy1fK3slF6KVUsEC+HnePb\/9VkXx2on7N7wMCy3e79EC9Hpf5rjQ9f+5561Xfe62H3lPIVFZr9P\/lwiJejuiq2Ti2dqxYP0oH4QoH4QoH4QoH4QoH4QoH4QoH4QoH4QoH4QoH4Qrkp8W\/fqurDa\/\/Xvf\/y4fDg2X2oqnvP7md9wafCc7f4271\/V9yU+ONF2PNj\/9e9v12hT4oPBB92f3g90j+ry\/7zlxf3\/n1xU+JbDo\/zX8f+08f\/I8X3uw\/\/4cjx\/e7Tx9ejc\/dFfLsS+MuL64D9cGvieyluNUO2dIm\/7H6+d\/9u2niPjnqkaUv5h8ZcCbxXbkz88\/3817X\/0E6O0+ZH6Y4fRrf7Xr9\/T+fwgBy\/AW\/fHqe\/7v2NS2y\/++jK8f1ut9g20sMD6vgNeO7y8rMzR097aPHD7oM3dNu2d\/ywnrvGAVr1YM9AfKFAfKFAfKFAfKFAfKFAfKFAfKFAfKFAfKFAfKFAfKEULF57jpJCO\/n2D2XXn9cw5xINxFNoT9OD+BtjEN9Npp7++lu\/TPP9vz9dcvxzdX\/61Cl\/bre0xzgXa10pxYvvl1BePg9PzvyjFX\/J8b34dqXWH0\/OtRlXTPHi+5VWveZ2ce2fmvhDt8WzWOtqKV78kON7zWqO\/\/ja5\/VLSeBcm3nFlCz+t6qt3vs6vhM\/1fHnz+2Tsvva\/d1fuzretd7qWilYvM21rJCVAOIV2oWUN1iqk0B8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8oUB8ofw\/AlrW9rO4IQIAAAAASUVORK5CYII=\" alt=\"plot of chunk unnamed-chunk-7\"\/> <\/p>\n<p>check the variables from imported shape file <\/p>\n<pre><code class=\"r\">str(FIN_adm2.pb)\n<\/code><\/pre>\n<pre><code>## Classes &#39;PolySet&#39; and &#39;data.frame&#39;:  136271 obs. of  5 variables:\n##  $ PID: int  1 1 1 1 1 1 1 1 1 1 ...\n##  $ SID: int  1 1 1 1 1 1 1 1 1 1 ...\n##  $ POS: int  1 2 3 4 5 6 7 8 9 10 ...\n##  $ X  : num  28.5 28.5 28.5 28.5 28.6 ...\n##  $ Y  : num  63.9 63.9 63.9 63.9 63.9 ...\n##  - attr(*, &quot;PolyData&quot;)=Classes &#39;PolyData&#39; and &#39;data.frame&#39;:  21 obs. of  12 variables:\n##   ..$ PID      : int  1 2 3 4 5 6 7 8 9 10 ...\n##   ..$ ID_0     : int  78 78 78 78 78 78 78 78 78 78 ...\n##   ..$ ISO      : Factor w\/ 1 level &quot;FIN&quot;: 1 1 1 1 1 1 1 1 1 1 ...\n##   ..$ NAME_0   : Factor w\/ 1 level &quot;Finland&quot;: 1 1 1 1 1 1 1 1 1 1 ...\n##   ..$ ID_1     : int  1 1 1 1 2 3 3 4 4 4 ...\n##   ..$ NAME_1   : Factor w\/ 5 levels &quot;Eastern Finland&quot;,..: 1 1 1 1 2 3 3 4 4 4 ...\n##   ..$ ID_2     : int  1 2 3 4 5 6 7 8 9 10 ...\n##   ..$ NAME_2   : Factor w\/ 19 levels &quot;Central Finland&quot;,..: 8 9 13 17 7 5 10 3 6 13 ...\n##   ..$ NL_NAME_2: Factor w\/ 0 levels: NA NA NA NA NA NA NA NA NA NA ...\n##   ..$ VARNAME_2: Factor w\/ 19 levels &quot;Etel\u00e4-Karjala|S\u00f6dra Karelen&quot;,..: 12 14 15 3 10 5 13 4 9 15 ...\n##   ..$ TYPE_2   : Factor w\/ 1 level &quot;Maakunta|landskap&quot;: 1 1 1 1 1 1 1 1 1 1 ...\n##   ..$ ENGTYPE_2: Factor w\/ 1 level &quot;Region&quot;: 1 1 1 1 1 1 1 1 1 1 ...\n##  - attr(*, &quot;parent.child&quot;)= num  1 1 1 1 0 1 0 0 0 0 ...\n##  - attr(*, &quot;shpType&quot;)= int 5\n##  - attr(*, &quot;prj&quot;)= __truncated__\n##  - attr(*, &quot;projection&quot;)= chr &quot;LL&quot;\n<\/code><\/pre>\n<p>To read dbf file we need foreign library<\/p>\n<pre><code class=\"r\">library(foreign)\n<\/code><\/pre>\n<p>reading shape file data from dbf file (ESRI files include normally: shp, dbf, prj, csv and shx data files)<\/p>\n<pre><code class=\"r\">labelData &lt;- read.dbf(&quot;FIN_adm2.dbf&quot;, as.is = TRUE)\n<\/code><\/pre>\n<p>Data checking<\/p>\n<pre><code class=\"r\">head(labelData)\n<\/code><\/pre>\n<pre><code>##   ID_0 ISO  NAME_0 ID_1          NAME_1 ID_2            NAME_2 NL_NAME_2\n## 1   78 FIN Finland    1 Eastern Finland    1     North Karelia      &lt;NA&gt;\n## 2   78 FIN Finland    1 Eastern Finland    2     North Savonia      &lt;NA&gt;\n## 3   78 FIN Finland    1 Eastern Finland    3 P\u00e4ij\u00e4nne Tavastia      &lt;NA&gt;\n## 4   78 FIN Finland    1 Eastern Finland    4  Southern Savonia      &lt;NA&gt;\n## 5   78 FIN Finland    2         Lapland    5           Lapland      &lt;NA&gt;\n## 6   78 FIN Finland    3            Oulu    6            Kainuu      &lt;NA&gt;\n##                         VARNAME_2            TYPE_2 ENGTYPE_2\n## 1   Pohjois-Karjala|Norra Karelen Maakunta|landskap    Region\n## 2      Pohjois-Savo|Norra Savolax Maakunta|landskap    Region\n## 3 P\u00e4ij\u00e4t-H\u00e4me|P\u00e4ij\u00e4nne Tavastland Maakunta|landskap    Region\n## 4        Etel\u00e4-Savo|S\u00f6dra Savolax Maakunta|landskap    Region\n## 5                  Lappi|Lappland Maakunta|landskap    Region\n## 6               Kainuu|Kajanaland Maakunta|landskap    Region\n<\/code><\/pre>\n<p>merging shape files area name and centroid points (x and y coordinates) labelData as data frame<\/p>\n<pre><code class=\"r\">labelData &lt;- as.data.frame(labelData)\n<\/code><\/pre>\n<p>merging labelData and centroid data table fields<\/p>\n<pre><code class=\"r\">label &lt;- data.frame(labelData$ID_2, FIN.centroids2$X, \n    FIN.centroids2$Y, labelData$NAME_2)\n<\/code><\/pre>\n<p>Event data function needs right column names<\/p>\n<pre><code class=\"r\">names(label) &lt;- c(&quot;EID&quot;, &quot;X&quot;, &quot;Y&quot;, &quot;label&quot;)\n<\/code><\/pre>\n<p>Note: Above routine used in eventdata function<\/p>\n<p>#######names(label)&lt;- c(&ldquo;PID&rdquo;, &ldquo;X&rdquo;, &ldquo;Y&rdquo;, &ldquo;label&rdquo;) &ndash;&gt; USED IN POLYDATA FUNCTION<\/p>\n<p>do the right data format into labelling <\/p>\n<pre><code class=\"r\">data &lt;- as.EventData(label, projection = NULL, \n    zone = NULL)\n<\/code><\/pre>\n<h2>Labelling into map<\/h2>\n<p>#col = 1 &ndash;&gt; font color 1=black, 2=&hellip;.<br \/>\n#font = 3 &ndash;&gt; font size <br \/>\n#cex = 0.5 &ndash;&gt;expansion factor<\/p>\n<pre><code class=\"r\">plotPolys(FIN_adm2.pb, projection = proj.abbr, \n    border = &quot;gray&quot;, xlab = &quot;Longitude&quot;, ylab = &quot;Latitude&quot;)\n\naddLabels(data, placement = &quot;DATA&quot;, polys = FIN_adm2.pb, \n    col = 1, font = 3, cex = 0.5)\n<\/code><\/pre>\n<pre><code>## Warning: The data&#39;s &#39;projection&#39; attribute (NULL) differs from the\n## projection of the plot region (LL).\n<\/code><\/pre>\n<p><img decoding=\"async\" src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAMAAACR9g9NAAACAVBMVEX9\/v0AAAAAABcAACAAACUAACsAADEAADkAAEwAAGUAF1EAKAAAK2sAOY8ATIcAWGYAZrUgAAAoUXsrKwArK2srTEwrTGsrTIcra6MxOWY2TGs5AAA5AGU5OQA5OTk5OWU5OY85UZw5ZmU5Zo85ZrU5j9pCTIdIOldMAABMACtMK0xMTABMTCtMTExMTGtMTIdMa2tMa6NMh4dMh75ROgBRURdRezlXOSBXOVFlAABlADllAGVlOQBlOWVlOY9lZgBlZjllZmVlZo9lZrVlj49lj9pltbVltdpltf1mFwBmZgBmZldmjzlrTABraytra0xra2tra4dra6Nrh2trh4drh75ro757UVF7ZlF7j1d7tXuBZlGBj4+B\/v2HTACHTCuHTEyHayuHa2uHh0yHh2uHh4eHh6OHvr6POQCPOTmPOWWPZgCPZmWPjzmPj2WPj3uPj4+Pj7WPj9qPtY+PtbWPtdqPtf2P29qP2\/2cUTmjayuja0yja2ujh0yjh4ejo2ujvqOjvr61ZgC1Zjm1ZmW1jzm1j4+1tWW1tYG1tY+1tbW1tdq124+1\/rW1\/tq1\/v27jyi7j4++h0y+h2u+o2u+voe+vqO+vr7ajyjajznaj2Xaj4\/atWXatY\/atbXa24\/a\/rXa\/tra\/v39nDn9tWX9tY\/9vFH924\/9\/rX9\/tr9\/v3Si3vGAAAAq3RSTlP\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/wC3Kj1sAAAACXBIWXMAAAsSAAALEgHS3X78AAAgAElEQVR4nO2djX8cx1nHu4CqCKwgXhNIG0xjQ1ogpZi3ONRtnLgVL4lV05h3A6FFTkyE4tJUJA6tLq6CG+hBYrOcFdr6OHGn\/SvZ2deZ3ZndmdlnX3Tz+9qf092+zO1nvje78\/Ls7EcC4CQf6fsAQD9AvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKNAvKOcIPEzj7HZ92EsCSdH\/PEre+Hr\/GPjvg9kOTg54hcXpuHr8auceO9E0VvOSTlB4s+XSnyneekzwj+2+0O8LeVrfJd56SfGIX4AdHj+9DPxtuYh3pL4VC\/QvfbA3jzEW7L47O9tFBa1mpe+n1zVRe\/WJ3uItySs1c88UX2beZm4zu1za6zUQ7wlUXNu5q11U6vP5ZY1l34KWkC8JXE7XqAT8fLV5ilCPB3t5WVdkYb4XmktL2tP5RYne4ino6281LFqrB7i6WhNvN5WZu4hno528tLAp4l6iKejjbw0PIXrF3uIp6MV8W3tAfF0tJCXFl0zmrtAPB3tnOrNd9HaCuLpGIp4rV0gno5hXOM194F4OoZxjdfcCeLpaEO81U4Q3y0tRF7ZRVlAfLcMpcRDfMe0UOBxqj8JDEa8zkYQT8dAxKMd3zXDqNuhy7ZzyPPSpsBjdK578ryUREHb0KJ3iCcky0vZfQ82tDY0F0A8JWleZrnf1Lzx7ojA6YUkL\/kb2xqlZ\/67MdgB4ulIxedLmhV5471NdoB4Osrim5o33R7ie6F8qm94sjcWb7AtxNPhJRV5siIP8ScDL7RcmpemiXmc6k8GXlTcy\/evWyeI5tzJwEsmoirkfgPxhruixPcDi8CRZL19kTfcE\/fO9QTLS1nmW5s3vXvKZHOIpyPKS6l5u\/SMTOI26f6IxUtWWBb5Nr1DPCFJiScp8objuubDwBBPR5yXMgPm1XPfcB9Zc6ISiKejNCybo3IiX25bG4T4fuACMcqT0Un3kBfsLrxDPCFcXuqKD2S\/kk5afxBPBy9er\/vOl63rptkP8XQIeVmyKal35xPUSpYagw6cvhDFS4ZqxJF6zrfwo+ikwEM8IWJeypt1acGXTUKdLOikwEM8IcW8VLbhFP3q6XT0dt+OLtveKOWlwry6k823j8ZHX31\/lPPSPITG+i4M9NX3h0S8qUf7phzE94c0L8273K0w3g\/i6ZDnJcHNkxpAfI8oJj+CeA1Otnj5Yl\/aa0cMxPdIVV5G6luUj8pdj9TkZbulHuL7ozYv2zzjQ3x\/6OQlxMtZfvFtlXl02faHXl62Zt5sc4inQzMvWzIP8b2hLb4V8xDfG9p52Y54BGL0Rb\/iEXPXGz2LN7uGQLwtM29j4q3ucUv0xbfWj6OfMsRbcvzqeLIRLC5M80UGJb7FHjzMXt0uofLZWaY\/X2SSl62O1+gkDvGWNCnxQcuD9DrmId6WBtf4YABDdRBPh2f0\/LGezUO8LbP1B9ftS3zQ7tke4lvj+JW9yVn7a3zQ7tm+PmmItyRUPtm0r9Uz2jzZQ3xr3B7PzwTz001KPO0BmSUO8bYsng4rc2tcgTfMy5bjbuuSh3g6DMW3dBSa6UM8HQMTX\/0FEE\/HsMTXmId4OqjFL24mb+ZvpG8uGX1D1VdAPB3U4kPfM291x1vd3Yn+7Ij9RRpfUfEdEE8Htfj7e8Fk9eLu+Ojia+zP4fPTw3HtTuJ3qL8E4umgFj\/y1j7nrd\/44usHN7z13enBt731ae1Out8C8XS0UrlbPPeO2fld81sgno5BdeCk36JaAfF0DGh0rv5LIJ6OAYpXn1Ygno4BnupRq+8Co7zse04kiKfDLPKqtcPgvwTiu2AwEXc6XwLxdBie6lufCgt99R1hnpc9RtdDPB1DE1+5FuLpGJZ4BGJ0xqDEI+auO4YlvmY9xNNhHlff4x0VEE\/HgOLqcSdNlwxokAb3znXJcEq8RtIQT8dgSjwmRuiWwZR4TIXSLca1+laOQvMXBfF0DOQ2ab1kIZ4O47zsczZjiKfDXHwb5jXThHg6jCJwIloQr\/tjgng6LPrq6c1jStPuschLevMQ3z02eQnxMe6JpzYP8d1jlZcQH+GgeGLzEN89dnlJLR7Nuc4ZRInXTRDi6bCr3PVUr4d4OgbRnMMgTfcMRLxemhBPx2DEYzy+W4YiHqFXHTMY8XgmTbcMSDyibLtkOOIRV98pgxiW1U0X4ukYRCBGmnLdBhBPh3HoFW23nZAaxHdIz+HVvqAed8t2R883VETJ+Rk1W0M8Hf2W+Cw139dJHOLp6DW6uiQa4jujX\/GGiUM8HRBvj1PiiRtzZqlDPB19V+4KyUF8V\/R861xJPKY07Yi+75Y1OtlDPB1999Ubteggng6r0TlC9eWKPearb4HZ+oPr4jM\/+72FSvYbwhMq6Dl+ZW9yNlhc4J792LN4s9Qh3pJQ+WQzOH6Ve9pr3yXeKHWIt+X2eH4mmJ8eTomXmYd4ehZPe563xj\/eue9752TilalDPB2e+exHAY159RA8xJNDVKunEq9IB+KpoarV24kvFHB1GhBPDVWt3sp8VMCbmYd4W4hq9YFNxd4v7ac0D\/HUENXqAwvxiWWtIg\/x7dOV+DysMl+CU32PWIuvHKopr+R918ZSq9ZBPB1NxKvkRQ214qLiNtWpyxdDPB324lUn63ipUUhVOQ35Yoino1le+mIlnTsLiJU409YfxLeOfV5G5Zq\/8U205fs6V\/OKtCVAPB2NxadiyQZuKn4qEE9H01N9qonKe+X5AeLpaJiXfpB0xJKV96qVEE8HRV7aXcdVaVWshHg6iPKyG\/MQTwfE2wPxEY3Mx1cLZWhGAsTTYRV5JaVRmU86AmtOHBBPB2FeNhCf+4b4jiDJy3KMhcG+YjmHeBJm3sZos3ILGvHWXThmN8hDvCaLZ\/7tzOKZcdUmBHnpC3+s9q1YwAPxmnQjPsGmxJf3gXgSZp7ndXCqZ1gNwhndJA3xlPTbjjeaDwPiKemrxMtH6muSgHgtoljqQjR1iV7EK3tqIJ4Gds+McMNUmT4qdxVx1dV9dxCvCavRD7BWbxVhy4B4XTqt1dNsiRLfDT0My9YUeIjvBKq81IvC0Qm6hfjmlO+RLNE0L7ko2xqjOmGZNYlAvAGL5ypXNw625Mq69bWbT64KiDeg3Vp9Glwff6jyplUHgHgaolP9RuUmDSNw8nN33dVbp8TjmTTd0bTE+\/69T16dv8He3bty7ya3iv3qWFNyES9cXPGjD\/HH4\/fkqbV6sNQMVzzrtRNmvCnTPC\/vfvnm0Y63etl7\/JG3vdUdb+1d71k2s9bRXli\/mHlrr535kC3Y8VbCrT4ffgjfru7uiJNvJUA8CaOor3695S7bo73DX\/v1\/aOL27fGh+PDz3mrF9nfcXB\/L5hfCkZn2YeD6dHF8QdPsNfw4274\/7Vws3JaEE9DTT89o3le3t6++8jEO7U1vXbuYHrwbW8jfA3\/sVP9+nTEYr8OphNv\/Yb3MHsNPxzshv9vSH+QEE\/B4sIP2m\/HZ9fl+v6b5rU7iKejeeVO8la6HUFDHuI1ab1yp9lFH98loyUeoVcEdFC50\/RukGDVuQHiNWm\/cte8FJukCvF0EE5+pNyILFWI12XUaq1e9wpvnrBiOcRrsnjm3y\/8oL3ROd2aHVnCEK\/J4pnvvvrd9mr1Wg1zu1vq5IshXpfJxqRudK5B6joRN3YJQ3zrtCzeLl1c4wlorwOnvo\/WtsCrVkC8Aa2Jr4+aVK5f3JyxoXr5oDzE09CW+Abeg\/k\/bASzeFB+Z\/Mo\/C8OzUM8BXXim8xbXr2BetVk4\/Z4snpxdzcem2ej91pJQ7wWrd40WX\/5rtji9nhx4Ya3vrs7PbjBRuqLwwkQ3z4tPn7Mfv4ziG+fNh84aG0e4tvHUrzeVhA\/XNp8mnSDc718V4ino9XnxzeYzly6I8TTYZGXRjYbFPryMoinw0p8WxvX7QnxdNiIb3FrYc\/SrhBPh2leUj1CTuuriksgng7DvCR7dqDWdxWXQHwTxOmsB1ziy7tCvCWS3nuLEm84haXZF1TuCvG2TDYalnjzmembiC\/uC\/HWzDYaiOcnOtLfyeAL6naFeHtm6zfsxefGDcQTDtJBfAPmjzUQL3tbt1ODbtviAoinwzMIweHFa8r09Z4rVvd1CRBvy2z9wXVx8hmzazz33nQPY\/fosqXj+JW9yVnxHlrtvCyUWou73Q3VY5CGjlD5ZFMMv6zPy6QFV9SmNXNteYm2egzLUnJ7PD8TzE8blXhfMZ+F5e0Uwtm\/aqVkX4i3pTyrsZ54mcTaOYtrV\/AJ+4WvQQROy9TlZVGIuKpqx4ok+b\/iF0B8V2iJVxVfO\/FpatIfU\/atsh0hno6avPR9vkQWdVROcFaZqmrn6h5hiKdDT7zqlG8pPlCLr5wXDeLp0LvGB7Gn0jnfWrzSb3RhUe0G8XRoi88+i2vVO2okq2gr4KbJDtA61ecfiwW1vu5euVqxDa7x7VOVl8fXo0cNhBqShw3cfcm\/96Tw+IF7V4JA\/viBRlPgKXaGeDqq8vKDq74\/Yw8eiJ4v8Ky3unXVLz9+4HK8mv0VHj\/QRLxiBcTTUZWXb4XiR9GDB9jcBZPVi1vbfunxA0c\/u88eP7B9a2tbePxAw0kv0Y5vmaq8DE\/ra69FDx5gzxcYeeu7U7\/0+IGN+PED+we7+y9zjx9oeiO1tNIP8XRU5KWs1i3TkVXExWqgDoY9QBBPR0UAjvbQTN70ztsAmuKr1iEQo030vVf02kta+s1LPO6kaZVyXvrZabvcS1fR3SZsqB1wUS2+uBbi6UjzUgiOiz4XS7d6\/KTYj2sQW1nbvyc92IGwLOKTi7PvSwdR0mu35mVf8+vrfiCieoinI89LbjAmjYouroveFlPQOAdUULshnxTE05HlpXBRjz9kpwAhJqqoSjG6qkn9hlxaEE8HL14o8eJnbg8d8fp32mhsCPFtkF3j\/bS5xhHIym6xxtVIvM6pIf9CiKdDFB8UtCub7Xy7XZasdq3eqGcX4unw+PjmUnlXiqntnCVpzhW\/DuLp4K7xmXuuTOtEx8vXa327diUQ4qkRxZdb62rzlRuw4I3ze\/M32PskTGPxtHgrB4va0O\/o8YWDHQbLIV7ocKsfZssr\/4r1wfyfbh6xR1B4j5\/5MIrPGG2yx1JcHN+PlsZRGxDfFywvSyMsQcWCfEV6UZCvDo72Dp99noVqHEb\/wmWjjcnqxTtfuxQFcCRRG0Z9uxBPBy++2GKPlqnF1AytsRs0Jyx4g8VsHEzD0u6dHXnr3z89jZYeJFEbEN8PaV6y2ynPhcU3jp5MYikXN+OP0sdF1YXRah6AZpGPNoN4OtK8ZEGUn37fe+jFR96OgyjZ\/8+feTu\/HBf3rB1e0TwCbfM+xBOS5uX9q\/78M\/61R\/MgyifC1+1b8XOiuCDKDK3AeR10q\/YQT0kaerV40js1veY9LAZRnguvzvnlWIRMvH6jH+LpiCp3spHYIKiTV3eu1z4G7Zo9xNNR6LkTV1pX3wxicILFJ6+KHT3x80cj5l\/nq5UQTwgvvlT0zAKjtHcUmX\/5itDRkzx\/dCesY67dEe\/NgXg6ShE4PHX+1OINzB\/t3RI6etLnj47ZTTxCtRLiCSkFYnDWaqtvlp07IrfHdx\/hO3rS54+GdcwNsVoJ8YTw4\/HZu3RZ9a6GQfHq7TTrA6jVUyJc4+M3FQG1ItXm9b7fIBYb7XhKcvHcG1nTTkZ1kdcR6itn1JJtCvF0JHkZ5v+9ZLycTYOwuMJZY9343FwIaTd++NdXdOPHSVZMZsNto42PvnpKUvEh8Xj5zrkP2DwHYXNq7V3Wa78XlOZCiLvxw22KcyEUqfcK8X2Ri2fj5VEv\/S02DcLh56IZEVj7qjwXwhNJ62t3W96Nn0FY5H2IJyWv1bPx8qiX\/hSbBuHg29GMCKx9pZoLYVqYC0ECYYnHsCwtufjiGx1qb32rT8Ak3BLi6eDa8clfk73r73msXY9AjH7gumy5V12alXjVvTrKlCCejlKJ190xmgQv2ilp5s3fiGoDD96L3sbcvVSRgtFJBuKJsRYfNvICYRK83b1I+FESqhX9X+FaepIhoHIwtwqEV1PjZe255K\/ujvf3fP8aPwnebtjw+7h3Ng3VCv88v3+YzXonGfP1y++UxFtAPB3FvnoN8fGW9570HnpRnAQvuB16TkO1dlmL8NR+2htf6sHlP9aPChQOdhgsk3its24SqaX1G0mNl0bh6u+4DUrrIZ4OYUaMumGX5FVvBIdPVkxCfJelWZFI8WCHwXKIT4bk1Fumt1HrDb9ke6k+Fs\/7ld9cOthBsDziK58Lwj9L2qyLXfa5dCbQieaBeDqKMXfq218FUdrpK8SbxPflayCejsKpXnluLlTJtdMv+VXVIdVnGohvAy4vpdU7cYKMfKlm8uUNVbFWOmcaiKeDz0tJ13lWoSug3ZqTLpTuKk9P2Bbi6SiJl1zMZUq0m\/HaRyLdVFwI8XQIeVkq8tLLfrpOuSbb2aDZJ0+s8GOAeDoqT\/WV3bjVjT+DekC2Q+0iiKejSnzNsGl1iTb2LkmrmAbE2zJjgdLHr3LRkeVafbFzrdJgVX+PKfVNPIi35PiVvWC0qRYfvWaftOIk2hRfWgDxliwuTIPgnW\/oic+XVyWp0++iiXrAPgXiLWElPnz5Bbl4oe3G9anXVd0VGzQu8uUEIN6W+cfGyUuKeI3PutX8gn81yhqeTZGXD+FmQDwdnlcYkRcLsX4gpGRxozIv2xvibZmtP7gu3utW7rmL3kmWqSEUX12zgHhLwmv85GxcxUsRxdf3ochQ9L2be6\/7wUG8JaHyibodH+j0mkqRi7+XzVOePl8+PQzJBOaFhOTfCvG2sImFg\/lpRYlXDcfUJqso8cl915tHhfuupROYi98l\/06It4XNcSCWNrE5Zyde2Z5L7rs+LNx3LZ\/APNknThDi26YgXmq+Lg3FBu+n912vF+67Vkxgnial7jqAeDoKp3qrIl8Mz9HdT5ZU5dAPxNOhI742kYL4rPfH+Gjymy+kQDwdhbyUxsjVJiJswc+YZ3gwabehaj3E06EhnlUJV7\/xXiA0ycRN+H4Xrsiai6\/eC+LpEEOvZOKj291v39nxnv3C+fT2Z89bjxpocStMiN3gu\/pNi7xfsxPE01GIuZOW+Psf9zbDJtlkdS+6Gzp8e\/Dg9fBt1grjxvCzneLQLJNDqZ8gA+Lp0BGf3v38g+dGcfPrweuLm6OooTYVdhJLvkmIbZCfJXQPtn+WS3wx4+sbc7K3cXE3Ea9zdYB4OorNubL4uhS4XWqG06uT0fg6iKejLL7YKK9NgkC8ZoAfxNPB52XYbHtoO\/CTJ8PEy3xlIy5HOrSiK16oDdTsBPF0eHxmjs757yfzE++EbbfNox1v5VeTRwJXINWmc6Yo9tbU\/Vggng4xEGP08JuswRY11fYOo2G1re2Jcppi1rOz972vfuur7\/3LV\/fzcjt\/zFsbcxIXN+XPKC1eWWpPEhBPB5+XbNTsGhtLY622tO1249xIOU0x69l59xOf2tq6\/DOfSvt22KDv\/FKY1ko63r722pn\/kUxuXu6drT1JQDwdpcpdcYOqYsh6dsLzwVYyr93uNht0ZxMfbgaTlafufO0SO1lE852XJzc3+ybZwfbOUokvrK7uhGE9O+H5YCuZ1253nw26h8sPx+HiU98\/PY16edh85y94p\/blg7fq76472N5ZIvGKmUtqeHOD21rStss68LNeXGkXofHB9s0SiS\/5MIyrl4fk57EV6ZRpklR1vgni6agZljW8k0aor\/mCeoOeXyUQT0dd5a7SR3bTlaIBz8fT1JjX6u6BeDrqAjGqfMVa\/cL8SGI553rlqg5Dr58P4unQaM6poqfL8Rfph\/RNelFPPynR7N+FeDpqx+MDlZZibaD0NjkhJJ8XTxbuq+GeE+\/f\/bryuYXKg+2fpRFvdqu7chivdFGPW3DXzvH31QjPiX\/oTuVzC6UHOwCWSbx8o9LixdPeytXvfefBd9773nfSWyH8pJM+PcWzKM24kz651Av31eTPid++VnxOvNbBDoClEa+8Dpd+EHdf8m+\/+4nP7u7ufOKzaSf9Q1kn\/VPjty57K2Inff48w\/i+muw58Vv7pefEax5s7yyNeEWJlwRnvBV30u+mnfTszriVpJN+9ak7f\/yZUie9svZmEq8B8XRoXOOL7fQk\/JIN3qWd9OzOuIeDtJP+mz+9fy3tpE9KsrLabhSoA\/F0FEu8JERSIj5hssEvFvYppqH6fqPIPIinQya+3FDTiI+qHnsjKfAQTwgXepV3wMpGXdLhFVWbTyzwuh2ARt4hnhB+zitJmRY73tMxVVlCdeIVtQe7gx0GyyKe85ANoiZL4ii6fMs0ii4KuksXZrtKzgo6tUaDgx0GyyG+MJZekBc20N\/PZ63JGuhR0B2b2GYnmeNmUzXiLlFsessFxFMinOqD\/ENBy+yc\/2Y0a43YQI9upxwfTI+e2j76uLf61HZy3pBqtonxUBzsMFgG8eViKmo53A7iWWtGQgM9up2SzWTzwhdfn3gPv8zi6tjgyzlh7+Sp4vfEWzPMvUM8IdnDpJXV9WyD4qI4qj7qsP\/XX\/rrq6nK\/\/sJb+0rV99KbqIPT\/8f\/cc\/jK4Fl3\/jfHQ5SPpw7Q92KCyDeGmNTKC8Lo6qjzvso775lctsjGby4788\/l\/vR373r\/75z0L33uof\/dyf\/pb3w7+9evG1lav+tUfD68JeYHWih3hKpM05GeVV97kO+63tW1vjW7\/CSvPtd79w4T9+6Of\/4s9\/8b8\/+jvej\/3kf\/7ln\/ym96N\/kNyj8fCbyU31jQ52GCyH+LrWlV969hDfYb+1f7C1\/\/LfxlH128e\/zzrsC\/WB7EsCW+8QT0hWuattVitq61GHfda3k1UW1GkJwfdWBzsUlkF8hEZ\/inLkNl2bt+PV4zKm\/TYZEE9H6SlUdRSH6rjKgWhbnZqtd4gnRNZzV+2lMFInDaxOP0p\/S9beIZ4QSa2+\/qaXvD++cjBW0VCw9g7xhPDi7z3pnXrw975\/941k2fwSNwyTk\/4+yuLLvbLlRfYFHuIJ4cXffcn3P7jsrWxdTkZdpMMwhZocj6o7Xrw2EBzsMFgO8aGTt37Ke\/QW64p5cYWFUD4\/PeSGYS6O2TAMi5zUj58TBvuSN0QHOwROtniuyP\/Ntu+\/vLX\/8tYL3qkofpIfhrnBhmE2Dqa8vPw6sLiSn\/cXzyXDMv5\/vRcnLNDkYBvs2wInWzz33g+48KqCoMVz7xQiLoJ42G0WD9A\/8ra39nfxI2fCpfejGbMe+splz\/NOvc368Vf3IH5Q8BE4gndllBy3IrwObCQD9Nsjb+Xi9i12i8zRXnB7d3xrPGLXjZe\/9aU3V556MdyoplvH6GCHwXKIT58KIpiXB0xlb0fxw2VYh\/y5aDg+ukVm\/pj36DXv1Avhgq1vfenelej9+rRZhb5wsINgOcQLxvl3xT2UsZZciGZhf7\/JyIzsYIfBUokvVsV0+1\/4MM3Cb6hpXV5ysMNgmcSXvKdt9iDfTpYM17YPOOf2A3HVBzsMlkm81Dtf8Bef9pMQOmGa0ly82JlPqR3iKZHcSBMUvAun6rsvJW21tTvCo2m4nwh\/ieCj7SkPdhAsm3hRvhiM538QtdUOWQtOeDRN4TxBW9LLBzsIlk58+Xqfb3\/3sbQFJz6aRlIxaPNgB8HSiS86FNpp5ST4Df3C3GdtHewgWA7xKutBdnVWPV8o3bRF48WDHQZLIb7cYyuUeLVX6Zmh5YMdCCdI\/Izd4pDOLheRiY9e5EU+OYNLU9QM2qEB4i05foWNsM0\/xk0sVn7unJ+d0Hnxii68rsp6DMRbsrjA7ms4flUu3s872tJXoQO2kFjVeaAlIN6SxfnqEh\/D1eGFs7+wiern0CoQb4v6Gs+R1+Y56Um9vbii9SMWgHg6vDT2StqcE+t5vl+uAnZ7sJ1+Wy0nSPxs\/cF1MVxa6MCpts\/1x5fGbroB4i0Ja\/WTs0kVL6Fwt6xSefFX0Yd3iLclVD7ZrKjVK7rvKk\/8XQLxttwez88E89PyEh+hJb7TY+aAeFsWT4eVuTV+YnjxNukg8BWn\/H7abwUgno7i7NW19HakAcRTIu2yFco3t5A6lMoUiKejXOLZXy+zz3fU58W9xkCz1er1EE+Hx\/9hxpOPcUvdixvvfrI6y3lPtrf2X+vNIZ4O70TRd26JnGjxwB6IdxSIdxSIdxSIdxSIdxSIdxSIdxSIdxSIdxSId5SlEn983fPWp\/Gren34ZnK2YvXE8zYqVo\/K8+MmjKLgb\/X6QbFU4mehsNFm\/KpezyL0peJrdk9Wr09n8p9VuH7xzFi9flgslXjGZDN\/Vayfn\/5QKj5ePXlWUeLj1fPT05lydSieRQI\/M1ZtMByWTXwsRa0mKZYq8eHq0Yb6d8PSnSmuIwE7y58NipHAQ2XJxI828lfV+gkbHJebT6UrfhjR6g3x\/j2RyVmU+B44vr6ZvarXByqx8eqZqsTHq9ViWaKTs7jG98AoKssjZYnO1sjFJ6snlXuzur3ihzWKKgeo1YMhA\/GOAvGOAvGOAvGOAvGOAvGOAvGOAvGOAvGOAvGOAvGO4rB4YR4lDjb4dsCteuckjLkYA\/EyhNn0IH7JSMRHg6nzxx+LwzTXPr8XlviRtzE\/EykfsSXR40j7PlxinBcfh1CG75OZM88z8WGJj8WzSK3ze8rYjBOM8+LjSKtYMwuufUcQP4mWVARrnVicF5+U+FgzX+JPT+OyHp4JlLGZJxiXxT\/msSdcuEQAAABxSURBVMt7fI2PxGfX+MXTbKbs+Oq++nh0jVfFW51UHBZf5qREyFIA8RwskHIJz+pSIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5RIN5R\/h9jpEOFiuNWrgAAAABJRU5ErkJggg==\" alt=\"plot of chunk unnamed-chunk-16\"\/> <\/p>\n<p>Have fun with mapping and RStudio<\/p>\n<p>Marko<\/p>\n<p>ps. This post is produced by using <a href=\"http:\/\/rstudio.org\/\" target=\"_blank\">RStudio<\/a> and <a href=\"http:\/\/yihui.name\/knitr\/demo\/rstudio\/\" target=\"_blank\">knitr<\/a>. Thanks to <a href=\"http:\/\/yihui.name\/\" target=\"_blank\">Yihui Xie<\/a> for amazing work with this issue.<\/body><\/p>\n<p><\/html><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Finding Centroid point from shape file and labelling map Finding Centroid point from shape file and map data labelling In this test I try to find centroid point in shp file. In additionally I will show how to use this centroid point information when to labelling area name into map. R Library PBSmapping offers wide <a class=\"read-more-excerpt\" href=\"https:\/\/science.ekqvist.fi\/blogi\/r_maps_spatial\/finding-centroid-from-shape-file\/\">[&#8230;] Read More<\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"gallery","meta":[],"categories":[22],"tags":[28,39,43,46,49,50],"_links":{"self":[{"href":"https:\/\/science.ekqvist.fi\/blogi\/wp-json\/wp\/v2\/posts\/262"}],"collection":[{"href":"https:\/\/science.ekqvist.fi\/blogi\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/science.ekqvist.fi\/blogi\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/science.ekqvist.fi\/blogi\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/science.ekqvist.fi\/blogi\/wp-json\/wp\/v2\/comments?post=262"}],"version-history":[{"count":0,"href":"https:\/\/science.ekqvist.fi\/blogi\/wp-json\/wp\/v2\/posts\/262\/revisions"}],"wp:attachment":[{"href":"https:\/\/science.ekqvist.fi\/blogi\/wp-json\/wp\/v2\/media?parent=262"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/science.ekqvist.fi\/blogi\/wp-json\/wp\/v2\/categories?post=262"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/science.ekqvist.fi\/blogi\/wp-json\/wp\/v2\/tags?post=262"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}