\u043f\u0435\u0440\u0432\u043e\u0439 \u0447\u0430\u0441\u0442\u0438<\/a><\/noindex>.<\/p>\n\u0418\u0441\u0445\u043e\u0434\u043d\u044b\u0439 \u043a\u043e\u0434 \u043f\u0430\u0440\u0441\u0435\u0440\u0430, \u043e\u043f\u0438\u0441\u0430\u043d\u043d\u044b\u0439 \u0432 \u0441\u0442\u0430\u0442\u044c\u0435, \u043f\u0440\u0438\u0432\u0435\u0434\u0435\u043d \u043f\u043e\u0434 \u043a\u0430\u0442\u043e\u043c. \u042d\u0442\u043e \u0442\u0435\u0441\u0442\u043e\u0432\u044b\u0439 \u043f\u0440\u0438\u043c\u0435\u0440, \u043d\u0435 \u043f\u0440\u0435\u0442\u0435\u043d\u0434\u0443\u044e\u0449\u0438\u0439 \u043d\u0430 production, \u043d\u043e \u043a\u043e\u0435-\u0447\u0442\u043e \u0432 \u043d\u0435\u043c \u0440\u0430\u0431\u043e\u0442\u0430\u0435\u0442, \u0438 \u043f\u0430\u0440\u0441\u0438\u0442\u044c \u0437\u0430\u043f\u0438\u0441\u0430\u043d\u043d\u044b\u0439 \u0432\u044b\u0448\u0435 \u0444\u0430\u0439\u043b, \u0438\u043c \u043c\u043e\u0436\u043d\u043e.\u0418\u0441\u0445\u043e\u0434\u043d\u044b\u0439 \u043a\u043e\u0434 (Python)<\/b><\/p>\nfrom __future__ import print_function\n\nfrom scipy.io import wavfile\nfrom scipy import signal\nimport matplotlib.pyplot as plt\nimport numpy as np\nimport math\nimport sys\n\n\ndef parse_message(data, start, bit_len):\n max_len = bit_len*128\n A = data[start:start + max_len]\n A = signal.resample(A, 10*max_len)\n bits = np.zeros(10*max_len)\n bit_len *= 10\n start_data = bit_len*8\n\n # Parse first 8 bits\n bits_str = \"\"\n for p in range(8):\n pos = start_data + bit_len*p\n p1, p2 = A[pos: pos + bit_len\/2], A[pos + bit_len\/2: pos + bit_len]\n avg1, avg2 = np.average(p1), np.average(p2)\n if avg1 < avg2:\n bits_str += \"0\"\n elif avg1 > avg2:\n bits_str += \"1\"\n\n df = int(bits_str[0:5], 2)\n\n # Aircraft address (db - https:\/\/junzis.com\/adb\/?q=3b1c5c )\n bits_str = \"\"\n for p in range(8, 32):\n pos = start_data + bit_len * p\n p1, p2 = A[pos: pos + bit_len \/ 2], A[pos + bit_len \/ 2: pos + bit_len]\n avg1, avg2 = np.average(p1), np.average(p2)\n if avg1 < avg2:\n bits_str += \"0\"\n elif avg1 > avg2:\n bits_str += \"1\"\n # print \"Aircraft address:\", bits_str, hex(int(bits_str, 2))\n address = hex(int(bits_str, 2))\n\n # Filter specific aircraft (optional)\n # if address != \"0x3c5ee2\":\n # return\n\n if df == 16 or df == 17 or df == 18 or df == 19 or df == 20 or df == 21:\n # print \"Pos:\", start, \"DF:\", msg_type\n\n # Data (56bit)\n bits_str = \"\"\n for p in range(32, 88):\n pos = start_data + bit_len*p\n p1, p2 = A[pos: pos + bit_len\/2], A[pos + bit_len\/2: pos + bit_len]\n avg1, avg2 = np.average(p1), np.average(p2)\n if avg1 < avg2:\n bits_str += \"0\"\n # bits[pos + bit_len \/ 2] = 50\n elif avg1 > avg2:\n bits_str += \"1\"\n\n # http:\/\/www.lll.lu\/~edward\/edward\/adsb\/DecodingADSBposition.html\n # print \"Data:\"\n # print bits_str[:8], bits_str[8:20], bits_str[20:22], bits_str[22:22+17], bits_str[39:39+17]\n # Type Code:\n tc, ec = int(bits_str[:5], 2), int(bits_str[5:8], 2)\n # print(\"DF:\", df, \"TC:\", tc)\n \n # 1 - 4 Aircraft identification\n # 5 - 8 Surface position\n # 9 - 18 Airborne position (w\/ Baro Altitude)\n # 19 Airborne velocities\n\n if tc >= 1 and tc <= 4: # and (df == 17 or df == 18):\n print(\"Aircraft address:\", address)\n print(\"Data:\")\n print(bits_str[:8], bits_str[8:14], bits_str[14:20], bits_str[20:26], bits_str[26:32], bits_str[32:38], bits_str[38:44])\n\n symbols = \"#ABCDEFGHIJKLMNOPQRSTUVWXYZ#####_###############0123456789######\"\n code_str = \"\"\n for p in range(8):\n c = int(bits_str[8 + 6*p:8 + 6*(p + 1)], 2)\n code_str += symbols[c]\n print(\"Aircraft Identification:\", code_str.replace('#', ''))\n print()\n if tc == 11:\n print(\"Aircraft address:\", address)\n print(\"Data: (11)\")\n print(bits_str[:8], bits_str[8:20], bits_str[20:22], bits_str[22:22+17], bits_str[39:39+17])\n\n # Bit 22 contains the F flag which indicates which CPR format is used (odd or even)\n # First frame has F flag = 0 so is even and the second frame has F flag = 1 so odd\n # f = bits_str[21:22]\n # print(\"F:\", int(f, 2))\n\n # Altitude\n alt1b = bits_str[8:20]\n if alt1b[-5] == '1':\n bits = alt1b[:-5] + alt1b[-4:]\n n = int(bits, 2)\n alt_ft = n*25 - 1000\n print(\"Alt (ft)\", alt_ft)\n\n # lat_dec = int(bits_str[22:22+17], 2)\n # lon_dec = int(bits_str[39:39+17], 2)\n # print(\"Lat\/Lon:\", lat_dec, lon_dec)\n\n # http:\/\/airmetar.main.jp\/radio\/ADS-B%20Decoding%20Guide.pdf\n print()\n if tc == 19:\n print(\"Aircraft address:\", address)\n print(\"Data:\")\n # print(bits_str)\n print(bits_str[:5], bits_str[5:8], bits_str[8:10], bits_str[10:13], bits_str[13] ,bits_str[14:24], bits_str[24], bits_str[25:35], bits_str[35:36], bits_str[36:65])\n\n subtype = int(bits_str[5:8], 2)\n # https:\/\/mode-s.org\/decode\/adsb\/airborne-velocity.html\n spd, hdg, rocd = -1, -1, -1\n if subtype == 1 or subtype == 2:\n print(\"Velocity Subtype 1: Ground speed\")\n \n v_ew_sign = int(bits_str[13], 2)\n v_ew = int(bits_str[14:24], 2) - 1 # east-west velocity\n \n v_ns_sign = int(bits_str[24], 2)\n v_ns = int(bits_str[25:35], 2) - 1 # north-south velocity\n \n v_we = -1*v_ew if v_ew_sign else v_ew\n v_sn = -1*v_ns if v_ns_sign else v_ns\n \n spd = math.sqrt(v_sn*v_sn + v_we*v_we) # unit in kts\n \n hdg = math.atan2(v_we, v_sn)\n hdg = math.degrees(hdg) # convert to degrees\n hdg = hdg if hdg >= 0 else hdg + 360 # no negative val\n if subtype == 3:\n print(\"Subtype Subtype 3: Airspeed\")\n hdg = int(bits_str[14:24], 2)\/1024.0*360.0\n spd = int(bits_str[25:35], 2)\n \n vr_sign = int(bits_str[36], 2)\n vr = int(bits_str[36:45], 2)\n rocd = -1*vr if vr_sign else vr # rate of climb\/descend\n print(\"Speed (kts):\", spd, \"Rate:\", rocd, \"Heading:\", hdg)\n print()\n\n # print()\n\ndef calc_coordinates():\n def _cprN(lat, is_odd):\n nl = _cprNL(lat) - is_odd\n return nl if nl > 1 else 1\n\n def _cprNL(lat):\n try:\n nz = 15\n a = 1 - math.cos(math.pi \/ (2 * nz))\n b = math.cos(math.pi \/ 180.0 * abs(lat)) ** 2\n nl = 2 * math.pi \/ (math.acos(1 - a\/b))\n return int(math.floor(nl))\n except:\n # happens when latitude is +\/-90 degree\n return 1\n \n def floor_(x):\n return int(math.floor(x))\n \n lat1b, lon1b, alt1b = \"10111000111010011\", \"10000110111111000\", \"000101111001\"\n lat2b, lon2b, alt2b = \"10010011101011100\", \"10000011000011011\", \"000101110111\"\n lat1, lon1, alt1 = int(lat1b, 2), int(lon1b, 2), int(alt1b, 2)\n lat2, lon2, alt2 = int(lat2b, 2), int(lon2b, 2), int(alt2b, 2)\n \n # 131072 is 2^17, since CPR lat and lon are 17 bits each\n cprlat_even, cprlon_even = lat1\/131072.0, lon1\/131072.0\n cprlat_odd, cprlon_odd = lat2\/131072.0, lon2\/131072.0\n print(cprlat_even, cprlon_even)\n\n j = floor_(59*cprlat_even - 60*cprlat_odd)\n print(j)\n\n air_d_lat_even = 360.0 \/ 60\n air_d_lat_odd = 360.0 \/ 59\n\n # Lat\n lat_even = float(air_d_lat_even * (j % 60 + cprlat_even))\n lat_odd = float(air_d_lat_odd * (j % 59 + cprlat_odd))\n if lat_even >= 270:\n lat_even = lat_even - 360\n if lat_odd >= 270:\n lat_odd = lat_odd - 360\n\n # Lon\n ni = _cprN(lat_even, 0)\n m = floor_(cprlon_even * (_cprNL(lat_even)-1) - cprlon_odd * _cprNL(lat_even) + 0.5)\n lon = (360.0 \/ ni) * (m % ni + cprlon_even)\n print(\"Lat\", lat_even, \"Lon\", lon)\n\n # Altitude\n # Q-bit (bit 48) indicates whether the altitude is encoded in multiples of 25 or 100 ft (0: 100 ft, 1: 25 ft)\n # The value can represent altitudes from -1000 to +50175 ft.\n if alt1b[-5] == '1':\n bits = alt1b[:-5] + alt1b[-4:]\n n = int(bits, 2)\n alt_ft = n*25 - 1000\n print(\"Alt (ft)\", alt_ft)\n\n\nfs, data = wavfile.read(\"adsb_20190311_191728Z_1090000kHz_RF.wav\")\nT = 1\/fs\n\nprint(\"Sample rate %f MS\/s\" % (fs \/ 1e6))\nprint(\"Cnt samples %d\" % len(data))\nprint(\"Duration: %f s\" % (T * len(data)))\n\ndata = data.astype(float)\n\ncnt = data.shape[0]\n# Processing only part on file (faster):\n# cnt = 10000000\n# data = data[:cnt]\nprint(\"Processing I\/Q...\")\nI, Q = data[:, 0], data[:, 1]\nA = np.sqrt(I*I + Q*Q)\n\nbits = np.zeros(cnt)\n\n# To see scope without any processing, uncomment\n# plt.plot(A)\n# plt.show()\n# sys.exit(0)\n\nprint(\"Extracting signals...\")\n\npos = 0\navg = 200\nmsg_start = 0\n# Find beginning of each signal\nwhile pos < cnt - 16*1024:\n # P1 - message start\n while pos < cnt - 16*1024:\n if A[pos] < avg and A[pos+1] > avg and pos - msg_start > 1000:\n msg_start = pos\n bits[pos] = 100\n pos += 4\n break\n pos += 1\n\n start1, start2, start3, start4 = msg_start, 0, 0, 0\n # P2\n while pos < cnt - 16*1024:\n if A[pos] < avg and A[pos+1] > avg:\n start2 = pos\n bits[pos] = 90\n pos += 1\n break\n pos += 1\n # P3\n while pos < cnt - 16*1024:\n if A[pos] < avg and A[pos+1] > avg:\n start3 = pos\n bits[pos] = 80\n pos += 1\n break\n pos += 1\n # P4\n while pos < cnt - 16*1024:\n if A[pos] < avg and A[pos+1] > avg:\n start4 = pos\n bits[pos] = 70\n pos += 1\n break\n pos += 1\n\n\n sig_diff = start4 - start1\n if 20 < sig_diff < 25:\n bits[msg_start] = 500\n bit_len = int((start4 - start1) \/ 4.5)\n # print(pos, start1, start4, ' - ', bit_len)\n # start = start1 + 8*bit_len\n parse_message(A, msg_start, bit_len)\n\n pos += 450\n\n# For debugging: check signal start\n# plt.plot(A)\n# plt.plot(bits)\n# plt.show()\n<\/code><\/pre>\n\u041d\u0430\u0434\u0435\u044e\u0441\u044c, \u043a\u043e\u043c\u0443-\u0442\u043e \u0431\u044b\u043b\u043e \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043d\u043e, \u0441\u043f\u0430\u0441\u0438\u0431\u043e \u0437\u0430 \u0432\u043d\u0438\u043c\u0430\u043d\u0438\u0435.<\/p>\n
\u0418\u0441\u0442\u043e\u0447\u043d\u0438\u043a: habr.com<\/a><\/p>\n","protected":false,"gt_translate_keys":[{"key":"rendered","format":"html"}]},"excerpt":{"rendered":"\u041f\u0440\u0438\u0432\u0435\u0442 \u0425\u0430\u0431\u0440. \u041d\u0430\u0432\u0435\u0440\u043d\u043e\u0435 \u043a\u0430\u0436\u0434\u044b\u0439, \u043a\u0442\u043e \u0445\u043e\u0442\u044c \u0440\u0430\u0437 \u0432\u0441\u0442\u0440\u0435\u0447\u0430\u043b \u0438\u043b\u0438 \u043f\u0440\u043e\u0432\u043e\u0436\u0430\u043b \u0440\u043e\u0434\u0441\u0442\u0432\u0435\u043d\u043d\u0438\u043a\u043e\u0432 \u0438\u043b\u0438 \u0434\u0440\u0443\u0437\u0435\u0439 \u043d\u0430 \u0441\u0430\u043c\u043e\u043b\u0435\u0442, \u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u043b\u0441\u044f \u0431\u0435\u0441\u043f\u043b\u0430\u0442\u043d\u044b\u043c \u0441\u0435\u0440\u0432\u0438\u0441\u043e\u043c Flightradar24. \u042d\u0442\u043e \u0432\u0435\u0441\u044c\u043c\u0430 \u0443\u0434\u043e\u0431\u043d\u044b\u0439 \u0441\u043f\u043e\u0441\u043e\u0431 \u043e\u0442\u0441\u043b\u0435\u0436\u0438\u0432\u0430\u043d\u0438\u044f \u043f\u043e\u043b\u043e\u0436\u0435\u043d\u0438\u044f \u0441\u0430\u043c\u043e\u043b\u0435\u0442\u0430 \u0432 \u0440\u0435\u0430\u043b\u044c\u043d\u043e\u043c \u0432\u0440\u0435\u043c\u0435\u043d\u0438. \u0412 \u043f\u0435\u0440\u0432\u043e\u0439 \u0447\u0430\u0441\u0442\u0438 \u0431\u044b\u043b \u043e\u043f\u0438\u0441\u0430\u043d \u043f\u0440\u0438\u043d\u0446\u0438\u043f \u0440\u0430\u0431\u043e\u0442\u044b \u0442\u0430\u043a\u043e\u0433\u043e \u043e\u043d\u043b\u0430\u0439\u043d-\u0441\u0435\u0440\u0432\u0438\u0441\u0430. \u0421\u0435\u0439\u0447\u0430\u0441 \u043c\u044b \u043f\u043e\u0439\u0434\u0435\u043c \u0434\u0430\u043b\u044c\u0448\u0435, \u0438 \u0432\u044b\u044f\u0441\u043d\u0438\u043c, \u043a\u0430\u043a\u0438\u0435 \u0434\u0430\u043d\u043d\u044b\u0435 \u043f\u0435\u0440\u0435\u0434\u0430\u044e\u0442\u0441\u044f \u0438 \u043f\u0440\u0438\u043d\u0438\u043c\u0430\u044e\u0442\u0441\u044f \u043e\u0442 \u0432\u043e\u0437\u0434\u0443\u0448\u043d\u043e\u0433\u043e \u0441\u0443\u0434\u043d\u0430 \u043a \u043f\u0440\u0438\u0435\u043c\u043d\u043e\u0439 \u0441\u0442\u0430\u043d\u0446\u0438\u0438, […]<\/p>\n","protected":false,"gt_translate_keys":[{"key":"rendered","format":"html"}]},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[702],"tags":[],"class_list":["post-42151","post","type-post","status-publish","format-standard","hentry","category-novosti-interneta"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t
![\"Flightradar24](\"\/wp-content\/uploads\/2019\/03\/f06537ac66ba8f63752cc36083cf5f26.png\")