Announcement

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

  • Tornado Map - spmap with line over USA state shapefile

    Goodmorning everyone,

  • #2
    I am willing to depict the path of tornado over a USA map. My data look like this:

    Code:
    * Example generated by -dataex-. For more info, type help dataex
    clear
    input long id int yr byte mag double(lat lon len) int wid byte expand_data float order double(lat_ip lon_ip _CX _CY) str2 STATEFP byte state str8 STATENS str11 AFFGEOID str2 STUSPS str44 state_name long _ID byte _merge
    56457 2012  2       .                   . 1.69   50 1 17  41.18641111111111 -100.86366666666666 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    19770 1979  0   41.05             -103.88  5.1   20 0  1              41.05             -103.88 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    40570 2000  1   41.63               -99.4   12  250 0  1              41.63               -99.4 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    47469 2005  0       .                   .  1.5   45 1  6  40.75555555555555              -98.85 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    65888 2020 -9       .                   .  .82   30 1  3             40.091 -101.86796666666666 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    49978 2007  0   42.78              -97.51  3.5  100 0  1              42.78              -97.51 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    38252 1998  0       .                   .  2.5   30 1 17  40.16555555555556  -97.49777777777777 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    16634 1975  0       .                   .   .2   10 1  9 41.120044444444446  -98.86995555555556 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    36144 1996  0       .                   .   .1   20 1 15              41.13              -96.97 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    49965 2007  1       .                   . 18.3  400 1 17 40.592222222222226  -97.09333333333332 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    11456 1969  1       .                   .   .1   10 1  3  41.40001111111111   -97.2499888888889 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     5357 1960  1       .                   .   .1   10 1 10           40.93005           -98.52995 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    24363 1984  1       .                   .    1   30 1  9  40.97004444444445 -101.06995555555555 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    35880 1996  0       .                   .   .1   20 1 18               41.9              -101.7 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    33629 1994  0       .                   .   .1   23 1  5  40.77002222222223  -97.76997777777778 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    16446 1975  0       .                   .    3   80 1  7 42.430033333333334  -96.74996666666667 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    25053 1985  1       .                   .   .5   30 1 18  40.23009444444445 -100.89990555555556 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    28986 1990  1       .                   .    2   50 1 13 41.480066666666666  -98.06993333333332 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     6577 1962  0       .                   .   .1   10 1  5  41.07002222222223  -97.76997777777778 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    10882 1968  0       .                   .   .1   10 1 11  41.40005555555556 -102.97994444444446 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    40353 1999  0       .                   .    1  100 1 14  41.40888888888889              -100.2 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    48985 2006  1       .                   .   .1   10 1 10               42.7             -102.45 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    33636 1994  0       .                   .   .1   23 1 12  40.53006111111112  -97.57993888888889 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    45868 2004  0       .                   .   .2   10 1  6               41.3             -100.05 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    63429 2018  0       .                   . 1.01   75 1  4           41.19795 -103.65298333333334 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    25818 1986  1       .                   .    2   50 1 12 40.492222222222225  -98.26166666666667 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    51668 2008  1       .                   . 3.79  440 1  4              41.07  -97.22833333333334 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     5343 1960  2       .                   .   .1   10 1 11  41.40005555555556 -102.34994444444445 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     5196 1960  1       .                   .  1.5   10 1 13  40.43006666666667  -99.19993333333333 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    50326 2007  0   41.97              -98.55   .1   10 0 19              41.97              -98.55 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    11943 1970  0       .                   .   .1   10 1 17   42.4500888888889  -98.47991111111112 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    45702 2004  0       .                   .    1   50 1  5              40.47  -97.39555555555556 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     5224 1960  1       .                   .  2.5   10 1  7 40.430033333333334  -99.34996666666666 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    31057 1992  0       .                   .   21   23 1 17 41.163333333333334  -99.96000000000001 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    24946 1985  1       .                   .   .5   50 1  4 41.870016666666665 -103.66998333333333 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    44667 2003  4    42.5              -97.32   13 1200 0  1               42.5              -97.32 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    18221 1977  0 41.9001  -98.09989999999999   .2   30 0 19            41.9001  -98.09989999999999 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     3433 1957  1       .                   .  6.2   10 1  3  40.23444444444444  -96.30666666666666 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    25388 1985  1       .                   .   .2   50 1 16  40.95008333333334  -99.26991666666666 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    30065 1991  0       .                   .   .5   30 1  6  40.47002777777778  -99.89997222222223 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    26001 1986  4       .                   . 52.5   73 1  8  42.48555555555556   -96.6588888888889 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     1031 1953  2       .                   .   31   10 1 16  42.37166666666667               -97.2 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    15409 1974  2       .                   .    1  233 1 10           40.03005  -98.09994999999999 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    40353 1999  0       .                   .    1  100 1  6 41.391111111111115              -100.2 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     3547 1957  1       .                   .   .1   10 1 12  41.38006111111112  -97.49993888888889 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    12673 1971  0       .                   .   .1   10 1 14  40.75007222222223 -100.74992777777778 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    24634 1984  1       .                   .    1  100 1 11 41.630055555555565  -98.21994444444445 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    25354 1985  0       .                   .   .1   40 1 14  40.77007222222223  -98.34992777777778 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    34927 1995  1       .                   .   .8  250 1 12  40.42006111111112   -96.3699388888889 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    18182 1977  2       .                   .  1.9   70 1  7              41.79 -103.48333333333333 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    28079 1989  0       .                   .   .1   30 1 16 40.320083333333336  -97.71991666666666 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    56652 2012  0       .                   .  .47   20 1 17  41.08594444444444            -100.798 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    48064 2005  0       .                   .   .4   50 1 13              41.67              -98.93 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    29054 1990  0       .                   .   .5   30 1 17   40.2800888888889  -98.81991111111111 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    36434 1996  0       .                   .   .2   50 1  6  40.26444444444444              -100.4 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    44960 2003  0       .                   .   .2   25 1 11              41.25             -103.12 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     6207 1962  2       .                   .    2   10 1 17 41.720088888888895  -98.77991111111112 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     7074 1963  0       .                   .   .1   10 1 17   42.9200888888889 -101.44991111111112 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    40546 2000  0       .                   .   .5   10 1 14  40.77722222222222  -98.11277777777778 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    31904 1992  0       .                   .   .1   10 1 10 41.020050000000005           -96.04995 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     4609 1959  1       .                   .   14   10 1  6 40.995000000000005  -99.56055555555555 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    16302 1975  0       .                   .   .2   10 1  4  42.45001666666667 -101.09998333333333 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    16302 1975  0       .                   .   .2   10 1  5  42.45002222222223 -101.09997777777778 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    33680 1994  1       .                   .   .1   23 1 10           41.93005           -96.74995 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    39999 1999  0       .                   .   .1   30 1 13              42.68             -103.87 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    10764 1968  1       .                   .   .1   10 1  5  42.45002222222223  -98.64997777777779 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    39539 1999  0       .                   .   .3   20 1  3              41.18              -98.47 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     1029 1953  1   40.32               -98.9   15   33 0 19              40.32               -98.9 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     7646 1964  0       .                   .   .1   10 1  2  40.90000555555556  -97.09999444444445 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    40217 1999  1       .                   .    2  200 1  8  42.90777777777778 -102.05444444444443 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    38726 1998  1       .                   .    2  500 1  5              41.45              -97.22 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    20752 1980  1       .                   .   .8  150 1  8   40.9500388888889  -98.36996111111112 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    16302 1975  0       .                   .   .2   10 1 11 42.450055555555565 -101.09994444444445 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    16298 1975  0 41.9801 -101.76989999999999   .2   10 0 19            41.9801 -101.76989999999999 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    47735 2005  0       .                   .   .1   25 1 18              41.88             -101.25 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    28102 1989  0       .                   .   .1   30 1  4  41.63001666666667  -99.21998333333333 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    31291 1992  0       .                   .   .1   73 1 14  40.92007222222223  -96.51992777777778 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     4509 1959  1       .                   .   .1   10 1  9  41.08004444444445  -97.77995555555556 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     9925 1967  2       .                   . 34.3   50 1 17  40.78555555555555              -98.49 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    47368 2005  1       .                   .    2   50 1 18  40.64888888888889 -100.59888888888888 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    43785 2003  0       .                   .   .1   40 1 13               40.8              -96.23 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
      166 1950  4       .                   .   10  440 1 14  41.74944444444444              -96.47 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    51508 2008  0       .                   .  .22   40 1  7              40.96              -96.44 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     3999 1958  1       .                   .  1.7   10 1  5 40.580022222222226  -96.54997777777778 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    41259 2000  0       .                   .   .5   20 1 16              40.65              -99.37 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    25156 1985  1       .                   .   .2   40 1 12 40.900061111111114   -97.4499388888889 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     4047 1958  2       .                   .    6  533 1  3 40.420011111111116   -99.3699888888889 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    16302 1975  0       .                   .   .2   10 1 17   42.4500888888889 -101.09991111111111 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    30188 1991  0       .                   .   .2   23 1  8 40.720038888888894 -101.71996111111112 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    34825 1995  0       .                   .   .1   10 1  9  42.28004444444445  -96.54995555555556 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     2475 1956  2       .                   . 12.3   10 1  9  40.64666666666666  -98.44222222222223 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     3585 1957  2       .                   . 34.2   10 1 13              40.13 -101.54666666666667 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    45680 2004  0       .                   .   .1   50 1 13               40.1              -99.83 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    43895 2003  0       .                   .    2   20 1 15  41.17777777777778  -99.50666666666666 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     5341 1960  0       .                   .   .1   10 1  7  41.23003333333333 -103.47996666666667 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    31665 1992  2       .                   .   .4   50 1  2  40.75000555555556 -100.26999444444445 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    42758 2001  0       .                   .    3  150 1  5  40.56555555555555  -97.97666666666667 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    27693 1989  1 40.2201            -99.6799   .2   70 0 19            40.2201            -99.6799 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
     2719 1956  1       .                   .  5.4   10 1 10              42.65             -99.875 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    37233 1997  1       .                   .    4  100 1 13 40.559999999999995  -96.67666666666668 -99.81080000561343 41.52714192771041 "31" 31 "01779792" "0400000US31" "NE" "Nebraska" 1 3
    end
    label values expand_data _merge
    label values _merge _merge
    label def _merge 0 "Original observation", modify
    label def _merge 1 "Duplicated observation", modify
    where lat_ip and lon_ip are ipolation of points between start and end point of each tornado path. I used geoinpoly in order to reference each point with a usa state shapefile (usaST_shp).
    However when I try to draw a map of all the path using:

    spmap using usaST_shp, id(_ID) osize(0.05 ..) point(data(tornado_path.dta) x(lon_ip) y(lat_ip))

    i got the following errore message:
    Variable _ID specified in option id() does not uniquely identify the observations

    Can you suggest me another way to depict tornado path over a USA state map?
    Thank you!!

    Comment


    • #3
      In another thread, your data seemed more suitable for this task. To depict a path you'll need starting and end lat/lon. Here's an example using your other data example


      Code:
      * Example generated by -dataex-. For more info, type help dataex
      clear
      input long om int yr byte(mo dy) str10 date str8 time byte tz str2 st byte stf int stn byte mag int(inj fat) double(loss closs slat slon elat elon len) int wid byte fc long _ID
      192 1950 10  1 "1950-10-01" "21:00:00" 3 "OK" 40 23 1  0  0 4 0 36.73 -102.52              36.88              -102.3 15.8   10 0   1
      193 1950 10  9 "1950-10-09" "02:15:00" 3 "NC" 37  9 3  3  0 5 0 34.17   -78.6 34.170100000000005  -78.59989999999999    2  880 0   2
      195 1950 11 20 "1950-11-20" "02:20:00" 3 "KY" 21  1 2  0  0 5 0 37.37   -87.2            37.3701            -87.1999   .1   10 0   3
      196 1950 11 20 "1950-11-20" "04:00:00" 3 "KY" 21  2 1  0  0 5 0  38.2   -84.5 38.200100000000006            -84.4999   .1   10 0   4
      197 1950 11 20 "1950-11-20" "07:30:00" 3 "MS" 28 14 1  3  0 4 0 32.42  -89.13 32.420100000000005  -89.12989999999999    2   37 0   5
      194 1950 11  4 "1950-11-04" "17:00:00" 3 "PA" 42  5 3  1  0 5 0  40.2  -76.12               40.4              -75.93 15.9  100 0   6
      198 1950 12  2 "1950-12-02" "15:00:00" 3 "IL" 17  7 2  3  1 4 0 38.97  -90.05              39.07              -89.72 18.8   50 0   7
      199 1950 12  2 "1950-12-02" "16:00:00" 3 "IL" 17  8 3 25  2 6 0 38.75  -89.67               38.9              -89.38   18  200 0   8
      200 1950 12  2 "1950-12-02" "16:25:00" 3 "AR"  5 12 3  0  0 1 0 36.12  -91.83              36.18              -91.72  7.8   10 0   9
      201 1950 12  2 "1950-12-02" "17:30:00" 3 "IL" 17  9 1  0  0 4 0 38.17  -89.78              38.22              -89.62  9.6   50 0  10
        4 1950  1 13 "1950-01-13" "05:25:00" 3 "AR"  5  1 3  1  1 3 0  34.4  -94.37            34.4001            -94.3699   .6   17 0  11
        5 1950  1 25 "1950-01-25" "19:30:00" 3 "MO" 29  2 2  5  0 5 0  37.6  -90.68              37.63              -90.65  2.3  300 0  12
        6 1950  1 25 "1950-01-25" "21:00:00" 3 "IL" 17  3 2  0  0 5 0 41.17  -87.33 41.170100000000005            -87.3299   .1  100 0  13
        7 1950  1 26 "1950-01-26" "18:00:00" 3 "TX" 48  1 2  2  0 0 0 26.88  -98.12              26.88              -98.05  4.7  133 0  14
        1 1950  1  3 "1950-01-03" "11:00:00" 3 "MO" 29  1 3  3  0 6 0 38.77  -90.22              38.83              -90.03  9.5  150 0  15
        2 1950  1  3 "1950-01-03" "11:55:00" 3 "IL" 17  2 3  3  0 5 0  39.1   -89.3              39.12              -89.23  3.6  130 0  16
        3 1950  1  3 "1950-01-03" "16:00:00" 3 "OH" 39  1 1  1  0 4 0 40.88  -84.58 40.880100000000006            -84.5799   .1   10 0  17
        8 1950  2 11 "1950-02-11" "13:10:00" 3 "TX" 48  2 2  0  0 4 0 29.42  -95.25              29.52              -95.13  9.9  400 0  18
        9 1950  2 11 "1950-02-11" "13:50:00" 3 "TX" 48  3 3 12  1 4 0 29.67  -95.05              29.83                 -95   12 1000 0  19
       10 1950  2 11 "1950-02-11" "21:00:00" 3 "TX" 48  4 2  5  0 5 0 32.35   -95.2              32.42               -95.2  4.6  100 0  20
       11 1950  2 11 "1950-02-11" "23:55:00" 3 "TX" 48  5 2  6  0 5 0 32.98  -94.63                 33               -94.7  4.5   67 0  21
       12 1950  2 12 "1950-02-12" "00:30:00" 3 "TX" 48  6 2  8  1 4 0 33.33  -94.42              33.45              -94.42    8  833 0  22
       13 1950  2 12 "1950-02-12" "01:15:00" 3 "TX" 48  7 1  0  0 4 0 32.08  -98.35               32.1              -98.33  2.3  233 0  23
       14 1950  2 12 "1950-02-12" "06:10:00" 3 "TX" 48  8 2  0  0 4 0 31.52  -96.55              31.57              -96.55  3.4  100 0  24
       15 1950  2 12 "1950-02-12" "11:57:00" 3 "TX" 48  9 1 32  0 5 0  31.8   -94.2              31.88              -94.12  7.7  100 0  25
       16 1950  2 12 "1950-02-12" "12:00:00" 3 "MS" 28  1 2  2  3 4 0  34.6  -89.12 34.600100000000005            -89.1199   .1   10 0  26
       17 1950  2 12 "1950-02-12" "12:00:00" 3 "MS" 28  2 1  0  0 0 0  34.6  -89.12 34.600100000000005            -89.1199    2   10 0  27
       18 1950  2 12 "1950-02-12" "12:00:00" 3 "TX" 48 10 3 15  3 5 0  31.8   -94.2               31.8              -94.18  1.9   50 0  28
       19 1950  2 12 "1950-02-12" "12:30:00" 3 "AR"  5  2 2  0  0 3 0 34.48   -92.4            34.4801            -92.3999   .1  100 0  29
       20 1950  2 12 "1950-02-12" "13:00:00" 3 "LA" 22  1 4 77 18 5 0 31.97     -94                 33               -93.3 82.6  100 0  30
       21 1950  2 12 "1950-02-12" "13:20:00" 3 "LA" 22  2 2 10  5 5 0  32.2  -93.58              32.97              -93.17 58.4  100 0  31
       22 1950  2 12 "1950-02-12" "14:00:00" 3 "LA" 22  4 3 25  5 5 0 31.63  -93.65              32.55              -93.03 74.5  100 0  32
       23 1950  2 12 "1950-02-12" "15:00:00" 3 "AR"  5  3 2  0  0 4 0 33.27  -92.95              33.35              -92.95  5.7  100 0  33
       24 1950  2 12 "1950-02-12" "23:00:00" 3 "LA" 22  5 1  0  0 4 0  32.6  -91.33 32.600100000000005            -91.3299   .5   33 0  34
       25 1950  2 13 "1950-02-13" "01:00:00" 3 "TN" 47  1 1  8  0 3 0 35.35  -89.77 35.350100000000005  -89.76989999999999   .2   10 0  35
       26 1950  2 13 "1950-02-13" "02:00:00" 3 "TN" 47  2 2  1  9 4 0 35.75  -89.48            35.7501            -89.4799   .2   10 0  36
       27 1950  2 27 "1950-02-27" "10:20:00" 3 "OK" 40  1 2  0  0 4 0 35.55   -97.6            35.5501  -97.59989999999999    2   50 0  37
       29 1950  3 16 "1950-03-16" "09:15:00" 3 "FL" 12  1 2  0  0 3 0 29.65  -81.22            29.6501            -81.2199  1.5  150 0  38
       30 1950  3 19 "1950-03-19" "07:30:00" 3 "LA" 22  6 1  2  0 4 0 30.45  -93.45            30.4501            -93.4499    2   33 0  39
       31 1950  3 19 "1950-03-19" "13:15:00" 3 "LA" 22  7 2  0  0 4 0  30.1     -91            30.1001            -90.9999    1   50 0  40
       32 1950  3 19 "1950-03-19" "13:15:00" 3 "LA" 22  8 0  0  0 4 0  29.7   -90.1              29.67               -89.8 18.1   27 0  41
       33 1950  3 26 "1950-03-26" "19:30:00" 3 "AR"  5  4 2  3  0 4 0 34.12  -93.07              34.32              -92.88 17.4  150 0  42
       34 1950  3 26 "1950-03-26" "19:31:00" 3 "AR"  5  5 3  1  0 5 0 36.15  -91.83               36.2              -91.75  5.7  200 0  43
       35 1950  3 26 "1950-03-26" "20:30:00" 3 "AR"  5  6 2  7  0 5 0  34.7  -92.35               34.8              -92.22 10.4  600 0  44
       36 1950  3 26 "1950-03-26" "21:15:00" 3 "AR"  5  7 3 20  0 5 0 34.98  -91.73              35.08               -91.5 14.9 1760 0  45
       37 1950  3 26 "1950-03-26" "23:30:00" 3 "AR"  5  8 2  2  0 0 0  35.1   -91.4              35.15              -91.33  5.4  833 0  46
       38 1950  3 27 "1950-03-27" "03:00:00" 3 "OK" 40  2 2  0  0 3 0 34.85  -95.75 34.850100000000005            -95.7499   .1   77 0  47
       39 1950  3 27 "1950-03-27" "05:00:00" 3 "MS" 28  6 2  2  1 0 0 33.17  -90.55 33.170100000000005            -90.5499   .1   10 0  48
       40 1950  3 27 "1950-03-27" "06:00:00" 3 "LA" 22  9 1  2  0 4 0 32.63  -91.28              32.73              -91.05 15.3   17 0  49
       41 1950  3 27 "1950-03-27" "07:30:00" 3 "MS" 28  7 2  0  0 0 0 33.13  -89.57              33.25              -89.42 11.9   10 0  50
       42 1950  3 27 "1950-03-27" "07:45:00" 3 "MS" 28  4 2  7  0 5 0 32.28  -90.27              32.98              -89.82 59.3   50 0  51
       43 1950  3 27 "1950-03-27" "07:45:00" 3 "MS" 28  5 2  6  0 5 0 32.28  -90.27              32.48              -89.78 31.9   50 0  52
       44 1950  3 27 "1950-03-27" "08:30:00" 3 "LA" 22 10 1  0  0 4 0 30.92  -91.13            30.9201  -91.12989999999999   .5   27 0  53
       45 1950  3 27 "1950-03-27" "08:45:00" 3 "MS" 28  9 2  0  0 0 0 33.67  -88.78              33.72              -88.73  4.7   10 0  54
       46 1950  3 27 "1950-03-27" "11:00:00" 3 "MS" 28  8 1  2  0 4 0 34.73  -88.27            34.7301  -88.26989999999999   .1   10 0  55
       47 1950  3 27 "1950-03-27" "15:00:00" 3 "TN" 47  3 1  0  0 3 0 35.68  -85.77            35.6801  -85.76989999999999   .2   10 0  56
       48 1950  3 27 "1950-03-27" "16:30:00" 3 "IL" 17  4 2  0  0 3 0 40.57   -88.6               40.6              -88.55    3   50 0  57
       28 1950  3  1 "1950-03-01" "02:30:00" 3 "MS" 28  3 1  0  0 0 0  32.5  -88.85            32.5001  -88.84989999999999   .1   10 0  58
       53 1950  4 18 "1950-04-18" "01:30:00" 3 "AL"  1  1 3 15  0 4 0 30.67   -88.2              30.85               -88.1   14  100 0  59
       54 1950  4 18 "1950-04-18" "01:45:00" 3 "AL"  1  2 2  0  0 3 0  30.7  -87.92            30.7001            -87.9199    2  150 0  60
       55 1950  4 24 "1950-04-24" "07:45:00" 3 "LA" 22 11 1  0  0 4 0  30.1  -90.03            30.1001            -90.0299    2  100 0  61
       56 1950  4 27 "1950-04-27" "16:00:00" 3 "SC" 45  1 1  0  0 4 0 33.07     -80 33.070100000000004            -79.9999   .7   37 0  62
       57 1950  4 28 "1950-04-28" "14:17:00" 3 "OK" 40  5 3  1  1 5 0 34.88  -99.28              35.17               -99.2 20.8  400 0  63
       58 1950  4 28 "1950-04-28" "17:00:00" 3 "KS" 20  1 1  0  0 4 0 38.12  -95.85              38.35              -95.63 19.8   10 0  64
       59 1950  4 28 "1950-04-28" "18:00:00" 3 "TX" 48 11 4  5  5 5 0 32.42   -99.5              32.42              -99.48  1.3  233 0  65
       60 1950  4 28 "1950-04-28" "19:05:00" 3 "OK" 40  6 4 32  5 5 0 35.08   -96.4              35.13              -96.35  4.5  200 0  66
       61 1950  4 28 "1950-04-28" "19:30:00" 3 "OK" 40  7 2  0  0 4 0 34.55   -96.2            34.5501            -96.1999   .8  100 0  67
       62 1950  4 29 "1950-04-29" "12:00:00" 3 "MS" 28 10 2  0  0 5 0 31.27  -89.83              31.38               -89.7 11.3  100 0  68
       63 1950  4 29 "1950-04-29" "15:30:00" 3 "TX" 48 12 1  0  0 4 0  31.9   -98.6              31.73               -98.6 11.5  200 0  69
       49 1950  4  2 "1950-04-02" "15:30:00" 3 "AR"  5  9 3  8  1 3 0 35.75  -91.22            35.7501            -91.2199  2.5   10 0  70
       50 1950  4  2 "1950-04-02" "18:00:00" 3 "OK" 40  3 2  0  0 4 0 35.82  -97.02 35.820100000000004  -97.01989999999999    1  100 0  71
       51 1950  4  2 "1950-04-02" "18:00:00" 3 "AR"  5 10 1  0  0 3 0 35.98  -89.92            35.9801            -89.9199   .8  100 0  72
       52 1950  4  3 "1950-04-03" "16:15:00" 3 "OK" 40  4 1  0  0 3 0 36.13  -95.83 36.130100000000006            -95.8299    1   33 0  73
       95 1950  5 10 "1950-05-10" "13:00:00" 3 "WY" 56  1 2  0  0 1 0 41.18 -110.42            41.1801           -110.4199    2   33 0  74
       96 1950  5 11 "1950-05-11" "18:30:00" 3 "OK" 40 15 1  0  0 2 0 36.82 -100.52 36.820100000000004 -100.51989999999999   .5   77 0  75
       97 1950  5 12 "1950-05-12" "14:00:00" 3 "NC" 37  1 1  0  0 4 0 35.75  -77.42            35.7501            -77.4199    2   10 0  76
       98 1950  5 12 "1950-05-12" "17:00:00" 3 "NC" 37  2 1  0  0 0 0 35.73   -78.4            35.7301            -78.3999    1  200 0  77
       99 1950  5 12 "1950-05-12" "17:00:00" 3 "NC" 37  3 0  0  0 0 0 35.73   -78.4            35.7301            -78.3999    1  200 0  78
      100 1950  5 13 "1950-05-13" "11:30:00" 3 "LA" 22 19 1  0  0 4 0 31.15  -91.95            31.1501            -91.9499    2   27 0  79
      101 1950  5 14 "1950-05-14" "16:00:00" 3 "TX" 48 14 1  0  0 3 0 29.53  -95.45            29.5301            -95.4499    1  100 0  80
      102 1950  5 14 "1950-05-14" "17:00:00" 3 "NC" 37  4 1  0  0 3 0  35.3   -79.8              35.33              -79.75  3.8   50 0  81
      103 1950  5 14 "1950-05-14" "17:30:00" 3 "NC" 37  5 2  5  0 0 0    35  -80.68            35.0001            -80.6799    2   10 0  82
      104 1950  5 15 "1950-05-15" "00:01:00" 3 "TX" 48 15 2  0  0 4 0 29.78  -98.83              29.67              -98.57   18   10 0  83
      105 1950  5 15 "1950-05-15" "11:00:00" 3 "FL" 12  3 1  0  0 4 0 28.58  -81.37 28.580099999999998            -81.3699   .1   10 0  84
      106 1950  5 15 "1950-05-15" "11:00:00" 3 "FL" 12  4 2  0  0 4 0  28.5  -81.37            28.5001            -81.3699   .1   10 0  85
      107 1950  5 16 "1950-05-16" "18:00:00" 3 "KS" 20 11 2  0  0 3 0 37.55  -98.42            37.5501            -98.4199   .2  200 0  86
      108 1950  5 16 "1950-05-16" "19:30:00" 3 "OK" 40 16 1  1  0 3 0 35.82  -99.12              35.87                 -99  7.3  100 0  87
      109 1950  5 18 "1950-05-18" "15:15:00" 3 "NE" 31  2 1  0  0 4 0 40.22 -100.63            40.2201 -100.62989999999999   .7  333 0  88
      110 1950  5 18 "1950-05-18" "19:00:00" 3 "KS" 20 12 3  2  0 5 0 38.12  -96.82              38.33              -96.82 14.9  100 0  89
      111 1950  5 18 "1950-05-18" "19:00:00" 3 "KS" 20 13 0  0  0 3 0 39.85   -98.3 39.850100000000005            -98.2999   .1   10 0  90
      112 1950  5 19 "1950-05-19" "01:30:00" 3 "KS" 20 14 0  0  0 3 0 37.82  -96.85 37.820100000000004  -96.84989999999999   .1   10 0  91
      113 1950  5 22 "1950-05-22" "12:39:00" 3 "OK" 40 17 1  0  0 0 0 35.67  -99.23 35.670100000000005            -99.2299  1.5  100 0  92
      114 1950  5 23 "1950-05-23" "18:00:00" 3 "KS" 20 15 0  0  0 4 0 39.75  -99.32            39.7501  -99.31989999999999   .2  250 0  93
      115 1950  5 24 "1950-05-24" "13:30:00" 3 "KS" 20 16 2  1  0 3 0 37.27  -99.48              37.28              -99.42  3.6   77 0  94
      116 1950  5 24 "1950-05-24" "17:30:00" 3 "KS" 20 17 1  0  0 4 0 39.38  -96.97              39.45              -96.78 10.7   40 0  95
      117 1950  5 24 "1950-05-24" "21:30:00" 3 "OK" 40 18 2  2  0 4 0 36.45  -99.03 36.450100000000006            -99.0299    1   33 0  96
      118 1950  5 25 "1950-05-25" "16:30:00" 3 "TX" 48 16 1  0  0 3 0 32.67 -101.88              32.65             -101.85    3  880 0  97
      119 1950  5 29 "1950-05-29" "14:48:00" 3 "OK" 40 19 1  0  0 2 0 36.38  -95.27 36.380100000000006  -95.26989999999999    1   33 0  98
      120 1950  5 29 "1950-05-29" "17:00:00" 3 "TX" 48 17 1  0  0 0 0 33.35  -96.92               33.3              -96.92  3.4  267 0  99
      121 1950  5 29 "1950-05-29" "23:00:00" 3 "LA" 22 20 1 13  0 4 0  32.1   -93.9 32.100100000000005            -93.8999    1   27 0 100
      end
      
      
      
      ren (slat    slon    elat    elon) (_Y1 _X1 _Y2 _X2)
      keep _ID _X1 _Y1 _X2 _Y2
      order _ID _X1 _Y1 _X2 _Y2
      
      save Tornado.dta, replace 
      
      copy "https://www2.census.gov/geo/tiger/GENZ2018/shp/cb_2018_us_county_500k.zip" "cb_2018_us_county_500k.zip"
      unzipfile cb_2018_us_county_500k.zip
      spshape2dta cb_2018_us_county_500k
      
      use cb_2018_us_county_500k, clear
      
      spmap using cb_2018_us_county_500k_shp if real(ST) <= 56 & !inlist(ST, "02", "15")  , id(_ID) osize(vthin   ..) ///
            ocolor(gs9 ..) arrow(data(Tornado.dta) direction(1) lcolor(red ..) hfcolor(red ..)  hsize(vtiny)    hbarbsize(vtiny) )
      
      exit

      Comment


      • #4
        of course it was better with the old dataset and with arrow command. Thank you a lot

        Comment

        Working...
        X