Announcement

Collapse
No announcement yet.
X
  • Filter
  • Time
  • Show
Clear All
new posts

  • Drawing coordinates on a map

    Hello Stata community, I have this kind of data to work with:

    Code:
    * Example generated by -dataex-. To install: ssc install dataex
    clear
    input int stopid str21 stopname float(stoplat stoplon)
      1 "Ain Drahem"             36.779 8.68833
      2 " Ain El Karma"         36.1752 8.39293
      3 " Ain Essobeh"          36.9579 8.91899
      4 "Aouana"                37.0302 9.33377
      5 "Babouch"               36.8021 8.67753
      6 " Bareket Essahel"      36.4041   10.56
      7 "Beja"                  36.7266 9.19202
      8 " Ben Aoun"             34.8561 9.14682
      9 "Ben Guerdene"          33.1555 11.2131
     10 "Bir Ali Ben Khelifa"   34.7347 10.0944
     11 "Bir Bouragba"          36.4318 10.5712
     12 "Bir El Hafey"          34.9325 9.19626
     13 "Bir El Mchergua"       36.5159 9.96051
     14 "Bizerte"               37.2695 9.87057
     15 "Borj El Aifa"          36.2186 8.85217
     16 "Bouarada"              36.3511 9.62311
     17 "Bouargoub"             36.5378 10.5472
     18 "Bouficha"               36.299 10.4603
     19 "Bouhajla"              35.3948  10.049
     20 "Boussalem"             36.6097 8.97289
     21 "Bouterfess"            36.9542 8.85778
     22 "Carrefour Haffouz"     35.5961 9.85787
     23 "Carrfour Sebitla"      35.1471 9.36316
     24 "Caserne Bouficha"      36.3279 10.5007
     25 "Dahmani"               35.9463 8.83067
     26 "Dhehiba"               32.0107 10.7001
     27 "Douz"                  33.4579 9.02284
     28 "El Ala"                35.6158 9.56252
     29 "El Batten"             35.7108 10.0007
     30 "El Faouar"             33.3562 8.67671
     31 "el hajeb"              35.3953 9.54222
     32 "El Houareb"            35.5675 9.76938
     33 "El Krib"               36.3256  9.1343
     34 "El Ksour"              35.8963 8.88281
     35 "El Oueslatia"          35.8485 9.59002
     36 "Enfidha"               36.1358 10.3814
     37 "Ennaima"               36.1099 8.60272
     38 "Errhima"               35.4784 9.66112
     39 "Feriana"               34.9475 8.56884
     40 "Fernenah"              36.6566 8.69617
     41 "Fondouk Jedid"         36.6596 10.4547
     42 "Fouchana"              36.6991 10.1696
     43 "Gabes"                 33.8823 10.0862
     44 "Gafsa"                 34.4164 8.78973
     49 "Ghomracen"               33.06 10.3396
     51 "Grombalia"             36.6016 10.4965
     53 "Haffouz"               35.6353 9.67362
     54 "Hajeb El Ayoun"        35.3959 9.54446
     55 "Hammam Lif"            36.7268 10.3372
     56 "Hazoua"                33.7345 7.59004
     57 "JERBA (houmet essouk)" 33.8721 10.8579
     58 "Jefna"                 37.0514 9.47382
     59 "Jelma"                 35.2719 9.41895
     60 "Jendouba"               36.495 8.77573
     61 "Jerissa"               35.8447   8.632
     62 "Jradou"                36.2777 8.44952
     63 "Kairouan"              35.6809 10.0866
     64 "Kalaat Senen"          35.7689 8.34217
     65 "Kasserine"             35.1776 8.84774
     66 "Kebili"                33.7261 8.97982
     67 "Kondar"                35.9352 10.3007
     68 "Koutine"               33.4392 10.3904
     69 "Ksar Lamssa"           36.0331 9.69247
     70 "Le Kef"                36.1629  8.7032
     71 "Lassouda"                 35.1 9.53519
     72 "Mateur"                37.0364 9.66717
     73 "Medenine"              33.3612 10.4694
     74 "Menzel Bourguiba"      37.1558 9.78861
     75 "Mareth"                 33.615 10.2903
     76 "Metbassta"             35.7937 10.1463
     77 "Metlaoui"              34.3159 8.40337
     78 "Midoun"                33.8052 10.9905
     79 "Medjez El Bab"         36.6399 9.61333
     80 "Mhamdia"               36.6738 10.1556
     81 "Nebeur"                36.2948  8.7669
     82 "Nefta"                   33.87 7.88268
     83 "Nefza"                 36.9755 9.08106
     84 "Ouechteta"             36.9616 9.00362
     85 "Oued Eremal"           36.1143 8.62922
     86 "Oued Ezzarka"          36.6645 9.44053
     87 "Oued Mellag"           36.1193 8.47501
     88 "Oum Labouab"            36.169 9.77139
     89 "Pont de Malleg"        36.4066 8.74816
     90 "Sejnene"                37.057 9.23965
     91 "Sakiet Sidi Youssef"   36.2219 8.35467
     92 "Sfax"                  34.7234 10.7578
     93 "Sidi Aouidet"          36.2217 9.74434
     94 "Sidi Bouzid"           35.0436 9.49667
     95 "Sidi Khiar"            36.3546 8.79179
     96 "Siliana"               36.0942 9.36984
     97 "Souk Lahad"            33.7786 8.84698
     98 "Sousse"                35.8141 10.6304
     99 "Tabarka"               36.9556 8.75575
    100 "Tajerouine"            35.8915 8.55263
    101 "Tamaghza"              34.3837 7.93326
    102 "Tamra"                 37.0587  9.1217
    103 "Tataouine"             32.9479 10.4572
    104 "Teboursouk"            36.4576 9.24982
    105 "Testour"               36.5516 9.43805
    106 "Thala"                 35.5727 8.67023
    end
    As you can see, those are some geographical coordinates of points. What I want to do is map those points on the Tunisia map using spmap package and so on. In other words, I want to have a Tunisia map that shows those points on it. I would also prefer to have the Tunisia map with governorate boundries, like, the administrative division by boundries.

    Any help please?

  • #2
    Use the -point()- option within spmap.

    Comment


    • #3
      Scott Merryman Thanks for the help, that worked. Still, I just had to not to use the labels on those points on the map because, with more than 200 points, it is going to be crowded as a map and it will be hard to even read the names of those bus stops.

      Thanks for the help!

      Comment

      Working...
      X