Advertisement
Guest User

Untitled

a guest
Mar 26th, 2019
153
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
SAS 2.01 KB | None | 0 0
  1. DATA quad;
  2. INPUT x y id;
  3. CARDS;
  4. 0 0 1
  5. 0 2 1
  6. 2 2 1
  7. 2 0 1
  8. 2 2 2
  9. 3 3 2
  10. 4 2 2
  11. ; RUN;
  12.  
  13. DATA a;
  14. id=2; v='Triângulo'; OUTPUT;
  15. id=1; v='Quadrado'; OUTPUT; RUN;
  16.  
  17. PROC GMAP MAP=quad DATA=a;
  18. ID id;
  19. CHORO v;
  20. RUN; QUIT;
  21.  
  22. * http://geoftp.ibge.gov.br/organizacao_do_territorio/malhas_territoriais/malhas_municipais/municipio_2018/UFs/GO/ ;
  23. PROC MAPIMPORT DATAFILE='C:\Users\User\Desktop\go\52MUE250GC_SIR.shp' OUT=goias;
  24. RUN;
  25.  
  26. DATA b;
  27. segment=2; v=2;
  28. RUN;
  29.  
  30. PROC GMAP MAP=goias DATA=b all;
  31. ID segment;
  32. CHORO v;
  33. RUN; QUIT;
  34.  
  35. DATA goias; SET goias;
  36. codigo=input(CD_GEOCMU, 7.); RUN;
  37.  
  38. DATA b;
  39. codigo=2; segment=1; v=2; RUN;
  40.  
  41. PROC GMAP MAP=goias DATA=b all;
  42. ID codigo;
  43. CHORO v;
  44. RUN; QUIT;
  45.  
  46. DATA b;
  47. codigo=5201108; v=1; OUTPUT;
  48. codigo=5208707; v=2; OUTPUT;
  49. codigo=5217302; v=3; OUTPUT; RUN;
  50.  
  51. PROC GMAP MAP=goias DATA=b all;
  52. ID codigo;
  53. CHORO v;
  54. RUN; QUIT;
  55.  
  56. PROC GMAP MAP=goias DATA=b all;
  57. ID codigo;
  58. PRISM v;
  59. RUN; QUIT;
  60.  
  61. PROC GMAP MAP=goias DATA=b all;
  62. ID codigo;
  63. BLOCK v;
  64. RUN; QUIT;
  65.  
  66. PROC GMAP MAP=goias DATA=b all;
  67. ID codigo;
  68. SURFACE v;
  69. RUN; QUIT;
  70.  
  71. DATA goias; SET goias;
  72. pop=ranpoi(2, 200);
  73. RUN;
  74.  
  75. PROC SORT DATA=goias(keep=codigo) OUT=goias1 NODUPKEY; BY codigo; RUN;
  76.  
  77. DATA goias1; SET goias1;
  78. pop=ranpoi(2, 200);
  79. RUN;
  80.  
  81. PROC GMAP MAP=goias DATA=goias1 all;
  82. ID codigo;
  83. CHORO pop;
  84. RUN; QUIT;
  85.  
  86. PROC SQL NOPRINT;
  87. SELECT min(x) INTO:minx FROM goias;
  88. SELECT max(x) INTO:maxx FROM goias;
  89. SELECT min(y) INTO:miny FROM goias;
  90. SELECT max(y) INTO:maxy FROM goias;
  91. QUIT;
  92. %put &minx &maxx &miny &maxy;
  93.  
  94. DATA grid;
  95. DO x=&minx TO &maxx BY 0.1;
  96. DO y=&miny TO &maxy BY 0.1;
  97. OUTPUT;
  98. END;
  99. END;
  100. RUN;
  101.  
  102. DATA grid; SET grid;
  103. RETAIN hsys '3' xsys ysys "2";
  104. function='label'; text='K'; color='red'; style='special'; size=1; when='a'; *Google : font special sas;
  105. RUN;
  106.  
  107. PROC GMAP MAP=goias DATA=b all ANNO=grid;
  108. ID codigo;
  109. CHORO v;
  110. RUN; QUIT;
  111.  
  112. PROC GINSIDE DATA=grid MAP=goias OUT=grid1 INSIDEONLY;
  113. ID codigo; RUN;
  114.  
  115. PROC GMAP MAP=goias DATA=b all ANNO=grid1;
  116. ID codigo;
  117. CHORO v;
  118. RUN; QUIT;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement