Announcement

Collapse
No announcement yet.
X
Collapse
  • Page of 1
  • Filter
  • Time
  • Show
Clear All
new posts
  • #1

    Counting Upwind and Nonupwind fires

    Hi, I have data on the cluster locations. with lat & log information and data on the fire location with lat & log information.
    Part of each data is shown below.
    I have to figure out the prevailing wind direction at the cluster location each day, and then count the number of fire events that happen on that day in the prevailing wind sector (defined as a 30 degrees circle around the prevailing wind direction).
    Fires occurring in the prevailing wind sector, are called upwind fires, they have pollutants blowing from the fire location towards the residential location (cluster location). Non-upwind fires are those that occur in all other sectors.
    This identification strategy is shown in the Figure attached.

    How I can count the upwind and non-upwind fire locations Please


    Click image for larger version Name: Figure 3.png Views: 1 Size: 38.1 KB ID: 1641511


    Last edited by Muhammad Ramzan; 18 Dec 2021, 22:56.
    Tags: None
  • #2
    Data of the FIRE locations


    Code:
    * Example generated by -dataex-. To install: ssc install dataex
clear
input long FID_F double(LATITUDE_F LONGITUDE_F BRIGHTNESS SCAN TRACK ACQ_DATE) str4 ACQ_TIME str5(SATELLITE INSTRUMENT) byte CONFIDENCE str4 VERSION double(BRIGHT_T31 FRP) str1 DAYNIGHT byte TYPE
 0 28.3755 69.9828 311.7 1.2 1.1 1.2939264e+12 "0607" "Terra" "MODIS" 66 "6.03" 291.6 12.3 "D" 0
 1 27.9079 69.1454   309 1.3 1.1 1.2939264e+12 "0607" "Terra" "MODIS" 60 "6.03" 293.5  9.8 "D" 0
 2 25.2932 69.0643 312.2 1.2 1.1 1.2939264e+12 "0608" "Terra" "MODIS" 57 "6.03" 296.6  9.6 "D" 0
 3 24.5419 69.1255 308.9 1.2 1.1 1.2939264e+12 "0608" "Terra" "MODIS" 37 "6.03" 295.5  6.3 "D" 0
 4 25.2143 68.8821   310 1.3 1.1 1.2939264e+12 "0608" "Terra" "MODIS" 45 "6.03" 294.9  9.5 "D" 0
 5 25.4614 69.2628 309.9 1.2 1.1 1.2939264e+12 "0608" "Terra" "MODIS" 40 "6.03" 296.5  6.4 "D" 0
 6 25.7585 68.8434 322.8 1.3 1.1 1.2939264e+12 "0608" "Terra" "MODIS" 77 "6.03" 294.4 25.6 "D" 0
 7 25.4621  69.269 310.2 1.2 1.1 1.2939264e+12 "0608" "Terra" "MODIS" 44 "6.03" 296.9  6.6 "D" 0
 8 24.9994 69.2118   311 1.2 1.1 1.2939264e+12 "0608" "Terra" "MODIS" 49 "6.03" 298.2  7.1 "D" 0
 9 25.7284 68.7147   319 2.2 1.4 1.2940128e+12 "0651" "Terra" "MODIS" 71 "6.03" 295.8 42.5 "D" 0
