【小花經濟學術】如何用Stata進行地圖數據可視化

【點石成金】 2017-03-12

展開全文

中國縣級地圖（Stata制作）

一、引言

在經濟學科學研究發(fā)展迅猛的時代背景下，地圖數據可視化在經濟學研究過程和研究成果展示中的重要性是不言而喻的了。Stata 是經濟學科研人員使用最多的軟件之一，它簡單的操作和較小的內存占用吸引了大量的用戶。習慣 Stata 軟件的粉絲級用戶總是希望用 Stata 來做更多的數據處理工作，小編作為 Stata 的愛好者之一，同樣有挖掘 Stata 在數據處理中無所不能的、無所不包的美好愿望。今天，給大家?guī)淼氖鞘褂肧tata制作地圖。

已經存在有不少的 Stata 作圖命令，但是往往由于Stata 地圖底圖數據的缺失，使 Stata 制作地圖的功能受到了很大的限制。今天的推送，能讓你在Stata上完成 GIS 軟件上制作成的任何地圖（小花之前也推送過如何用 ArcGIS 軟件進行地圖數據可視化）。

而且在一般的地圖數據可視化上，小編認為 Stata 能做到比一般的 GIS 軟件（e.g. ArcGIS）更簡單。下面就來開始制作你的 Stata 地圖吧。

二、數據準備

一、坐標數據準備

其實我們可以很形象的把地圖中的邊界想象成一個一個坐標點連接而成的線條。那么關鍵在于我們如何獲得這些構成邊界的點坐標呢？地理信息學科已經為我們提供了大量的 GIS 數據了，我們可以找到各個地區(qū)的 shapefile 和 dbf 文件。（點擊 “閱讀原文” 獲?。?/p>

在 Stata 中已經有命令 shp2dta 幫助我們將shapefile 和 dbf 文件轉成 stata 所需的地圖坐標數據。比如，筆者有整個世界城市級別的 GIS 數據。

shp2dta -- Converts shape boundary files to Stata datasets

clear

set more off

cd c:\users\admin\desktop\Stata_map\

ssc install shp2dta, replace

help shp2dta

shp2dta using world_city_level, database(world_city_level_data) coordinates(world_city_level_co)

這樣你就可以看到兩個 dta 文件，world_city_level_data.dta 和 world_city_level_co.dta。其中，后面這個dta文件是坐標文件，包含_ID, _Y, _X。

三、畫圖

use world_city_level_co.dta , clear

scatter _Y _X ,msize(vtiny)

同理，我們可以將中國縣級的 shapefile 轉成 stata 文件。

cd c:\users\admin\desktop\Stata_map\中國縣級數據\

shp2dta using china_county, database(china_county_data) coordinates(china_county_co)

use china_county_co.dta

scatter _Y _X ,msize(vtiny)

把 shapefile 和 dbf 中包含的 GIS 數據庫與坐標文件 china_county_co.dta 合并

merge m:m _ID using china_county_data.dta

/* 畫出廣東省的地圖 */

scatter _Y _X if NAME_1=="Guangdong",msize(vtiny)

spmap face using Guangdong.dta , id(_ID) title("Stata展示暨南大學各校區(qū)位置圖") point(data("Jinan_Point.dta") xcoord(yc) ycoord(xc) fcolor(black) size(Large) legenda(off) legcount leglabel(暨南大學各校區(qū))) label(label(po_name) xcoord(yc) ycoord(xc)) fcolor(Oranges) legenda(off)

/* 畫出暨南大學各校區(qū)的點圖 */