Announcement

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

  • Drawing bus stops on a country's map with points

    Hello Stata users;

    I have two sorts of data available to me. The first one is of a country's governorates coordinates (basically, I draw the map using the spmap command):
    Code:
    * Example generated by -dataex-. To install: ssc install dataex
    clear
    input float(Shape_Leng Shape_Area) str12 ADM2_EN str19 ADM2_AR str4 ADM2_PCODE str9 ADM2_REF str25 ADM1_AR str3 ADM1_PCODE str7 ADM0_EN str8 ADM0_AR byte idgouv
    1.6094854 .02651095 "Tunis"        "تونس"            "TN11" ""          "الشمال الشرقي" "TN1" "Tunisia" "تونس"  1
     1.665691 .04211232 "Ariana"       "أريانة"        "TN12" ""          "الشمال الشرقي" "TN1" "Tunisia" "تونس"  2
    1.7179574 .06661047 "Ben Arous"    "بن عروس"       "TN13" ""          "الشمال الشرقي" "TN1" "Tunisia" "تونس"  3
    2.2618155  .1153453 "Manubah"      "منوبة"          "TN14" ""          "الشمال الشرقي" "TN1" "Tunisia" "تونس"  4
      3.30774  .2840685 "Nabeul"       "نابل"            "TN15" ""          "الشمال الشرقي" "TN1" "Tunisia" "تونس"  5
    3.2403774 .28479052 "Zaghouan"     "زغوان"          "TN16" ""          "الشمال الشرقي" "TN1" "Tunisia" "تونس"  6
     5.343521  .3661806 "Bizerte"      "بنزرت"          "TN17" ""          "الشمال الشرقي" "TN1" "Tunisia" "تونس"  7
    4.3867736  .3681952 "Béja"        "باجة"            "TN21" "Beja"      "الشمال الغربي" "TN2" "Tunisia" "تونس"  8
     3.647018  .3133408 "Jendouba"     "جندوبة"        "TN22" ""          "الشمال الغربي" "TN2" "Tunisia" "تونس"  9
    4.1811414 .51084507 "Le Kef"       "الكاف"          "TN23" ""          "الشمال الغربي" "TN2" "Tunisia" "تونس" 10
     5.552421  .4684452 "Siliana"      "سليانـــة"  "TN24" ""          "الشمال الغربي" "TN2" "Tunisia" "تونس" 11
     3.156125  .2278667 "Sousse"       "سوسة"            "TN31" ""          "الوسط الشرقي"   "TN3" "Tunisia" "تونس" 12
    2.2377114 .10218735 "Monastir"     "المنستير"    "TN32" ""          "الوسط الشرقي"   "TN3" "Tunisia" "تونس" 13
     4.073021 .28655812 "Mahdia"       "المهدية"      "TN33" ""          "الوسط الشرقي"   "TN3" "Tunisia" "تونس" 14
     7.548694  .6909879 "Sfax"         "صفاقس"          "TN34" ""          "الوسط الشرقي"   "TN3" "Tunisia" "تونس" 15
     4.946872  .6479417 "Kairouan"     "القيروان"    "TN41" ""          "الوسط الغربي"   "TN4" "Tunisia" "تونس" 16
     4.989357  .8122157 "Kassérine"   "القصرين"      "TN42" "Kasserine" "الوسط الغربي"   "TN4" "Tunisia" "تونس" 17
     5.466782  .7327916 "Sidi Bou Zid" "سيدي بوزيد" "TN43" ""          "الوسط الغربي"   "TN4" "Tunisia" "تونس" 18
    4.5108953  .7290807 "Gabès"       "قابس"            "TN51" "Gabes"     "الجنوب الشرقي" "TN5" "Tunisia" "تونس" 19
     9.757171  .8966198 "Médenine"    "مدنين"          "TN52" "Medenine"  "الجنوب الشرقي" "TN5" "Tunisia" "تونس" 20
    11.073194 3.6204395 "Tataouine"    "تطاوين"        "TN53" ""          "الجنوب الشرقي" "TN5" "Tunisia" "تونس" 21
     4.548769  .7402925 "Gafsa"        "قفصة"            "TN61" ""          "الجنوب الغربي" "TN6" "Tunisia" "تونس" 22
     3.798989  .5907553 "Tozeur"       "توزر"            "TN62" ""          "الجنوب الغربي" "TN6" "Tunisia" "تونس" 23
     7.504424 2.1741836 "Kebili"       "قبلي"            "TN63" ""          "الجنوب الغربي" "TN6" "Tunisia" "تونس" 24
    end
    This first data is about a country's map (I've used shapefiles to get this. It is possible to get the map of the country just by using spmap).

    The sedond database that I have is about the coordinates of bus stops:
    Code:
    * Example generated by -dataex-. To install: ssc install dataex
    clear
    input int stopid str24 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 "Babbouche"               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 "Alaa"                    35.6158 9.56252
     29 "El Batten"               35.7108 10.0007
     30 "Fawar"                   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 "Ghomrassen"                33.06 10.3396
     51 "Grombalia"               36.6016 10.4965
     52 "Gare routiere Tunis Sud" 36.7874 10.1798
     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 Houmt Souk"        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 Senan"            35.7689 8.34217
     65 "Gasserine"               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 "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 Ezzarga"            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
    end
    What I wanna do is get those bus stops coordinates on the map of the country. Since I have the coordinates of each bus stop, I wanna show them on the map with point (a point represent a bus stop)
    Should I already merge the two datasets? There is no common variable between the two that I could use for merging already. Yet, I do think it is possible to do so.

    Any help, please?

  • #2
    Any help please ??

    Comment

    Working...
    X