10 25.7276   68.72 322.8 2.2 1.4 1.2940128e+12 "0651" "Terra" "MODIS" 77 "6.03" 295.7   52 "D" 0
11 24.5625 67.8861 311.6   2 1.4 1.2940128e+12 "0651" "Terra" "MODIS" 53 "6.03" 300.1 15.3 "D" 0
12  25.832 68.8779 315.6 2.1 1.4 1.2940992e+12 "0555" "Terra" "MODIS" 64 "6.03" 296.7 28.2 "D" 0
13 24.4254 68.1614 310.5 2.2 1.5 1.2940992e+12 "0556" "Terra" "MODIS" 54 "6.03" 298.1 16.6 "D" 0
14 25.3055 68.4361 310.5 2.2 1.4 1.2940992e+12 "0556" "Terra" "MODIS" 36 "6.03" 294.6 19.6 "D" 0
15 25.4869 68.8397 310.3 2.1 1.4 1.2940992e+12 "0556" "Terra" "MODIS" 26 "6.03" 295.4 15.5 "D" 0
16 32.8975 73.0531 300.8 1.9 1.3 1.2941856e+12 "0636" "Terra" "MODIS" 42 "6.03" 289.4  8.6 "D" 0
17 25.4181 68.6151 313.8 1.4 1.2 1.2941856e+12 "0638" "Terra" "MODIS" 58 "6.03" 298.5 13.3 "D" 0
18  25.529  68.975 313.1 1.4 1.2 1.2941856e+12 "0638" "Terra" "MODIS" 38 "6.03" 300.8  9.4 "D" 0
19 25.5221 68.9815 312.6 1.4 1.2 1.2941856e+12 "0638" "Terra" "MODIS" 37 "6.03"   302  8.9 "D" 2
20 25.8422 69.1473 313.4 1.4 1.2 1.2941856e+12 "0638" "Terra" "MODIS" 57 "6.03" 300.3 10.5 "D" 0
21 26.3668 68.7947 313.6 1.3 1.2 1.2941856e+12 "0638" "Terra" "MODIS" 46 "6.03"   301  9.5 "D" 0
22 25.0426 68.1704 310.9 1.3 1.1 1.2941856e+12 "0639" "Terra" "MODIS" 26 "6.03" 299.9  5.9 "D" 0
23 24.5502 68.2785 313.6 1.3 1.1 1.2941856e+12 "0639" "Terra" "MODIS" 61 "6.03" 299.8 11.8 "D" 0
24 25.5185 68.8008 311.6 3.5 1.8  1.294272e+12 "0543" "Terra" "MODIS" 25 "6.03" 295.3 42.4 "D" 0
25 25.5122 68.8078   316 3.5 1.8  1.294272e+12 "0543" "Terra" "MODIS" 67 "6.03" 295.6 71.3 "D" 0
26 32.9467 72.5605 311.6 1.2 1.1 1.2943584e+12 "0624" "Terra" "MODIS" 71 "6.03" 283.5 15.3 "D" 0
27  25.578 68.5578 315.1   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 63 "6.03" 297.2 10.3 "D" 0
28 25.5834 68.8838 316.6   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 65 "6.03" 302.1  9.8 "D" 0
29 25.5481 68.5728 310.7   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 41 "6.03" 295.5  7.3 "D" 0
30 25.6577 68.8208 327.8   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 81 "6.03" 300.9 21.6 "D" 0
31 25.5849 68.8738   317   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 66 "6.03" 302.1  9.4 "D" 0
32 25.0266  68.299 314.5   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 58 "6.03" 298.2    7 "D" 0
33 25.5659 68.6374 313.4   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 57 "6.03" 299.9  8.5 "D" 0
34 25.6593 68.8108 322.6   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 75 "6.03" 300.7 15.9 "D" 0
35 25.0296 68.2793 335.5   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 87 "6.03" 300.6 30.2 "D" 0
36 27.4382 68.5207 309.6   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 50 "6.03" 294.2  7.4 "D" 0
37 25.5344 68.6624 315.6   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 64 "6.03" 295.8 10.4 "D" 0
38 26.0908 68.1532 315.2   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 64 "6.03" 297.9  9.5 "D" 0
39 24.3432 68.0574 314.6   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 53 "6.03" 300.6  7.4 "D" 0
40 25.7842  68.353 311.5   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 48 "6.03" 298.1    5 "D" 0
41 26.0923 68.1433 310.5   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 38 "6.03" 297.9  5.9 "D" 0
42 25.1073 69.0352 319.7   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 69 "6.03"   301   12 "D" 0
43 24.6484 68.1765 315.4   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 61 "6.03" 300.7  9.1 "D" 0
44 24.9399 68.8046 312.8   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 46 "6.03" 299.9  7.2 "D" 0
45  25.555 68.7689 310.6   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 30 "6.03" 297.6  6.2 "D" 0
46  25.103 68.6452 311.6   1   1 1.2943584e+12 "0626" "Terra" "MODIS" 33 "6.03" 299.8  5.4 "D" 0
47 32.8428 74.0425 301.4   1   1 1.2945312e+12 "0612" "Terra" "MODIS" 47 "6.03" 289.4  3.6 "D" 0
48 32.8445 74.0316   301   1   1 1.2945312e+12 "0612" "Terra" "MODIS" 44 "6.03" 289.4  3.6 "D" 0
49 24.5014 67.8054 307.3 1.2 1.1 1.2945312e+12 "0614" "Terra" "MODIS" 37 "6.03" 287.2  8.3 "D" 0
50 25.1945 68.8834 307.9 1.1   1 1.2945312e+12 "0614" "Terra" "MODIS" 62 "6.03" 290.5    7 "D" 0
51 28.5047 70.0123 304.5   3 1.6 1.2946176e+12 "0656" "Terra" "MODIS" 38 "6.03" 290.7 24.4 "D" 0
52 28.5119 69.9885 308.8   3 1.6 1.2946176e+12 "0656" "Terra" "MODIS" 58 "6.03" 290.6 33.9 "D" 0
53  28.511 69.9829 305.2   3 1.6 1.2946176e+12 "0656" "Terra" "MODIS" 43 "6.03" 290.5 25.9 "D" 0
54 25.2322 69.3618 312.7 3.2 1.7 1.2946176e+12 "0657" "Terra" "MODIS" 60 "6.03" 297.3 42.2 "D" 0
55 25.2388  69.331 310.9 3.2 1.7 1.2946176e+12 "0657" "Terra" "MODIS" 54 "6.03" 298.3 34.4 "D" 0
56 25.2349 69.3556 315.9 3.2 1.7 1.2946176e+12 "0657" "Terra" "MODIS" 67 "6.03" 297.7 55.9 "D" 0
57 33.7988 71.8645 306.4 1.4 1.2  1.294704e+12 "0559" "Terra" "MODIS" 63 "6.03" 293.7  9.9 "D" 0
58 25.5115 68.7621 318.7 1.6 1.3  1.294704e+12 "0601" "Terra" "MODIS" 71 "6.03" 295.3   29 "D" 0
59 25.5963 69.2253 310.7 1.5 1.2  1.294704e+12 "0601" "Terra" "MODIS" 45 "6.03" 297.4 10.6 "D" 0
60 25.5097 68.7783 309.6 1.6 1.3  1.294704e+12 "0601" "Terra" "MODIS" 31 "6.03" 295.7 13.2 "D" 0
61 24.7311 69.2197 310.2 1.5 1.2  1.294704e+12 "0602" "Terra" "MODIS" 24 "6.03" 297.2  9.2 "D" 0
62 24.7361 69.1754 311.7 1.5 1.2  1.294704e+12 "0602" "Terra" "MODIS" 50 "6.03" 298.1 11.7 "D" 0
63 24.8765 68.9395 312.4 1.6 1.2  1.294704e+12 "0602" "Terra" "MODIS" 53 "6.03" 296.9 14.2 "D" 0
64  33.945 72.7319 303.7   2 1.4  1.294704e+12 "1701" "Terra" "MODIS" 57 "6.03"   279 26.6 "N" 0
65 33.9481  72.725 309.8   2 1.4  1.294704e+12 "1701" "Terra" "MODIS" 78 "6.03" 279.1 37.5 "N" 0
66 26.8429 68.0396 308.7 1.5 1.2 1.2947904e+12 "0644" "Terra" "MODIS" 33 "6.03" 294.8 10.7 "D" 0
67 25.8817  68.855 310.2 1.7 1.3 1.2947904e+12 "0644" "Terra" "MODIS" 43 "6.03" 297.6 11.6 "D" 0
68 33.7155 73.5648 300.4 1.6 1.3 1.2948768e+12 "0547" "Terra" "MODIS" 35 "6.03" 287.6  9.3 "D" 0
69  33.717 73.5505 304.8 1.6 1.3 1.2948768e+12 "0547" "Terra" "MODIS" 60 "6.03" 287.7 14.5 "D" 0
70  32.962 72.5819 303.6 1.4 1.2 1.2949632e+12 "0630" "Terra" "MODIS" 56 "6.03" 291.2  6.8 "D" 0
71 32.9646 72.5675 307.8 1.4 1.2 1.2949632e+12 "0630" "Terra" "MODIS" 66 "6.03" 291.6   11 "D" 0
72 24.2588 67.8528   314 1.1   1 1.2949632e+12 "0632" "Terra" "MODIS" 56 "6.03" 303.8  5.7 "D" 0
73 24.8524 68.3444 333.6 1.1 1.1 1.2949632e+12 "0632" "Terra" "MODIS" 85 "6.03" 300.9 32.7 "D" 0
74 25.5166 68.8074 311.9 1.1 1.1 1.2949632e+12 "0632" "Terra" "MODIS" 39 "6.03" 299.8  7.7 "D" 0
75 26.3151 68.3934 311.3 1.1   1 1.2949632e+12 "0632" "Terra" "MODIS" 33 "6.03" 295.4  7.2 "D" 0
76 26.1548 68.2405 309.9 1.1   1 1.2949632e+12 "0632" "Terra" "MODIS" 39 "6.03" 296.1    7 "D" 0
77 25.0206 68.2934   312 1.1   1 1.2949632e+12 "0632" "Terra" "MODIS" 34 "6.03" 299.3  5.4 "D" 0
78 25.1174  68.779   315 1.2 1.1 1.2949632e+12 "0632" "Terra" "MODIS" 56 "6.03" 300.5  8.8 "D" 0
79 27.0304 67.7323 312.6   1   1 1.2949632e+12 "0632" "Terra" "MODIS" 51 "6.03" 295.4  9.3 "D" 0
80 24.9115  68.457 313.1 1.1 1.1 1.2949632e+12 "0632" "Terra" "MODIS" 47 "6.03" 301.3  7.1 "D" 0
81 25.0263 68.4126 312.1 1.1 1.1 1.2949632e+12 "0632" "Terra" "MODIS" 38 "6.03" 300.3  6.5 "D" 0
82 25.3589 68.9391 312.3 1.2 1.1 1.2949632e+12 "0632" "Terra" "MODIS" 37 "6.03" 301.2  7.1 "D" 0
83 33.6932 73.0225 315.1   1   1  1.295136e+12 "0617" "Terra" "MODIS" 75 "6.03" 289.8 13.9 "D" 2
84 33.8535 70.7187 305.5   1   1  1.295136e+12 "0617" "Terra" "MODIS" 47 "6.03" 284.6    8 "D" 0
85 34.5231 73.6189 302.8 1.1   1  1.295136e+12 "0617" "Terra" "MODIS" 24 "6.03" 290.6  7.6 "D" 0
86 34.8225  72.905 303.6   1   1  1.295136e+12 "0617" "Terra" "MODIS" 30 "6.03" 286.7  5.9 "D" 0
87 29.8186 72.7161 307.9 1.1   1  1.295136e+12 "0618" "Terra" "MODIS" 64 "6.03" 294.1    7 "D" 0
88 32.1207 69.6685 311.5   1   1  1.295136e+12 "0618" "Terra" "MODIS" 70 "6.03" 294.5  7.9 "D" 0
89 29.2967 71.5168 307.4   1   1  1.295136e+12 "0618" "Terra" "MODIS" 48 "6.03" 295.7  5.7 "D" 0
90 26.0844 68.1442 310.1   1   1  1.295136e+12 "0619" "Terra" "MODIS" 43 "6.03" 298.8  5.8 "D" 0
91 25.8544 69.0566 313.3   1   1  1.295136e+12 "0619" "Terra" "MODIS" 54 "6.03" 299.9    7 "D" 0
92  27.988 69.0124 311.4   1   1  1.295136e+12 "0619" "Terra" "MODIS" 54 "6.03" 298.8    6 "D" 0
93 24.9739 68.5816 319.6   1   1  1.295136e+12 "0620" "Terra" "MODIS" 68 "6.03" 303.3 10.6 "D" 0
94  25.029 68.6478 313.9   1   1  1.295136e+12 "0620" "Terra" "MODIS" 49 "6.03" 301.7  6.3 "D" 0
95 25.4271 68.9523 312.7   1   1  1.295136e+12 "0620" "Terra" "MODIS" 38 "6.03" 301.5  5.2 "D" 0
96 24.9708 68.4086 319.2   1   1  1.295136e+12 "0620" "Terra" "MODIS" 67 "6.03" 300.9 11.9 "D" 0
97 24.9725 68.5915 314.4   1   1  1.295136e+12 "0620" "Terra" "MODIS" 51 "6.03" 302.6    6 "D" 0
98 25.5034 69.1292   312   1   1  1.295136e+12 "0620" "Terra" "MODIS" 30 "6.03"   299  5.4 "D" 0
99 25.4838 68.8164 317.3   1   1  1.295136e+12 "0620" "Terra" "MODIS" 65 "6.03" 299.5 10.1 "D" 0
end
format %tc ACQ_DATE

    Comment

    • #3
      Data of the Cluster locations


      Code:
      * Example generated by -dataex-. To install: ssc install dataex
