210 lines
1.5 MiB
Plaintext
210 lines
1.5 MiB
Plaintext
|
{
|
||
|
|
"cells": [
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 1,
|
||
|
|
"id": "f39816ec-dfd0-4d31-b068-699f97c12fd9",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"# 哨兵相干性计算说明\n",
|
||
|
|
"from osgeo import gdal\n",
|
||
|
|
"import numpy as np\n",
|
||
|
|
"import os\n",
|
||
|
|
"import re\n",
|
||
|
|
"from matplotlib import pyplot as plt"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 2,
|
||
|
|
"id": "462aed17-93a7-4c5f-aac1-19b2ca82d017",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"def read_tiff_dataset(tiff_data,band_idx=0):\n",
|
||
|
|
" dataset=gdal.Open(tiff_data,gdal.GA_ReadOnly)\n",
|
||
|
|
" data=dataset.ReadAsArray(0,0,dataset.RasterXSize,dataset.RasterYSize)\n",
|
||
|
|
" del dataset\n",
|
||
|
|
" return data"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 18,
|
||
|
|
"id": "116ea506-ce66-46c7-b630-efe9490f1679",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"master shape:\t (13644, 24731)\n",
|
||
|
|
"salver shape:\t (13644, 24731)\n"
|
||
|
|
]
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"# 哨兵数据\n",
|
||
|
|
"master_path=r\"F:\\Djiuzhaigou\\S1A_IW_SLC__1SDV_20170818T230347_20170818T230414_017984_01E2E5_8AA2.SAFE\\measurement\\s1a-iw3-slc-vv-20170818t230347-20170818t230412-017984-01e2e5-006.tiff\"\n",
|
||
|
|
"salver_path=r\"F:\\Djiuzhaigou\\S1A_IW_SLC__1SDV_20170806T230347_20170806T230414_017809_01DD96_76EA.SAFE\\measurement\\s1a-iw3-slc-vv-20170806t230347-20170806t230412-017809-01dd96-006.tiff\"\n",
|
||
|
|
"master_data=read_tiff_dataset(master_path)\n",
|
||
|
|
"salver_data=read_tiff_dataset(salver_path)\n",
|
||
|
|
"print(\"master shape:\\t\",master_data.shape)\n",
|
||
|
|
"print(\"salver shape:\\t\",salver_data.shape)"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 19,
|
||
|
|
"id": "1a00e820-dd5f-4670-be14-df83ee5ad9a3",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"(13644, 24731) complex64\n",
|
||
|
|
"(13644, 24731) complex64\n",
|
||
|
|
"(13644, 24731) float32\n"
|
||
|
|
]
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"# 做干涉\n",
|
||
|
|
"master_salver=master_data*np.conj(salver_data)\n",
|
||
|
|
"# 输出二进制图像\n",
|
||
|
|
"print(master_salver.shape,master_salver.dtype)\n",
|
||
|
|
"writer = open(r\"D:\\forestHeight\\Result\\master_salver.bin\",\"wb\")\n",
|
||
|
|
"writer.write(master_salver)\n",
|
||
|
|
"writer.close()\n",
|
||
|
|
"print(master_salver.shape,master_salver.dtype)\n",
|
||
|
|
"master_salver_phi=np.arctan2(master_salver.imag,master_salver.real)\n",
|
||
|
|
"writer = open(r\"D:\\forestHeight\\Result\\master_salver_phi.bin\",\"wb\")\n",
|
||
|
|
"writer.write(master_salver_phi)\n",
|
||
|
|
"writer.close()\n",
|
||
|
|
"print(master_salver_phi.shape,master_salver_phi.dtype)\n",
|
||
|
|
"# 数据建议在 ENVI 中查看看\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 39,
|
||
|
|
"id": "06678ed3-5e3e-4c91-8699-35e690baf55f",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"ename": "AttributeError",
|
||
|
|
"evalue": "'TarInfo' object has no attribute 'find'",
|
||
|
|
"output_type": "error",
|
||
|
|
"traceback": [
|
||
|
|
"\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
|
||
|
|
"\u001b[1;31mAttributeError\u001b[0m Traceback (most recent call last)",
|
||
|
|
"Cell \u001b[1;32mIn [39], line 6\u001b[0m\n\u001b[0;32m 4\u001b[0m \u001b[38;5;28;01mwith\u001b[39;00m tarfile\u001b[38;5;241m.\u001b[39mopen(tar_path, \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mr\u001b[39m\u001b[38;5;124m\"\u001b[39m) \u001b[38;5;28;01mas\u001b[39;00m file:\n\u001b[0;32m 5\u001b[0m \u001b[38;5;28;01mfor\u001b[39;00m i \u001b[38;5;129;01min\u001b[39;00m file\u001b[38;5;241m.\u001b[39mgetmembers():\n\u001b[1;32m----> 6\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m \u001b[43mi\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mfind\u001b[49m(\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mmeta.xml\u001b[39m\u001b[38;5;124m\"\u001b[39m)\u001b[38;5;241m>\u001b[39m\u001b[38;5;241m0\u001b[39m:\n\u001b[0;32m 7\u001b[0m xml_text\u001b[38;5;241m=\u001b[39m\u001b[38;5;28;01mNone\u001b[39;00m\n\u001b[0;32m 8\u001b[0m \u001b[38;5;28;01mwith\u001b[39;00m TarIO\u001b[38;5;241m.\u001b[39mTarIO(tar_path, i) \u001b[38;5;28;01mas\u001b[39;00m fp:\n\u001b[0;32m 9\u001b[0m \u001b[38;5;66;03m# 解析坐标范围\u001b[39;00m\n",
|
||
|
|
"\u001b[1;31mAttributeError\u001b[0m: 'TarInfo' object has no attribute 'find'"
|
||
|
|
]
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
" "
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 11,
|
||
|
|
"id": "29d59868-768e-45da-bd3f-566e10eb10bf",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAsUAAAGOCAYAAAB/vH4qAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOy9eXxU9b3//zznzGSSEEQUQTKTjR1RCQmEJKxuJCQhbIpaahetbbXaa7Xf7vdq29+99naxq0oXq1elIoqREAIpiiCQhISQKCI7ZJkJuCACIdvMOZ/fH59zPjOjXbxX1Krn9XjkAZmcmTnns74/7/fr/XprQgiBCxcuXLhw4cKFCxefYugf9Q24cOHChQsXLly4cPFRwzWKXbhw4cKFCxcuXHzq4RrFLly4cOHChQsXLj71cI1iFy5cuHDhwoULF596uEaxCxcuXLhw4cKFi089XKPYhQsXLly4cOHCxacerlHswoULFy5cuHDh4lMP1yh24cKFCxcuXLhw8amHaxS7cOHChQsXLly4+NTDNYpduHDxL4XOzk7uueceWlpaPupbOSvQNI177rnno76N941HHnkETdPYsWPHP7129uzZzJ49+4O/KRcuXLg4i/B81DfgwoULF7Ho7Ozkhz/8IZmZmWRnZ3/Ut+Pi/4AHHnjgo74FFy5cuPhfwzWKXbhw8alAT08PiYmJaJr2Ud/KB45wOIymaXg8H80Sf9FFF30k3+vChQsX7wcufcKFCxcA3HPPPWiaxssvv8w111zDoEGDOO+887jzzjuJRCLs27eP4uJiBg4cSGZmJj/96U/Ve3t7e7nrrrvIzs5W7ysoKGD16tXv+p6nnnqKqVOnMmjQIJKTkxkxYgQ33ngjAJs2bWLKlCkAfPGLX0TTtHfRD3bs2EF5eTnnnXceiYmJTJo0iZUrV8Z9hxPq/+tf/8qNN97IBRdcQHJyMn19fe+pLR588EEmTpxISkoKAwcOZNy4cXzve99Tf3/jjTe49dZbueiii0hJSWHo0KFcfvnlbNmy5R9+7ksvvYSmaTz00EPv+tu6devQNI3Kykr12oEDB/jMZz7D0KFD8fl8jB8/nvvvvz/ufZs2bULTNB577DHuuusu/H4/Pp+PgwcPvqdn1TSN2267jd///veMGTMGn8/HRRddxIoVK/7m9adPn+aWW25hyJAhnH/++SxatIjOzs64a1z6hAsXLj6OcI1iFy5cxGHJkiVMnDiRVatWcfPNN/PLX/6Sb3zjGyxYsIDS0lIqKiq4/PLL+fa3v80zzzwDQF9fH2+99Rbf/OY3efbZZ3niiSeYPn06ixYt4tFHH1WfXVdXx7XXXsuIESNYsWIFa9eu5T/+4z+IRCIA5OTk8PDDDwPwgx/8gLq6Ourq6vjSl74EwAsvvMC0adN4++23WbZsGatXryY7O5trr72WRx555F3PcuONN+L1ennsscd4+umn8Xq9//T5V6xYwa233sqsWbOoqKjg2Wef5Rvf+AZnzpxR17z11lsA3H333axdu5aHH36YESNGMHv2bDZt2vR3P3vixIlMmjRJPWMsHnnkEYYOHUpJSQkAr776KlOmTOGVV17hF7/4BVVVVZSWlvL1r3+dH/7wh+96/3e/+13a29tZtmwZa9asYejQof/0WR1UVlbym9/8hh/96Ec8/fTTZGRkcP311/P000+/69ovfelLeL1e/vKXv/DTn/6UTZs28dnPfvY9f5cLFy5c/MtCuHDhwoUQ4u677xaA+MUvfhH3enZ2tgDEM888o14Lh8PiggsuEIsWLfqbnxWJREQ4HBY33XSTmDRpknr95z//uQDE22+//Xfvo7GxUQDi4Ycfftffxo0bJyZNmiTC4XDc62VlZWL48OHCNE0hhBAPP/ywAMTnPve5f/rc78Rtt90mzj333P/Ve5znveKKK8TChQvj/gaIu+++W/3+m9/8RgBi37596rW33npL+Hw+cdddd6nXioqKRCAQECdPnnzX/SUmJoq33npLCCHECy+8IAAxc+bM/9U9x95fUlKSOHbsWNzzjBs3TowaNUq95rTprbfeGvf+n/70pwIQR48eVa/NmjVLzJo16/90Py5cuPjg0NPTI06ePHlWfnp6ej7qxznrcD3FLly4iENZWVnc7+PHj0fTNObOnate83g8jBo1ira2NvXaU089xbRp00hJScHj8eD1ennooYfYs2ePusahRixZsoSVK1cSCoXe830dPHiQvXv3snTpUgAikYj6KSkp4ejRo+zbty/uPYsXL37vD24jLy+Pt99+m+uvv57Vq1fz5ptv/s3rli1bRk5ODomJiep5n3/++bjn/VtYunQpPp8vzrP9xBNP0NfXxxe/+EVA0lGef/55Fi5cSHJy8ruetbe3l/r6+vf9rA6uuOIKhg0bpn43DINrr72WgwcPEgwG464tLy+P+/3SSy8FiBsLLly4+NdDb28vWRkpDBo06Kz8ZGVl0dvb+1E/1lmFaxS7cOEiDuedd17c7wkJCSQnJ5OYmPiu150F8ZlnnmHJkiX4/X4ef/xx6urqaGxs5MYbb4xbNGfOnMmzzz5LJBLhc5/7HIFAgIsvvpgnnnjin97Xa6+9BsA3v/lNvF5v3M+tt94K8C4Ddvjw4f/r57/hhhv485//TFtbG4sXL2bo0KFMnTqVDRs2qGvuu+8+brnlFqZOncqqVauor6+nsbGR4uJienp6/uHnn3feeZSXl/Poo49imiYgqRN5eXlMmDABgOPHjxOJRPjtb3/7rmd16BVn41kdXHjhhX/3tePHj8e9fv7558f97vP5AP7pc7tw4eKjRX9/P8deNznSlMHx/Vnv6+dIUwbHjh2jv7//o36sswpXfcKFCxfvG48//jhZWVk8+eSTceoOfyuxbf78+cyfP5++vj7q6+u59957+cxnPkNmZiYFBQV/9zuGDBkCSO7sokWL/uY1Y8eOjfv9/6o08cUvfpEvfvGLnDlzhhdffJG7776bsrIy9u/fT0ZGBo8//jizZ8/mwQcfjHvf6dOn3/PnP/XUU2zYsIH09HQaGxvjPmvw4MEYhsENN9zA1772tb/5GVlZWXG/vx9VjWPHjv3d195pBLtw4eLjjXMG6pwz0PWJ/i24RrELFy7eNzRNIyEhIc4wO3bs2N9Un3Dg8/mYNWsW5557LjU1NTQ3N1NQUPB3PY9jx45l9OjRvPTSS/zXf/3XB/Mg78CAAQOYO3cu/f39LFiwgN27d5ORkYGmaeo+Hbz88svU1dWRlpb2Tz93zpw5+P1+Hn74YdLT00lMTOT6669Xf09OTuayyy6jubmZSy+9lISEhLP+bLF4/vnnee211xSFwjRNnnzySUaOHEkgEPhAv9uFCxcfLkxhYYr3/xmfRLhGsQsXLt43ysrKeOaZZ7j11lu5+uqr6ejo4Mc//jHDhw/nwIED6rr/+I//IBgMcsUVVxAIBHj77bf59a9/jdfrZdasWQCMHDmSpKQkli9fzvjx40lJSSE1NZXU1FR+//vfM3fuXIqKivjCF76A3+/nrbfeYs+ePezcuZOnnnrqfT/LzTffTFJSEtOmTWP48OEcO3aMe++9l0GDBilOdFlZGT/+8Y+5++67mTVrFvv27eNHP/oRWVlZSknjH8EwDD73uc9x3333cc4557Bo0SIGDRoUd82vf/1rpk+fzowZM7jlllvIzMzk9OnTHDx4kDVr1rBx48b3/awOhgwZwuWXX86///u/M2DAAB544AH27t37d2XZXLhw8fGFhcDi/VnF7/f9/6pw/ecuXLh43/jiF7/IT37yE9atW0dJSQn//d//zXe+8x0+85nPxF03depUjh07xre//W3mzJnDl7/8ZZKSkti4caPi0yYnJ/PnP/+Z48ePM2fOHKZMmcIf/vAHAC677DIaGho499xzueOOO7jyyiu55ZZbeO6557jyyivPyrPMmDGDV155hX/7t3/jqquu4hvf+AZjxoxhy5YtXHDBBQB8//vf56677uKhhx6itLSUP/3pTyxbtozp06f/r9qsr6+PN954QyXYxeKiiy5i586dXHzxxfzgBz9gzpw53HTTTTz
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 870.044x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"# 展示强度图--master\n",
|
||
|
|
"plt.matshow(master_salver_phi)\n",
|
||
|
|
"plt.colorbar()\n",
|
||
|
|
"plt.title(\"master_salver_phi\")\n",
|
||
|
|
"plt.show()"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 7,
|
||
|
|
"id": "d37c87bc-f965-44c4-9b02-0aeab694a3bb",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAArUAAAGXCAYAAACtLKjhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOy9e1xNafv4/2427ansKUoRjVJEpKjRiChFI0JEI5JDRHLIIDnkfIjImYicaRxCRJQiyp4pZaIRRWlETZndpNixn/37Yzf79/V9ns/nmfnO0cx6v17rVXvte611r3XttdZ1X/d10FAqlUoEBAQEBAQEBAQE3mM++LM7ICAgICAgICAgIPBrEZRaAQEBAQEBAQGB9x5BqRUQEBAQEBAQEHjvEZRaAQEBAQEBAQGB9x5BqRUQEBAQEBAQEHjvEZRaAQEBAQEBAQGB9x5BqRUQEBAQEBAQEHjvEZRaAQEBAQEBAQGB9x5BqRUQEBAQEBAQEHjvEZRaAQEBgd+IsrIyli5dSm5u7p/dFQEBAYF/HIJSKyAgIPAbUVZWxrJlywSlVkBAQOBPQFBqBQQEBP7ivHr1CqVS+Wd3Q0BAQOAvjaDUCggIvPcsXboUDQ0NvvnmG0aMGIGuri7NmjVj9uzZvH37loKCAj777DMkEgmmpqasW7dOve3r16/54osvsLW1VW/Xo0cPzp49+2/HOXHiBA4ODujq6qKtrU3btm2ZMGECAGlpaXzyyScAjB8/Hg0NDTQ0NFi6dKl6+6ysLAYPHkyzZs348MMP6dq1K19++eU7x9i/fz8aGhpcvnyZCRMm0Lx5c7S1tZHL5b/DlRMQEBD4+yAotQICAn8bRo4ciY2NDadOnWLSpElERUUREhLC0KFDGThwIPHx8fTt25fQ0FBOnz4NgFwu58WLF8yZM4czZ85w7NgxevXqxbBhwzh48KB635mZmfj4+NC2bVuOHz/OhQsXCA8P5+3btwB069aN2NhYABYtWkRmZiaZmZkEBAQAkJqaSs+ePZHJZOzatYuzZ89ia2uLj48P+/fv/7dzmTBhAo0bN+bQoUOcPHmSxo0b/85XT0BAQOA9RykgICDwnrNkyRIloNywYcM7621tbZWA8vTp0+p1b968UTZv3lw5bNiw/7ivt2/fKt+8eaOcOHGismvXrur1kZGRSkApk8n+x358/fXXSkAZGxv7b9916NBB2bVrV+WbN2/eWT9o0CBly5YtlQqFQqlUKpWxsbFKQDl27Nj/et4CAgICAv8/jf5UjVpAQEDgN2TQoEHvfO7YsSN37txhwIAB6nWNGjXCwsKCkpIS9boTJ06wadMm7ty5Q21trXr9hx9+qP7/J9eCkSNHMnHiRHr27EmrVq1+Vr8KCwu5f/8+kZGRAGrrLoCHhwfnz5+noKCAjh07qtcPHz78Z+1bQEBA4Ofy+vVr6uvrf5N9aWpqvvOM/CsguB8ICAj8bWjWrNk7nzU1NdHW1v63B6+mpiavX78G4PTp04wcOZJWrVpx+PBhMjMz+frrr5kwYYK6DUDv3r05c+YMb9++ZezYsbRu3ZrOnTtz7Nix/9qv8vJyAObMmUPjxo3fWYKCggCorKx8Z5uWLVv+8gsgICAg8D/w+vVrmmtpoaur+5ssZmZm7zwj/woIlloBAYF/NIcPH8bMzIy4uDg0NDTU6/9TYNaQIUMYMmQIcrmcW7dusWbNGnx9fTE1NaVHjx7/4zEMDAwACAsLY9iwYf+xjaWl5Tuf/8++CAgICPxa6uvreQnMBcS/cl9yYP3z59TX1/+lrLWCUisgIPCPRkNDA01NzXeUyOfPn//H7Ac/IRaL6dOnD3p6eiQlJZGTk0OPHj0Qi1WvilevXr3T3tLSknbt2nHnzh1Wr179+5yIgICAwM+gCfBr1dC/atiqoNQKCAj8oxk0aBCnT58mKCgIb29vSktLWbFiBS1btuThw4fqduHh4Xz33Xe4urrSunVrZDIZmzdvpnHjxvTp0wcAc3NztLS0OHLkCB07dqRJkyYYGxtjbGxMdHQ0AwYMwN3dnXHjxtGqVStevHjBt99+y+3btzlx4sSfdQkEBAQE/hYIPrUCAgL/aMaPH8/atWu5ePEiHh4eREREMH/+fHx9fd9p5+DgwPPnzwkNDaV///5MnjwZLS0trl69SqdOnQDQ1tZm3759VFVV0b9/fz755BN2794NgIuLC1999RV6enrMmjULNzc3pk6dSnJyMm5ubn/4eQsICPwzafQbLX9FNJRKoUyNgICAgICAgMDfmR9//BFdXV0iAK1fua9XQChQXV3NRx999Os79xshWGoFBAQEBAQEBATee/6qFmQBAQEBAQEBAYHfmN/CfeCvqjz+VfslICAgICAgICDwG9OYX5+94O1/b/KnILgfCAgICAgICAgIvPcIlloBAQEBAQEBgX8IgvuBgICAgICAgIDAe08jfr37wZvfoiO/A4L7gYCAgICAgICAwO/Gzp076dKlCx999BEfffQRPXr04OLFi+rvx40bh4aGxjvLp59++ouPI1hqBQQEBAQEBAT+IfwZ7getW7dm7dq1WFhYAHDgwAGGDBlCTk6OunjNZ599RmxsrHobTU3NX9yvf7SldseOHZiZmfHhhx9iZ2dHenr6n92lvz1Lly79t9FYixYt1N8rlUqWLl2KsbExWlpaODs7c+/evXf2IZfLmT59OgYGBujo6DB48GC+++67d9r88MMP+Pn5oauri66uLn5+fshksj/iFP8WXL9+HU9PT4yNjdHQ0ODMmTPvfP9HyunJkyd4enqio6ODgYEBM2bMoL6+/vc47fea/yazn2MJEWT2x7FmzRo++eQTJBIJhoaGDB06lIKCgnfaCPeZwO9B499o+SV4enri4eFB+/btad++PatWraJJkybcunVL3UYsFtOiRQv10qxZs198bv9YpTYuLo5Zs2axcOFCcnJycHJyYsCAATx58uTP7trfnk6dOvHs2TP1kpeXp/5u3bp1bNy4kW3btvH111/TokUL+vXrR01NjbrNrFmziI+P5/jx49y4cYOXL18yaNAgFAqFuo2vry+5ublcunSJS5cukZubi5+f3x96nu8ztbW12NjYsG3btv/4/R8lJ4VCwcCBA6mtreXGjRscP36cU6dO8cUXX/x+J/+e8t9kBipLyP957yUmJr7zvSCzP45r164xbdo0bt26xZUrV3j79i39+/entrZW3Ua4zwR+D35LpfbHH398Z5HL5f/1+AqFguPHj1NbW0uPHj3U69PS0jA0NKR9+/ZMmjSJioqKX35yyn8o3bt3V06ZMuWddR06dFDOnz//T+rRP4MlS5YobWxs/uN3//rXv5QtWrRQrl27Vr3u9evXSl1dXeWuXbuUSqVSKZPJlI0bN1YeP35c3ebp06fKDz74QHnp0iWlUqlU5ufnKwHlrVu31G0yMzOVgPL+/fu/w1n9vQGU8fHx6s9/pJwSExOVH3zwgfLp06fqNseOHVOKxWJldXX173K+fwf+b5kplUqlv7+/csiQIf/jNoLM/lwqKiqUgPLatWtKpVK4zwR+e6qrq5WA8jgoz/3K5Tgo+Q/LkiVL/sfjf/PNN0odHR2lSCRS6urqKi9cuKD+7vjx48rz588r8/LylOfOnVPa2NgoO3XqpHz9+vUvOsd/pKW2vr6e7Oxs+vfv/876/v37k5GR8Sf16p/Dw4cPMTY2xszMjM8//5xHjx4B8PjxY54/f/6OXMRiMX369FHLJTs7mzdv3rzTxtjYmM6dO6vbZGZmoquri4ODg7rNp59+iq6uriDf34A/Uk6ZmZl07twZY2NjdRt3d3fkcjnZ2dm/63n+HfnfLCGCzP5cqqurAdRTrsJ9JvB70eg3WgBKS0uprq5WL2FhYf/jcS0tLcnNzeXWrVtMnToVf39/8vPzAfDx8WHgwIF07twZT09PLl68yIMHD7hw4cIvPrd/HJWVlSgUCoyMjN5Zb2RkxPPnz/+kXv0zcHBw4ODBg7Rv357y8nJWrlyJo6Mj9+7dU1/7/ySXkpISAJ4/f46mpiZNmzb9tzY/bf/8+XMMDQ3/7diGhoa
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 870.044x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAArUAAAGXCAYAAACtLKjhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOy9e1xNefv//2w27ak0NaLYaWRCJho1NfaIFEUjQkREcohIDhE5n88iZxo5JKIRRSOiFJHZphQRoZGJqJGpSTU79r1/f+ys38fnvu/PPfMdM8M96/l4rAet/d5rv9d+77XW9b7e1/W6tNRqtRoRERERERERERGRd5j3/uoOiIiIiIiIiIiIiPxeRKNWRERERERERETknUc0akVERERERERERN55RKNWRERERERERETknUc0akVERERERERERN55RKNWRERERERERETknUc0akVERERERERERN55RKNWRERERERERETknUc0akVERERERERERN55RKNWRERE5A1RUlLC4sWLyc3N/au7IiIiIvK3QzRqRURERN4QJSUlLFmyRDRqRURERP4CRKNWRERE5C2ntrYWtVr9V3dDRERE5K1GNGpFRETeeRYvXoyWlhbXr19n8ODBGBgY0LhxY6ZPn87Lly8pKCjgyy+/RF9fH3Nzc9auXSu895dffmHGjBnY2NgI7+vcuTPHjx//p885cuQIcrkcAwMDdHV1+fjjjxkzZgwA6enpfP755wCMHj0aLS0ttLS0WLx4sfD+rKws+vXrR+PGjXn//fextbXl66+/fu0z9u3bh5aWFmfOnGHMmDE0bdoUXV1dlErlH/DNiYiIiPz3IBq1IiIi/zUMGTKEjh07cvToUcaNG0d4eDjBwcEMGDCAPn36EB8fT48ePQgNDeXYsWMAKJVKnj17RkhICAkJCRw6dIiuXbsycOBA9u/fLxz78uXLeHt78/HHH3P48GFOnjzJwoULefnyJQCfffYZe/fuBWD+/PlcvnyZy5cv4+/vD0BaWhpdunShoqKCnTt3cvz4cWxsbPD29mbfvn3/dC5jxoyhYcOGREdHExcXR8OGDf/gb09ERETkHUctIiIi8o6zaNEiNaBev379a/ttbGzUgPrYsWPCvhcvXqibNm2qHjhw4L881suXL9UvXrxQjx07Vm1rayvsDwsLUwPqioqKf9uP7777Tg2o9+7d+0+vtWvXTm1ra6t+8eLFa/v79u2rbt68uVqlUqnVarV67969akA9cuTI/3jeIiIiIiL/Pw3+UotaRERE5A3St2/f1/7+5JNPuHbtGr179xb2NWjQgNatW/PgwQNh35EjR9i4cSPXrl2jurpa2P/+++8L/38VWjBkyBDGjh1Lly5dMDU1/VX9unfvHrdv3yYsLAxA8O4CuLu7880331BQUMAnn3wi7B80aNCvOraIiIjIr+WXX36hrq7ujRxLW1v7tXvk24AYfiAiIvJfQ+PGjV/7W1tbG11d3X+68Wpra/PLL78AcOzYMYYMGYKpqSkHDhzg8uXLfPfdd4wZM0ZoA9CtWzcSEhJ4+fIlI0eOpEWLFnTo0IFDhw79x36VlpYCEBISQsOGDV/bAgMDAXj69Olr72nevPlv/wJERERE/g2//PILTXV0MDAweCNbq1atXrtHvg2InloREZG/NQcOHKBVq1bExsaipaUl7P9XiVn9+/enf//+KJVKvv32W1atWoWPjw/m5uZ07tz5335GkyZNAJgzZw4DBw78l20sLS1f+/t/9kVERETk91JXV8dzYCYg/Z3HUgLrnjyhrq7urfLWikatiIjI3xotLS20tbVfMyKfPHnyL9UPXiGVSnFycsLQ0JDk5GRycnLo3LkzUqnmUVFbW/tae0tLS9q0acO1a9dYuXLlH3MiIiIiIr+CRsDvNUPf1rRV0agVERH5W9O3b1+OHTtGYGAgXl5eFBcXs2zZMpo3b87du3eFdgsXLuThw4e4uLjQokULKioq2LRpEw0bNsTJyQkACwsLdHR0OHjwIJ988gmNGjVCJpMhk8mIiIigd+/euLm5MWrUKExNTXn27Bm3bt3i6tWrHDly5K/6CkRERET+KxBjakVERP7WjB49mtWrV3Pq1Cnc3d1Zs2YNs2fPxsfH57V2crmcJ0+eEBoaSq9evRg/fjw6OjqcO3eO9u3bA6Crq8uePXsoLy+nV69efP7553z11VcAdO/enStXrmBoaMi0adNwdXVl4sSJpKSk4Orq+qeft4iIyN+TBm9oexvRUqvFMjUiIiIiIiIiIv/N/PzzzxgYGLAG0Pmdx6oFQoHKyko++OCD39+5N4ToqRUREREREREREXnneVs9yCIiIiIiIiIiIm+YNxE+8LYaj29rv0RERERERERERN4wDfn96gUv/3OTvwQx/EBEREREREREROSdR/TUioiIiIiIiIj8TRDDD0RERERERERERN55GvD7ww9evImO/AGI4QciIiIiIiIiIiLvPKJRKyIiIiIiIiLyN+GvKL6wY8cOPv30Uz744AM++OADOnfuzKlTp4TXR40ahZaW1mvbF1988ZvP7W9t1G7fvp1WrVrx/vvvY2dnR0ZGxl/dpf96Fi9e/E8/3GbNmgmvq9VqFi9ejEwmQ0dHB2dnZ27evPnaMZRKJZMnT6ZJkybo6enRr18/Hj58+Fqbn376CV9fXwwMDDAwMMDX15eKioo/4xT/K7hw4QIeHh7IZDK0tLRISEh47fU/c5x++OEHPDw80NPTo0mTJkyZMoW6uro/4rTfaf7TmP2ah4Y4Zn8eq1at4vPPP0dfXx9jY2MGDBhAQUHBa23E60zkj6DhG9p+Cy1atGD16tVkZWWRlZVFjx496N+//2u/5y+//JLHjx8LW1JS0m8+t7+tURsbG8u0adOYN28eOTk5ODo60rt3b3744Ye/umv/9bRv3/61H25eXp7w2tq1a9mwYQNbt27lu+++o1mzZvTs2ZOqqiqhzbRp04iPj+fw4cNcvHiR58+f07dvX1QqldDGx8eH3NxcTp8+zenTp8nNzcXX1/dPPc93merqajp27MjWrVv/5et/1jipVCr69OlDdXU1Fy9e5PDhwxw9epQZM2b8cSf/jvKfxgz+80NDHLM/j/PnzzNp0iS+/fZbzp49y8uXL+nVqxfV1dVCG/E6E/kj+CuMWg8PD9zd3Wnbti1t27ZlxYoVNGrUiG+//VZoI5VKadasmbA1btz4t5+c+m9Kp06d1BMmTHhtX7t27dSzZ8/+i3r092DRokXqjh07/svX/vGPf6ibNWumXr16tbDvl19+URsYGKh37typVqvV6oqKCnXDhg3Vhw8fFto8evRI/d5776lPnz6tVqvV6vz8fDWg/vbbb4U2ly9fVgPq27dv/wFn9d8NoI6Pjxf+/jPHKSkpSf3ee++pHz16JLQ5dOiQWiqVqisrK/+Q8/1v4H+PmVqtVvv5+an79+//b98jjtlfS1lZmRpQnz9/Xq1Wi9eZyJunsrJSDagPg/rE79wOgxr4f/p9vHz5Un3o0CG1tra2+ubNm2q1WnN/MjAwUDdt2lTdpk0btb+/v7q0tPQ3H/tv6amtq6sjOzubXr16vba/V69eZGZm/kW9+vtw9+5dZDIZrVq1YujQoXz//fcA3L9/nydPnrw2LlKpFCcnJ2FcsrOzefHixWttZDIZHTp0ENpcvnwZAwMD5HK50OaLL77AwMBAHN83wJ85TpcvX6ZDhw7IZDKhjZubG0qlkuzs7D/0PP8bSU9Px9jYmLZt2zJu3DjKysqE18Qx+2uprKwEELxT4nUm8kfxJmNqf/7559c2pVL5bz83Ly+PRo0aIZVKmTBhAvHx8VhZWQHQu3dvDh48yLlz51i/fj3fffcdPXr0+D+P96/4Wxq1T58+RaVSYWJi8tp+ExMTnjx58hf16u+BXC5n//79JCcns2vXLp48eYKDgwPl5eXCd/9/jcuTJ0/Q1tbmww8//D/bGBsb/9NnGxsbi+P7Bvgzx+nJkyf/9Dkffvgh2tra4lj+Rv7TQ0Mcs78OtVrN9OnT6dq1Kx06dADE60zkj+OVpNfv2V4ZtWZ
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 870.044x480 with 2 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"# 展示强度图--master\n",
|
||
|
|
"master_power_data=10*np.log10((master_data.real**2+master_data.imag**2)**0.5+1)\n",
|
||
|
|
"plt.matshow(master_power_data,cmap='jet')\n",
|
||
|
|
"plt.colorbar()\n",
|
||
|
|
"plt.title(\"master\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"\n",
|
||
|
|
"salver_power_data=10*np.log10((salver_data.real**2+salver_data.imag**2)**0.5+1)\n",
|
||
|
|
"plt.matshow(salver_power_data,cmap='jet')\n",
|
||
|
|
"plt.colorbar()\n",
|
||
|
|
"plt.title(\"master\")\n",
|
||
|
|
"plt.show()\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "c3709429-57b6-4053-acaf-823be7a78b30",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"metadata": {
|
||
|
|
"kernelspec": {
|
||
|
|
"display_name": "Python 3 (ipykernel)",
|
||
|
|
"language": "python",
|
||
|
|
"name": "python3"
|
||
|
|
},
|
||
|
|
"language_info": {
|
||
|
|
"codemirror_mode": {
|
||
|
|
"name": "ipython",
|
||
|
|
"version": 3
|
||
|
|
},
|
||
|
|
"file_extension": ".py",
|
||
|
|
"mimetype": "text/x-python",
|
||
|
|
"name": "python",
|
||
|
|
"nbconvert_exporter": "python",
|
||
|
|
"pygments_lexer": "ipython3",
|
||
|
|
"version": "3.8.13"
|
||
|
|
}
|
||
|
|
},
|
||
|
|
"nbformat": 4,
|
||
|
|
"nbformat_minor": 5
|
||
|
|
}
|