python如何导入arcgis数据_新晋克拉克奖得主Dell地理断点论文技术分析III：ArcGIS地理数据导入stata形成计算数据等...

Almost 100%抹平与世界前沿的技术处理差距，希望后续者在此基础上创新超越！可惜中文依然无法主导世界前沿。1982年中国青年足球队2:0战胜英格兰队，那是一次实质性提升，而可惜没有延续。现在的中国人在社会科学研究方法上几乎不存在瓶颈，但在英语主导的学科主流，依然不能推陈出新，达到世界主流的认可。始终处在追随者状态。这是需要深思和反思的！Dell的很多研究都使用断点回归，最有名的是mita这篇，地理断点回归的方法论贡献，她使用的数据背后所凝结的艰辛(唯有Dave Donaldson可堪一比)和她实证策略的神乎其技确实是每一位regression runner高山仰止的典范。在用Python操控ArcGIS后形成地理断点回归所需要的dbf等数据，然后用stata读入地理数据，再转换为stata格式数据，需要安装AccessDatabaseEngine_X64插件，然后在stata平台输入:

odbc list//出现：Data Source Name                   Driver-------------------------------------------------------------------------------dBASE Files                        Microsoft Access dBASE Driver (*.dbf, *.ndxExcel Files                        Microsoft Excel Driver (*.xls, *.xlsx, *.xlMS Access Database                 Microsoft Access Driver (*.mdb, *.accdb)ttest                              Microsoft Access dBASE Driver (*.dbf, *.ndx-------------------------------------------------------------------------------r; t=0.00 18:04:16. odbc query "dBASE Files", dialog(complete)点击dBASE FilesDataSource: dBASE FilesPath      : G:\mita\data随后操作下框命令形成stata格式的gis数据：

*This file prepares the GIS data by district and grid cell***********************************************************clearset mem 400000capture log closecd G:\mita\datalog using prepare_gis, replace/*Data Files Used  dbnd.csv  elevdis.dbf  sldist.dbf  d2_potd.csv  d2_huand.csv  dis_xy.dta  inmita.dbf  outmita.dbf  gr_coor.dbf  e3.dbf  s3.dbf*Data Files Created as Final Product    gis_dist  gis_grid*Data Files Created as Intermediate Product  temp*/***********************************Geographic controls by district*************************************----distance to mita boundary*gives distance from district capitals to the nearest point on the mita boundary, as well as the*geographic location of that pointclearinsheet using dbnd.csvgen d_bnd=near_dist/1000label var d_bnd "distance to  mita boundary (km)"**4 bfegen temp=.replace temp=1 if (near_y<=8400000)recode temp (.=0)tab tempgen bfe4_1=.replace bfe4_1=1 if (temp==1 & near_x<=824796.228990515)recode bfe4_1 (.=0)gen bfe4_2=.replace bfe4_2=1 if (temp==1 & bfe4_1 !=1)recode bfe4_2 (.=0)gen bfe4_3=.replace bfe4_3=1 if (temp==0 & near_x<=808697.919841043)recode bfe4_3 (.=0)tab bfe4_1 tab bfe4_2 tab bfe4_3keep d_bnd bfe4* ubigeo pothuan_mita border cusco near_x near_ysort ubigeosave   gis_dist, replace***--------------