clear
input int FID_C str14 DHSID int DHSCLUST str2 CCFIPS byte DHSREGCO str11 DHSREGNA str1 URBAN_RURA double(LATNUM_C LONGNUM_C)
 0 "PK201700000001"   1 "PK" 3 "kpk"    "R" 36.4499182254 72.5715578674
 1 "PK201700000002"   2 "PK" 3 "kpk"    "U" 35.8919142536 71.7268731073
 2 "PK201700000003"   3 "PK" 3 "kpk"    "R" 35.1695664132 71.8344582034
 3 "PK201700000004"   4 "PK" 3 "kpk"    "R" 35.4247289555 72.1639311741
 4 "PK201700000005"   5 "PK" 3 "kpk"    "R" 35.0056955272 71.7764776856
 5 "PK201700000006"   6 "PK" 3 "kpk"    "R"  34.944432422 71.6551094148
 6 "PK201700000007"   7 "PK" 3 "kpk"    "R" 34.7355764204 71.8375579917
 7 "PK201700000008"   8 "PK" 3 "kpk"    "R"  34.689715398 72.1238682541
 8 "PK201700000009"   9 "PK" 3 "kpk"    "R" 34.8138638781 72.3664785607
 9 "PK201700000010"  10 "PK" 3 "kpk"    "U" 34.7486443502 72.4169229558
10 "PK201700000011"  11 "PK" 3 "kpk"    "U" 34.7435888145 72.4244339897
11 "PK201700000012"  12 "PK" 3 "kpk"    "U" 34.9987162105  72.540403036
12 "PK201700000013"  13 "PK" 3 "kpk"    "U" 34.8614502539 72.2755876121
13 "PK201700000014"  14 "PK" 3 "kpk"    "R" 34.8487858784 72.8823515247
14 "PK201700000015"  15 "PK" 3 "kpk"    "R" 34.4607255601 72.4609016637
15 "PK201700000016"  16 "PK" 3 "kpk"    "R" 34.4672412698 71.8554040871
16 "PK201700000017"  17 "PK" 3 "kpk"    "U" 34.5657254918 71.9057528327
17 "PK201700000022"  22 "PK" 3 "kpk"    "R" 34.4484616298 73.3836077956
18 "PK201700000023"  23 "PK" 3 "kpk"    "R" 34.5177441604  73.265701841
19 "PK201700000024"  24 "PK" 3 "kpk"    "U" 34.4309865058 73.2175568765
20 "PK201700000025"  25 "PK" 3 "kpk"    "R" 34.6745435574 73.0526476984
21 "PK201700000026"  26 "PK" 3 "kpk"    "R"  34.089194344  73.383530112
22 "PK201700000027"  27 "PK" 3 "kpk"    "R" 34.0123092115 73.2395696716
23 "PK201700000028"  28 "PK" 3 "kpk"    "U"  34.138836295 73.2200455137
24 "PK201700000029"  29 "PK" 3 "kpk"    "U" 33.9656834462 73.2177445651
25 "PK201700000030"  30 "PK" 3 "kpk"    "R"  34.134098912 72.9400355457
26 "PK201700000031"  31 "PK" 3 "kpk"    "R"  33.809855724 73.1174442725
27 "PK201700000032"  32 "PK" 3 "kpk"    "U" 34.0005178374 72.9282441939
28 "PK201700000033"  33 "PK" 3 "kpk"    "R" 34.5448600486 72.8395490711
29 "PK201700000034"  34 "PK" 3 "kpk"    "R" 34.4866438792 72.0836053773
30 "PK201700000035"  35 "PK" 3 "kpk"    "R" 34.1774184331 72.1705950662
31 "PK201700000036"  36 "PK" 3 "kpk"    "R" 34.4187948321 72.3278790408
32 "PK201700000037"  37 "PK" 3 "kpk"    "U" 34.2801182307 72.1756213181
33 "PK201700000038"  38 "PK" 3 "kpk"    "U"  34.272924996 72.1737534586
34 "PK201700000039"  39 "PK" 3 "kpk"    "U" 34.2766122933 71.9418244891
35 "PK201700000040"  40 "PK" 3 "kpk"    "R" 34.2839028716 72.3767401211
36 "PK201700000041"  41 "PK" 3 "kpk"    "R" 34.1860364687 72.3979275223
37 "PK201700000042"  42 "PK" 3 "kpk"    "U" 34.1277885467  72.475738516
38 "PK201700000043"  43 "PK" 3 "kpk"    "U" 34.1172697963 72.5035644427
39 "PK201700000044"  44 "PK" 3 "kpk"    "R" 34.1342646468 71.6908071196
40 "PK201700000045"  45 "PK" 3 "kpk"    "R" 34.3032010258 71.6118536887
41 "PK201700000046"  46 "PK" 3 "kpk"    "U" 34.1654929177 71.7598850094
42 "PK201700000047"  47 "PK" 3 "kpk"    "U" 34.1927121501 71.7638071487
43 "PK201700000048"  48 "PK" 3 "kpk"    "R" 34.1033784498 71.6019087575
44 "PK201700000049"  49 "PK" 3 "kpk"    "R" 33.9583254271 71.6146843403
45 "PK201700000050"  50 "PK" 3 "kpk"    "R"  33.945303642 71.6865737739
46 "PK201700000051"  51 "PK" 3 "kpk"    "U" 33.9903730476 71.5375839682
47 "PK201700000052"  52 "PK" 3 "kpk"    "U" 34.0113739915 71.5212273276
48 "PK201700000053"  53 "PK" 3 "kpk"    "U" 34.0086142429 71.5288680618
49 "PK201700000054"  54 "PK" 3 "kpk"    "U" 34.0100262472 71.5353768295
50 "PK201700000055"  55 "PK" 3 "kpk"    "U" 34.0007967437 71.5339475817
51 "PK201700000056"  56 "PK" 3 "kpk"    "U" 34.0004557253 71.5777313753
52 "PK201700000057"  57 "PK" 3 "kpk"    "U" 33.9872542054 71.5971479859
53 "PK201700000058"  58 "PK" 3 "kpk"    "U" 33.9979765861  71.552242676
54 "PK201700000059"  59 "PK" 3 "kpk"    "U" 33.9931561414 71.5269282239
55 "PK201700000060"  60 "PK" 3 "kpk"    "U" 33.9972634656  71.538230657
56 "PK201700000061"  61 "PK" 3 "kpk"    "U" 34.1566380784 71.4504940907
57 "PK201700000062"  62 "PK" 3 "kpk"    "U" 34.0089296333 71.5420730915
58 "PK201700000063"  63 "PK" 3 "kpk"    "U" 33.9959747569 71.5266948364
59 "PK201700000064"  64 "PK" 3 "kpk"    "R"   33.91341027 72.0176179786
60 "PK201700000065"  65 "PK" 3 "kpk"    "R"  34.005454073  71.779015226
61 "PK201700000066"  66 "PK" 3 "kpk"    "U" 34.1227926199  71.960556304
62 "PK201700000067"  67 "PK" 3 "kpk"    "U" 33.9249335029 72.0390322529
63 "PK201700000068"  68 "PK" 3 "kpk"    "U"  33.938676527 71.7990767178
64 "PK201700000069"  69 "PK" 3 "kpk"    "R" 33.4840898942 71.3995560447
65 "PK201700000070"  70 "PK" 3 "kpk"    "U" 33.5800871476 71.4282835193
66 "PK201700000071"  71 "PK" 3 "kpk"    "U" 33.5742304548 71.4584730272
67 "PK201700000072"  72 "PK" 3 "kpk"    "U" 33.3626754137 70.5539741567
68 "PK201700000073"  73 "PK" 3 "kpk"    "R" 33.1748642707 71.1865990396
69 "PK201700000074"  74 "PK" 3 "kpk"    "R" 32.9790426464 70.6790446519
70 "PK201700000075"  75 "PK" 3 "kpk"    "U" 32.9899955976 70.6147092986
71 "PK201700000076"  76 "PK" 3 "kpk"    "R" 32.5176163005  70.921225542
72 "PK201700000077"  77 "PK" 3 "kpk"    "U" 32.6152064682 70.9164198863
73 "PK201700000078"  78 "PK" 3 "kpk"    "R" 31.7568504356 70.7176522238
74 "PK201700000079"  79 "PK" 3 "kpk"    "R" 31.9540009492  70.347534046
75 "PK201700000080"  80 "PK" 3 "kpk"    "U" 31.8439080198 70.8969695143
76 "PK201700000081"  81 "PK" 3 "kpk"    "U" 31.7958158111 70.8967685943
77 "PK201700000082"  82 "PK" 3 "kpk"    "U" 31.5949781671 70.8234721928
78 "PK201700000083"  83 "PK" 3 "kpk"    "R" 32.1718730258 70.2509650811
79 "PK201700000084"  84 "PK" 3 "kpk"    "U" 32.2237092477 70.3756924373
80 "PK201700000085"  85 "PK" 1 "punjab" "R" 33.8684398977 72.7646453356
81 "PK201700000086"  86 "PK" 1 "punjab" "R" 33.1642327843 72.1126681286
82 "PK201700000087"  87 "PK" 1 "punjab" "U" 33.4969445401 71.9532672364
83 "PK201700000088"  88 "PK" 1 "punjab" "R" 33.3056905393 73.2131831488
84 "PK201700000089"  89 "PK" 1 "punjab" "R" 33.4829668086  73.477666433
85 "PK201700000090"  90 "PK" 1 "punjab" "R" 33.8861385011 73.4286611655
86 "PK201700000091"  91 "PK" 1 "punjab" "U" 33.6533953758  73.410579747
87 "PK201700000092"  92 "PK" 1 "punjab" "U" 33.7355745855 72.7908112076
88 "PK201700000093"  93 "PK" 1 "punjab" "U" 33.4804065011  73.419872597
89 "PK201700000094"  94 "PK" 1 "punjab" "U" 33.9046198718 73.3889983098
90 "PK201700000095"  95 "PK" 1 "punjab" "R" 32.9495356629 73.5862306661
91 "PK201700000096"  96 "PK" 1 "punjab" "U" 32.6370501442 73.1454932241
92 "PK201700000097"  97 "PK" 1 "punjab" "R" 32.9951480988 72.8431774034
93 "PK201700000098"  98 "PK" 1 "punjab" "R" 32.7557417861 71.9793290391
94 "PK201700000099"  99 "PK" 1 "punjab" "U" 32.9278043706   72.87326944
95 "PK201700000100" 100 "PK" 1 "punjab" "U"  32.866188617 72.7629057467
96 "PK201700000101" 101 "PK" 1 "punjab" "R" 32.3227490272 72.9209185559
97 "PK201700000102" 102 "PK" 1 "punjab" "R" 32.1487643466  72.911008241
98 "PK201700000103" 103 "PK" 1 "punjab" "R" 32.2318100144 72.4440213163
99 "PK201700000104" 104 "PK" 1 "punjab" "U"  32.035856394 72.7546317476
end

      Comment

      • #4
        Hi Stata Users,

        Is there anyone who can help with this question? I have exactly the same request.

        Soodeh

        Comment

        • #5
          You can calculate the angle of the points within some radius of the cluster location and then find all points within reference line plus or minus _pi/6.

          Code:
          clear
set seed 2457
set obs 1
gen lat = 0
gen lon = 0
//50 random points in the unit cirlc
set obs 50
gen r = sqrt(runiform())
gen theta = runiform()*2*_pi
gen y = r*sin(theta)
gen x = r*cos(theta)
gen angle = atan2(y, x)
gen radius = sqrt(x^2 + y^2)

//Tag points +/- 30 degrees and less than unit radius
gen tag = cond(angle < _pi/6 & angle > 0 & radius <=1 ,1, /// 
    cond(angle< 0 & angle > -_pi/6 & radius <=1 , 1,0))


scatter lat lon /// 
    || function y = sqrt(1 - x^2), range(-1 1) lc(black)  /// 
    || function y = -sqrt(1 - x^2), range(-1 1) lc(black) /// 
    || scatteri 0 0 .5  `=sqrt(3)/2' , recast(line) lc(gs4) /// 
    || scatteri 0 -1.05 0 1.05 , recast(line)  /// 
    || scatteri 0 1.05 "UpWind" , ms(A) ms(p) mlabcolor(black) /// 
    || scatteri 0 0  -.5  `=sqrt(3)/2', recast(line)  lc(gs4) /// 
    || scatter y x if tag == 1, mlabel(tag) mcolor(red) mlabcolor(red) /// 
    || scatter y x if tag == 0,  mcolor(gs10) ms(Oh) /// 
    || ,  aspect(1) legend(off)
          Click image for larger version Name: wind.png Views: 1 Size: 64.5 KB ID: 1724881

          Comment

          • #6
            Thanks Scott Merryman. This is very helpful.

            Comment

            Previous Next
            Working...
            X